Adaptive numerical simulations with Trixi.jl: A case study of Julia for scientific computing

Abstract

We present Trixi.jl, a Julia package for adaptive high-order numerical simulations of hyperbolic partial differential equations. Utilizing Julia's strengths, Trixi.jl is extensible, easy to use, and fast. We describe the main design choices that enable these features and compare Trixi.jl with a mature open source Fortran code that uses the same numerical methods. We conclude with an assessment of Julia for simulation-focused scientific computing, an area that is still dominated by traditional high-performance computing languages such as C, C++, and Fortran.

References in corpus (6)

• Forward-Mode Automatic Differentiation in Julia
• A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics
• Preventing pressure oscillations does not fix local linear stability issues of entropy-based split-form high-order schemes
• Efficient implementation of modern entropy stable and kinetic energy preserving discontinuous Galerkin methods for conservation laws
• On a linear stability issue of split form schemes for compressible flows
• Opening the Blackbox: Accelerating Neural Differential Equations by Regularizing Internal Solver Heuristics

Cited by in corpus (23)

• HORSES3D: a high-order discontinuous Galerkin solver for flow simulations and multi-physics applications
• Subcell limiting strategies for discontinuous Galerkin spectral element methods
• Evaluating performance and portability of high-level programming models: Julia, Python/Numba, and Kokkos on exascale nodes
• Efficient implementation of modern entropy stable and kinetic energy preserving discontinuous Galerkin methods for conservation laws
• GALÆXI: Solving complex compressible flows with high-order discontinuous Galerkin methods on accelerator-based systems
• Lax-Wendroff flux reconstruction method for hyperbolic conservation laws
• On error-based step size control for discontinuous Galerkin methods for compressible fluid dynamics
• On the robustness of high-order upwind summation-by-parts methods for nonlinear conservation laws
• Computing with B-series
• Julia as a unifying end-to-end workflow language on the Frontier exascale system
• Many-Stage Optimal Stabilized Runge-Kutta Methods for Hyperbolic Partial Differential Equations
• Generalized upwind summation-by-parts operators and their application to nodal discontinuous Galerkin methods
• Scientific Machine Learning of Flow Resistance Using Universal Shallow Water Equations with Differentiable Programming
• Multirate Time-Integration based on Dynamic ODE Partitioning through Adaptively Refined Meshes for Compressible Fluid Dynamics
• An entropy stable discontinuous Galerkin method for the two-layer shallow water equations on curvilinear meshes
• A discontinuous Galerkin discretization of elliptic problems with improved convergence properties using summation by parts operators
• Structure-Preserving Numerical Methods for Two Nonlinear Systems of Dispersive Wave Equations
• Compact Runge-Kutta flux reconstruction methods for non-conservative hyperbolic equations
• Fourth-Order Paired-Explicit Runge-Kutta Methods
• Multiderivative time integration methods preserving nonlinear functionals via relaxation
• Stability of step size control based on a posteriori error estimates
• A Note on Numerical Fluxes Conserving a Member of Harten's One-Parameter Family of Entropies for the Compressible Euler Equations
• Mimetic Metrics for the DGSEM