Related papers: GAMERA-OP: A three-dimensional finite-volume MHD s…
This paper introduces a novel geometric multigrid solver for unstructured curved surfaces. Multigrid methods are highly efficient iterative methods for solving systems of linear equations. Despite the success in solving problems defined on…
We propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modelling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved…
This paper describes the main features of a pioneering unsteady solver for simulating ideal two-fluid plasmas on unstructured grids, taking profit of GPGPU (General-purpose computing on graphics processing units). The code, which has been…
Relativistic macroscopic plasma dynamics can be described by general-relativistic magnetohydrodynamics. In many high-energy astrophysical settings, such as the interior dynamics of magnetized stars, the ideal GRMHD approximation, in which…
We present a novel spectral solver for general relativistic magnetohydrodynamics on dynamical spacetimes. By combining a high order discontinuous spectral method on mapped Chebyshev Fourier grids, our scheme attains exponential convergence.…
The phase-field method has become a useful tool for the simulation of classical metallurgical phase transformations as well as other phenomena related to materials science. The thermodynamic consistency that forms the basis of these…
Tokamak fusion reactors are actively studied as a means of realizing energy production from plasma fusion. However, due to the substantial cost and time required to construct fusion reactors and run physical experiments, numerical…
Following the previous work of Ferretti and Yang on the role of magnetic fields in the theory of conformal turbulence, we show that non-unitary minimal model solutions to 2-dimensional magnetohydrodynamics (MHD) obtained by dimensional…
A new fully discrete linearized $H^1$-conforming Lagrange finite element method is proposed for solving the two-dimensional magneto-hydrodynamics equations based on a magnetic potential formulation. The proposed method yields numerical…
This paper presents a new computer code to solve the general relativistic magnetohydrodynamics (GRMHD) equations using distributed parallel adaptive mesh refinement (AMR). The fluid equations are solved using a finite difference Convex ENO…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with…
The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an…
Equations of ideal magnetohydrodynamics (MHD) play an important role in the studies of turbulence, astrophysics, and plasma physics. These equations possess remarkable geometric structures and symmetries. Indeed, they admit a geodesic…
In this paper we propose a novel thermodynamically compatible finite volume scheme for the numerical solution of the equations of magnetohydrodynamics (MHD) in one and two space dimensions. As shown by Godunov in 1972, the MHD system can be…
We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an…
The JOREK extended magneto-hydrodynamic (MHD) code is a widely used simulation code for studying the non-linear dynamics of large-scale instabilities in divertor tokamak plasmas. Due to the large scale-separation intrinsic to these…
We describe the construction of stepped-pressure equilibria as extrema of a multi-region, relaxed magnetohydrodynamic (MHD) energy functional that combines elements of ideal MHD and Taylor relaxation, and which we call MRXMHD. The model is…
Multigrid solvers are the standard in modern scientific computing simulations. Domain Decomposition Aggregation-Based Algebraic Multigrid, also known as the DD-$\alpha$AMG solver, is a successful realization of an algebraic multigrid solver…
A numerical scheme is developed for systems of conservation laws on manifolds which arise in high speed aerodynamics and magneto-aerodynamics. The systems are presented in an arbitrary coordinate system on the manifold and involve source…
The transport of energy through radiation is very important in many astrophysical phenomena. In dynamical problems the time-dependent equations of radiation hydrodynamics have to be solved. We present a newly developed…