Search NASA⌕ Search

SEARCH · Search NASA

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Three-dimensional structure of coronal mass ejections

The 3D structure of coronal mass ejections (CMEs) are examined on the basis of numerical solutions of the MHD equations. The numerical solution in a meridional plane centered over the region of newly emerged flux has features similar to those computed in previous 2D simulations. The CME is found to be shaped more like an arcade with an axis parallel to the Equator than like a bubble. This characteristic arcade shape occurs both in the density distribution in meridional and azimuthal planes and in the integrated white-light brightness in meridional and polar views. Examination of the brightness change due to the CME shows that the CME in a meridional projection as viewed from the earth has the same loop shape even if the CME originates almost directly in front of or behind the sun.

Steinolfson, R. S.↗
NASA NTRS2022

Global Driving of Auroral Precipitation: 1. Balance of Sources

The accurate determination of auroral precipitation in global models has remained a daunting and rather inexplicable obstacle. Understanding the calculation and balance of multiple sources that constitute the aurora, and their eventual conversion into ionospheric electrical conductance, is critical for improved prediction of space weather events. In this study, we present a semi-physical global modeling approach that characterizes contributions by four types of precipitation—monoenergetic, broadband, electron, and ion diffuse—to ionospheric electrodynamics. The model uses a combination of adiabatic kinetic theory and loss parameters derived from historical energy flux patterns to estimate auroral precipitation from magnetohydrodynamic (MHD) quantities. It then converts them into ionospheric conductance that is used to compute the ionospheric feedback to the magnetosphere. The model has been employed to simulate the 5–7 April 2010 Galaxy15 space weather event. Comparison of auroral fluxes show good agreement with observational data sets like NOAA-DMSP and OVATION Prime. The study shows a dominant contribution by electron diffuse precipitation, accounting for ∼74% of the auroral energy flux. However, contributions by monoenergetic and broadband sources dominate during times of active upstream solar conditions, providing for up to 61% of the total hemispheric power. The study also finds a greater role played by broadband precipitation in ionospheric electrodynamics which accounts for ∼31% of the Pedersen conductance.

Agnit Mukhopadhyay↗
NASA NTRS2024

Improved Modelling of SEP Event Onset Within the WSA–Enlil–SEPMOD Framework

Multi-spacecraft observations of solar energetic particle (SEP) events not only enable a deeper understanding and development of particle acceleration and transport theories but also provide important constraints for model validation efforts. However, because of computational limitations, a given physics-based SEP model is usually best suited to capture a particular phase of an SEP event, rather than its whole development from onset through decay. For example, magnetohydrodynamic (MHD) models of the heliosphere often incorporate solar transients only at the outer boundary of their so-called coronal domain – usually set at a heliocentric distance of 20–30 R⊙. This means that particle acceleration at coronal mass ejection (CME)-driven shocks is also computed from this boundary onwards, leading to simulated SEP event onsets that can be many hours later than observed, since shock waves can form much lower in the solar corona. In this work, we aim to improve the modelled onset of SEP events by inserting a “fixed source” of particle injection at the outer boundary of the coronal domain of the coupled WSA–Enlil 3D MHD model of the heliosphere. The SEP model that we employ for this effort is Solar Energetic Particle MODel (SEPMOD), a physics-based test-particle code based on a field line tracer and adiabatic invariant conservation. We apply our initial tests and results of SEPMOD’s fixed-source option to the 2021 October 9 SEP event, which was detected at five well-separated locations in the inner heliosphere – Parker Solar Probe, STEREO-A, Solar Orbiter, BepiColombo, and near-Earth spacecraft.

Solar eruptions↗
NASA NTRS1989

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Chung, T. J.↗
NASA NTRS1994

Transition to turbulence in solar surges

Transition to turbulence in magnetohydrodynamic (MHD) tearing jets has been invoked as a mechanism underlying some of the complex behavior observed in solar surges, including deceleration of the upflowing plasma and temporal correlations with types I and III radio bursts. In this paper we investigate a possible mechanism for this transition: three-dimensional secondary instabilities on two-dimensional saturated states. We find through linear analysis that these MHD configurations -- in particular, the tearing jet -- are secondarily unstable, with the dominant energy transfer from the one-dimensional field into the 3-dimensional fields. Using nonlinear simulations, we also investigate the system evolution after the secondary modes attain finite amplitude. When the tearing jet transitions to turbulence, the total kinetic energy drops rapidly corresponding to the deceleration of the jet. The electric field grows rapidly as the primary mode saturates and the three-dimensional secondary mode develops, and then decays quickly as the tearing jet becomes turbulent, providing a possible explanation for the finite duration of the associated meter-wave bursts. The electric field decays as the magnetic and velocity fields both decay. The system is dominated at late times by spanwise modes, which strongly resemble the magnetic field-aligned filamentary flows characteristic of many surges.

Dahlburg, Russell B.↗
NASA NTRS1998

The Paradox of Filamented Coronal Hole Flow but Uniform High Speed Wind

Plumes and rays in coronal holes are nearly radially aligned density striations that follow the ambient magnetic field. They have long been known, but have gained new interest with growing awareness that coronal hole flow is inherently filamentary. In retrospect, filamentary flow should have been no surprise. This is because,Beta much less than 1 in coronal holes inside approximately 10 Solar radius, allowing the flow to be filamentary down to the smallest scale of photospheric magnetic activity. While the magnetic field itself is locally smooth across any height above ca. 50,000 km, SOHO/MDI has shown that the photospheric magnetic field is a complex array of rapidly evolving small bipoles that are constantly emerging, evolving, and cancelling. The resulting activity is manifested in microflares, concentrated in the magnetic network, that produce Impulsive injections at the footpoints of coronal field lines. The uneven distribution of this activity in space and time is the source of coronal hole filamentation. What is surprising is that the radial flow speed also exhibits filamentary structure. It is not well described as smooth, spherically symmetric, diverging flow, but instead ranges from 300 to over 1000 km/s at 5.5 Solar radius among field-aligned filaments like those seen in plumes and rays [Feldman et al., JGR, Dec. 1997]. This is completely unlike the constant high speed solar wind reported beyond 0.3 AU. Consequently, plumes and filamentary structure must be strongly mixed, and the mixing must be far along by 0.3 AU to be consistent with Helios observations. The paradox is what causes the mixing? Existing models of coronal heating and solar wind acceleration hardly address this issue. One possibility we are investigating is the MHD Kelvin-Helmholtz instability, to which the shear between plumes and interplume corona is expected to become unstable at 5-10 Solar radius. This instability can be simulated and followed far into the nonlinear regime and may lead to Alfvenic fluctuations like those seen at 1 AU.

Suess, Steven T.↗
NASA NTRSDate not supplied

Impact Of Rotation on the Sun’s Subsurface Structure

It is well known that the dynamical evolution of subsurface flows of the Sun and magnetic fields is strongly linked to rotation. Observational investigation of sub surface structuring and dynamics can be performed only with relatively low resolution using helioseismic inversion techniques. To confidently model the coupling of turbulent convection, magnetic fields, and rotation in the solar interior, we perform advanced radiative hydrodynamic and MHD modeling of local scales. To explore the coupling effects in the upper layers of the convection zone, we present results of analysis from a series of hydrodynamic simulations performed for a 25-Mm-deep computational domain at different latitudes obtained for 1) the Carrington rotation rate at all latitudes and 2) the realistic rotation rate from the global helioseismic inversions during the solar minimum. Using additional 3D radiative hydrodynamic models, we reveal how physical properties of the solar convection are affected by rotation.

Heliophysics↗
NASA NTRS2016

Impacts of Spontaneous Hot Flow Anomalies on the Magnetosheath and Magnetopause

Spacecraft observations and global hybrid (kinetic ions and fluid electrons) simulations have demonstrated that ion dissipation processes at the quasi-parallel bow shock are associated with the formation of structures called spontaneous hot flow anomalies (SHFAs). Previous simulations and recent spacecraft observations have also established that SHFAs result in the formation of magnetosheath filamentary structures(MFS). In this paper we demonstrate that in addition to MFS, SHFAs also result in the formation of magnetos heath cavities that are associated with decreases in density, velocity, and magnetic field and enhancements in temperature. We use the results of a global MHD run to determine the change in the magnetosheath properties associated with cavities due to ion kinetic effects. The results also show the formation of regions of high flow speed called magnetosheath jets whose properties as a function of solar wind Mach number are described in this study. Comparing the properties of the simulated magnetosheath cavities and jets to past spacecraft observations provides good agreement in both cases. We also demonstrate that pressure variations associated with cavities and SHFAs in the sheath result in a continuous sunward and anti sunward magnetopause motion. This result is consistent with previous suggestions that SHFAs may be responsible for the generation of ion cyclotron waves and precipitation of ring current protons in the outer magnetosphere.

Omidi, N.↗
NASA NTRS1994

Transition to whistler mediated magnetic reconnection

The transition in the magnetic reconnection rate from the resistive magnetohydrodynamic (MHD) regime where the Alfen wave controls reconnection to a regime in which the ions become unmagnetized and the whistler wave mediates reconnection is explored with 2-D hybrid simulations. In the whistler regime the electrons carry the currents while the ions provide a neutralizing background. A simple physical picture is presented illustrating the role of the whistler mediated reconnection is calculated analytically. The development of an out-of-plane component of the magnetic field is an observable signature of whistler driven reconnection.

Mandt, M. E.↗
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 NTRS1982

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

This paper describes a computer code for calculating the flow dynamics of constant density flow in the second stage trumpet shaped nozzle section of a two stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite difference computer code has been developed to allow the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for efficient solution of the equations. Turbulence simulation is by way of a two-equation Kappa-epsilon model. The code produces as output the flowfield map of the nondimensional stream function, axial, and swirl velocity. Good agreement was obtained between the theoretical predictions and the qualitative experimental results. The best seed injector location for uniform seed distribution at combustor exit is with injector located centrally on the combustor axis at entrance to the second stage combustor.

Gupta, A. K.↗
NASA NTRS1986

Filament formation due to photospheric shear

Recently, Wu et al. (1983, 1986) developed a three-dimensional, time-dependent axisymmetric ideal magnetohydrodynamic (MHD) model which has been used to study lower solar atmospheric response due to photospheric motion. They have shown that various observed solar features can be formed from these numerical simulation experiments. Some work currently in progress which has resulted from a specific numerical simulation experiment is reported. These results may elucidate the possible physical mechanism of formation of filament.

Wu, S. T.↗
NASA NTRS2004

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made substantial progress toward the goal of producing a unified model of the basic physical processes responsible for solar wind acceleration. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a 1D model of plasma heating and acceleration. The accomplishments in Year 2 are divided into these two categories: 1a. Focused Study of Kinetic Magnetohydrodynamic (MHD) Turbulence. lb. Focused Study of Non - WKB Alfven Wave Rejection. and 2. The Unified Model Code. We have continued the development of the computational model of a time-study open flux tube in the extended corona. The proton-electron Monte Carlo model is being tested, and collisionless wave-particle interactions are being included. In order to better understand how to easily incorporate various kinds of wave-particle processes into the code, the PI performed a detailed study of the so-called "Ito Calculus", i.e., the mathematical theory of how to update the positions of particles in a probabilistic manner when their motions are governed by diffusion in velocity space.

Cranmer, Steven R.↗
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 NTRS2007

Community Coordinated Modeling Center Support of Operations: Real-Time Simulations and V & V.

In support of Operations Community Coordinated Modeling Center (CCMC) performing validation and verification of space weather models. To identify suitable metrics the CCMC focus on parameters most useful to operations that CCMC resident models can provide. The real time simulations carried out at CCMC are an essential tool to test model performance and stability by using input conditions that may occur in nature at any time. Since 2001, the magnetospheric MHD model BATSRUS has been run in real time using ACE real time data. CCMC staff developed an experimental real-time system that controls uploading of the real-time ACE data, monitors continuous model execution, initiates automatic recovery procedure in case of data gaps or hardware failures, synchronizes BATSRUS and FRC runs, and periodically runs IDL based visualization software.

Kuznetsova, M.↗
NASA NTRS2010

Modeling of the Convection and Interaction of Ring Current, Plasmaspheric and Plasma Sheet Plasmas in the Inner Magnetosphere

Distinctive sources of ions reside in the plasmasphere, plasmasheet, and ring current regions at discrete energies constitute the major plasma populations in the inner/middle magnetosphere. They contribute to the electrodynamics of the ionosphere-magnetosphere system as important carriers of the global current system, in triggering; geomagnetic storm and substorms, as well as critical components of plasma instabilities such as reconnection and Kelvin-Helmholtz instability at the magnetospheric boundaries. Our preliminary analysis of in-situ measurements shoves the complexity of the plasmas pitch angle distributions at particularly the cold and warm plasmas, vary dramatically at different local times and radial distances from the Earth in response to changes in solar wind condition and Dst index. Using an MHD-ring current coupled code, we model the convection and interaction of cold, warm and energetic ions of plasmaspheric, plasmasheet, and ring current origins in the inner magnetosphere. We compare our simulation results with in-situ and remotely sensed measurements from recent instrumentation on Geotail, Cluster, THEMIS, and TWINS spacecraft.

Fok, Mei-Ching↗
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.↗