Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD simulation”

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 1,027 records · Page 57

Yield and PDV Measurements of the Godiva-IV Critical Assembly

A high dynamic range, scinillator based detector, the MHD-240, was used to measure the gamma and neutron flux of several prompt critical bursts of the Godiva-IV assembly. The measurements provided the alpha and FWHM for several different burst temperatures. Calibrations of the MHD240 gamma and neutron sensitivity along with MCNP based simulation of the Godiva-IV assembly allowed us to also determine the total fission yield of the burst integrated over approximately 5 ms. In addition a Photo Doppler Velocimetry measurement was conducted to determine the surface motion vibrational frequency of Godiva-IV assembly during the bursts. The combined data were collected to be used as input to multi-physics simulation of the Godiva-IV assembly. This is a preliminary report, there is an ongoing effort to update the MCNP simulations used in the analysis and develop a correction for the room return component of the measured signal.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The electrostatic effect for the collisionless tearing mode

Electron dynamics has not been self-consistently considered in collisionless tearing mode theories to date because of the mathematical complexity of the Vlasov-Maxwell equations. It has been found, using computer simulations, that electrostatic fields play an important role in the tearing mode. Vlasov theory, including the electrostatic field, is investigated for topologies with both antiparallel and nonantiparallel magnetic field lines. The electrostatic field influences the resonant current in the neutral sheet which is a non-MHD effect, and modifies the linear growth rate. At the magnetopause, where the field lines are not antiparallel, the electrostatic effect acts to raise the linear growth rate of the tearing mode. On the other hand, in the magnetotail, where magnetic field lines are antiparallel, the electrostatic effect reduces the tearing mode growth rate.

Hoshino, Masahiro↗

Coronal heating by the resonant absorption of Alfven waves: The effect of viscous stress tensor

The time-dependent linearized magnetohydrodynamics (MHD) equations for a fully compressible, low-beta, viscoresistive plasma are solved numerically using an implicit integration scheme. The full viscosity stress tensor (Braginskii 1965) is included with the five parameters eta(sub i) i = 0 to 4. In agreement with previous studies, the numerical simulations demonstrate that the dissipation on inhomogeneities in the background Alfven speed occurs in a narrow resonant layer. For an active region in the solar corona the values of eta(sub i) are eta(sub o) = 0.65 g/cm/s, eta(sub 1) = 3.7 x 10(exp -12) g/cm/s, eta(sub 2) = 4 eta(sub 1), eta(sub 3) = 1.4 x 10(exp -6) g/cm/s, eta(sub 4) = 2 eta(sub 3), with n = 10(exp 10)/cu cm, T = 2 x 10(exp 6) K, and B = 100 G. When the Lundquist number S = 10(exp 4) and R(sub 1) much greater than S (where R(sub 1) is the dimensionless shear viscous number) the width of the resistive dissipation layer d(sub r) is 0.22a (where a is the density gradient length scale) and d(sub r) approximately S(exp -1/3). When S much greater than R(sub 1) the shear viscous dissipation layer width d(sub r) scales as R(sub 1)(exp -1/3). The shear viscous and the resistive dissipation occurs in an overlapping narrow region, and the total heating rate is independent of the value of the dissipation parameters in agreement with previous studies. Consequently, the maximum values of the perpendicular velocity and perpendicular magnetic field scale as R(sub 1)(exp -1/3). It is evident from the simulations that for solar parameters the heating due to the compressive viscosity (R(sub 0) = 560) is negligible compared to the resistive and the shear viscous (R(sub 1)) dissipation and it occurs in a broad layer of order a in width. In the solar corona with S approximately equals 10(exp 4) and R(sub 1) approximately equals 10(exp 14) (as calculated from the Braginskii expressions), the shear viscous resonant heating is of comparable magnitude to the resistive resonant heating.

Ofman, L.↗

Nested Pebble Bed Blanket (NesPeB)

Recent advances in magnetic confinement fusion technology have attracted billions of dollars of investments in startups from venture capitals and corporations, resulting in the development of devices aiming to demonstrate net energy gain in a self-heated burning plasma, such as SPARC (under construction) and others. However, future fusion power plants must operate in regimes that will require technologies far beyond current experience. According to a National Academies of Science, Engineering, and Medicine report, to have nuclear fusion power plants contributing in a timely manner to the planned reduction of atmospheric carbon dioxide, a pilot plant should be built by 2035, and it should demonstrate fusion power production and the performance of the tritium fuel system (requiring a high enough tritium breeding) by 2040. A recognized key technology gap by [26] is the fusion first wall and blanket since no current blanket concept is considered satisfactory or has been built and proven. The first wall and blanket in magnetic fusion reactors form a vital and complex system, as it must satisfy different functions such as power extraction, tritium breeding, plasma containment, radiation shielding, and safety. The list of design requirements is even longer: high enough tritium production for fusion self-sufficiency, low material activation, decay heat and shutdown dose rates, high thermal efficiency, high-capacity factor, high magnets-divertor-vacuum vessel-first wall life, low corrosion, low cost, and intrinsically safe (requiring minimal licensing). Despite fifty-plus years of research, the first wall and blanket concepts proposed suffer from fundamental technical problems and immaturity (TRL=2-3) that jeopardize the timely delivery of a commercial fusion power plant. A fusion first-wall blanket has never been built nor tested, and a "winning", practical functioning design requires enough engineering margins (high enough tritium breeding considering the uncertainty, etc.), manufacturing simplicity, ease of continuous operation, maintenance, and low cost. A new, groundbreaking blanket concept called "Nested Pebble Bed Blanket" (NesPeB) was developed at ORNL under the successful ARPA-E GAMOW FERMI project (patent application allowed by the USPTO). The NesPeB blanket concept addresses current blanket concepts' shortcomings and technical immaturity, paving the way for accelerated delivery of fusion power plants. NesPeB is based on nested pebbles, which are binary-sized lithium-ceramic pebbles enclosed in "Beryllide" perforated and coated spherical shells, which are also binary-sized, stacked on top of each other, forming a "bed" and cooled by Nitrogen gas also "sweeping" the Helium and Tritium generated by the neutron irradiation of Lithium; the vacuum vessel plasma facing material is Molybdenum-96 and -97 with the first wall cooled by Helium while the divertor armor is made of Tungsten. The simulations of the NesPeB blanket using Fusion Reactors Models Integrator (FERMI) are encouraging as they estimate a tritium breeding ratio (TBR) greater than 1.2 using natural Lithium, acceptable pressure drop, and excellent heat transfer properties. Furthermore, the NesPeB blanket is not limited by magneto-hydro-dynamics (MHD) effects, is designed for online refueling, relies on existing tritium extraction technologies, has a simple construction, and limits the corrosion and chemical reactivity problems. NesPeB has the potential to be transformational and disruptive since it can solve all the main, challenging technical problems of fusion device blankets and accelerate a pilot plant delivery for 10 or more years.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Substorm Cycle as Reproduced by Global MHD Models

Recently, Gordeev et al. (2015) suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field, and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized2 h north then 2 h south interplanetary magnetic field (IMF) Bz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly available simulations which is related to post processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

Gordeev, E.↗

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↗

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stability and Control of Burning Tokamak Plasmas with Resistive Walls (Final Report)

This research has focused on quantitative prediction of the stability, control, and equilibrium state evolution in toroidal burning plasmas. The stability of long pulse burning plasmas is highly sensitive to the physics of resonant layers in the plasma, sources of momentum and flow, kinetic effects of energetic particles, and boundary conditions at the wall, including feedback control and error fields. In ITER in particular, the low toroidal flow equilibrium state, sustained primarily by energetic alpha particles from fusion reactions, will require the consideration of all of these key elements to predict quantitatively the stability and evolution. The principal investigators on this proposal are leading experts in the relevant theoretical and computational areas, and aimed to perform computations guided by analytic modeling, to address this physics in realistic configurations. The overall goal is to understand the key physics mechanisms that describe resistive toroidal burning plasmas, surrounded by a resistive wall, under active feedback control. With the physics of the energetic ions, resonant layers, resistive wall, and toroidal momentum transport included, this study will extend from recent publications in theory and simulation of individual effects and move toward predictive modeling for burning plasmas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Theoretical analysis of energetic-ion-driven resistive interchange mode stabilization strategies using a Landau closure model

The aim of the present study is to perform a theoretical analysis of different strategies to stabilize energetic-ion-driven resistive interchange mode (EIC) in Large Helical Device (LHD) plasma. In this work, we use a reduced MHD for the thermal plasma coupled with a gyrofluid model for the energetic particles (EP) species. The hellically trapped EP component is introduced through a modification of the drift frequency to include their precessional drift. The stabilization trends of the 1/1 EIC observed experimentally with respect to the thermal plasma density and temperature are reproduced by the simulations, showing a reasonable agreement with the data. The LHD operation scenarios with stable 1/1 EIC are identified, leading to the stabilization of the 1/1 EIC if the thermal plasma density and temperature are above a given threshold. Additionally, the 1/1 EIC are also stabilized if the rotational transform is modified in a way that the 1/1 rational surface is located further away than 0.9 times the normalized radius, or the magnetic shear in the plasma periphery is enhanced. Also, LHD discharges with large magnetic fields show a higher EIC destabilization threshold with respect to the thermal plasma density. If the perpendicular NBI deposition region is moved further inward than 0.875 times the normalized radius the 1/1 EIC are also stabilized. Furthermore, increasing the perpendicular NBI voltage such that the EP energy is higher than 30 keV stabilizes the 1/1 EIC. Moreover, deuterium plasmas show a higher stability threshold for the 1/1 EIC than hydrogen plasmas. The experimental data shows a larger time interval between EIC events as the power of the tangential NBI is increased providing that the perpendicular NBI power is at least 13 MW. This implies a stabilizing effect of the tangential NBI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗

Analysis of fusion alphas interaction with RF waves in D-T plasma at JET

This work studies the influence of RF waves in ICRH range of frequency on fusion alphas during the recent JET D-T campaign. Fusion alphas from D-T reactions are born with energies of about 3.5MeV and therefore have significant Doppler shift enabling synergistic interaction between them and RF waves at broad range of frequencies including the ones foreseen for future fusion machines ITER and SPARC. Resonant interaction between RF waves and alphas, also called synergistic effects, will modify the alpha distribution and ultimately will have an impact on alpha orbit losses and heating. Data from JET 3.43T/2.3MA pulses based on the hybrid scenario during the DTE2 campaign were used for the analysis in this study. The impact of synergistic effects on alpha orbit losses and alpha heating is assessed. Conclusions are based on analysis of experimental data for fast alphas losses, i.e. measurements from neutral particle analyser, fast ion losses scintillator detector, Faraday cups, and TRANSP simulations. Experimental data and TRANSP analysis indicate that there are indeed changes in the alphas' distribution function due to interaction with RF waves. Data from the scintillator detector and the Faraday cups were compared for pulses with and without ICRH power and versus cases with enhanced alpha losses due to MHD activities. The trends from these diagnostics consistently show no additional alpha losses due to interaction with RF waves. TRANSP predictions for the impact of the synergistic effects on alpha heating show up to 42% increase in alpha electron heating and up to 25% increase in alpha ion heating. These effects however become negligibly small, less than 1%, when alpha heating is compared to the total auxiliary hearting power in the investigated JET pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

MHD Turbulence Sheared in Fixed and Rotating Frames

We consider homogeneous turbulence in a conducting fluid that is exposed to a uniform external magnetic field while being sheared in fixed and rotating frames. We take both the frame-rotation axis and the applied magnetic field to be aligned in the direction normal to the plane of the mean shear. Here a systematic parametric study is carried out in a series of Direct Numerical Simulations (DNS) in order to clarify the main effects determining the structural anisotropy and stability of the flow.

Kassinos, S. C.↗

Three-dimensional, time-dependent, MHD model of a solar flare-generated interplanetary shock wave

A three-dimensional time-dependent MHD model of the propagation of an interplanetary shock wave into an ambient three-dimensional heliospheric solar wind is initialized with a peak velocity of 1000 km/s at the center of a right circular cone of 18 deg included angle at 18 solar radii. Differences from a previous 2-1/2 simulation (Wu et al., 1983; Gislason et al., 1984; Dryer et al., 1984) include diminuation of the solar peak velocity and concentration of the peak density at each radius. The IMF magnitude starts with high-latitude peaks, and helical-like IMF rotation is noted due to a large-amplitude nonlinear Alfven wave in the shocked plasma.

Dryer, M.↗

The role of magnetohydrodynamics in heliospheric space plasma physics research

Magnetohydrodynamics (MHD) is a fairly recent extension of the field of fluid mechanics. While much remains to be done, it has successfully been applied to the contemporary field of heliospheric space plasma research to evaluate the 'macroscopic picture' of some vital topics via the use of conducting fluid equations and numerical modeling and simulations. Some representative examples from solar and interplanetary physics are described to demonstrate that the continuum approach to global problems (while keeping in mind the assumptions and limitations therein) can be very successful in providing insight and large scale interpretations of otherwise intractable problems in space physics.

Dryer, Murray↗

A Two-Fluid, MHD Coronal Model

We describe first results from a numerical two-fluid MHD model of the global structure of the solar corona. The model is two-fluid in the sense that it accounts for the collisional energy exchange between protons and electrons. As in our single-fluid model, volumetric heat and momentum sources are required to produce high speed wind from coronal holes, low speed wind above streamers, and mass fluxes similar to the empirical solar wind. By specifying different proton and electron heating functions we obtain a high proton temperature in the coronal hole and a relatively low proton temperature in the streamer (in comparison with the electron temperature). This is consistent with inferences from SOHO/UVCS, and with the Ulysses/SWOOPS proton and electron temperature measurements which we show from the fast latitude scan. The density in the coronal hole between 2 solar radii and 5 solar radii (2RS and 5RS) is similar to the density reported from SPARTAN 201-01 measurements by Fisher and Guhathakurta. The proton mass flux scaled to 1 AU is 2.4 x 10(exp 8)/sq cm s, which is consistent with Ulysses observations. Inside the closed field region, the density is sufficiently high so that the simulation gives equal proton and electron temperatures due to the high collision rate. In open field regions (in the coronal hole and above the streamer) the proton and electron temperatures differ by varying amounts. In the streamer, the temperature and density are similar to those reported empirically by Li et al and the plasma beta is larger than unity everywhere above approx. 1.5 R(sub s), as it is in all other MHD coronal streamer models.

Suess, Steven T.↗

Thermomagnetic instability of plasma composition gradients

In this paper, we show that, under Braginskii magneto-hydrodynamics, anti-parallel gradients in an average ion charge state and electron temperature can be unstable to the growth of self-generated magnetic fields. The instability is analogous to the field-generating thermomagnetic instability, although it is driven by the collisional thermal force magnetic source term rather than the Biermann battery term. The gradient in ion charge state causes a gradient in collisionality, which couples with temperature perturbations to create a self-generated magnetic field. This magnetic field deflects the electron heat flux in a way that reinforces the temperature perturbation. The derived linearized growth rate, typically on hydrodynamic timescales, includes resistive and thermal smoothing. It increases with large ion composition gradients and electron heat flux, conditions typical of the hohlraum walls or contaminant mix jets in inertial confinement fusion implosions. However, extended magneto-hydrodynamic simulations indicate that the instability is usually dominated and stabilized by nonlinear Nernst advection, in a similar manner to the standard thermomagnetic instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗