3,108 research outputs found
Hyperbolic systems of conservation laws in one space dimension
Aim of this paper is to review some basic ideas and recent developments in
the theory of strictly hyperbolic systems of conservation laws in one space
dimension. The main focus will be on the uniqueness and stability of entropy
weak solutions and on the convergence of vanishing viscosity approximations
A Jacobian-free Newton-Krylov method for time-implicit multidimensional hydrodynamics
This work is a continuation of our efforts to develop an efficient implicit
solver for multidimensional hydrodynamics for the purpose of studying important
physical processes in stellar interiors, such as turbulent convection and
overshooting. We present an implicit solver that results from the combination
of a Jacobian-Free Newton-Krylov method and a preconditioning technique
tailored to the inviscid, compressible equations of stellar hydrodynamics. We
assess the accuracy and performance of the solver for both 2D and 3D problems
for Mach numbers down to . Although our applications concern flows in
stellar interiors, the method can be applied to general advection and/or
diffusion-dominated flows. The method presented in this paper opens up new
avenues in 3D modeling of realistic stellar interiors allowing the study of
important problems in stellar structure and evolution.Comment: Accepted for publication in A&
Fornax: a Flexible Code for Multiphysics Astrophysical Simulations
This paper describes the design and implementation of our new multi-group,
multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a
suite of code tests to validate its application in a wide range of physical
regimes. Instead of focusing exclusively on tests of neutrino radiation
hydrodynamics relevant to the core-collapse supernova problem for which Fornax
is primarily intended, we present here classical and rigorous demonstrations of
code performance relevant to a broad range of multi-dimensional hydrodynamic
and multi-group radiation hydrodynamic problems. Our code solves the
comoving-frame radiation moment equations using the M1 closure, utilizes
conservative high-order reconstruction, employs semi-explicit matter and
radiation transport via a high-order time stepping scheme, and is suitable for
application to a wide range of astrophysical problems. To this end, we first
describe the philosophy, algorithms, and methodologies of Fornax and then
perform numerous stringent code tests, that collectively and vigorously
exercise the code, demonstrate the excellent numerical fidelity with which it
captures the many physical effects of radiation hydrodynamics, and show
excellent strong scaling well above 100k MPI tasks.Comment: Accepted to the Astrophysical Journal Supplement Series; A few more
textual and reference updates; As before, one additional code test include
High-order conservative reconstruction schemes for finite volume methods in cylindrical and spherical coordinates
High-order reconstruction schemes for the solution of hyperbolic conservation
laws in orthogonal curvilinear coordinates are revised in the finite volume
approach. The formulation employs a piecewise polynomial approximation to the
zone-average values to reconstruct left and right interface states from within
a computational zone to arbitrary order of accuracy by inverting a
Vandermonde-like linear system of equations with spatially varying
coefficients. The approach is general and can be used on uniform and
non-uniform meshes although explicit expressions are derived for polynomials
from second to fifth degree in cylindrical and spherical geometries with
uniform grid spacing. It is shown that, in regions of large curvature, the
resulting expressions differ considerably from their Cartesian counterparts and
that the lack of such corrections can severely degrade the accuracy of the
solution close to the coordinate origin. Limiting techniques and monotonicity
constraints are revised for conventional reconstruction schemes, namely, the
piecewise linear method (PLM), third-order weighted essentially non-oscillatory
(WENO) scheme and the piecewise parabolic method (PPM).
The performance of the improved reconstruction schemes is investigated in a
number of selected numerical benchmarks involving the solution of both scalar
and systems of nonlinear equations (such as the equations of gas dynamics and
magnetohydrodynamics) in cylindrical and spherical geometries in one and two
dimensions. Results confirm that the proposed approach yields considerably
smaller errors, higher convergence rates and it avoid spurious numerical
effects at a symmetry axis.Comment: 37 pages, 12 Figures. Accepted for publication in Journal of
Compuational Physic
A dynamically adaptive multigrid algorithm for the incompressible Navier-Stokes equations: Validation and model problems
An algorithm is described for the solution of the laminar, incompressible Navier-Stokes equations. The basic algorithm is a multigrid based on a robust, box-based smoothing step. Its most important feature is the incorporation of automatic, dynamic mesh refinement. This algorithm supports generalized simple domains. The program is based on a standard staggered-grid formulation of the Navier-Stokes equations for robustness and efficiency. Special grid transfer operators were introduced at grid interfaces in the multigrid algorithm to ensure discrete mass conservation. Results are presented for three models: the driven-cavity, a backward-facing step, and a sudden expansion/contraction
Energy preserving model order reduction of the nonlinear Schr\"odinger equation
An energy preserving reduced order model is developed for two dimensional
nonlinear Schr\"odinger equation (NLSE) with plane wave solutions and with an
external potential. The NLSE is discretized in space by the symmetric interior
penalty discontinuous Galerkin (SIPG) method. The resulting system of
Hamiltonian ordinary differential equations are integrated in time by the
energy preserving average vector field (AVF) method. The mass and energy
preserving reduced order model (ROM) is constructed by proper orthogonal
decomposition (POD) Galerkin projection. The nonlinearities are computed for
the ROM efficiently by discrete empirical interpolation method (DEIM) and
dynamic mode decomposition (DMD). Preservation of the semi-discrete energy and
mass are shown for the full order model (FOM) and for the ROM which ensures the
long term stability of the solutions. Numerical simulations illustrate the
preservation of the energy and mass in the reduced order model for the two
dimensional NLSE with and without the external potential. The POD-DMD makes a
remarkable improvement in computational speed-up over the POD-DEIM. Both
methods approximate accurately the FOM, whereas POD-DEIM is more accurate than
the POD-DMD
- …