Geometry-aware numerics
High-order finite-volume reconstruction on Cartesian, cylindrical, spherical, and other orthogonal curvilinear grids.
GAMERA-OP
GAMERA-OP (Orthogonal-Plus) is a three-dimensional finite-volume magnetohydrodynamics solver for orthogonal curvilinear geometries. Rewritten in C with a modular design, it provides a practical framework for planetary, space, and astrophysical plasma simulations that require high-order accuracy and robust treatment of curved coordinates.
Core capabilities
Divergence control
Constrained transport preserves the solenoidal magnetic-field condition to machine precision.
Rotating systems
Angular momentum is conserved to round-off in cylindrical and spherical coordinates, supporting rapidly rotating planetary systems.
Extended physics
Options include ring averaging near the axis, semi-relativistic correction, background-field splitting, and anisotropic MHD.
Numerical design
The solver combines geometry-consistent high-order reconstruction with flexible numerical fluxes and time integrators. Its architecture is designed to simplify case setup, addition of new physics, and coupling to other first-principles models. Standard benchmarks across multiple geometries demonstrate accuracy, low numerical diffusion, and robust treatment of coordinate singularities and rotating flows.
Reference
H. Luo, B. Zhang, J. Tian, J. Cai, J. Chen, E. Feng, Z. Zheng, S. Xi, and J. G. Lyon (2026), “GAMERA-OP: A Three-dimensional Finite-volume Magnetohydrodynamic Solver for Orthogonal Curvilinear Geometries,” The Astrophysical Journal Supplement Series, 285(1), 15. https://doi.org/10.3847/1538-4365/ae7344