Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD solvers”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

The Adaptively-Refined MHD Solver (ARMS)

ARMS is a robust code for solving the single-fluid, fully compressible equations of magnetohydrodynamics. In this presentation I will describe how to access and run ARMS and how to use HelioSpace to analyze the results. The main features of the code are the following: 1. Equations are solved in either Cartesian or Spherical coordinates. For the latter, either a linear or a Log(r) coordinate can be used. 2. The latest flux-corrected-transport schemes developed at NRZ, are used to advance the convective equations. All waves and shocks are accurately captured. 3. A fully adaptive grid can be used. The code incorporates the PARAMESH technology for block-adaptation. Arbitrary adaptation criteria can be specified. 4. The code maintains a divergence-free magnetic field to machine accuracy throughout all refinement and de-refinement. 5. A uniform resistivity and thermal conduction can be included. Solution is either explicit or implicit. 6. Arbitrary boundary and initial conditions can be used in either Cartesian or spherical coordinates. 7. An accompanying visualization package, HelioSpace, is available for viewing and analyzing the native ARMS output, (without need for interpolating onto a structured grid).

Antichos, S. K.↗
DOE OSTI2020

Multidimensional Tests of a Finite-Volume Solver for MHD With a Real-Gas Equation of State

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Gas-Kinetic Theory Based Flux Splitting Method for Ideal Magnetohydrodynamics

A gas-kinetic solver is developed for the ideal magnetohydrodynamics (MHD) equations. The new scheme is based on the direct splitting of the flux function of the MHD equations with the inclusion of "particle" collisions in the transport process. Consequently, the artificial dissipation in the new scheme is much reduced in comparison with the MHD Flux Vector Splitting Scheme. At the same time, the new scheme is compared with the well-developed Roe-type MHD solver. It is concluded that the kinetic MHD scheme is more robust and efficient than the Roe- type method, and the accuracy is competitive. In this paper the general principle of splitting the macroscopic flux function based on the gas-kinetic theory is presented. The flux construction strategy may shed some light on the possible modification of AUSM- and CUSP-type schemes for the compressible Euler equations, as well as to the development of new schemes for a non-strictly hyperbolic system.

Xu, Kun↗
NASA NTRS2013

A Simple GPU-Accelerated Two-Dimensional MUSCL-Hancock Solver for Ideal Magnetohydrodynamics

We describe our experience using NVIDIA's CUDA (Compute Unified Device Architecture) C programming environment to implement a two-dimensional second-order MUSCL-Hancock ideal magnetohydrodynamics (MHD) solver on a GTX 480 Graphics Processing Unit (GPU). Taking a simple approach in which the MHD variables are stored exclusively in the global memory of the GTX 480 and accessed in a cache-friendly manner (without further optimizing memory access by, for example, staging data in the GPU's faster shared memory), we achieved a maximum speed-up of approx. = 126 for a sq 1024 grid relative to the sequential C code running on a single Intel Nehalem (2.8 GHz) core. This speedup is consistent with simple estimates based on the known floating point performance, memory throughput and parallel processing capacity of the GTX 480.

graphics processing units↗
NASA NTRS2012

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2025

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

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 alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Investigation of Anchor Link Magnetohydrodynamic Effects in the Toroidally Symmetric Lead-Lithium (TSLL) Blanket Concept

This paper presents a comprehensive study on the magnetohydrodynamic (MHD) flow in a slotted channel with a cylindrical anchor link, which is an imperative component of the toroidally symmetric lead-lithium (TSLL) blanket concept. Following the validation of the MHD solver implemented in COMSOL and a proposed subtraction approach to compute the anchor link pressure drop, the effects of computational domain size on the pressure drop and velocity distribution are examined. The results show that the pressure drop associated with the anchor link and the maximum flow velocity follow an asymptotic trend as the domain width increases, with wall-induced pressure drop being more dominant than that of the anchor link. The velocity distribution analysis revealed the formation of an internal boundary layer extending along the magnetic field direction, which is a unique feature of the investigated MHD flow. An estimation of the total MHD pressure drop associated with the array of anchor links under the TSLL blanket conditions suggests ~0.3 MPa, which is significantly lower than the recommended maximum allowable blanket pressure drop of 2 MPa. In conclusion, this work offers valuable insights into the anchor link–associated MHD phenomena and serves as a foundation for further development of the TSLL blanket concept for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Long-term 3D MHD simulations of black hole accretion discs formed in neutron star mergers

ABSTRACT We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in flash 4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disc simulations that vary the initial magnetic field geometry and disc compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modelling the BH’s gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively inefficient state in the disc. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disc compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively inefficient state at late times. The electron fraction distribution of the disc outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding ${\sim} 50{{\ \rm per\ cent}}$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

Fahlman, Steven (ORCID:0000000348759940)↗
DOE OSTI2021

Toroidal modeling of runaway avalanche in DIII-D discharges

Abstract A toroidal modeling tool is developed to study the runaway electron (RE) avalanche production process in tokamak plasmas, by coupling the Rosenbluth–Putvinski avalanche model (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355) with an n = 0 magneto-hydrodynamic (MHD) solver. Initial value numerical simulations are carried out for two DIII-D discharges with different plasma shapes (one near circular, and the other with high elongation). It is found that, assuming the same level of about 1% seed current level, the Rosenbluth–Putvinski model somewhat underestimates the RE plateau current for the circular-shaped plasma, as compared with that measured in DIII-D experiments. For an elongated, higher current plasma, simulations find strong runaway current avalanche production despite the lack of measured plateau RE current in experiments. A possible reason for this discrepancy is a lack of additional RE dissipation physics in the present two-dimensional model. Systematic scans of the plasma boundary shape, at fixed pre-disruption plasma current, find that the plasma elongation helps to reduce the RE avalanche production, confirming recent results obtained with an analytic model (Fülöp et al 2020 J. Plasma Phys. 86 474860101). The effect of the plasma triangularity (either positive or negative), on the other hand, has a minor effect. On the physics side, the avalanche process involves two competing mechanisms associated with the electric field. On the one hand, a stronger electric field produces a higher instantaneous avalanche growth rate. On the other hand, a fast growing RE current quickly reduces the fraction of the conduction current together with the electric field, and hence a faster dissipation of the poloidal flux. As a final result of these two dynamic processes, the runaway plateau current is not always the largest with the strongest initial electric field. These results lay the foundation for future self-consistent inclusion of the MHD dynamics and the RE amplification process.

Physics↗
NASA NTRS2023

Interchange Reconnection Dynamics in A Solar Coronal Pseudo-Streamer

Context. The generation of the slow solar wind remains an open problem in heliophysics. One of the current theories among those aimed at explaining the injection of coronal plasma in the interplanetary medium is based on interchange reconnection. It assumes that the exchange of magnetic connectivity between closed and open fields allows the injection of coronal plasma in the interplanetary medium to travel along the newly reconnected open field. However, the exact mechanism underlying this effect is still poorly understood. Aims. Our objective is to study this scenario in a particular magnetic structure of the solar corona: a pseudo-streamer. This topological structure lies at the interface between open and closed magnetic field and is thought to be involved in the generation of the slow solar wind. Methods. We performed innovative 3D magnetohydrodynamic (MHD) simulations of the solar corona with a pseudo-streamer, using the Adaptively Refined MHD Solver (ARMS). By perturbing the quasi-steady ambient state with a simple photospheric, large-scale velocity flow, we were able to generate a complex dynamics of the open-and-closed boundary of the pseudo-streamer. We studied the evolution of the connectivity of numerous field lines to understand its precise dynamics. Results. We witnessed different scenarios of opening of the magnetic field initially closed under the pseudo-streamer: one-step interchange reconnection dynamics, along with more complex scenarios, including a coupling between pseudo-streamer and helmet streamer, as well as back-and-forth reconnections between open and closed connectivity domains. Finally, our analysis revealed large-scale motions of a newly opened magnetic field high in the corona that may be explained by slipping reconnection. Conclusions. By introducing a new analysis method for the magnetic connectivity evolution based on distinct closed-field domains, this study provides an understanding of the precise dynamics underway during the opening of a closed field, which enables the injection of closed-field, coronal plasma in the interplanetary medium. Further studies shall provide synthetic observations for these diverse outgoing flows, which could be measured by Parker Solar Probe and Solar Orbiter.

magnetic fields↗
NASA NTRS2012

The Effect of Magnetic Topology on the Escape of Flare Particles

Magnetic reconnection in the solar atmosphere is believed to be the driver of most solar explosive phenomena. Therefore, the topology of the coronal magnetic field is central to understanding the solar drivers of space weather. Of particular importance to space weather are the impulsive Solar Energetic particles that are associated with some CME/eruptive flare events. Observationally, the magnetic configuration of active regions where solar eruptions originate appears to agree with the standard eruptive flare model. According to this model, particles accelerated at the flare reconnection site should remain trapped in the corona and the ejected plasmoid. However, flare-accelerated particles frequently reach the Earth long before the CME does. We present a model that may account for the injection of energetic particles onto open magnetic flux tubes connecting to the Earth. Our model is based on the well-known 2.5D breakout topology, which has a coronal null point (null line) and a four-flux system. A key new addition, however, is that we include an isothermal solar wind with open-flux regions. Depending on the location of the open flux with respect to the null point, we find that the flare reconnection can consist of two distinct phases. At first, the flare reconnection involves only closed field, but if the eruption occurs close to the open field, we find a second phase involving interchange reconnection between open and closed. We argue that this second reconnection episode is responsible for the injection of flare-accelerated particles into the interplanetary medium. We will report on our recent work toward understanding how flare particles escape to the heliosphere. This work uses high-resolution 2.5D MHD numerical simulations performed with the Adaptively Refined MHD Solver (ARMS).

Antiochos, S. K.↗
NASA NTRS2012

A Model for SEP Escape

Magnetic reconnection in the solar atmosphere is believed to be the driver of most solar explosive phenomena. Therefore, the topology of the coronal magnetic field is central to understanding the solar drivers of space weather. Of particular importance to space weather are the impulsive Solar Energetic particles that are associated with some CME/eruptive flare events. Observationally, the magnetic configuration of active regions where solar eruptions originate appears to agree with the standard eruptive flare model. According to this model, particles accelerated at the flare reconnection site should remain trapped in the corona and the ejected plasmoid. However, flare-accelerated particles frequently reach the Earth long before the CME does. We present a model that may account for the injection of energetic particles onto open magnetic flux tubes connecting to the Earth. Our model is based on the well-known 2.5D breakout topology, which has a coronal null point (null line) and a four-flux system. A key new addition, however, is that we include an isothermal solar wind with open-flux regions. Depending on the location of the open flux with respect to the null point, we find that the flare reconnection can consist of two distinct phases. At first, the flare reconnection involves only closed field, but if the eruption occurs close to the open field, we find a second phase involving interchange reconnection between open and closed. We argue that this second reconnection episode is responsible for the injection of flare-accelerated particles into the interplanetary medium. We will report on our recent work toward understanding how flare particles escape to the heliosphere. This work uses high-resolution 2.5D MHD numerical simulations performed with the Adaptively Refined MHD Solver (ARMS).

Antiochos, S. K.↗
NASA NTRS2013

Computational Studies of Magnetic Nozzle Performance

An extensive literature review of magnetic nozzle research has been performed, examining previous work, as well as a review of fundamental principles. This has allow us to catalog all basic physical mechanisms which we believe underlie the thrust generation process. Energy conversion mechanisms include the approximate conservation of the magnetic moment adiabatic invariant, generalized hall and thermoelectric acceleration, swirl acceleration, thermal energy transformation into directed kinetic energy, and Joule heating. Momentum transfer results from the interaction of the applied magnetic field with currents induced in the plasma plume., while plasma detachment mechanisms include resistive diffusion, recombination and charge exchange collisions, magnetic reconnection, loss of adiabaticity, inertial forces, current closure, and self-field detachment. We have performed a preliminary study of Hall effects on magnetic nozzle jets with weak guiding magnetic fields and weak expansions (p(sub jet) approx. = P(sub background)). The conclusion from this study is that the Hall effect creates an azimuthal rotation of the plasma jet and, more generally, creates helical structures in the induced current, velocity field, and magnetic fields. We have studied plasma jet expansion to near vacuum without a guiding magnetic field, and are presently including a guiding magnetic field using a resistive MHD solver. This research is progressing toward the implementation of a full generalized Ohm's law solver. In our paper, we will summarize the basic principle, as well as the literature survey and briefly review our previous results. Our most recent results at the time of submittal will also be included. Efforts are currently underway to construct an experiment at the University of Michigan Plasmadynamics and Electric Propulsion Laboratory (PEPL) to study magnetic nozzle physics for a RF-thruster. Our computational study will work directly with this experiment to validate the numerical model, in order to study magnetic nozzle physics and optimize magnetic nozzle design. Preliminary results from the PEPL experiment will also be presented.

Ebersohn, Frans H.↗
DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension)

An approximate Riemann solver is developed for the governing equations of ideal magnetohydrodynamics (MHD). The Riemann solver has an eight-wave structure, where seven of the waves are those used in previous work on upwind schemes for MHD, and the eighth wave is related to the divergence of the magnetic field. The structure of the eighth wave is not immediately obvious from the governing equations as they are usually written, but arises from a modification of the equations that is presented in this paper. The addition of the eighth wave allows multidimensional MHD problems to be solved without the use of staggered grids or a projection scheme, one or the other of which was necessary in previous work on upwind schemes for MHD. A test problem made up of a shock tube with rotated initial conditions is solved to show that the two-dimensional code yields answers consistent with the one-dimensional methods developed previously.

Powell, Kenneth G.↗
NASA NTRS2002

The Delta x B = 0 Constraint Versus Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. It has previously been shown by Powell et al. (1995) that an 'almost' equivalent MHD system in non-conservative form can be derived. This non-conservative system has a better conditioned eigensystem. Aside from Powell et al., the MHD equations can be derived from basic principles in either conservative or non-conservative form. The Delta x B = 0 constraint of the MHD equations is only an initial condition constraint, it is very different from the incompressible Navier-Stokes equations in which the divergence condition is needed to close the system (i.e., to have the same number of equations and the same number of unknown). In the MHD formulations, if Delta x B = 0 initially, all one needs is to construct appropriate numerical schemes that preserve this constraint at later time evolutions. In other words, one does not need the Delta x B condition to close the MHD system. We formulate our new scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the divergence of the magnetic condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗