Version 3 (modified by 14 years ago) ( diff ) | ,
---|
solver.data
Explained
The solver.data
file controls the options used by AstroBEAR's hyperbolic equation integrators. Some of the options in this file are a bit esoteric, so wherever possible the default options will be clearly labeled.
SolverData
This short namelist governs the numerical solvers used by AstroBEAR, as well as the integration scheme used to implement those solvers.
- iScheme: Selects the numerical scheme to run. Currently, we only have option 1 (the sweep method) implemented in AstroBEAR 2.0.
- iSolver:: Selects the solver option to use for this problem. AstroBEAR currently has one exact Riemann solver and five HLL-based approximate solvers.
0 -- Default
: With this option, AstroBEAR selects the appropriate solver option based on the equation of state and whether or not the problem is MHD. If you aren't certain which solver to use, then select this option.1 -- iExactRS
: An exact Riemann solver. Such solvers produce the most accurate solutions to the Riemann problem, but are also the most computationally expensive.2 -- HLLC
: A three-wave version of the HLL solver. A good choice for hydrodynamic problems.3 -- HLL
: The standard HLL solver from Chapter 10 of Toro. HLL uses a two-wave approximation of the intermediate region of the Riemann problem, sacrificing accuracy for speed.4 -- HLLC_ISO
: A variant of the HLLC solver option to be used in isothermal conditions (i.e., wheniEOS = 4
).6 -- HLLD
: A three-wave HLL-derived solver that can resolve discontinuities in MHD. Consequently, it is a good solver to use when magnetic fields are present (lMHD = true
).8 -- HLLD_ISO
: An isothermal variant of HLLD to use theniEOS = 4
.
When in doubt, the HLLC and HLLD solvers are the best choices, as they are faster than the exact solver and more robust than the standard HLL solver. The appropriate HLLx solver choice is determined by two factors: whether or not the problem is MHD, and whether or not the problem is isothermal. To select the appropriate solver based on the physics of your problem, consult the table below:
lMHD iEOS solver F 0 HLLC (2) T 0 HLLD (6) F 4 HLLC_ISO (4) T 4 HLLD_ISO (8)