SEARCH · Engineering Papers
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.
Medium- and heavy-duty electric truck charging simulation results for the Houston and Dallas megaregion along the I-45 corridor under 96 scenarios
Explore the source record for details and available documents.
ABSTRACT In this work, we find empirical evidence that the scale-dependent statistical properties of solar wind and magnetohydrodynamic (MHD) turbulence can be described in terms of a family of parametric probability distribution functions (PDFs) known as Normal Inverse Gaussian (NIG). Understanding these PDFs is one of the most important goals in turbulence theory, as they are inherently connected to the intermittent properties of solar wind turbulence. We investigate the properties of PDFs of Elsasser increments based on a large statistical sample from solar wind observations and high-resolution numerical simulations of MHD turbulence. In order to measure the PDFs and their corresponding properties, three experiments are presented: fast and slow solar wind for experimental data and a simulation of reduced MHD (RMHD) turbulence. Conditional statistics on a 23-yr-long sample of WIND data near 1 au and high-resolution pseudo-spectral simulation of steadily driven RMHD turbulence on a $2048^3$ mesh are used to construct scale-dependent PDFs. The empirical PDFs are fitted to NIG distributions, which depend on four free parameters. Our analysis shows that NIG distributions accurately capture the evolution of the PDFs, with scale-dependent parameters, from large scales characterized by a Gaussian distribution, turning to exponential tails within the inertial range and stretched exponentials at dissipative scales. We also show that empirically-measured NIG parameters exhibit well-defined scaling properties that are similar across the three empirical data sets, which may be indicative of universal behaviour.
On May 23, 1995, the Comprehensive Plasma Instrumentation (CPI) onboard the Geotail spacecraft observed a complex and structured ion distribution function near the magnetotail midplane at x approximately -10 R(sub E). On the same day, the Wind spacecraft observed a very high density (approximately 40/cubic cm) solar wind and an interplanetary magnetic field (IMF) that was predominantly northward but had several southward turnings. We have inferred the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time using time-dependent electric and magnetic fields obtained from a global MHD (magnetohydrodynamic) simulation. Wind data were used as input for the MHD model. We found that three sources contributed to this distribution: the ionosphere, the plasma mantle which had near-Earth and distant tail components, and the low latitude boundary layer (LLBL). Moreover, distinct structures in the low energy part of the distribution function were found to be associated with individual sources. Structures near 0 deg pitch angle were made up of either ionospheric or plasma mantle ions, while structures near 90 deg pitch angle were dominated by ions from the LLBL source. Particles that underwent nonadiabatic acceleration were numerous in the higher energy part of the ion distribution function, whereas ionospheric and LLBL ions were mostly adiabatic. A large proportion of the near-Earth mantle ions underwent adiabatic acceleration, while most of the distant mantle ions experienced nonadiabatic acceleration.
We report here the observations of a plasma jet evolving through a macro- to micro-scale progression sequence. This leads to a fast magnetic reconnection that results in the jet breaking off from its originating electrode and forming a force-free state. A sausage-like pinching occurs first and squeezes an initially fat, short magnetized jet so that it becomes thin. The thin jet then becomes kink unstable. The lengthening of the jet by the kinking thins the jet even more since the kink is an incompressible instability. When the jet radius becomes comparable to the ion-skin depth, Hall and electron inertial physics become important and establish the environment for a fast magnetic reconnection. This fast reconnection occurs, disrupting the jet and establishing a force-free state. X-ray bursts and whistler waves, evidence of a magnetic reconnection, are observed when the plasma jet breaks off from the electrode. This experimentally confirmed sequence of successive thinning from pinching followed by kinking is reproduced in a three-dimensional ideal Magnetohydrodynamic (MHD) numerical simulation. The results of the experiment and the numerical simulation, together demonstrate a viable path from macro-scale MHD physics to micro-scale non-MHD physics where fast reconnection occurs.
This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.
ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.
Abstract Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at 2048 3 cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.
A two-dimensional MHD model simulating the formation of Kippenhahn-Schluter (1957) quiescent prominence (QP) magnetic field configurations is used to explore symmetric mass injection into a dipole magnetic field. An optimum magnetic field strength for QP formation by mass injection is obtained. It is found that a weaker magnetic field strength is more favorable for the condensation of the injected plasma but that a stronger field is more favorable for supporting the condensed plasma against gravity.
Recent developments in the field of numerical simulation models for the study of shock wave propagation in the corona are presented. These models are based on gasdynamic (GD) and ideal (that is, dissipationless, except at shocks) magnetohydrodynamic (MHD) theories. The characteristics and physical interpretations of the results derived from these models are discussed in some detail. The most significant physical results obtained to date are provided by the two-dimensional non-planar, time-dependent, MHD numerical simulation model. In this model, the non-linear interaction among the three essential MHD waves, i.e., fast-, slow-, and Alfven waves are demonstrated. Finally, the physical relevance of these numerical simulation models in relation to observed solar activity is presented.
We present numerical results from long-term CPU and GPU simulations of rotating, homogeneous, magnetohydrodynamic (MHD) turbulence, and discuss their connection to the spherically bounded case. We compare our numerical results with a statistical theory of geodynamo action that has evolved from the absolute equilibrium ensemble theory of ideal MHD turbulence, which is based on the ideal MHD invariants are energy, cross helicity and magnetic helicity. However, for rotating MHD turbulence, the cross helicity is no longer an exact invariant, although rms cross helicity becomes quasistationary during an ideal MHD simulation. This and the anisotropy imposed by rotation suggests an ansatz in which an effective, nonzero value of cross helicity is assigned to axisymmetric modes and zero cross helicity to non-axisymmetric modes. This hybrid statistics predicts a large-scale quasistationary magnetic field due to broken ergodicity , as well as dipole vector alignment with the rotation axis, both of which are observed numerically. We find that only a relatively small value of effective cross helicity leads to the prediction of a dipole moment vector that is closely aligned (less than 10 degrees) with the rotation axis. We also discuss the effect of initial conditions, dissipation and grid size on the numerical simulations and statistical theory.
We present the numerical simulations of winds from evolved giant stars using a fully non-linear, time dependent 2.5-dimensional magnetohydrodynamic (MHD) code. This study extends our previous fully non-linear MHD wind simulations to include a broadband frequency spectrum of Alfven waves that drive winds from red giant stars. We calculated four Alfven wind models that cover the whole range of Alfven wave frequency spectrum to characterize the role of freely propagated and reflected Alfven waves in the gravitationally stratified atmosphere of a late-type giant star. Our simulations demonstrate that, unlike linear Alfven wave-driven wind models, a stellar wind model based on plasma acceleration due to broadband non-linear Alfven waves, can consistently reproduce the wide range of observed radial velocity profiles of the winds, their terminal velocities and the observed mass loss rates. Comparison of the calculated mass loss rates with the empirically determined mass loss rate for alpha Tau suggests an anisotropic and time-dependent nature of stellar winds from evolved giants.
We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.
Global MHD (magnetohydrodynamic) simulations have shown a remarkable ability to describe the global dynamics of geospace. The limitations of the physical approximations underlying MHD would seem to limit the effectiveness of these codes, since kinetic and hybrid effects should manifest themselves by cross-scale coupling from microscales to mesoscales to global scales. However three effects appear to allow the codes to operate much more successfully than one would at first believe. They are:(l) the globally self-consistent nature of the codes with very well defined exterior boundary conditions (the solar wind) which allows the proper intercommunication between magnetospheric regions on MHD scales, (2) the control by global dynamics of the boundary layer locations where micro and meso scale processes operate, and (3) the critical role of numerical diffusion and with a sufficiently high resolution grid, the use of an empirical resistivity term, which if set at a level where the major magnetosphere boundaries properly calibrate against their observed locations, appear to well represent the effects of kinetic and hybrid processes on the global dynamics. The effectiveness of the global MHD codes, which have been developed under the ISTP mission, in describing Wind, Polar and Geotail observations, as well as ground-based observations are described. Particular emphasis is placed upon the Polar imaging data which when combined with ground-based data and global MHD-based synthetic aurora and convection patterns provide a powerful tool in understanding the final link in the solar-terrestrial chain: coupling into the atmosphere and ionosphere.
Calculations of the variations of interplanetary scintillation (IPS) from a disturbance simulated by a 3-D magnetohydrodynamical (MHD) model of the solar wind are presented. The simulated maps are compared with observations and it is found that the MHD model reproduces the qualitative features of observed disturbances. The disturbance produced by the MHD simulation is found to correspond in strength with the weakest disturbance which can be reliably detected by existing single station IPS observations.
The effect of adding a By magnetic component to a solar wind already carrying a Bz magnetic component is studied using three-dimensional MHD computer simulations. It is found that the tail structure and the Bx lobe fields are rotated and become asymmetric. In addition, a By component is generated in the neutral sheet. These results are found to correlate with observational data.
The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.
Analysis is provided in this report of using two head-on magnetohydrodynamic (MHD) shocks to achieve supercritical nuclear fission in an axially elongated cylinder filled with UF4 gas as an energy source for deep space missions. The motivation for each aspect of the design is explained and supported by theory and numerical simulations. A subsequent report will provide detail on relevant experimental work to validate the concept. Here the focus is on the theory of and simulations for the proposed gas core reactor conceptual design from the onset of shock generations to the supercritical state achieved when the shocks collide. The MHD model is coupled to a standard nuclear code (MCNP) to observe the neutron flux and fission power attributed to the supercritical state brought about by the shock collisions. Throughout the modeling, realistic parameters are used for the initial ambient gaseous state and currents to ensure a resulting supercritical state upon shock collisions.