Search NASA⌕ Search

SEARCH · Search NASA

Results for “implicit magnetohydrodynamics”

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
DOE OSTI2010

Toward textbook multigrid efficiency for fully implicit resistive magnetohydrodynamics

Multigrid methods can solve some classes of elliptic and parabolic equations to accuracy below the truncation error with a work-cost equivalent to a few residual calculations - so-called 'textbook' multigrid efficiency. We investigate methods to solve the system of equations that arise in time dependent magnetohydrodynamics (MHD) simulations with textbook multigrid efficiency. We apply multigrid techniques such as geometric interpolation, full approximate storage, Gauss-Seidel smoothers, and defect correction for fully implicit, nonlinear, second-order finite volume discretizations of MHD. We apply these methods to a standard resistive MHD benchmark problem, the GEM reconnection problem, and add a strong magnetic guide field, which is a critical characteristic of magnetically confined fusion plasmas. We show that our multigrid methods can achieve near textbook efficiency on fully implicit resistive MHD simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Toward textbook multigrid efficiency for fully implicit resistive magnetohydrodynamics

Multigrid methods can solve some classes of elliptic and parabolic equations to accuracy below the truncation error with a work-cost equivalent to a few residual calculations – so-called “textbook” multigrid efficiency. We investigate methods to solve the system of equations that arise in time dependent magnetohydrodynamics (MHD) simulations with textbook multigrid efficiency. We apply multigrid techniques such as geometric interpolation, full approximate storage, Gauss-Seidel smoothers, and defect correction for fully implicit, nonlinear, second-order finite volume discretizations of MHD. We apply these methods to a standard resistive MHD benchmark problem, the GEM reconnection problem, and add a strong magnetic guide field, which is a critical characteristic of magnetically confined fusion plasmas. We show that our multigrid methods can achieve near textbook efficiency on fully implicit resistive MHD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Implicit Methods for the Magnetohydrodynamic Description of Magnetically Confined Plasmas

Implicit algorithms are essential for predicting the slow growth and saturation of global instabilities in today’s magnetically confined fusion plasma experiments. Present day algorithms for obtaining implicit solutions to the magnetohydrodynamic (MHD) equations for highly magnetized plasma have their roots in algorithms used in the 1960s and 1970s. However, today’s computers and modern linear and non-linear solver techniques make practical much more comprehensive implicit algorithms than were previously possible. Combining these advanced implicit algorithms with highly accurate spatial representations of the vector fields describing the plasma flow and magnetic fields and with improved methods of calculating anisotropic thermal conduction now makes possible simulations of fusion experiments using realistic values of plasma parameters and actual configuration geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

A full-implicit-continuous-Eulerian (FICE) scheme for multidimensional transient magnetohydrodynamic (MHD) flows

A full implicit continuous Eulerian (FICE) scheme is developed for solving multidimensional transient MHD flow problems. The physical system under consideration is a general case of a transient MHD flow in which an initial steady state is subject to a finite amplitude disturbance. The governing equations are described, their finite difference formulation is presented and the FICE algorithm is given. The boundary conditions are treated by classifying them into physical and computational ones. The usefulness of the FICE algorithm is demonstrated using a physical example concerning the dynamical response of the static solar atmosphere due to a representative photospheric disturbance.

Wu, S. T.↗
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↗
DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

FLIP MHD - A particle-in-cell method for magnetohydrodynamics

The fluid-implicit-particle, or 'FLIP' method presently extended to 2D and 3D MHD flow incorporates a Lagrangian field representation and yields a grid magnetic Reynolds number of up to 16 while preserving contact continuities that retain the Galilean invariance of the MHD flow equations. Analytical arguments and numerical examples demonstrate the conservation of mass, momentum, magnetic flux, and energy; 2D calculation results for the illustrative cases of contact discontinuity convection, Rayleigh-Taylor unstable flow.

Brackbill, J. U.↗
DOE OSTI2022

General Relativistic Implicit Monte Carlo Radiation-hydrodynamics

Abstract We report on a new capability added to our general relativistic radiation-magnetohydrodynamics code, Cosmos++ : an implicit Monte Carlo (IMC) treatment for radiation transport. The method is based on a Fleck-type implicit discretization of the radiation-hydrodynamics equations, but generalized for both Newtonian and relativistic regimes. A multiple reference frame approach is used to geodesically transport photon packets (and solve the hydrodynamics equations) in the coordinate frame, while radiation–matter interactions are handled either in the fluid or electron frames then communicated via Lorentz boosts and orthonormal tetrad bases attached to the fluid. We describe a method for constructing estimators of radiation moments using path-weighting that generalizes to arbitrary coordinate systems in flat or curved spacetime. Absorption, emission, scattering, and relativistic Comptonization are among the matter interactions considered in this report. We discuss our formulations and numerical methods, and validate our models against a suite of radiation and coupled radiation-hydrodynamics test problems in both flat and curved spacetimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

ON THE MODE OF DYNAMO ACTION IN A GLOBAL LARGE-EDDY SIMULATION OF SOLAR CONVECTION

In this paper, we examine the mode of dynamo action in the implicit large-eddy magnetohydrodynamical simulation of solar convection reported upon in Ghizaru et al. Motivated by the presence of a strong and well-defined large-scale axisymmetric magnetic component undergoing regular polarity reversals, we define the fluctuating component of the magnetic field as the difference between the total field and its zonal average. The subsequent analysis follows the physical logic and mathematical formulation of mean-field electrodynamics, whereby a turbulent electromotive force (EMF) is computed by the suitable averaging of cross-correlations between fluctuating flow and field components and expressed in terms of the mean field via a linear truncated tensorial expansion. We use singular value decomposition to perform a linear least-squares fit of the temporal variation of the EMF to that of the large-scale magnetic component, which yields the components of the full {alpha}-tensor. Its antisymmetric component, describing general turbulent pumping, is also extracted. The {alpha}-tensor so calculated reproduces a number of features already identified in local, Cartesian simulations of magnetohydrodynamical rotating convection, including an {alpha}{sub {phi}{phi}} component positive in the northern solar hemisphere, peaking at high latitudes, and reversing sign near the bottom of the convection zone; downward turbulent pumping throughout the convecting layer; and significant equatorward turbulent pumping at mid latitudes, and poleward at high latitudes in subsurface layers. We also find that the EMF contributes significantly to the regeneration of the large-scale toroidal magnetic component, which from the point of view of mean-field dynamo models would imply that the simulation operates as an {alpha}{sup 2}{Omega} dynamo. We find little significant evidence of {alpha}-quenching by the large-scale magnetic field. The amplitude of the magnetic cycle appears instead to be regulated primarily by a magnetically driven reduction of the differential rotation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Final Report: Towards Optimal Petascale Simulations (TOPS), ER25785

Multiscale, multirate scientific and engineering applications in the SciDAC portfolio possess resolution requirements that are practically inexhaustible and demand execution on the highest-capability computers available, which will soon reach the petascale. While the variety of applications is enormous, their needs for mathematical software infrastructure are surprisingly coincident; moreover the chief bottleneck is often the solver. At their current scalability limits, many applications spend a vast majority of their operations in solvers, due to solver algorithmic complexity that is superlinear in the problem size, whereas other phases scale linearly. Furthermore, the solver may be the phase of the simulation with the poorest parallel scalability, due to intrinsic global dependencies. This project brings together the providers of some of the world's most widely distributed, freely available, scalable solver software and focuses them on relieving this bottleneck for many specific applications within SciDAC, which are representative of many others outside. Solver software directly supported under TOPS includes: hypre, PETSc, SUNDIALS, SuperLU, TAO, and Trilinos. Transparent access is also provided to other solver software through the TOPS interface. The primary goals of TOPS are the development, testing, and dissemination of solver software, especially for systems governed by PDEs. Upon discretization, these systems possess mathematical structure that must be exploited for optimal scalability; therefore, application-targeted algorithmic research is included. TOPS software development includes attention to high performance as well as interoperability among the solver components. Support for integration of TOPS solvers into SciDAC applications is also directly supported by this proposal. The role of the UCSD PI in this overall CET, is one of direct interaction between the TOPS software partners and various DOE applications scientists' specifically toward magnetohydrodynamics (MHD) simulations with the Center for Extended Magnetohydrodynamic Modeling (CEMM) SciDAC and Applied Partial Differential Equations Center (APDEC) SciDAC, and toward core-collapse supernova simulations with the previous Terascale Supernova Initiative (TSI) SciDAC and in continued work on INCITE projects headed by Doug Swesty, SUNY Stony Brook. In addition to these DOE applications scientists, the UCSD PI works to bring leading-edge DOE solver technology to applications scientists in cosmology and large-scale galactic structure formation. Unfortunately, the funding for this grant ended after only two years of its five-year duration, in August 2008, due to difficulties at DOE in transferring the grant to the PI's new faculty position at Southern Methodist University. Therefore, this report only describes two years' worth of effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

EXPLICIT-IMPLICIT SCHEME FOR RELATIVISTIC RADIATION HYDRODYNAMICS

We propose an explicit-implicit scheme for numerically solving special relativistic radiation hydrodynamic equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, zeroth and first moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion approximation. For an hyperbolic term, of which the timescale is the light crossing time when the flow velocity is comparable to the speed of light, is explicitly solved using an approximate Riemann solver. Source terms describing an exchange of energy and momentum between the matter and the radiation via the gas-radiation interaction are implicitly integrated using an iteration method. The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in optically thick media, where heating/cooling and scattering timescales could be much shorter than the dynamical timescale. We show that our numerical code can pass test problems of one- and two-dimensional radiation energy transport, and one-dimensional radiation hydrodynamics. Our newly developed scheme could be useful for a number of relativistic astrophysical problems. We also discuss how to extend our explicit-implicit scheme to the relativistic radiation magnetohydrodynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

General-relativistic Resistive Magnetohydrodynamics with Robust Primitive-variable Recovery for Accretion Disk Simulations

Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black holeat the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting “hot spot”nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the developmentof novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related tosuch hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin currentlayers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magneticreconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) codeBHACand present the implementationof an Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithmis tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transportmethod to keep the magnetic field solenoidal. Several novel methods for primitive variable recovery, a key part inrelativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency.We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magneticpressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron starmagnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accountingfor dissipative effects as reconnection

(magnetohydrodynamics:) MHD↗
DOE OSTI2016

On the performance of exponential integrators for problems in magnetohydrodynamics

Exponential integrators have been introduced as an efficient alternative to explicit and implicit methods for integrating large stiff systems of differential equations. Over the past decades these methods have been studied theoretically and their performance was evaluated using a range of test problems. While the results of these investigations showed that exponential integrators can provide significant computational savings, the research on validating this hypothesis for large scale systems and understanding what classes of problems can particularly benefit from the use of the new techniques is in its initial stages. Resistive magnetohydrodynamic (MHD) modeling is widely used in studying large scale behavior of laboratory and astrophysical plasmas. In many problems numerical solution of MHD equations is a challenging task due to the temporal stiffness of this system in the parameter regimes of interest. In this paper we evaluate the performance of exponential integrators on large MHD problems and compare them to a state-of-the-art implicit time integrator. Both the variable and constant time step exponential methods of EPIRK-type are used to simulate magnetic reconnection and the Kevin–Helmholtz instability in plasma. Performance of these methods, which are part of the EPIC software package, is compared to the variable time step variable order BDF scheme included in the CVODE (part of SUNDIALS) library. We study performance of the methods on parallel architectures and with respect to magnitudes of important parameters such as Reynolds, Lundquist, and Prandtl numbers. We find that the exponential integrators provide superior or equal performance in most circumstances and conclude that further development of exponential methods for MHD problems is warranted and can lead to significant computational advantages for large scale stiff systems of differential equations such as MHD.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

On performance of Krylov smoothing for fully-coupled AMG preconditioners for VMS resistive MHD

Here, this study explores the performance and scaling of a GMRES Krylov method employed as a smoother for an algebraic multigrid (AMG) preconditioned Newton- Krylov solution approach applied to a fully-implicit variational multiscale (VMS) nite element (FE) resistive magnetohydrodynamics (MHD) formulation. In this context a Newton iteration is used for the nonlinear system and a Krylov (GMRES) method is employed for the linear subsystems. The efficiency of this approach is critically dependent on the scalability and performance of the AMG preconditioner for the linear solutions and the performance of the smoothers play a critical role. Krylov smoothers are considered in an attempt to reduce the time and memory requirements of existing robust smoothers based on additive Schwarz domain decomposition (DD) with incomplete LU factorization solves on each subdomain. Three time dependent resistive MHD test cases are considered to evaluate the method. The results demonstrate that the GMRES smoother can be faster due to a decrease in the preconditioner setup time and a reduction in outer GMRESR solver iterations, and requires less memory (typically 35% less memory for global GMRES smoother) than the DD ILU smoother.

97 MATHEMATICS AND COMPUTING↗