Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field theoretic 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 73 records · Page 4

Homogeneous buoyancy-generated turbulence

Using a theoretical analysis of fundamental equations and a numerical simulation of the flow field, the statistically homogeneous motion that is generated by buoyancy forces after the creation of homogeneous random fluctuations in the density of infinite fluid at an initial instant is examined. It is shown that analytical results together with numerical results provide a comprehensive description of the 'birth, life, and death' of buoyancy-generated turbulence. Results of numerical simulations yielded the mean-square density mean-square velocity fluctuations and the associated spectra as functions of time for various initial conditions, and the time required for the mean-square density fluctuation to fall to a specified small value was estimated.

Batchelor, G. K.↗

Statistics Analysis of the Uncertainties in Cloud Optical Depth Retrievals Caused by Three-Dimensional Radiative Effects

This paper presents a simple approach to estimate the uncertainties that arise in satellite retrievals of cloud optical depth when the retrievals use one-dimensional radiative transfer theory for heterogeneous clouds that have variations in all three dimensions. For the first time, preliminary error bounds are set to estimate the uncertainty of cloud optical depth retrievals. These estimates can help us better understand the nature of uncertainties that three-dimensional effects can introduce into retrievals of this important product of the MODIS instrument. The probability distribution of resulting retrieval errors is examined through theoretical simulations of shortwave cloud reflection for a wide variety of cloud fields. The results are used to illustrate how retrieval uncertainties change with observable and known parameters, such as solar elevation or cloud brightness. Furthermore, the results indicate that a tendency observed in an earlier study, clouds appearing thicker for oblique sun, is indeed caused by three-dimensional radiative effects.

Varnai, Tamas↗

Hybrid Simulations of FRC Merging and Compression

An improved understanding of field-reversed configuration (FRC) merging and stability in high acceleration and compression magnetic fields is needed to speed up the development of the pulsed fusion concept developed at Helion Energy. All previous theoretical and simulation work on FRC merging and compression was performed using two-dimensional (2D) magnetohydrodynamic (MHD) models. The results of novel 2D hybrid simulations (fluid electrons and full-orbit kinetic ions) of FRC merging and compression are presented. Results of kinetic and MHD simulations, computed using the HYM code, are compared and analyzed. In cases without axial magnetic compression, both the MHD and hybrid simulations show a high sensitivity to the initial parameters (i.e. FRC separation, velocity, normalized separatrix radius, and plasma viscosity), showing that FRCs with large elongation and separatrix radius either do not merge or merge partially, forming a doublet FRC. In conclusion, application of a mirror coil field at the FRC ends with increasing strength is shown to lead to fast and complete merging of the FRCs in MHD and kinetic simulations.

FRC↗

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.↗

Mixing, Noise and Thrust Benefits Using Corrugated Designs

These projects are directed toward the analysis of several concepts for nozzle and inlet performance improvement and noise reduction from jet exhausts. Currently. The FM&AL also initiates new joint research between the HU/FM&AL, the Hyper-X Program Team at the LaRC, and the Central Institute of Aviation Motors (CIAM), Moscow, Russia in the field of optimization of fuel injection and mixing in air-breathing propulsion systems. The main results of theoretical, numerical simulation and experimental tests obtained in the previous research are in the papers and patents. The goals of the 14U/FM&AL programs are twofold: 1) to improve the working efficiency of the HU/FM&AL team in generating new innovative ideas and in conducting research in the field of fluid dynamics and acoustics, basically for improvement of supersonic and subsonic aircraft engines, and 2) to attract promising minority students to this research and training and, in cooperation with other HU departments, to teach them basic knowledge in Aerodynamics, Gas Dynamics, and Theoretical and Experimental Methods in Aeroacoustics and Computational Fluid Dynamics (CFD). The research at the HU/FM&AL supports reduction schemes associated with the emission of en 'ne pollutants for commercial aircraft and concepts for reduction of 91 observables for military aircraft. These research endeavors relate to the goals of the NASA Strategic Enterprise in Aeronautics concerning the development of environmentally acceptable aircraft. It is in this precise area, where the US aircraft industry, academia, and Government are in great need of trained professionals and which is a high priority goal of the Minority University Research and Education (MLTREP) Program, that the HU/FM&AL can make its most important contribution.

Morgan, Morris H., III↗

Observational evidence for groundwater influence on crop yields in the United States

As climate change shifts crop exposure to dry and wet extremes, a better understanding of factors governing crop response is needed. Recent studies identified shallow groundwater—groundwater within or near the crop rooting zone—as influential, yet existing evidence is largely based on theoretical crop model simulations, indirect or static groundwater data, or small-scale field studies. Here, we use observational satellite yield data and dynamic water table simulations from 1999 to 2018 to provide field-scale evidence for shallow groundwater effects on maize yields across the United States Corn Belt. We identify three lines of evidence supporting groundwater influence: 1) crop model simulations better match observed yields after improvements in groundwater representation; 2) machine learning analysis of observed yields and modeled groundwater levels reveals a subsidy zone between 1.1 and 2.5 m depths, with yield penalties at shallower depths and no effect at deeper depths; and 3) locations with groundwater typically in the subsidy zone display higher yield stability across time. We estimate an average 3.4% yield increase when groundwater levels are at optimum depth, and this effect roughly doubles in dry conditions. Groundwater yield subsidies occur ~35% of years on average across locations, with 75% of the region benefitting in at least 10% of years. Overall, we estimate that groundwater-yield interactions had a net monetary contribution of approximately $10 billion from 1999 to 2018. This study provides empirical evidence for region-wide groundwater yield impacts and further underlines the need for better quantification of groundwater levels and their dynamic responses to short- and long-term weather conditions.

60 APPLIED LIFE SCIENCES↗

Acceleration of energetic protons by interplanetary shocks

The University of Iowa instrument aboard Pioneer 11 detected 69 energetic proton events (EPE) (in the 0.6-3.4 MeV energy range) during 1973-1974 in the heliocentric radial range 1-5 AU. Sixty percent of the EPE peak within plus or minus 5 hours of a corotating interaction region (CIR) boundary, while 19% peak inside and 21% peak outside the interaction regions. Of the CIR boundaries at which an EPE peaks with plus or minus 5 hours, 80% have associated shocks. The observed intensities and pitch angle distributions of protons near shock fronts are consistent with a theoretical simulation of the acceleration of protons by a drift in the electric field at the shock front.

Pesses, M. E.↗

Data analysis study and performance evaluation of the scanning laser Doppler system

A simulation program which provided information on theoretically expected vortex spectra, evaluations of potential algorithms, and expected location accuracies for given scan patterns is presented. Field tests using an aircraft engine flow field and aircraft vortices during flyby tests were compared to the results of the simulation. From these studies, a vortex location algorithm was developed which provided vortex location for one or two vortices as a function of time. Results of this algorithm used on data from flyby tests were used to study vortex transport, to evaluate system performance, and to provide suggestions for real-time vortex location algorithms. The results of real-time analysis were compared to those which were expected based on theoretical considerations.

Sonnenschein, C. M.↗

Statistical Performance of Forced Oscillation Detectors in the Presence of Missing Measurements

In bulk power systems, measurement-based monitoring for large oscillations can help maintain system reliability. One of the challenges encountered in a recent field demonstration was the unavailability of measurements due to underlying measurement quality or communication problems. During the demonstration, the oscillation detector ignored a measurement location if even 10 seconds of data was missing. To extend the detector's ability to operate in these conditions, this paper evaluates the impact of three methods for addressing missing data. The strengths and weaknesses of each approach are evaluated using theoretical expressions for the probability of detection along with results from simulated data and publicly available field measurements. Based on these results, a suitable approach is identified that can extend the oscillation detector's performance when large segments of data are missing.

Follum, James D.↗

Thermalization from quantum entanglement: Jet simulations in the massive Schwinger model

We investigate the emergence of thermalization in a quantum-field-theoretic model mimicking the production of jets in QCD: the massive lattice Schwinger model coupled to external sources. Specifically, we compute the expectation values of local operators as functions of time and compare them to their thermal counterparts, quantify the overlap between the evolving density matrix and the thermal one, and compare the dynamics of the energy-momentum tensor to predictions from relativistic hydrodynamics. Through these studies, we find that the system approaches thermalization at late times and elucidate the mechanisms by which quantum entanglement drives thermalization in closed field-theoretic systems. Our results show how thermodynamic behavior emerges in real time from unitary quantum dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

FAST/Polar Conjunction Study of Field-Aligned Auroral Acceleration and Corresponding Magnetotail Drivers

The discrete aurora results when energized electrons bombard the Earth's atmosphere at high latitudes. This paper examines the physical processes that can cause field-aligned acceleration of plasma particles in the auroral region. A data and theoretical study has been carried out to examine the acceleration mechanisms that operate in the auroral zone and to identi@ the magnetospheric drivers of these acceleration mechanisms. The observations used in the study were collected by the Fast Auroral Snapshot (FAST) and Polar satellites when the two satellites were in approximate magnetic conjunction in the auroral region. During these events FAST was in the middle of the auroral zone and Polar was above the auroral zone in the near-Earth plasma sheet. Polar data were used to determine the conditions in the magnetotail at the time field-aligned acceleration was measured by FAST in the auroral zone. For each of the magnetotail drivers identified in the data study, the physics of field-aligned acceleration in the auroral region was examined using existing theoretical efforts and/or a long-system particle in cell simulation to model the magnetically connected region between the two satellites. Results from the study indicate that there are three main drivers of auroral acceleration: (1) field-aligned currents that lead to quasistatic parallel potential drops (parallel electric fields), (2) earthward flow of high-energy plasma beams from the magnetotail into the auroral zone that lead to quasistatic parallel potential drops, and (3) large-amplitude Alfven waves that propagate into the auroral region from the magnetotail. The events examined thus far confm the previously established invariant latitudinal dependence of the drivers and show a strong dependence on magnetic activity. Alfven waves tend to occur primarily at the poleward edge of the auroral region during more magnetically active times and are correlated with intense electron precipitation. At lower latitudes away from the poleward edge of the auroral zone is the primary field-aligned current region which results in the classical field- aligned acceleration associated with the auroral zone (electrons earthward and ion beams tailward). During times of high magnetic activity, high-energy ion beams originating from the magnetotail are observed within, and overlapping, the regions of primary and return field-aligned current. Along the field lines where the high-energy magnetotail ion beams are located, field-aligned acceleration can occur in the auroral zone leading to precipitating electrons and upwelling ionospheric ion beams. Field-aligned currents are present during both quiet and active times, while the Alfven waves and magnetotail ion beams were observed only during more magnetically active events.

Schriver, D.↗

Planar laser-induced fluorescence system for spatiotemporal ion velocity distribution function measurements

In this work, we present a planar laser-induced fluorescence (PLIF) system for measuring two-dimensional (2D), spatiotemporally resolved ion velocity distribution functions (IVDFs). A continuous-wave tunable diode laser produces a laser sheet that irradiates the plasma, and the resulting fluorescence is captured by an intensified CCD (ICCD) camera. Fluorescence images recorded at varying laser wavelengths are converted into 2D IVDFs using the Doppler shift principle. The developed diagnostic is implemented in an electron beam generated E × B plasma with a bulk plasma density of $\sim\!\!{10^{10}}{\text{c}}{{\text{m}}^{ - 3}}$ . The developed diagnostic is validated against a conventional single-point LIF method using photomultiplier tube-based detection, while significantly reducing the total measurement time by the number of spatial positions measured. The time-resolving capability of this diagnostic is tested by oscillating the plasma between two nominal operating modes with different density profiles and triggering the ICCD camera by the externally driven plasma oscillation. The measured 2D IVDF maps reveal several signatures of ion dynamics in this plasma source, including radially outflowing ions across the electric field and anomalous ion heating at the periphery, consistent with recent kinetic simulations and theoretical studies. A possible correlation between these ion kinetic features and rotating spoke structures is discussed.

E×B plasma↗

Time evolution from linear to nonlinear stages in magnetohydrodynamic parametric instabilities

The nonlinear evolution of the magnetohydrodynamic (MHD) parametric instability of wave fluctuations propagating along an unperturbed magnetic field is investigated. Both a magnetohydrodynamic perturbation-theoretical approach and a nonlinear MHD simulation are used. It is shown that high harmonic waves are rapidly excited by wave-wave coupling, and that the wave spectrum evolves from a state containing a small number of degrees of freedom in k space to one which contains a large number of degrees of freedom. It is found that the spectral evolution prior to nonlinear saturation is well described by the prturbation theory. During this stage, the ratio of the growth rate of the nth harmonic wave to the linear growth rate of the fundamental wave is n. The nonlinear saturation stage is characterized by a frequency shift of the fundamental wave that destroys the wave-wave resonance condition which, in turn, causes the wave amplitude to cease its growth.

Hoshino, M.↗

A simulation study of multiple ion wave generation downstream of low Mach number quasiperpendicular shocks

The behavior of minor ions just downstream of a low Mach number quasi-perpendicular shock is investigated both theoretically and by computer simulations. Because all ions see the same cross shock electric field their deceleration depends on their charge to mass ratio, yielding different downstream velocities. It is shown that these differences in velocity can lead to coherent wave structures in the downstream region of quasi-perpendicular shocks with a narrow transition layer. These waves are shown to be multi ion hybrid waves in contrast to mirror waves and ion cyclotron waves. Under favorable conditions these waves should be observable both at interplanetary shocks and at planetary bowshocks.

Motschmann, Uwe↗

The Onset of Magnetic Reconnection in Tail-Like Equilibria

Magnetic reconnection is a fundamental mode of dynamics in the magnetotail, and is recognized as the basic mechanisms converting stored magnetic energy into kinetic energy of plasma particles. The effects of the reconnection process are well documented by spacecraft observations of plasmoids in the distant magnetotail, or bursty bulk flows, and magnetic field dipolarizations in the near Earth region. Theoretical and numerical analyses have, in recent years, shed new light on the way reconnection operates, and, in particular, which microscopic mechanism supports the dissipative electric field in the associated diffusion region. Despite this progress, however. the question of how magnetic reconnection initiates in a tail-like magnetic field with finite flux threading the current i.sheet remains unanswered. Instead, theoretical studies supported by numerical simulations support the point-of-view that such plasma and current sheets are stable with respect to collisionless tearing mode. In this paper, we will further investigate this conclusion, with emphasis on the question whether it remains valid in plasma sheets with embedded thin current sheets. For this purpose, we perform particle-in-cell simulations of the driven formation of thin current sheets, and their subsequent evolution either to equilibrium or to instability of a tearing-type mode. In the latter case we will pay particular attention to the nature of the electric field contribution which unmagnetizes the electrons.

Hesse, Michael↗

The time-dependent dynamical response of the thermosphere to major geomagnetic disturbances

Several recent observations of thermospheric dynamics, made in the polar regions during extremely disturbed geomagnetic periods are reviewed. In general, the magnitude and the variability of winds in the thermospheric polar regions increases with magnetic activity, as measured by any of the conventional indices. However, none of the conventional indices is a particularly good aid to predicting wind magnitudes. In very general terms, two major factors may be considered in describing the wind system. The magnitude of the Interplanetary Magnetic Field (IMF) and, in particular, its southward component, determine the size of the auroral oval, and the strength of the cross-polar cap potential. This determines the size of the auroral oval, the magnitude of the sunward winds in the auroral oval and of those blowing anti-sunward over the polar cap, and is probably the major factor in determining the rate of geomagnetic energy deposition in the thermosphere. Superimposed on this enhanced polar circulation system are the effects of discrete auroral substorms. From a global view point, the effect of substorms is to generate a series of strong disturbances which propagate from their source region, usually near magnetic midnight in the auroral oval. The energy associated with discrete substorms is, however, usually a rather small proportion of the total global geomagnetic input during disturbed periods. These observations of thermospheric wind disturbances will be evaluated by comparison with global simulations of the thermospheric response to theoretical and semi-empirical models of the polar electric field, and of the effects of magnetospheric particle precipitation.

Rees, D.↗

The excavation stage of basin formation - A qualitative model

One of the most complex problems in planetary geology and geophysics is the determination of the nature of the impact cratering processes at scales of tens to thousands of kilometers that produce the complex morphological structures of multiring basins. The cratering process is frequently considered to be divided into three stages, including a short high-pressure stage of initial contact between the projectile and the planetary crust, a longer excavation or cratering flow stage culminating in the formation of a transient crater, and a still longer modification stage during which the transient crater is modified into the observed final geologic form. The transient crater may be considered as the initial boundary condition of the modification stage. In the present investigation, the nature of the transient crater is indicated by the cratering flow field determined from numerical simulations of the excavation stage. Attention is given to empirical and theoretical scaling.

Croft, S. K.↗