SEARCH · Search NASA

Results for “Numerical magnetohydrodynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS2004

A Massively Parallel Hybrid Dusty-Gasdynamics and Kinetic Direct Simulation Monte Carlo Model for Planetary Applications

In order to understand the global structure, dynamics, and physical and chemical processes occurring in the upper atmospheres, exospheres, and ionospheres of the Earth, the other planets, comets and planetary satellites and their interactions with their outer particles and fields environs, it is often necessary to address the fundamentally non-equilibrium aspects of the physical environment. These are regions where complex chemistry, energetics, and electromagnetic field influences are important. Traditional approaches are based largely on hydrodynamic or magnetohydrodynamic (MHD) formulations and are very important and highly useful. However, these methods often have limitations in rarefied physical regimes where the molecular collision rates and ion gyrofrequencies are small and where interactions with ionospheres and upper neutral atmospheres are important. At the University of Michigan we have an established base of experience and expertise in numerical simulations based on particle codes which address these physical regimes. The Principal Investigator, Dr. Michael Combi, has over 20 years of experience in the development of particle-kinetic and hybrid kinetichydrodynamics models and their direct use in data analysis. He has also worked in ground-based and space-based remote observational work and on spacecraft instrument teams. His research has involved studies of cometary atmospheres and ionospheres and their interaction with the solar wind, the neutral gas clouds escaping from Jupiter s moon Io, the interaction of the atmospheres/ionospheres of Io and Europa with Jupiter s corotating magnetosphere, as well as Earth s ionosphere. This report describes our progress during the year. The contained in section 2 of this report will serve as the basis of a paper describing the method and its application to the cometary coma that will be continued under a research and analysis grant that supports various applications of theoretical comet models to understanding the inner comae of comets (grant NAGS- 13239 from the Planetary Atmospheres program).

Combi, Michael R.↗

NASA NTRS2012

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗

NASA NTRS1997

The Geometric Spreading of Coronal Plumes and Coronal Holes

The geometric spreading in plumes and in the interplume region in coronal holes is calculated, using analytic and numerical theoretical models, between 1.0 and 5.0 solar radius. We apply two scale approximation that permits the rapid local spreading at the base of plumes to be evaluated separately from the global spreading imposed by coronal hole geometry. We show that fl can be computed from a potential field model and fg can be computed from global magnetohydrodynamic simulations of coronal structure. The approximations are valid when the plasma beta is small with respect to unity and for a plume separation small with respect to a solar radius.

Suess, S. T.↗

NASA NTRS1976

Magnetohydrodynamic and gasdynamic aspects of solar-wind flow around terrestrial planets: A critical review

Problems in space physics are discussed whose models, in simplified form, reduce to a supersonic flow scheme with a detached shock wave, namely: (1) solar wind interaction with an intrinsic planetary magnetic field; (2) solar wind interaction with the ionized component of the atmosphere of a comet; and (3) solar wind interaction with the ionosphere of a planet which does not possess its own magnetic field. The numerical study of the above problems is performed with the use of magnetogasdynamic equations for an ideal single fluid model. From the physical veiwpoint, the problems are solved in terms of as simple phenomena as possible; the principal objective is to make recently developed methods of numerical analysis of mixed flows applicable to space physics problems.

Spreiter, J. R.↗

NASA NTRS1983

Propagation of solar disturbances - Theories and models

Recent theoretical developments and construction of several models for the propagation of solar disturbances from the sun and their continuation throughout heliospheric space are discussed. Emphasis centers on physical mechanisms as well as mathematical techniques (i.e., analytical and numerical methods). This outline will lead to a discussion of the state-of-the-art of theoretically based modeling efforts in this area. It is shown that the fundamental theory for the study of propagation of disturbances in heliospheric space is centered around the self-consistent analysis of wave and mass motion within the context of magnetohydrodynamics in which the small scale structures will be modified by kinetic effects. Finally, brief mention is made of some interesting problems for which attention is needed for advancement of the understanding of the physics of large scale propagation of solar disturbances in heliospheric space.

Wu, S. T.↗

NASA NTRS1991

A study of physical mechanisms for filament eruption and coronal mass ejection via numerical simulation

It is well-known from both observation and theory that photospheric shear motion has played a key role in physical processes of the energy build-up and release for the solar flare. In order to further our understanding of the occurrence of solar flares it is necessary to investigate the triggering mechanism. One popular scenario for the onset of eruptive solar flares is that in response to photospheric shear motions the magnetic field evolves slowly through a series of magnetohydrodynamic-equilibria until a threshold is reached where magnetohydrodynamic (MHD) non-equilibrium sets in. Thus, a magnetic eruption occurs, causing the solar flare. To substantiate this claim we have employed our newly developed three-dimensional time-dependent MHD code with gravity to simulate the evolution of the coronal field. We use plasma beta = 0.1 to closely approximate the condition in the actual corona. Some preliminary results are presented.

Wu, S. T.↗

NASA NTRS2022

The Dynamic Coupling of Streamers and Pseudostreamers to the Heliosphere

The slow solar wind is generally believed to result from the interaction of open and closed coronal magnetic flux at streamers and pseudostreamers. We use three-dimensional magnetohydrodynamic simulations to determine the detailed structure and dynamics of open-closed interactions that are driven by photospheric convective flows. The photospheric magnetic field model includes a global dipole giving rise to a streamer together with a large parasitic polarity region giving rise to a pseudostreamer that separates a satellite coronal hole from the main polar hole. Our numerical domain extends out to 30Rꙩ and includes an isothermal solar wind, so that the coupling between the corona and heliosphere can be calculated rigorously. This system is driven by imposing a large set of quasi-random surface flows that capture the driving of coronal flux in the vicinity of streamer and pseudostreamer boundaries by the supergranular motions. We describe the resulting structures and dynamics. Interchange reconnection dominates the evolution at both streamer and pseudostreamer boundaries, but the details of the resulting structures are clearly different from one another. Additionally, we calculate in situ signatures of the reconnection and determine the dynamic mapping from the inner heliosphere back to the Sun for a test spacecraft orbit. We discuss the implications of our results for interpreting observations from inner heliospheric missions, such as Parker Solar Probe and Solar Orbiter, and for space weather modeling of the slow solar wind.

V. Aslanyan↗

NASA NTRS2023

Understanding the Deflection of the “Cartwheel CME”: Data Analysis and Modeling

We study the low corona evolution of the “Cartwheel” coronal mass ejection (CME; 2008 April 9) by- reconstructing its three-dimensional path and modeling it with magnetohydrodynamic simulations. This event exhibited a double deflection that has been reported and analyzed in previous works but whose underlying cause remained unclear. The Cartwheel CME traveled toward a coronal hole (CH) and against the magnetic gradients. Using a high-cadence, full-trajectory reconstruction, we accurately determine the location of the magnetic flux rope (MFR) and, consequently, the magnetic environment in which it is immersed. We find a pseudostreamer (PS) structure whose null point may be responsible for the complex evolution of the MFR at the initial phase. From the preeruptive magnetic field reconstruction, we estimate the dynamic forces acting on the MFR and provide a new physical insight into the motion exhibited by the 2008 April 9 event. By setting up a similar magnetic configuration in a 2.5D numerical simulation we are able to reproduce the observed behavior, confirming the importance of the PS null point. We find that the magnetic forces directed toward the null point cause the first deflection, directing the MFR toward the CH. Later, the magnetic pressure gradient of the CH produces the reversal motion of the MFR.

Solar prominences↗

NASA NTRS2015

Deep Chandra Observation and Numerical Studies of the Nearest Cluster Cold Front in the Sky

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are approximately 10 per cent brighter than the surrounding gas and are separated by approximately 15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches approximately 5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

hydrodynamics↗

NASA NTRS1991

MHD shocks in coronal mass ejections

The primary objective of this research program is the study of the magnetohydrodynamic (MHD) shocks and nonlinear simple waves produced as a result of the interaction of ejected lower coronal plasma with the ambient corona. The types of shocks and nonlinear simple waves produced for representative coronal conditions and disturbance velocities were determined. The wave system and the interactions between the ejecta and ambient corona were studied using both analytic theory and numerical solutions of the time-dependent, nonlinear MHD equations. Observations from the SMM coronagraph/polarimeter provided both guidance and motivation and are used extensively in evaluating the results. As a natural consequence of the comparisons with the data, the simulations assisted in better understanding the physical interactions in coronal mass ejections (CME's).

Steinolfson, R. S.↗

NASA NTRS1995

Radial energy transport by magnetospheric ULF waves: Effects of magnetic curvature and plasma pressure

The 'radial' transport of energy by internal ULF waves, stimulated by dayside magnetospheric boundary oscillations, is analyzed in the framework of one-fluid magnetohydrodynamics. (the term radial is used here to denote the direction orthogonal to geomagnetic flux surfaces.) The model for the inhomogeneous magnetospheric plasma and background magnetic field is axisymmetric and includes radial and parallel variations in the magnetic field, magnetic curvature, plasma density, and low but finite plasma pressure. The radial mode structure of the coupled fast and intermediate MHD waves is determined by numerical solution of the inhomogeneous wave equation; the parallel mode structure is characterized by a Wentzel-Kramer-Brillouin (WKB) approximation. Ionospheric dissipation is modeled by allowing the parallel wave number to be complex. For boudnary oscillations with frequencies in the range from 10 to 48 mHz, and using a dipole model for the background magnetic field, the combined effects of magnetic curvature and finite plasma pressure are shown to (1) enhance the amplitude of field line resonances by as much as a factor of 2 relative to values obtained in a cold plasma or box-model approximation for the dayside magnetosphere; (2) increase the energy flux delivered to a given resonance by a factor of 2-4; and (3) broaden the spectral width of the resonance by a factor of 2-3. The effects are attributed to the existence of an 'Alfven buoyancy oscillation,' which approaches the usual shear mode Alfven wave at resonance, but unlike the shear Alfven mode, it is dispersive at short perpendicular wavelengths. The form of dispersion is analogous to that of an internal atmospheric gravity wave, with the magnetic tension of the curved background field providing the restoring force and allowing radial propagation of the mode. For nominal dayside parameters, the propagation band of the Alfven buoyancy wave occurs between the location of its (field line) resonance and that of the fast mode cutoff that exists at larger radial distances.

Kouznetsov, Igor↗

NASA NTRS2011

Geomagnetic Variability and Predictability

The Earth's magnetic field varies on time scales from seconds to millions of years. These variations can be measured directly by ground observatories and satellites, and indirectly through laboratories. They are from different electromagnetic processes in magnetosphere and ionosphere (external), ocean and crust (surface), and deep in the Earth's fluid outer core (internal). Of the measured magnetic signals, 95% is from the core, and is called the core field, or internal field. Variation of the core field, often called the geomagnetic secular variation (SV), is the manifestation of the magnetohydrodynamic processes in the core. Therefore, SV provides rich information on the core dynamical state that is then critical for understanding the geophysical mechanisms of the SV, the core-mantle interactions, and contribution of mass transport inside the core on Earth's gravity variation. Traditionally, the core state is probed via two independent methodologies: core flow inversion from the observed SV, and pure numerical dynamo simulation. The former is driven by observations, but does not include necessary interactions among physical variables; the latter is dynamically consistent, but is not constrained by any observation. Therefore, both are very limited in providing appropriate estimates of the core state. To avoid these limitations, a new approach, geomagnetic data assimilation, appears recently. In this approach, simulation results (model forecasts) are constantly corrected with observations, the corrected solutions (analysis), are then used as the initial conditions to make more accurate forecasts of future. Recent studies have demonstrated successfully that assimilation solutions are different from those of free running models. These new solutions are better estimates of the core state, and have been used to provide accurate prediction of SV for the period from 2010 to 2015. This prediction is part of the IGRF field model for international community applications.

Kuang, Weijia↗

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 NTRS1989

Helical, dissipative, magnetohydrodynamic states with flow

It is shown that for an axially periodic column of magnetofluid driven by an applied axial electric field, the total rate of energy dissipation (Ohmic plus viscous) can be lowered by permitting a helical component with vortical flow in the solution. The principle of minimum energy-dissipation rate suggests that this partially helical state will be preferred to the axisymmetric one that exists for the same parameters. The result is consistent with the repeated appearance of such partially-helical states in several fully three-dimensional numerical computations and is not inconsistent with the data from some confinement experiments.

Montgomery, David↗

NASA NTRS1995

Small scales formation via Alfven wave propagation in compressible nonuniform media

In weakly dissipative media governed by the magnetohydrodynamics (MHD) equations, any efficient mechanism of energy dissipation requires the formation of small scales. The possibility to produce small scales has been studied by Malara et al. in the case of MHD disturbances propagating in an incompressible and inhomogeneous medium, for a strictly 2D geometry. We extend the work of Malara et al. to include both compressibility and the third component for vector quantities. Using numerical simulations we show that, when an Alfven wave propagates in a compressible nonuniform medium, the two dynamical effects responsible for the small scales formation in the incompressible case are still at work: energy pinching and phase-mixing. Moreover, the interaction between the initial Alfven wave and the inhomogeneity gives rise to the formation of compressible perturbations (fast and slow waves or a static entropy wave). Some of these compressive fluctuations are subject to the steepening of the wave front and become shock waves, which are extremely efficient in dissipating their energy, their dissipation being independent of the Reynolds number. A rough estimate of the typical times which the various dynamical processes take to produce small scales and then to dissipate the energy show that these times are consistent with those required to dissipate inside the solar corona the energy of Alfven waves of photospheric origin.

Malara, F.↗

NASA NTRS1976

A two-region model of the solar wind including azimuthal velocity

The two-region model of the solar wind divides the interplanetary space into two regions: it assumes that the solar wind is one-fluid in an inner region within 0.4 AU and two-fluid in an outer region beyond 0.4 AU. This paper includes the angular motion of the solar wind in the two-region model. The flow in the one-fluid region is governed by the one-fluid magnetohydrodynamic equations. The second and third moment equations of the Vlasov equation together with other conservation equations are used to describe the solar-wind flow in the two-fluid region. The predicted azimuthal velocity at 1 AU is less than 2 km/s. All other macroscopic and microscopic properties from this model are in good agreement with experimental quiet-time observations at 1 AU. The numerical results also confirm that when the azimuthal velocity is included in the analysis, the amount of magnetic-field energy converted into kinetic energy in the solar wind is only a small fraction of the total expansion energy flux and has little effect upon the final radial expansion velocity.

Acuna, M. H.↗

NASA NTRS1980

Selective decay hypothesis at high mechanical and magnetic Reynolds numbers

Implications of certain applications of turbulence theory to two-dimensional turbulence and magnetohydrodynamic flow are discussed. It is shown that the use of the Navier-Stokes equation (NSE) for measurements of turbulent fluctuations has been effective only for three-dimensional flows. For two-dimensional flows, used for the study of large-scale motions in the atmosphere or ocean, enstrophy is cascaded to high wave numbers and dissipated at a finite rate even at infinite Re. MHD flows are numerically calculated for the two-dimensional case and analytically for the three-dimensional case, for which discrepancies in the relative rates of energy and cross helicity decay lead to a recommendation that numerical calculations for the three-dimensional case be carried out to determine the precise decayed states.

Matthaeus, W. H.↗

NASA NTRS1998

The Geometric Spreading of Coronal Plumes and Coronal Holes

The geometric spreading in plumes and in the interplume region in coronal holes is calculated, using analytic and numerical theoretical models, between 1.0 and 5.0 solar radius. We apply a two-scale approximation that permits the rapid local spreading at the base of plumes (f(sub t)) to be evaluated separately from the global spreading (f(sub g)) imposed by coronal hole geometry. We show that f(sub t) can be computed from a potential-field model and f(sub g) can be computed from global magnetohydrodynamic simulations of coronal structure. The approximations are valid when the plasma beta is mail with respect to unity and for a plume separation small with respect to a solar radius.

Suess, S. T.↗