Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetohydrodynamical simulations”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

iharm3D: Vectorized General Relativistic Magnetohydrodynamics

iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM (“High Accuracy Relativistic Magnetohydrodynamics”) algorithm outlined in Gammie et al. (2003) with updates as outlined in McKinney & Gammie (2004) and Noble et al. (2006). The code is most directly derived from Ryan et al. (2015) but with radiative transfer portions removed. HARM is a conservative finite-volume scheme for solving the equations of ideal GRMHD, a hyperbolic system of partial differential equations, on a logically Cartesian mesh in arbitrary coordinates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Enhanced plasma current spike formation due to onset of 1/1 kink-tearing reconnection during a massive gas injection process

The formation of the plasma current spike at the end of the thermal quench (TQ) phase is studied systematically, which is found to strongly correlate with the onset of the 1/1 kink-tearing reconnection only in the simulation results presented here. Here, the magnetohydrodynamic activity on the q = 1 surface plays a critical role in the spike formation and the disruption process. In particular, when the safety factor in the magnetic axis q 0 exceeds 1, the plasma major disruption transits into successive minor disruptions, and the start of the TQ phase is delayed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Whole device modeling of the fuze sheared-flow-stabilized Z pinch

Abstract The FuZE sheared-flow-stabilized Z pinch at Zap Energy is simulated using whole-device modeling employing an axisymmetric resistive magnetohydrodynamic formulation implemented within the discontinuous Galerkin WARPXM framework. Simulations show formation of Z pinches with densities of approximately 10 22 m −3 and total DD fusion neutron rate of 10 7 per µ s for approximately 2 µ s. Simulation-derived synthetic diagnostics show peak currents and voltages within 10% and total yield within approximately 30% of experiment for similar plasma mass. The simulations provide insight into the plasma dynamics in the experiment and enable a predictive capability for exploring design changes on devices built at Zap Energy.

Physics↗
DOE OSTI2022

The structure of 3-D collisional magnetized bow shocks in pulsed-power-driven plasma flows

Here we investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ($M_S$ $\sim 7$), super-Alfvénic ($M_A>1$) magnetized flows with frozen-in magnetic flux $(R_M\gg1)$. These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation of Chemistry in the Athena++ Code

Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magnetohydrodynamic code Athena++ . We implement several chemical networks and heating and cooling processes suitable for simulating the interstellar medium (ISM). A general chemical network framework in the KIDA format is also included, allowing users to easily implement their own chemistry. Radiation transfer and cosmic-ray ionization are coupled with chemistry and solved with the simple six-ray approximation. The chemical and thermal processes are evolved as a system of coupled ordinary differential equations with an implicit solver from the CVODE library. We perform and present a series of tests to ensure the numerical accuracy and convergence of the code. Many tests combine chemistry with gas dynamics, including comparisons with analytic solutions, 1D problems of the photodissociation regions and shocks, and realistic 3D simulations of the turbulent ISM. We release the code with the new public version of Athena++ , aiming to provide a robust and flexible code for the astrochemical simulation community.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Self-organized confinement in whole-device modeling of laboratory magnetospheres

Turbulent self-organization driven by global ideal interchange modes in a dipole-confined plasma is explored with self-consistent, whole-plasma simulations using a flux-tube averaged magnetohydrodynamic model in dipole magnetic geometry. We show the existence of robust particle pinch driven by ideal interchange-mode fluctuations, in which the particles are transported up the density gradient. It is found that the plasma profiles in a dipole field spontaneously relax to a marginally stable state as centrally peaked pressure and density are created by global interchange-mode transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Variety in the Variability of Accreting Supermassive Binary Black Holes

Accreting supermassive binary black holes are key multi-messenger sources for LISA, yet are challenging to simulate realistically as solving the radiation magnetohydrodynamics (MHD) equations over the full dynamic spatio-temporal range of the problem is computationally infeasible at present. We report on our collaboration's progress simulating these systems using general relativistic MHD simulations and dynamic GR. In order to cover a larger temporal range in one set of simulations, we constrain our view to the circumbinary disk region and measure the effects binary mass ratio, accretion disk size, and black hole spin have on the structure and variability of the accretion flow. We particularly emphasize how these parameters influence the overdensity feature, which orbits the binary near the edge of the cavity, since it is responsible for most of the electromagnetic emission's variability and variability is a key signature of a system being a binary. Extending to smaller length scales, we will report on simulations following accretion all the down to the event horizons so that we can begin to investigate how black hole spin affects mini-disk dynamics, accretion rate, and jet power.

Scott C Noble↗
NASA NTRS2008

Recovering Photospheric Velocities from Vector Magnetograms by using a Three-Dimensional, Fully Magnetohydrodynamic (MHD) Model

We introduce a numerical simulation method for recovering the photospheric velocity field from the vector magnetograms. The traditional method is local correlation tracking (LCT) which is based on measuring the relative displacements of features in blocks of pixels between successive white-light images or magnetograms. Within this method, there are a variety of implementations. One of recently developed implementations is induction local correlation tracking (ILCT) as described by Welsch et al. (2004). They employ the normal component of magnetic induction equation as a constraint to assure consistent solutions. Our numerical method uses the fully three-dimensional MHD equations to recover the photospheric velocity field with individual vector magnetograms. We compare our method to the ILCT method using NOAA AR8210 as an example. The differences and similarities are discussed in detail.

Wang, A. H.↗
DOE OSTI2021

Viscous effects on plasmoid formation from nonlinear resistive tearing growth in a Harris sheet

Abstract In this work, the evolution of a highly unstable m = 1 resistive tearing mode, leading to plasmoid formation in a Harris sheet, is studied in the framework of full MHD model using the Non-Ideal Magnetohydrodynamics with Rotation, Open Discussion simulation. Following the initial nonlinear growth of the primary m = 1 island, the X-point develops into a secondary elongated current sheet that eventually breaks into plasmoids. Two distinctive viscous regimes are found for the plasmoid formation and saturation. In the low viscosity regime (i.e. P r ≲ 1), the plasmoid width increases sharply with viscosity, whereas in the viscosity dominant regime (i.e. P r ≳ 1), the plasmoid size gradually decreases with viscosity. Such a finding quantifies the role of viscosity in modulating the plasmoid formation process through its effects on the plasma flow and the reconnection itself.

Physics↗
NASA NTRS1987

Plasma physics abstracts, 1 January - 31 December, 1986

Topics addressed include: ion-cyclotron waves; plasma waves; solar wind lithium releases; bow shock; Pi2 wave bursts; auroral kilometric radiation; ion energization; magnetic field corrections; electric fields; magnetospheric processes; electron acceleration; inner heliosphere; nightside auroral zone; computerized simulation; plasma wave turbulence; and magnetohydrodynamic waves in plasma sheets.

Gurnett, D. A.↗
DOE OSTI2026

Modeling Material Interfaces with the Six-Equation Model in FLEXO

FLEXO is a multiphysics code developed at Sandia National Laboratories for predictive simulation of pulsed power target physics with extended magnetohydrodynamics modeling. Given this application space, FLEXO must be able to perform high-fidelity simulations of pulsed-power systems that include multiple materials. This work presents the augmentation of FLEXO with a six-equation model for multimaterial flow. We design and analyze discretizations of the six-equation model, emphasizing desirable properties for robust multimaterial simulation. We also present solution limiting, bounds preservation, and pressure equilibration/relaxation techniques that enhance the robustness of the overall multimaterial scheme. We present a suite of test problems to document these new capabilities of the FLEXO code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Magnetohydrodynamic (MHD) analyses of various forms of activity and their propagation through helio spheric space

Theoretical and numerical modeling of solar activity and its effects on the solar atmosphere within the context of magnetohydrodynamics were examined. Specifically, the scientific objectives were concerned with the physical mechanisms for the flare energy build-up and subsequent release. In addition, transport of this energy to the corona and solar wind was also investigated. Well-posed, physically self-consistent, numerical simulation models that are based upon magnetohydrodynamics were sought. A systematic investigation of the basic processes that determine the macroscopic dynamic behavior of solar and heliospheric phenomena was conducted. A total of twenty-three articles were accepted and published in major journals. The major achievements are summarized.

Wu, S. T.↗
DOE OSTI2023

Magnetohydrodynamics with physics informed neural operators

Abstract The modeling of multi-scale and multi-physics complex systems typically involves the use of scientific software that can optimally leverage extreme scale computing. Despite major developments in recent years, these simulations continue to be computationally intensive and time consuming. Here we explore the use of AI to accelerate the modeling of complex systems at a fraction of the computational cost of classical methods, and present the first application of physics informed neural operators (NOs) (PINOs) to model 2D incompressible magnetohydrodynamics (MHD) simulations. Our AI models incorporate tensor Fourier NOs as their backbone, which we implemented with the TensorLY package. Our results indicate that PINOs can accurately capture the physics of MHD simulations that describe laminar flows with Reynolds numbers Re ⩽ 250 . We also explore the applicability of our AI surrogates for turbulent flows, and discuss a variety of methodologies that may be incorporated in future work to create AI models that provide a computationally efficient and high fidelity description of MHD simulations for a broad range of Reynolds numbers. The scientific software developed in this project is released with this manuscript.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2021

Variety in the Variability of Accreting Supermassive Binary Black Holes

Accreting supermassive binary black holes are key multi-messenger sources for LISA, yet are challenging to simulate realistically as solving the radiation magnetohydrodynamics (MHD) equations over the full dynamic spatio-temporal range of the problem is computationally infeasible at present. We will provide a brief summary of the the progress made in the field to understand these systems theoretically and what new directions groups are pursuing. We will also report on our collaboration's progress to simulate these systems using general relativistic MHD and dynamic GR. In order to cover a larger temporal range in one set of simulations, we constrain our view to the circumbinary disk region and measure how the binary mass ratio, accretion disk size, and black hole spin have on the structure and variability of the accretion flow. We particularly emphasize how these parameters influence the overdensity feature, which orbits the binary near the edge of the cavity, since it is responsible for most of the electromagnetic emission's variability and variability is a key signature of a system being a binary. Extending to smaller length scales, we will report on simulations following accretion all the down to the event horizons so that we can begin to investigate how black hole spin affects mini-disk dynamics, accretion rate, and jet power. The novel computational methods enabling inclusion of the black holes in the domain, including multi-patch methods, will be described.

Scott C Noble↗
DOE OSTI2025

Hybrid Simulations of FRC Merging and Compression

An improved understanding of field-reversed configuration (FRC) merging and stability in high acceleration and compression magnetic fields is needed to speed up the development of the pulsed fusion concept developed at Helion Energy. All previous theoretical and simulation work on FRC merging and compression was performed using two-dimensional (2D) magnetohydrodynamic (MHD) models. The results of novel 2D hybrid simulations (fluid electrons and full-orbit kinetic ions) of FRC merging and compression are presented. Results of kinetic and MHD simulations, computed using the HYM code, are compared and analyzed. In cases without axial magnetic compression, both the MHD and hybrid simulations show a high sensitivity to the initial parameters (i.e. FRC separation, velocity, normalized separatrix radius, and plasma viscosity), showing that FRCs with large elongation and separatrix radius either do not merge or merge partially, forming a doublet FRC. In conclusion, application of a mirror coil field at the FRC ends with increasing strength is shown to lead to fast and complete merging of the FRCs in MHD and kinetic simulations.

FRC↗
NASA NTRS1988

Helical magnetohydrodynamic turbulence and the coronal heating problem

Numerical simulations are used to investigate the relaxation of an unconfined, helically turbulent, fully three-dimensional magnetofluid, with conditions similar to those which are thought to result in the heating of the solar corona. In these simulations, the system evolves through a succession of force free states. After a relatively quiescent period of Ohmic decay, a phase of accelerated magnetic energy dissipation occurs. Some magnetic energy is transformed into kinetic energy, and the magnitude of entrophy created is a nontrival fraction of the mean square electric current. Concentrated vorticity structures are seen to play almost as important a role as electric current sheets in the heating process. Coincident with this accelerated dissipation process, a reorganization of the magnetic fields occurs, with transfer of magnetic energy to both shorter and longer wavelength modes than are initially present. The ratio of the magnetic field to the electric current density, alpha does not in general tend to assume a constant value in the force free regions during the evolution of the magnetofluid.

Dahlburg, R. B.↗
DOE OSTI2022

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗