Search NASA⌕ Search

SEARCH · Search NASA

Results for “magnetohydrodynamical simulations”

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 307 records · Page 17

General Relativistic Simulations of Magnetized Plasmas Around Merging Supermassive Black Holes

Coalescing supermassive black hole binaries are produced by the mergers of galaxies and are the most powerful sources of gravitational waves accessible to space-based gravitational observatories. Some such mergers may occur in the presence of matter and magnetic fields and hence generate an electromagnetic counterpart. In this paper we present the first general relativistic simulations of magnetized plasma around merging supermassive black holes using the general relativistic magnetohydrodynamic code Whisky. By considering different magnetic field strengths, going from non-magnetically dominated to magnetically dominated regimes, we explore how magnetic fields affect the dynamics of the plasma and the possible emission of electromagnetic signals. In particular we observe, total amplification of the magnetic field of approx 2 orders of magnitude which is driven by the accretion onto the binary and that leads to stronger electromagnetic signals than in the force-free regime where such amplifications are not possible.

Giacomazzo, Bruno↗

Initial Results from the GEM Challenge on the Spacecraft Surface Charging Environment

Spacecraft surface charging during geomagnetically disturbed times is one of the most important causes of satellite anomalies. Predicting the surface charging environment is one prevalent task of the geospace environment models. Therefore, the Geospace Environment Modeling (GEM) Focus Group "Inner Magnetosphere Crossenergy/Population Interactions" initiated a communitywide challenge study to assess the capability of several inner magnetosphere ring current models in determining surface charging environment for the Van Allen Probes orbits during the 17 March 2013 storm event. The integrated electron flux between 10 and 50 keV is used as the metrics. Various skill scores are applied to quantitatively measure the modeling performance against observations. Results indicate that no model consistently perform the best in all of the skill scores or for both satellites. We find that from these simulations the ring current model with observational flux boundary condition and Weimer electric potential driver generally reproduces the most realistic flux level around the spacecraft. A simple and weaker VollandStern electric field is not capable of effectively transporting the same plasma at the boundary toward the Earth. On the other hand, if the ring current model solves the electric field selfconsistently and obtains similar strength and pattern in the equatorial plane as the Weimer model, the boundary condition plays another crucial role in determining the electron flux level in the inner region. When the boundary flux spectra based on magnetohydrodynamics (MHD) model/empirical model deviate from the shape or magnitude of the observed distribution function, the simulation produces poor skill scores along Van Allen Probes orbits.

Yu, Yiqun↗

Research and advanced concepts

Gaseous breakdown of simulated planetary atmospheres and illustration of Greenhouse effect, thermal radiation from ionized argon, liquid magnetohydrodynamic power conversion

GASEOUS DIFFUSION↗

LLNL contributions to MPD thrusters for SEI

Some of the topics covered with respect to the Lawrence Livermore National Laboratory's (LLNL's) contributions to Magnetoplasmadynamic (MPD) Thrusters for the Space Exploration Initiative (SEI) include: an IR camera, plasma-induced erosion/redeposition, the Mirror Fusion Test Facility-B (MFTF-B), the Thruster Lifetime Test Facility, the RACE Compact Torus Accelerator Facility, and a RACE program summary. Some of the other topics addressed include: flux contours for HAM simulation, comparison of RACE data of plasma ring formation with the HAM 2-D magnetohydrodynamic code, and the 2-D Ring Acceleration Code (TRAC).

Hooper, Edwin Bickford↗

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2-D, two-temperature, single fluid magnetohydrodynamic code which incorporates classical plasma transport coefficients and Hall effects was developed to predict steady-state, self-field magnetoplasmadynamic (MPD) thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects, the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries, assuming a fully ionized argon propellant. A maximum predicted specific impulse of 1680 s is obtained for a thruster with an anode radius of 2.5 cm, a cathode radius of 0.5 cm, and equal electrode lengths of 2.5 cm. A scaling relation is developed to predict, within limits, the onset of cylindrical, self-field thruster instability as a function of geometry and operating condition.

Lapointe, M.↗

Magnetic effects in Venus/solar wind interaction

The overall objective of this research program is to better understand the interaction of a magnetized solar wind with the Venus atmosphere through the use of numerical solutions of the time-dependent, 2-D and 3-D magnetohydrodynamic (MHD) equations. Due to the more modest CPU requirements for the 2-D simulations, they were used for studies in which useful information not dependent on the third dimension could be obtained. The 2-D simulations served several purposes in addition to providing useful physical insight. They were used to determine the numerical parameters required in the 3-D studies, such as the grid spacing required to resolve particular features, and the damping that must be included to remove high-frequency oscillations. Among the specific studies performed with support from this grant that are discussed in this report are the following: comparison with other available models for purposes of testing the code and obtaining a baseline with which to evaluate the effects of additional physical processes, effects of a finite planet conductivity, bow shock standoff distance, and formation of the magnetic barrier and slippage of the magnetic field around the planet. A brief description of the methodology is presented before discussing the results.

Steinolfson, Richard S.↗

Magnetic Damping of g-Jitter Induced Double-Diffusive Convection

This paper describes a numerical study of the g-jitter driven double diffusive convective flows, thermal and concentration distributions in binary alloy melt systems subject to an external magnetic field. The study is based on the finite element solution of transient magnetohydrodynamic equations governing the momentum, thermal and solutal transport in the melt pool. Numerical simulations are conducted using the synthesized single- and multi- frequency g-jitter as well as the real g-jitter data taken during space flights with or without an applied magnetic field. It is found that for the conditions studied, the main melt flow follows approximately a lineal- superposition of velocity components induced by individual g-jitter components, regardless of whether a magnetic field exists or not. The flow field is characterized by a recirculating double diffusive convection loop oscillating in time with a defined frequency equal to that of the driving g-jitter force. An applied magnetic field has little effect on the oscillating recirculating pattern, except around the moment in time when the flow reverses its direction. The field has no effect on the oscillation period, but it changes the phase angle. It is very effective in suppressing the flow intensity and produces a notable reduction of the solutal striation and time fluctuations in the melt. For a given magnetic field strength, the magnetic damping effect is more pronounced on the velocity associated with the largest g-jitter component present and/or the g-jitter spiking peaks. A stronger magnetic field is more effective in suppressing the melt convection and also is more helpful in bringing the convection in phase with the g-jitter driving force. The applied field is particularly useful in suppressing the effect of real g-jitter spikes on both flow and solutal distributions. With appropriately selected magnetic fields, the convective flows caused by g-jitter can be reduced sufficiently and diffusion dominant. solutal transport in the melt is possible.

Shu, Y.↗

Annual Research Briefs, 2004: Center for Turbulence Research

This report contains the 2004 annual progress reports of the Research Fellows and students of the Center for Turbulence Research in its eighteenth year of operation. Since its inception in 1987, the objective of the CTR has been to advance the physical understanding of turbulent flows and development of physics based predictive tools for engineering analysis and turbulence control. Turbulence is ubiquitous in nature and in engineering devices. The studies at CTR have been motivated by applications where turbulence effects are significant; these include a broad range of technical areas such as planetary boundary layers, formation of planets, solar convection, magnetohydrodynamics, environmental and eco systems, aerodynamic noise, propulsion systems and high speed transportation. Numerical simulation has been the predominant research tool at CTR which has required a critical mass of researchers in numerical analysis and computer science in addition to core disciplines such as applied mathematics, chemical kinetics and fluid mechanics. Maintaining and promoting this interdisciplinary culture has been a hallmark of CTR and has been responsible for the realization of the results of its basic research in applications. The first group of reports in this volume are directed towards development, analysis and application of novel numerical methods for ow simulations. Development of methods for large eddy simulation of complex flows has been a central theme in this group. The second group is concerned with turbulent combustion, scalar transport and multi-phase ows. The nal group is devoted to geophysical turbulence where the problem of solar convection has been a new focus of considerable attention recently at CTR.

Moin, Parviz↗

Development of Numerical Tools for the Investigation of Plasma Detachment from Magnetic Nozzles

A multidimensional numerical simulation framework aimed at investigating the process of plasma detachment from a magnetic nozzle is introduced. An existing numerical code based on a magnetohydrodynamic formulation of the plasma flow equations that accounts for various dispersive and dissipative processes in plasmas was significantly enhanced to allow for the modeling of axisymmetric domains containing three.dimensiunai momentum and magnetic flux vectors. A separate magnetostatic solver was used to simulate the applied magnetic field topologies found in various nozzle experiments. Numerical results from a magnetic diffusion test problem in which all three components of the magnetic field were present exhibit excellent quantitative agreement with the analytical solution, and the lack of numerical instabilities due to fluctuations in the value of del(raised dot)B indicate that the conservative MHD framework with dissipative effects is well-suited for multi-dimensional analysis of magnetic nozzles. Further studies will focus on modeling literature experiments both for the purpose of code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗

Hybrid-Vlasov Simulation of Soft X-Ray Emissions at the Earth’S Dayside Magnetospheric Boundaries

Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause, the magnetosheath and the polar cusps by remote sensing techniques. The Solar wind–Magnetosphere–Ionosphere Link Explorer (SMILE) and Lunar Environment heliospheric X-ray Imager (LEXI) missions aim to obtain soft X-ray images of near-Earth space thanks to their Soft X-ray Imager (SXI) instruments. While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission, the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma. To investigate the possible signatures of ion-kinetic-scale processes in soft X-ray images, we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model. The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field. We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations, providing meridional and equatorial views. We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirror mode wave structures in the magnetosheath and flux transfer events (FTEs) at the magnetopause. Our results suggest that, although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers, mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images. For instance, a local increase by 30% in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12% enhancement in the line-of-sight- and time integrated soft X-ray emissivity originating from this region. Likewise, a proton density increase by 14% in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%. These are likely conservative estimates, given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind. These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale, transient phenomena occurring on the dayside.

magnetosphere↗

On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby

The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.

Kabin, K.↗

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

In this study, we simulate the Solar Wind Charge Exchange (SWCX) soft X-ray emissions at dayside magnetosheath and cusps by using magnetohydrodynamic (MHD) and LAtmos TEst Particle (LaTeP) models. MHD models are unable to resolve the particle kinetic effects, such as the different behaviors of ions with different q/m, or distinguish the magnetospheric plasma from the solar wind plasma. We investigate these effects with the LaTeP model. As the LaTeP model does not self-compute magnetic and electric field, the magnetic and electric field data obtained from Open Geospace General Circulation Model (OpenGGCM) and Lagrangian version of the piecewise parabolic method (PPMLR) MHD model are used as the input to LaTeP model. The soft X-ray emissivity maps simulated from pure OpenGGCM and PPMLR MHD approaches and from LaTeP-OpenGGCM and LaTeP-PPMLR approaches are presented and compared. The results indicate that the LaTeP model can well resolve the kinetic effects and can be used to investigate the individual spectral characteristics. Therefore, the LaTeP model is a complementary approach for simulating the X-ray emissions near the dayside magnetopause. We also calculate the ratio of integrated OVII/OVIII line intensities, produced by charge exchange of O7+ ions and O8+ ions, respectively. We find a relatively higher ratio at the bow shock compared to the surrounding areas, suggesting that this ratio can be an effective parameter to identify the bow shock location.

Qiuyu Xu↗

Collisionless Shock Waves in the Solar Terrestrial Environment

Natural shocks, magnetohydrodynamics, shock extensions, global asymmetries, shock profiles, quasiperpendicular supercritical process; ion dynamics and field structure, numerical simulations, ion reflection machinery, electrons, instabilities, phenomenology, and ion acceleration were described.

Source record↗

Self-excitation of auroral arcs in a three-dimensionally coupled magnetosphere-ionosphere system

This paper presents the first full three-dimensional dynamic simulation of auroral arc formation. The magnetospheric and ionospheric dynamics are represented by one-fluid magnetohydrodynamic equations and two-fluid weakly ionized plasma equations, respectively. The feedback coupling between magnetospheric Alfven waves and ionospheric density waves results in a spontaneous generation of longitudinally elongated striations of field-aligned currents and ionospheric electron densities, which compare very well with many features of quiet auroral arcs.

Watanabe, Kunihiko↗

Multispacecraft Observations and Modeling of the 22/23 June 2015 Geomagnetic Storm

The magnetic storm of 22-23 June 2015 was one of the largest in the current solar cycle. We present in situ observations from the Magnetospheric Multiscale Mission (MMS) and the Van Allen Probes (VAP) in the magnetotail, field-aligned currents from AMPERE (Active Magnetosphere and Planetary Electrodynamics Response), and ionospheric flow data from Defense Meteorological Satellite Program (DMSP). Our real-time space weather alert system sent out a "red alert," correctly predicting Kp indices greater than 8. We show strong outflow of ionospheric oxygen, dipolarizations in the MMS magnetometer data, and dropouts in the particle fluxes seen by the MMS Fast Plasma Instrument suite. At ionospheric altitudes, the AMPERE data show highly variable currents exceeding 20 MA. We present numerical simulations with the Block Adaptive Tree-Solarwind - Roe - Upwind Scheme (BATS-R-US) global magnetohydrodynamic model linked with the Rice Convection Model. The model predicted the magnitude of the dipolarizations, and varying polar cap convection patterns, which were confirmed by DMSP measurements.

Reiff, P. H.↗

Global simulation of the time-dependent magnetosphere

The paper presents preliminary results from time-dependent two-dimensional numerical modelling of the magnetohydrodynamic interaction of the solar wind with the magnetosphere. A southward solar wind-field produces a magnetospheric topology consistent with Dungey's (1961) model. The interaction appears to be fundamentally unsteady; the shock, magnetosheath, and magnetopause are highly turbulent. A 'substorm' is modelled as the passage of a rotational discontinuity over the magnetosphere; the onset of enhanced reconnection in the magnetospheric tail produces a closed magnetic island which convects downstream.

Leboeuf, J. N.↗

On the Existence of Fast Modes in Compressible Magnetohydrodynamic Turbulence

We study the existence and properties of fast magnetosonic modes in 3D compressible MHD turbulence by carrying out a number of simulations with compressible and incompressible driving conditions. We use two approaches to determine the presence of fast modes: mode decomposition based on spatial variations only and spatio-temporal 4D fast Fourier transform (4D FFT) analysis of all fluctuations. The latter method enables us to quantify fluctuations that satisfy the dispersion relation of fast modes with finite frequency. Overall, we find that the fraction of fast modes identified via the spatio-temporal 4D FFT approach in total fluctuation power is either tiny with nearly incompressible driving or ∼2% with highly compressible driving. We discuss the implications of our results for understanding the compressible fluctuations in space and astrophysical plasmas.

Zhaoming Gan↗

Numerical Simulation of Turbulent MHD Flows Using an Iterative PNS Algorithm

A new parabolized Navier-Stokes (PNS) algorithm has been developed to efficiently compute magnetohydrodynamic (MHD) flows in the low magnetic Reynolds number regime. In this regime, the electrical conductivity is low and the induced magnetic field is negligible compared to the applied magnetic field. The MHD effects are modeled by introducing source terms into the PNS equation which can then be solved in a very efficient manner. To account for upstream (elliptic) effects, the flowfields are computed using multiple streamwise sweeps with an iterated PNS algorithm. Turbulence has been included by modifying the Baldwin-Lomax turbulence model to account for MHD effects. The new algorithm has been used to compute both laminar and turbulent, supersonic, MHD flows over flat plates and supersonic viscous flows in a rectangular MHD accelerator. The present results are in excellent agreement with previous complete Navier-Stokes calculations.

Kato, Hiromasa↗