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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nature and Scalings of Density Fluctuations of Compressible Magnetohydrodynamic Turbulence with Applications to the Solar Wind

Abstract The solar wind is a magnetized and turbulent plasma. Its turbulence is often dominated by Alfvénic fluctuations and often deemed as nearly incompressible far away from the Sun, as shown by in situ measurements near 1 au. However, for solar wind closer to the Sun, the plasma β decreases (often lower than unity) while the turbulent Mach number M t increases (can approach unity, e.g., transonic fluctuations). These conditions could produce significantly more compressible effects, characterized by enhanced density fluctuations, as seen by several space missions. In this paper, a series of 3D MHD simulations of turbulence are carried out to understand the properties of compressible turbulence, particularly the generation of density fluctuations. We find that, over a broad range of parameter space in plasma β , cross helicity, and polytropic index, the turbulent density fluctuations scale linearly as a function of M t , with the scaling coefficients showing weak dependence on parameters. Furthermore, through detailed spatiotemporal analysis, we show that the density fluctuations are dominated by low-frequency nonlinear structures, rather than compressible MHD eigenwaves. These results could be important for understanding how compressible turbulence contributes to solar wind heating near the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Simulations of coronal disconnection events

The lack of evidence for magnetic disconnection of coronal mass ejections (CMEs) from the sun has long been a puzzle, as it implies a buildup of the interplanetary magnetic field (IMF) magnitude over time. Such a buildup is ruled out by observations. Magnetic reconnection above helmet streamer configurations could provide a mechanism for maintaining the observed relative constancy of the IMF (McComas et al., 1989). McComas et al. (1991) showed observational evidence of reconnection above a streamer. This interpretation is investigated using time-dependent MHD simulations. The opening of new magnetic flux on the sun (as might occur in a CME or other transient event) is modeled as an increase in magnetic flux at the poles of a simulated corona. It is found that this perturbation can in fact cause reconnection above an equatorial helmet streamer, and the resultant density signature is similar to the observations of McComas et al. (1991).

Linker, Jon A.↗
NASA NTRS2023

Is the Global MHD Modeling of the Magnetosphere Adequate for GIC Prediction: the May 27–28, 2017 Storm

Practical steps taken by the international community to reduce the damage to technological systems from space weather include the development of numerical models capable of real-time predictions of electromagnetic disturbances at the Earth’s surface. Here we examine the feasibility of a version of the Space Weather Modeling Framework (SWMF) global MHD simulation code similar to that used by the NOAA Space Weather Prediction Center to predict the level of geomagnetic field variability, and consequently geomagnetically induced currents (GICs). We consider the contribution of geomagnetic disturbances to the bursts of GIC in the electric power line of the Kola Peninsula during the May 27–28, 2017 storm and compare the observations with results of the global MHD model. During the maximal disturbance magnetic field variations at East Scandinavian stations become more irregular, as intense Pi3 pulsations are superposed on the magnetic bay. These pulsations are not quasi-sinusoidal waves like typical Pc5 pulsations, but they are rather a quasi-periodic sequence of magnetic impulses with time scales ~5–15 min. During this period with elevated Pi3 activity very high values of GIC were recorded (variations >100 A) in the electric power transmission line. The SWMF modeling reasonably well reproduces the global magnetospheric parameters, such as SYM-H index or cross-polar potential. However, the magnetic field variability dB/dt in the East Scandinavia predicted by the modeling has turned out to be more than order of magnitude less than that observed. Thus, the version of SWMF with the grid used by NOAA SWPC still cannot adequately predict for the May 27–28 event the fine structure of the storm/substorm—Pi3 geomagnetic disturbances, and consequently the magnitude of the GIC that they drive.

geomagnetically induced currents↗
NASA NTRS1998

The Distant Tail at 200 R(sub E): Comparison Between Geotail Observations and the Results from a Global Magnetohydrodynamic Simulation

This paper reports a comparison between Geotail observations of plasmas and magnetic fields at 200 R(sub E) in the Earth's magnetotail with results from a time-dependent, global magnetohydrodynamic simulation of the interaction of the solar wind with the magnetosphere. The study focuses on observations from July 7, 1993, during which the Geotail spacecraft crossed the distant tail magnetospheric boundary several times while the interplanetary magnetic field (IMF) was predominantly northward and was marked by slow rotations of its clock angle. Simultaneous IMP 8 observations of solar wind ions and the IMF were used as driving input for the MHD simulation, and the resulting time series were compared directly with those from the Geotail spacecraft. The very good agreement found provided the basis for an investigation of the response of the distant tail associated with the clock angle of the IMF. Results from the simulation show that the stresses imposed by the draping of magnetosheath field lines and the asymmetric removal of magnetic flux tailward of the cusps altered considerably the shape of the distant tail as the solar wind discontinuities convected downstream of Earth. As a result, the cross section of the distant tail was considerably flattened along the direction perpendicular to the IMF clock angle, the direction of the neutral sheet following that of the IMF. The simulation also revealed that the combined action of magnetic reconnection and the slow rotation of the IMF clock angle led to a braiding of the distant tail's magnetic field lines along the axis of the tail, with the plane of the braid lying in the direction of the IMF.

Berchem, J.↗
NASA NTRS1996

Magnetic Field Orientation Effects on the Standoff Distance of Earth's Bow Shock

Three-dimensional, global MHD simulations of solar wind flow onto a prescribed magnetopause obstacle are used to show that a bow shock's nose location a(sub s), and the relative subsolar magnetosheath thickness Delta(sub ms)/a(sub mp) are strong functions of the IMF cone angle theta (between v(sub sw) and B(sub sw)) and the Alfven Mach number M(sub A). For a given M(sub A) the shock is more distant for higher theta (restricted to the interval 0-90deg by symmetries), while a(sub s)/a(sub mp) and Delta(sub ms/a(sub mp) increase with decreasing M(sub A) for theta greater than or approximately 20deg but decrease with decreasing M(sub A) for theta approximately Odeg. Large differences in Delta(sub ms/a(sub mp) are predicted between theta = Odeg and 90deg at low M(sub A), with smaller differences remaining even at M(sub A) approximately 10. The theta = Odeg results confirm and extend the previous work of Spreiter and Rizzi [1974]. The simulations show that successful models for the subsolar shock location cannot subsume the dependences on M(sub A) and theta into a sole dependence on M(ms). Instead, they confirm a recent prediction [Cairns and Grabbe, 1994] that a(sub s)/a(sub mp) and Delta(sub ms)/a(sub mp) should depend strongly on theta and M(sub A) for M(sub A) less than or approximately 10 (as well as on other MHD variables). Detailed comparisons between theory and data remain to be done. However, preliminary comparisons show good agreement, with distant shock locations found for low M(sub A) and large theta greater than or approximately 45deg and closer locations found for theta less than or approximately 20deg even at M/A approximately 8.

Cairns, Iver H.↗
DOE OSTI2024

On the Geometry of the Near-core Magnetic Field in Massive Stars

It is well-known that the cores of massive stars sustain a stellar dynamo with a complex magnetic field configuration. However, the same cannot be said for the field's strength and geometry at the convective–radiative boundary, which are crucial when performing asteroseismic inference. In this Letter, we present 3D magnetohydrodynamic (MHD) simulations of a 7 M ⊙ mid-main-sequence star, with particular attention given to the convective–radiative boundary in the near-core region. Our simulations reveal that the toroidal magnetic field is significantly stronger than the poloidal field in this region, contrary to recent assumptions. Moreover, the rotational shear layer, also important for asteroseismic inference, is specifically confined within the extent of the Brunt–Väisälä frequency peak. These results, which are based on the inferred properties of HD 43317, have widespread implications for asteroseismic studies of rotation, mixing, and magnetism in stars. While we expect our results to be broadly applicable across stars with similar Brunt–Väisälä frequency profiles and stellar masses, we also expect the MHD parameters (e.g., Re m ) and the initial stellar rotation rate to impact the geometry of the field and differential rotation at the convective–radiative interface.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

The Production of Jets From Magnetic Accretion Disks: Simulation of the Blandford-Payne Mechanism

We have performed magnetohydrodynamic (MRD) simulations of the production of jets from magnetized accretion disks with a factor of 5 greater extent in space and time, and with more models, than any study published so far. We find that jets are produced by such disks in a broad range of parameter space, and by at least two different mechanisms. We also are able to follow the propagation of the jet well beyond the accretion disk into the region of hydrodynamic collimation. The code used is our MHD simulation code FLOW (K. Lind, D. Payne, D. Meier, and R. Blandford, 1989), converted to run on Caltech's massively parallel Intel Touchstone Delta supercomputer. Some of these models may be directly applicable to observed radio sources.

Blandford-Payne↗
NASA NTRS1992

Observation and simulation of the radial evolution and stream structure of solar wind turbulence

This paper presents a brief overview of the observed evolution in a variety of quantities describing the turbulent evolution of the interplanetary plasma and describes simulation results consistent with many features of the evolution. The turbulence is manifested through a dissipation at small scales in the inner heliosphere with a corresponding evolution in the breakpoint between a relatively flat and a Kolmogoroff spectrum; an evolution from kinetically to (slightly) magnetically dominated energy of the plasma fluctuations; a general decrease in the cross helicity or 'Alfvenicity'; changes in the anisotropy of the fluctuations; and the increasing predominance of quasi-pressure-balanced structures in the compressive component of the fluctuations. MHD simulations with shear layers either side of a central current sheet show that even in the absence of compressibility the lack of a mean field along the direction of the main flow in the current sheet leads to rapid nonlinear evolution and the observed characteristics of 'Elsasser spectra' of the fields in the inner heliosphere. Adding compressibility to the simulations does not greatly change the 'incompressive' quantities but leads in addition to observed correlations between a measure of compression and other quantities.

Roberts, D. A.↗
NASA NTRS2022

How a Realistic Magnetosphere Alters the Polarizations of Surface, Fast Magnetosonic, and Alfvén Waves

System-scale magnetohydrodynamic (MHD) waves within Earth's magnetosphere are often understood theoretically using box models. While these have been highly instructive in understanding many fundamental features of the various wave modes present, they neglect the complexities of geospace such as the inhomogeneities and curvilinear geometries present. Here, we show global MHD simulations of resonant waves impulsively excited by a solar wind pressure pulse. Although many aspects of the surface, fast magnetosonic (cavity/waveguide), and Alfvén modes present agree with the box and axially symmetric dipole models, we find some predictions for large-scale waves are significantly altered in a realistic magnetosphere. The radial ordering of fast mode turning points and Alfvén resonant locations may be reversed even with monotonic wave speeds. Additional nodes along field lines that are not present in the displacement/velocity occur in both the perpendicular and compressional components of the magnetic field. Close to the magnetopause, the perpendicular oscillations of the magnetic field have the opposite handedness to the velocity. Finally, widely used detection techniques for standing waves, both across and along the field, can fail to identify their presence. We explain how all these features arise from the MHD equations when accounting for a non-uniform background field and propose modified methods that might be applied to spacecraft observations.

ULF Waves↗
NASA NTRS1995

The tilts of corotating interaction regions at mid heliographic latitudes

Ulysses observations of corotating interaction regions (CIRs) at mid heliographic latitudes have shown that the flow downstream of the forward shock (or wave) on the leading edge of a CIR generally turns northward and westward, while the flow downstream of the reverse wave on the trailing edge generally turns southward and eastward. These systematic flow deflections are a natural consequence of large scale pressure gradients associated with the CIRs, and indicate that the forward waves tend to propagate toward and across the equator with increasing heliocentric distance, while the reverse waves tend to propagate toward the pole. Recent determinations of CIR shock normals using the Ulysses magnetic field data appear to confirm these plasma results (Burton, private communication). Numerical simulations indicate that these effects (which imply that CIRs are systematically tilted in the north-south direction at mid latitudes) are a natural consequence of the tilt of the solar magnetic dipole axis relative to the solar rotation axis. The present work utilizes a variety of techniques to analyze the flow deflections observed within CIRs from which we can infer the overall orientations of the CIRs and the speeds and directions of propagation of the waves. Where possible, the observations are quantitatively compared with the results of 3-dimensional MHD simulations.

Riley, Pete↗
DOE OSTI2024

Simulation of DIII-D disruption with argon pellet injection and runaway electron beam

Abstract The next generation of large tokamaks, including ITER, will be equipped with a disruption mitigation system (DMS) that can be activated if a disruption is deemed to be imminent. Introducing impurities by pellet (large or shattered) or massive gas injection has been shown to be an effective mitigation mechanism on many tokamaks. The goal of the mitigation is to lessen the thermal and electromagnetic loads from the disruption without generating enough high-energy (runaway) electrons to damage the device. Variations of this mitigation process with impurity injection are presently being tested on many experiments. We have modeled one such impurity injection experiment on DIII-D using the M3D-C1 nonlinear 3D extended MHD code (Jardin et al 2012 Comput. Sci. Discovery 6 014002), The model includes an argon large pellet injection and ablation model, impurity ionization, recombination, and radiation, and runaway electron formation and subsequent evolution, including both Dreicer and avalanche sources. We obtain reasonable agreement with the experimental results for the timescale of the thermal and current quench and for the magnitude of the runaway electron plateau formed during the mitigation. This is the first 3D full MHD simulation with pellets and REs to simulate the disruption process and it also provides a partial validation of the M3D-C1 DMS model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Improved melt model for power flow

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm 2 in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO ® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO ® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe + ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

43 PARTICLE ACCELERATORS↗
NASA NTRS2004

Evolutionary Models of Cold, Magnetized, Interstellar Clouds

We modeled the long-term and small-scale evolution of molecular clouds using direct 2D and 3D magnetohydrodynamic (MHD) simulations. This work followed up on previous research by our group under auspices of the ATP in which we studied the energetics of turbulent, magnetized clouds and their internal structure on intermediate scales. Our new work focused on both global and smallscale aspects of the evolution of turbulent, magnetized clouds, and in particular studied the response of turbulent proto-cloud material to passage through the Galactic spiral potential, and the dynamical collapse of turbulent, magnetized (supercritical) clouds into fragments to initiate the formation of a stellar cluster. Technical advances under this program include developing an adaptive-mesh MHD code as a successor to ZEUS (ATHENA) in order to follow cloud fragmentation, developing a shearing-sheet MHD code which includes self-gravity and externally-imposed gravity to follow the evolution of clouds in the Galactic potential, and developing radiative transfer models to evaluate the internal ionization of clumpy clouds exposed to external photoionizing UV and CR radiation. Gammie's work at UIUC focused on the radiative transfer aspects of this program.

Gammie, Charles F.↗
DOE OSTI2024

Stellar wind impact on early atmospheres around unmagnetized Earth-like planets

ABSTRACT Stellar rotation at early ages plays a crucial role in the survival of primordial atmospheres around Earth-mass exoplanets. Earth-like planets orbiting fast-rotating stars may undergo complete photoevaporation within the first few hundred Myr driven by the enhanced stellar XUV [X-rays and extreme ultraviolet (EUV)] radiation, while planets orbiting slow-rotating stars are expected to experience difficulty in losing their primordial envelopes. Besides the action of stellar radiation, stellar winds induce additional erosion on these primordial atmospheres, altering their morphology, extent, and causing supplementary atmospheric losses. In this paper, we study the impact of activity-dependent stellar winds on primordial atmospheres to evaluate the extent to which the action of these winds can be significant in the whole planetary evolution at early evolutionary stages. We performed 3D magnetohydrodynamical (MHD) simulations of the interaction of photoevaporating atmospheres around unmagnetized Earth-mass planets in the time span between 50 and 500 Myr, analysing the joint evolution of stellar winds and atmospheres for both fast- and slow-rotating stars. Our results reveal substantial changes in the evolution of primordial atmospheres when influenced by fast-rotating stars, with a significant reduction in extent at early ages. In contrast, atmospheres embedded in the stellar winds from slow-rotating stars remain largely unaltered. The interaction of the magnetized stellar winds with the ionized upper atmospheres of these planets allows us to evaluate the formation and evolution of different MHD structures, such as double bow shocks and induced magnetospheres. This work will shed light on the first evolutionary stages of Earth-like exoplanets, which are of crucial relevance in terms of planet habitability.

Astronomy & Astrophysics↗
DOE OSTI2022

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counterpropagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scattering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. Finally, we compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counter- propagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scat- tering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. We compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Shear-induced opening of the coronal magnetic field

This work describes the evolution of a model solar corona in response to motions of the footpoints of its magnetic field. The mathematics involved is semianalytic, with the only numerical solution being that of an ordinary differential equation. This approach, while lacking the flexibility and physical details of full MHD simulations, allows for very rapid computation along with complete and rigorous exploration of the model's implications. We find that the model coronal field bulges upward, at first slowly and then more dramatically, in response to footpoint displacements. The energy in the field rises monotonically from that of the initial potential state, and the field configuration and energy appraoch asymptotically that of a fully open field. Concurrently, electric currents develop and concentrate into a current sheet as the limiting case of the open field is approached. Examination of the equations shows rigorously that in the asymptotic limit of the fully open field, the current layer becomes a true ideal MHD singularity.

Wolfson, Richard↗
NASA NTRS1986

Turbulent magnetic reconnection

The results of two-dimensional spectral-method MHD simulations of the effect of turbulence on magnetic-reconnection processes are presented in extensive graphs and tables and characterized in detail. In a periodic sheet-pinch configuration, early-onset, persistent turbulence is found to occur when a low level of broadband fluctuations is introduced into the initial data. Features observed include small-scale unsteady coherent electric-current and vorticity structures, enhanced viscous and resistive dissipation, unsteady and spatially asymmetric fluid flow, large-scale and small-scale magnetic islands, multiple X points, large-amplitude electric-field fluctuations, and spectral characteristics indicating that the turbulence is broadband in both spatial directions and similar to homogeneous MHD turbulence. The importance of reconnection for studies of solar, magnetospheric, and laboratory plasmas in inhomogeneous magnetic fields is indicated.

Matthaeus, W. H.↗