Search NASA⌕ Search

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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Models for application of radiation boundary condition for MHD waves in collapse calculations

The problem of reflection of magnetohydrodynamic (MHD) waves at the boundary of a numerical grid has to be resolved in order to obtain reliable results for the end state of the (isothermal) collapse of a rotating, magnetic protostellar cloud. Since the goal of investigating magnetic braking in collapse simulations is to see if the transport of angular momentum via alfven waves is large enough to solve the angular momentum problem an approximation that artificially suppresses large amplitudes in the MHD waves can be self-defeating. For this reason, four alternate methods of handling reflected waves where no assumptions are made regarding the amplitudes of the waves were investigated. In order to study this problem (of reflection) without interference from other effects these methods were tried on two simpler cases. The four methods are discussed.

Vanajakshi, C. T.↗
NASA NTRS1988

Dynamical evolution of a solar coronal magnetic field arcade

Calculations of the long-term dynamical evolution of a solar coronal magnetic field arcade which is subjected to shearing photospheric flows are presented. The evolution is obtained by numerical solution of a subset of the resistive magnetohydrodynamic equations. For a simplified model of the bipolar magnetic field observed in the solar corona, it is found that photospheric flow produces a slow evolution of the magnetic field, with a buildup of magnetic energy. For certain photospheric shear profiles, the field configuration produced is linearly unstable to an ideal magnetohydrodynamic mode when the shear exceeds a critical value. The nonlinear evolution of this instability shows the spontaneous formation of current sheets. Reconnection of the magnetic field produces a rapid release of magnetic energy. The major fraction of the energy is dissipated resistively, while a small fraction is converted into kinetic energy of an ejected plasmoid. The relevance of these results to two-ribbon flares is discussed.

Mikic, Z.↗
NASA NTRS1991

Numerical simulation of self-field MPD thrusters

A fully two dimensional magnetohydrodynamics code was developed to predict self-field, steady-state MPD thruster performance. The governing equations are outlined, and methods of solution are presented. Model predictions are compared with experimental data for two thruster geometries over a range of discharge currents and mass flow rates. Model limitations are evaluated, and issues concerning quasi-steady versus steady-state thruster comparisons are discussed.

Lapointe, Michael R.↗
NASA NTRS2000

Residual Distribution Schemes for Conservation Laws Via Adaptive Quadrature

This paper considers a family of nonconservative numerical discretizations for conservation laws which retains the correct weak solution behavior in the limit of mesh refinement whenever sufficient order numerical quadrature is used. Our analysis of 2-D discretizations in nonconservative form follows the 1-D analysis of Hou and Le Floch. For a specific family of nonconservative discretizations, it is shown under mild assumptions that the error arising from non-conservation is strictly smaller than the discretization error in the scheme. In the limit of mesh refinement under the same assumptions, solutions are shown to satisfy an entropy inequality. Using results from this analysis, a variant of the "N" (Narrow) residual distribution scheme of van der Weide and Deconinck is developed for first-order systems of conservation laws. The modified form of the N-scheme supplants the usual exact single-state mean-value linearization of flux divergence, typically used for the Euler equations of gasdynamics, by an equivalent integral form on simplex interiors. This integral form is then numerically approximated using an adaptive quadrature procedure. This renders the scheme nonconservative in the sense described earlier so that correct weak solutions are still obtained in the limit of mesh refinement. Consequently, we then show that the modified form of the N-scheme can be easily applied to general (non-simplicial) element shapes and general systems of first-order conservation laws equipped with an entropy inequality where exact mean-value linearization of the flux divergence is not readily obtained, e.g. magnetohydrodynamics, the Euler equations with certain forms of chemistry, etc. Numerical examples of subsonic, transonic and supersonic flows containing discontinuities together with multi-level mesh refinement are provided to verify the analysis.

Barth, Timothy↗
NASA NTRS1989

Intensive MHD-structures penetration in the middle atmosphere initiated in the ionospheric cusp under quiet geomagnetic conditions

In connection with the recently detected quasiperiodical magnetic disturbances in the ionospheric cusp, the penetration of compressional surface magnetohydrodynamic (MHD) waves through the middle atmosphere is modelled numerically. For the COSPAR International Reference Atmosphere (CIRA) 72 model the respective energy density flux of the disturbances in the middle atmosphere is determined. On the basis of the developed model certain conclusions are reached about the height distribution of the structures (energy losses, currents, etc.) initiated by intensive magnetic cusp disturbances.

Mateev, L. N.↗
NASA NTRS1990

Magnetic reconnection driven by current repulsion

The evolution of an equilibrium consisting of two magnetic islands with oppositely directed currents embedded in a strong magnetic field is investigated, using numerical simulation methods. The rapid development of an ideal magnetohydrodynamic instability is observed, which first rotates and then expels the islands. The growth rate is on the order of the inverse of the Alfven transit time and is much higher than that for magnetic island coalescence. In the nonlinear stage, resistivity becomes important as the reconnection process ensues and dissipates the magnetic energy. The growth rate of the instability is a weak function of the plasma beta and other plasma parameters such as S, the magnetic Reynolds number. An energy principle analysis, based on eigenfunctions obtained from the simulation, confirms the existence of the instability.

Richard, R. L.↗
NASA NTRS1996

MHD Modeling of the Interaction of the Solar Wind With Venus

The primary objective of this research program is to improve our understanding of the physical processes occurring in the interaction of the solar wind with Venus. This will be accomplished through the use of numerical solutions of the two- and three-dimensional magnetohydrodynamic (MHD) equations and through comparisons of the computed results with available observations. A large portion of this effort involves the study of processes due to the presence of the magnetic field and the effects of mass loading. Published papers are included in the appendix.

Steinolfson, R. S.↗
NASA NTRS2002

Extension of Low Dissipative High Order Hydrodynamics Schemes for MHD Equations

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD (magnetohydrodynamic) equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls to limit the amount and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvi-linear grids. The three features of the present MHD scheme over existing schemes in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion magnetized flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike existing schemes for the conservative MHD equations which suffer from ill-conditioned eigen-decompositions, the present scheme makes use of a well-conditioned eigen-decomposition to solve the conservative form of the MHD equations. This is due to, partly. the fact that the divergence of the magnetic field condition is a different type of constraint from its incompressible Navier-Stokes cousin. Third, a new approach to minimize the numerical error of the divergence free magnetic condition for high order scheme is introduced.

Yee, H. C.↗
NASA NTRS2003

Pulsed Inductive Thruster (PIT): Modeling and Validation Using the MACH2 Code

Numerical modeling of the Pulsed Inductive Thruster exercising the magnetohydrodynamics code, MACH2 aims to provide bilateral validation of the thruster's measured performance and the code's capability of capturing the pertinent physical processes. Computed impulse values for helium and argon propellants demonstrate excellent correlation to the experimental data for a range of energy levels and propellant-mass values. The effects of the vacuum tank wall and massinjection scheme were investigated to show trivial changes in the overall performance. An idealized model for these energy levels and propellants deduces that the energy expended to the internal energy modes and plasma dissipation processes is independent of the propellant type, mass, and energy level.

Schneider, Steven↗
NASA NTRS1991

Numerical simulation of self-field MPD thrusters

A fully two-dimensional magnetohydrodynamics code has been developed to predict self-field, steady-state MPD thruster performance. The governing equations are outlined, and methods of solution are presented. Model predictions are compared with experimental data for two thruster geometries over a range of discharge currents and mass flow rates. Model limitations are evaluated, and issues concerning quasi-steady versus steady-state thruster comparisons are discussed.

La Pointe, Michael R.↗
NASA NTRS2001

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.↗
NASA NTRS1992

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.↗
NASA NTRS1993

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.↗
NASA NTRS2023

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↗
NASA NTRS1994

Anisotropic Alfven-ballooning modes in Earth's magnetosphere

We have carried out a theoretical analysis of the stability and parallel structure of coupled shear Alfven and slow magnetosonic waves in Earth's inner magnetopause (i.e., at equatorial distances between about five and ten Earth radii) including effects of finite anisotropic Grad-Shafranov equation yields an approximate self-consistent magnetohydrodynamic (MHD) equilibrium. This MHD equilibrium is used in the numerical solution of a set of eigenmode equations which describe the field line eigenfrequency, linear stability, and parallel eigenmode structure. We call these modes anisotropic Alfven-ballooning modes. The main results are: (1) The field line eigenfrequency can be significantly lowered by finite pressure effects. (2) The parallel mode structure of the transverse wave components is fairly insensitive to changes in the plasma pressure, but the compressional magnetic component can become highly peaked near the magnetic equator as a result of increased pressure, especially when P(sub perpendicular to) is greater than P(sub parallel) (here P(sub perpendicular to) and P(sub parallel) are the perpendicular and parallel plasma pressure). (3) For the isotropic (P(sub parallel) = P(sub perpendicular to) = P) case ballooning instability can occur when the ratio of the plasma presure to the magnetic pressure, beta = P/(B squared/8 pi), exceeds a critical value beta(sup B)(sub 0) is approximately equal to 3.5 at the equator. (4) Compared to the isotropic case the critical beta value is lowered by anisotropy, either due to decreased field line bending stabilization when P(sub parallel) is greater than P(sub perpendicular to) or due to increased ballooning-mirror destabilization when P(sub perpendicular to) is greater than P(sub parallel). (5) We use a beta-delta stability diagram to display the regions of instability with respect to the equatorial values of the parameters bar beta and delta, where bar beta = (1/3)(beta(sub parallel) + 2 beta(sub perpendicular to)) is an average beta value and delta = 1 - P(sub parallel)/P(sub perpendicular to) is a measure of the plasma anisotropy. The diagram is divided into regions corresponding to the firehose, mirror and ballooning instabilities. It appears that observed values of the plasma pressure are below the critical value for the isotropic ballooning instability but it may be possible to approach a ballooning-mirror instability when P(sub perpendicular to)/P(sub parallel) is greater than or approximately 2.

Chan, Anthony A.↗