5,333 research outputs found
ORB5: a global electromagnetic gyrokinetic code using the PIC approach in toroidal geometry
This paper presents the current state of the global gyrokinetic code ORB5 as
an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177
409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of
equations using a PIC scheme and also includes collisions and strong flows. The
code assumes multiple gyrokinetic ion species at all wavelengths for the
polarization density and drift-kinetic electrons. Variants of the physical
model can be selected for electrons such as assuming an adiabatic response or a
``hybrid'' model in which passing electrons are assumed adiabatic and trapped
electrons are drift-kinetic. A Fourier filter as well as various control
variates and noise reduction techniques enable simulations with good
signal-to-noise ratios at a limited numerical cost. They are completed with
different momentum and zonal flow-conserving heat sources allowing for
temperature-gradient and flux-driven simulations. The code, which runs on both
CPUs and GPUs, is well benchmarked against other similar codes and analytical
predictions, and shows good scalability up to thousands of nodes
Thermoacoustic instability - a dynamical system and time domain analysis
This study focuses on the Rijke tube problem, which includes features
relevant to the modeling of thermoacoustic coupling in reactive flows: a
compact acoustic source, an empirical model for the heat source, and
nonlinearities. This thermo-acoustic system features a complex dynamical
behavior. In order to synthesize accurate time-series, we tackle this problem
from a numerical point-of-view, and start by proposing a dedicated solver
designed for dealing with the underlying stiffness, in particular, the retarded
time and the discontinuity at the location of the heat source. Stability
analysis is performed on the limit of low-amplitude disturbances by means of
the projection method proposed by Jarlebring (2008), which alleviates the
linearization with respect to the retarded time. The results are then compared
to the analytical solution of the undamped system, and to Galerkin projection
methods commonly used in this setting. This analysis provides insight into the
consequences of the various assumptions and simplifications that justify the
use of Galerkin expansions based on the eigenmodes of the unheated resonator.
We illustrate that due to the presence of a discontinuity in the spatial
domain, the eigenmodes in the heated case, predicted by using Galerkin
expansion, show spurious oscillations resulting from the Gibbs phenomenon. By
comparing the modes of the linear to that of the nonlinear regime, we are able
to illustrate the mean-flow modulation and frequency switching. Finally,
time-series in the fully nonlinear regime, where a limit cycle is established,
are analyzed and dominant modes are extracted. The analysis of the saturated
limit cycles shows the presence of higher frequency modes, which are linearly
stable but become significant through nonlinear growth of the signal. This
bimodal effect is not captured when the coupling between different frequencies
is not accounted for.Comment: Submitted to Journal of Fluid Mechanic
A note on adjoint error estimation for one-dimensional stationary balance laws with shocks
We consider one-dimensional steady-state balance laws with discontinuous
solutions. Giles and Pierce realized that a shock leads to a new term in the
adjoint error representation for target functionals.This term disappears if and
only if the adjoint solution satisfies an internal boundary condition.
Curiously, most computer codes implementing adjoint error estimation ignore the
new term in the functional, as well as the internal adjoint boundary condition.
The purpose of this note is to justify this omission as follows: if one
represents the exact forward and adjoint solutions as vanishing viscosity
limits of the corresponding viscous problems, then the internal boundary
condition is naturally satisfied in the limit
- …