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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

Interchange Reconnection Dynamics in A Solar Coronal Pseudo-Streamer

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

magnetic fields↗
NASA NTRS2014

Realistic Modeling of Multi-Scale MHD Dynamics of the Solar Atmosphere

Realistic 3D radiative MHD simulations open new perspectives for understanding the turbulent dynamics of the solar surface, its coupling to the atmosphere, and the physical mechanisms of generation and transport of non-thermal energy. Traditionally, plasma eruptions and wave phenomena in the solar atmosphere are modeled by prescribing artificial driving mechanisms using magnetic or gas pressure forces that might arise from magnetic field emergence or reconnection instabilities. In contrast, our 'ab initio' simulations provide a realistic description of solar dynamics naturally driven by solar energy flow. By simulating the upper convection zone and the solar atmosphere, we can investigate in detail the physical processes of turbulent magnetoconvection, generation and amplification of magnetic fields, excitation of MHD waves, and plasma eruptions. We present recent simulation results of the multi-scale dynamics of quiet-Sun regions, and energetic effects in the atmosphere and compare with observations. For the comparisons we calculate synthetic spectro-polarimetric data to model observational data of SDO, Hinode, and New Solar Telescope.

Realistic Modeling↗
NASA NTRS2004

Plasma Sail Concept Fundamentals

The mini-magnetospheric plasma propulsion (M2P2) device, originally proposed by Winglee et al., predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to unprecedented speeds of 50 C80 km/s after an acceleration period of 3 mo. Such velocities will enable travel out of the solar system in period of 7 yr almost an order of magnitude improvement over present chemical-based propulsion systems. However, for the parameters of the simulation of Winglee et al., a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the magnetohydrodynamic (MHD) fluid model, normally applied to planetary magnetospheres, that the characteristic scale size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or comparable to the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. A two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail has been adopted. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the forces delivered to the innermost regions of the plasma sail are considerably ( 10 times) smaller than the MHD counterpart, are dominated by the magnetic field pressure gradient, and are directed primarily in the transverse direction.

Khazanov, G. V.↗
NASA NTRS1991

Forced magnetic reconnection in a plasma sheet with localized resistivity profile excited by lower hybrid drift type instability

A forced magnetic reconnection process with a temporal evolution of resistivity is studied for a plasma sheet with a nonuniform resistivity profile based on the nonlocal mode structure of the lower hybrid drift type instability. The growth rate of the mode found is almost independent of the resistivity at the neutral sheet, but depends on the resistivity of the region of maximum density gradient away from the neutral sheet. This is studied by using both a nonlinear numerical MHD simulation and a linear theory. The mode may be relevant to the prevalent theoretical concept of MHD reconnection and the localized anomalous resistivity profile based on the lower hybrid drift instability.

Hoshino, M.↗
NASA NTRS1990

Simulation of period doubling of recurrent solar wind structures

Based on satellite observations of a recurrent solar wind structure conducted in 1974, an MHD simulation model, and input functions generated from plasma and magnetic field data, the continuing evolution of the solar wind structure outside 5 AU is studied. The model uses the Rankine-Hugoniot relations to describe the jumps in flow properties across the shocks, and it treats shocks as surfaces of discontinuity with zero thickness. Two interaction processes (the collision and the merging of shocks) play important roles in restructuring the solar wind in the outer heliosphere. The simulation result shows that period doubling occurs between 5 and 10 AU. The recurrent solar wind appears to be a persistent new structure between 10 and 20 AU, and it consists of one merged interaction region per solar rotation.

Whang, Y. C.↗
NASA NTRS1990

Parallel electric fields in a simulation of magnetotail reconnection and plasmoid evolution

Properties of the electric field component parallel to the magnetic field are investigate in a 3D MHD simulation of plasmoid formation and evolution in the magnetotail, in the presence of a net dawn-dusk magnetic field component. The spatial localization of E-parallel, and the concept of a diffusion zone and the role of E-parallel in accelerating electrons are discussed. A localization of the region of enhanced E-parallel in all space directions is found, with a strong concentration in the z direction. This region is identified as the diffusion zone, which plays a crucial role in reconnection theory through the local break-down of magnetic flux conservation.

Hesse, M.↗
NASA NTRS1987

Magnetic field draping in the Comet Halley coma - Comparison of Vega observations with computer simulations

During the Vega-1 encounter with Comet Halley, the magnetometer observed draping and compression of the interplanetary magnetic field. These are reproduced well by a three-dimensional MHD simulation of the cometary interaction. Rotations in the magnetic field similar to those at closest approach are also observed 2.75 hours earlier. It is suggested that both rotations correspond to the same IMF interval and that the spacecraft had overtaken the plasma and encountered 'older' magnetic field as it penetrated the coma. Analysis of the MHD model indicates that it should take about 3 to 5 hours for a solar wind parcel to pass from the unperturbed solar wind to Vega-1 at closest approach. A simulated magnetic field profile composed of nested sections for different IMF orientations closely resembles the observations. This result supports the hypothesis of layered magnetic orientations in the coma.

Schwingenschuh, K.↗
DOE OSTI2023

Magnetic Reconnection in Solar Flares [Slides]

Magnetic reconnection is the underlying mechanism in solar flares. Magnetic reconnection is an astrophysical process in an electrically conducting plasma where neighboring magnetic field lines running anti-parallel break and reconnection therefore exchanging the direction of their paths. The result of this process is the magnetic field energy converting to kinetic and thermal energy which is explosively released. The scale of a solar is extremely large compared to Earth so kinetic simulations are unable to portray an accurate model. In this project, we are using an astrophysical magnetohydrodynamic (MHD) simulation code known as Athena++ with the reconnection problem generator file created by Li et al. 2018. To further study this simulation, we began testing a heating/cooling term. The results of the original file will be compared with the modified file with an implemented heating/cooling term to better understand how solar flares function.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Generation of slow-mode waves in front of the dayside magnetopause

Slow-mode waves have been observed to appear frequently in front of the dayside magnetopause. It is found based on two-dimensional global magnetohydrodynamic (MHD) simulations that slow-mode waves are generated through the interaction between the bow shock and various MHD waves (fast-mode, Alfven-mode, or slow-mode waves) in the upstream solar wind. The generated slow-mode waves stay in front of the magnetopause for a long time (over 15 minutes) before the wave energy is convected away tailward. Since various waves are often present in the solar wind, this mechanism may lead to the frequent appearance of slow-mode waves in front of the magnetopause.

Yan, M.↗
NASA NTRS1987

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗
DOE OSTI2024

Three-dimensional simulations of magneto-inertial Magnetic-Direct-Drive targets

For the cylindrically symmetric targets that are normally fielded on the Z machine, two dimensional axisymmetric MHD simulations provide the backbone of our target design capability. These simulations capture the essential operation of the target and allow for a wide range of physics to be addressed at a substantially lower computational cost than 3D simulations. This approach, however, makes some approximations that may impact its ability to accurately provide insight into target operation. As an example, in 2D simulations, targets are able to stagnate directly to the axis in a way that is not entirely physical, leading to uncertainty in the impact of the dynamical instabilities that are an important source of degradation for ICF concepts. In this report, we have performed a series of 3D calculations in order to assess the importance of this higher fidelity treatment on MagLIF target performance.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Parallel electric fields in a simulation of magnetotail reconnection and plasmoid evolution

Properties of the electric field component parallel to the magnetic field (E sub parallel) in a three-dimensional MHD simulation of plasmoid formation and evolution in the magnetotail in the presence of a net dawn-dusk magnetic field component were observed. Particularly emphasized was the spatial location of E(sub parallel), the concept of a diffusion zone and the role of E(sub parallel) in accelerating electrons. A localization of the region of enhanced E(sub parallel) in all space directions with a strong concentration in the z direction was found. This region was identified as the diffusion zone, which plays a crucial role in reconnection theory through the local break-down of magnetic flux conservation. The presence of B(sub y) implies a north-south asymmetry of the injection of accelerated particles into the near-earth region, if the net B(sub y) field is strong enough to force particles to follow field lines through the diffusion region. It is estimated that for a typical net B(sub y) field this should affect the injection of electrons into the near-earth dawn region, so that precipitation into the Northern (Southern) Hemisphere should dominate for duskward (dawnward) net B(sub y). In addition, a spatial clottiness of the expected injection of adiabatic particles which could be related to the appearance bright spots in auroras was observed.

Hesse, Michael↗
NASA NTRS1993

Three-dimensional magnetohydrodynamics of the emerging magnetic flux in the solar atmosphere

The nonlinear evolution of an emerging magnetic flux tube or sheet in the solar atmosphere is studied through 3D MHD simulations. In the initial state, a horizontal magnetic flux sheet or tube is assumed to be embedded at the bottom of MHD two isothermal gas layers, which approximate the solar photosphere/chromosphere and the corona. The magnetic flux sheet or tube is unstable against the undular mode of the magnetic buoyancy instability. The magnetic loop rises due to the linear and then later nonlinear instabilities caused by the buoyancy enhanced by precipitating the gas along magnetic field lines. We find by 3D simulation that during the ascendance of loops the bundle of flux tubes or even the flux sheet develops into dense gas filaments pinched between magnetic loops. The interchange modes help produce a fine fiber flux structure perpendicular to the magnetic field direction in the linear stage, while the undular modes determine the overall buoyant loop structure. The expansion of such a bundle of magnetic loops follows the self-similar behavior observed in 2D cases studied earlier. Our study finds the threshold flux for arch filament system (AFS) formation to be about 0.3 x 10 exp 20 Mx.

Matsumoto, R.↗
NASA NTRS2003

High Performance Parallel Methods for Space Weather Simulations

This is the final report of our NASA AISRP grant entitled 'High Performance Parallel Methods for Space Weather Simulations'. The main thrust of the proposal was to achieve significant progress towards new high-performance methods which would greatly accelerate global MHD simulations and eventually make it possible to develop first-principles based space weather simulations which run much faster than real time. We are pleased to report that with the help of this award we made major progress in this direction and developed the first parallel implicit global MHD code with adaptive mesh refinement. The main limitation of all earlier global space physics MHD codes was the explicit time stepping algorithm. Explicit time steps are limited by the Courant-Friedrichs-Lewy (CFL) condition, which essentially ensures that no information travels more than a cell size during a time step. This condition represents a non-linear penalty for highly resolved calculations, since finer grid resolution (and consequently smaller computational cells) not only results in more computational cells, but also in smaller time steps.

Hunter, Paul↗
NASA NTRS1989

Computer simulation of Alfven resonance in a cylindrical, axially bounded flux tube

The resonant absorption of Alfven waves in an axially bounded cylindrical flux tube is investigated in a dissipative MHD simulation. It is found that in an axially bounded flux tube, in contrast to an infinite periodic model, the resonant frequency is nearly independent of the poloidal component of the magnetic field. This is a consequence of the 'ballooning' structure of the resonant Alfven waves. The scaling with resistivity and viscosity of the width of the resonance layer, the dissipation rate, and the time for steady state absorption to occur, are all in agreement with theory.

Strauss, H. R.↗
NASA NTRS2010

3D Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability: Instability of a Static Column - 1

We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.

Mizuno, Yosuke↗
NASA NTRS2010

Jovian Plasma Torus Interaction with Europa: 3D Hybrid Kinetic Simulation. First results

The hybrid kinetic model supports comprehensive simulation of the interaction between different spatial and energetic elements of the Europa-moon-magnetosphere system with respect to variable upstream magnetic field and flux or density distributions of plasma and energetic ions, electrons, and neutral atoms. This capability is critical for improving the interpretation of the existing Europa flyby measurements from the Galileo orbiter mission, and for planning flyby and orbital measurements, (including the surface and atmospheric compositions) for future missions. The simulations are based on recent models of the atmosphere of Europa (Cassidy etal.,2007;Shematovichetal.,2005). In contrast to previous approaches with MHD simulations, the hybrid model allows us to fully take into account the finite gyro radius effect and electron pressure, and to correctly estimate the ion velocity distribution and the fluxes along the magnetic field (assuming an initial Maxwellian velocity distribution for upstream background ions).Non-thermal distributions of upstream plasma will be addressed in future work. Photoionization,electron-impact ionization, charge exchange and collisions between the ions and neutrals are also included in our model. We consider two models for background plasma:(a) with O(++) ions; (b) with O(++) and S(++) ions. The majority of O2 atmosphere is thermal with an extended cold population (Cassidyetal.,2007). A few first simulations already include an induced magnetic dipole; however, several important effects of induced magnetic fields arising from oceanic shell conductivity will be addressed in later work.

Lipatov, A. S.↗
NASA NTRS2016

Do We Know the Actual Magnetopause Position for Typical Solar Wind Conditions?

We compare predicted magnetopause positions at the subsolar point and four reference points in the terminator plane obtained from several empirical and numerical MHD (magnetohydrodynamics) models. Empirical models using various sets of magnetopause crossings and making different assumptions about the magnetopause shape predict significantly different magnetopause positions (with a scatter greater than 1 Earth radius (R (sub E)) even at the subsolar point. Axisymmetric magnetopause models cannot reproduce the cusp indentations or the changes related to the dipole tilt effect, and most of them predict the magnetopause closer to the Earth than non axisymmetric models for typical solar wind conditions and zero tilt angle. Predictions of two global non axisymmetric models do not match each other, and the models need additional verification. MHD models often predict the magnetopause closer to the Earth than the non axisymmetric empirical models, but the predictions of MHD simulations may need corrections for the ring current effect and decreases of the solar wind pressure that occur in the foreshock. Comparing MHD models in which the ring current magnetic field is taken into account with the empirical Lin et al. model, we find that the differences in the reference point positions predicted by these models are relatively small for B (sub z) equals 0 (note: B (sub z) is when the Earth's magnetic field points north versus Sun's magnetic field pointing south). Therefore, we assume that these predictions indicate the actual magnetopause position, but future investigations are still needed.

Samsonov, A. A.↗