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.
Report discusses simplified recursive solutions for equations of class of fading-memory filters. Under appropriate assumptions regarding noise and exponential decay of effects of previous states, filter equations reduced to analytic formulas in closed form.
Solution of Boltzmann, and rate equations for bound electronic states in nonequilibrium magnetohydrodynamic plasma
When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.
The modelling of dense plasmas using finite-temperature density functional theory has proven very successful in determining transport properties and the equation of state of systems where quantum many-body effects and correlations play a key role in their structure. Here we show how excited-state projector augmented- wave potentials can be used to extend these calculations to explicitly model core-hole states, allowing for the calculation of the electronic structure of a range of integer charge configurations embedded in a dense plasma environment. Our excited-state potentials show good agreement with all-electron calculations at finite-temperatures, motivating their use as an efficient approach in modelling from first principles both the structure of strongly-coupled non-equilibrium plasmas and their interaction with intense X-rays.
The nonlinear interfacial instability of a liquid jet in a coflowing compressible airstream is studied numerically. A high-resolution scheme which has second-order accuracy in space and time is coupled with a Lagrangian marker particle algorithm to visualize the large-scale motion of the interfaces in compressible flow. A numerical algorithm based on an approximate equation of state of a compressible liquid is developed to allow this two-fluid system to be governed by the nonlinear unsteady Euler equations in conservative form. The initial growth of small disturbances given by the simulations agrees well with linear theory. The process of jet disruption in compressible flow is demonstrated to consist of the formation of liquid spikes, interweaving of the gas and liquid and stretching and detachment of the liquid main center core.
In order to provide a rational background for the analysis of experimental observations of blast wave phenomena, the conservation equations governing their nonsteady flow field are formulated in a general manner, without the usual restrictions imposed by an equation of state, and with proper account taken, by means of source terms, of other effects which, besides the inertial terms that conventionally dominate these equations, can affect the flow. Taking advantage of the fact that a blast wave can be generally considered as a spatially one-dimensional flow field whose nonsteady behavior can be regarded, consequently, as a function of just two independent variables, two generalized blast wave coordinates are introduced, one associated with the front of the blast wave and the other with its flow field. The conservation equations are accordingly transformed into this coordinate system, acquiring thereby a comprehensive character, in that they refer then to any frame of reference, being applicable, in particular, to problems involving either space or time profiles of the gas-dynamic parameters in the Eulerian system, or time profiles in the Lagrangian system.
The trace anomaly Δ ≡ 1/3 −𝑃/𝜖 =1/3 −𝜙 quantifies the possibly broken conformal symmetry in supradense matter under pressure 𝑃 at energy density 𝜖. Perturbative QCD (pQCD) predicts a vanishing Δ at extremely high energy or baryon densities when the conformal symmetry is realized but its behavior at intermediate densities reachable in neutron stars (NSs) is still very uncertain. The extraction of Δ from NS observations strongly depends on the employed model for nuclear equation of state (EOS). Using the IPAD-TOV method based on an intrinsic and perturbative analysis of the dimensionless (IPAD) Tolman-Oppenheimer-Volkoff (TOV) equations that are further verified numerically by using 10 5 EOSs generated randomly with a metamodel in a very broad EOS parameter space constrained by terrestrial nuclear experiments and astrophysical observations, here we first show that the compactness 𝜉 ≡ 𝐺𝑀 NS /𝑅𝑐 2 ≡ 𝑀 NS /𝑅 of a NS with mass 𝑀 NS and radius 𝑅 scales very accurately with $\bar{Π}$ c ≡ $Π$ c · (1 +18X/25) ≡ X/(1 +3X 2 +4X) · (1 +18X/25) where X ≡ 𝜙 c = 𝑃 c /𝜖 c is the ratio of pressure over energy density at NS centers. The scaling of NS compactness thus enables one to readily read off the central trace anomaly Δ c = 1/3 −X directly from the observational data of either the mass-radius or red-shift measurements. Finally, we then demonstrate indeed that the available NS data themselves from recent X-ray and gravitational wave observations can determine model insensitively the trace anomaly as a function of energy density in NS cores, providing a stringent test of existing NS models and a clear guidance in a new direction for further understanding the nature and EOS of supradense matter.
The paper concentrates on a linear approximation method for predicting the changes occurring in steady-state numerical solutions of the Euler equations as a consequence of small changes in the independent variables which control the problem. The importance of proper boundary-condition treatment and other issues concerning the problem are covered along with the importance of proper algorithm selection for a fully supersonic inviscid flow. The method is applied to a subsonic nozzle involving variation of the pressure on the outflow boundary and to a supersonic inlet involving variation of the inflow Mach number. In the subsonic test case, the comparisons between the predicted and conventional numerical solutions are shown to be good, while in the supersonic test case, the agreement between the approximation method and conventional numerical solution starts out well but rapidly degenerates at some point in the flowfield as the perturbation of the boundary conditions is increased.
Steady and unsteady state problems and shock wave structure using Krook model of Boltzmann equation
High-pressure neutron powder diffraction data from PbNCN were collected on the high-pressure diffraction beamline SNAP located at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (Tennessee, USA). The diffraction data were analyzed using the novel method of multidimensional (two dimensions for now, potentially more in the future) Rietveld refinement and, for comparison, employing the conventional Rietveld method. To achieve two-dimensional analysis, a detailed description of the SNAP instrument characteristics was created, serving as an instrument parameter file, and then yielding both cell and spatial parameters as refined under pressure for the first time for solid-state cyanamides/carbodiimides. The bulk modulus B 0 = 25.1 (15) GPa and its derivative B ′ 0 = 11.1 (8) were extracted for PbNCN following the Vinet equation of state. Surprisingly, an internal transition was observed beyond 2.0 (2) GPa, resulting from switching the bond multiplicities (and bending direction) of the NCN 2− complex anion. The results were corroborated using electronic structure calculation from first principles, highlighting both local structural and chemical bonding details.
Iron is an important metal, scientifically and technologically. It is a common metal on Earth, forming the main constituent of the planet's inner core, where it is believed to be in solid state at high pressure and high temperature. It is also the main component of many important structural materials used in quasistatic and dynamic conditions. Laser-driven Rayleigh–Taylor instability provides a means of probing material strength at high pressure and high temperature. The unavoidable phase transition in iron at relatively low pressure induces microstructural changes that ultimately affect its strength in this extreme regime. This inevitable progression can make it difficult to design experiments and understand their results. Here, we address this challenge with the introduction of a new approach: a direct-drive design for Rayleigh–Taylor strength experiments capable of reaching up to 400 GPa over a broad range of temperatures. We use 1D and 2D hydrodynamic simulations to optimize target components and laser pulse shape to induce the phase transition and compress the iron to high pressure and high temperature. At the simulated pressure–temperature state of 350 GPa and 4000 K, we predict a ripple growth factor of 3–10 depending on the strength with minimal sensitivity to the equation of state model used. The growth factor is the primary observable, and the measured value will be compared to simulations to enable the extraction of the strength under these conditions. These experiments conducted at high-energy laser facilities will provide a unique way to study an important metal.
A general analysis is presented of the structure of a steady state, plane-parallel shock wave in which both thermal conduction and radiative cooling are important. The fluid is assumed to have a perfect-gas equation of state, with radiative cooling a function only of its temperature and density. Conduction in both diffusive and saturated regimes is treated. For the case of a strong shock, with conductivity and cooling function varying as power laws in temperature, approximate analytic solutions describing the shock wave are derived. For a plasma of solar composition, conduction is found to have a significant effect on the shock temperature and overall thickness of the postshock layer only for shock velocities greater than about 30,000 km/s, corresponding to shock temperatures greater than about 10 to the 10th K, but it affects the local structure of parts of the shock wave at much lower velocities. The effects of conduction are greatly enhanced if the heavy-element abundance is increased.
Abstract Conventional phase diagram generation involves experimentation to provide an initial estimate of the set of thermodynamically accessible phases and their boundaries, followed by use of phenomenological models to interpolate between the available experimental data points and extrapolate to experimentally inaccessible regions. Such an approach, combined with high throughput first-principles calculations and data-mining techniques, has led to exhaustive thermodynamic databases (e.g. compatible with the CALPHAD method), albeit focused on the reduced set of phases observed at distinct thermodynamic equilibria. In contrast, materials during their synthesis, operation, or processing, may not reach their thermodynamic equilibrium state but, instead, remain trapped in a local (metastable) free energy minimum, which may exhibit desirable properties. Here, we introduce an automated workflow that integrates first-principles physics and atomistic simulations with machine learning (ML), and high-performance computing to allow rapid exploration of the metastable phases to construct “metastable” phase diagrams for materials far-from-equilibrium. Using carbon as a prototypical system, we demonstrate automated metastable phase diagram construction to map hundreds of metastable states ranging from near equilibrium to far-from-equilibrium (400 meV/atom). We incorporate the free energy calculations into a neural-network-based learning of the equations of state that allows for efficient construction of metastable phase diagrams. We use the metastable phase diagram and identify domains of relative stability and synthesizability of metastable materials. High temperature high pressure experiments using a diamond anvil cell on graphite sample coupled with high-resolution transmission electron microscopy (HRTEM) confirm our metastable phase predictions. In particular, we identify the previously ambiguous structure of n -diamond as a cubic-analog of diaphite-like lonsdaelite phase.
Accurate pressure calibration is fundamental to quantitative high-pressure science, yet inconsistencies among widely used secondary pressure scales persist. In this work, we report simultaneous synchrotron x-ray diffraction meatogether in diamond anvil cells to 140 GPa. This simultaneous-measurement strategy minimizes transitive errors and run-to-run inconsistencies inherent to paired-measurement approaches, yielding direct experimental volume-volume (V-V) relations among the calibrants. These V-V relations link the compression state of each phase to that of every other phase measured in the same experiment. By anchoring these V-V relations to the reduced 300 K equation of state (EOS) of copper derived from ramp-compression measurements, we derive an internally consistent set of Cu-referenced pressure scales via Vinet EOS fits. Because the primary experimental result is the V-V dataset, it provides an EOS-independent constraint that can be rereferenced if primary standards are revised. The dataset further enables a direct experimental cross-check of consistency between the adopted reduced 300 K Cu EOS and the corresponding reduced 300 K ramp-derived EOSs for Pt, Au, Ta, and hcp-Fe under a fixed V-V constraint. The residuals quantify small, material-dependent mismatches between each phase EOS and the Cu reference EOS.
This paper presents the application of the Generalized Fluid System Simulation Program (GFSSP) to model the time-dependent flow in a complex secondary flow circuit of the turbopump of the Fastrac engine currently under development at Marshall Space Flight Center. GFSSP is a general purpose computer program for analyzing steady-state and time-dependant flowrates, pressures, temperatures, and concentrations in a complex flow network. The program employs a finite volume formulation of mass, momentum and energy conservation equations in conjunction with the thermodynamic equation of state of real fluids. GFSSP was used to calculate the axial thrust and internal flow distribution of the Fastrac engine turbopump during the start and shut down transients. The models discussed in this paper use boundary conditions that were extracted from turbopump test data. The GFSSP predicted turbopump secondary flow passage pressures and temperatures were compared with actual measured values.
The application of an integrated rotorcraft identification method to the linear modeling of rotor system dynamics is studied. Two approaches used to describe the rotor flapping parameters in the rotor state dynamic equations and the period coefficients of the blade flapping equation of motion are presented. In the first approach, the parameters are identified in the nonrotating reference frame; in the second, the blade equivalent damping and spring periodic coefficients as well as other periodic coefficients are identified in the rotating reference frame.
Using deep convolutional neural network (CNN), the nature of the QCD transition can be identified from the final-state pion spectra from hybrid model simulations of heavy-ion collisions that combines a viscous hydrodynamic model with a hadronic cascade “after-burner”. Two different types of equations of state (EoS) of the medium are used in the hydrodynamic evolution. The resulting spectra in transverse momentum and azimuthal angle are used as the input data to train the neural network to distinguish different EoS. Different scenarios for the input data are studied and compared in a systematic way. A clear hierarchy is observed in the prediction accuracy when using the event-by-event, cascade-coarse-grained and event-fine-averaged spectra as input for the network, which are about 80%, 90% and 99%, respectively. A comparison with the prediction performance by deep neural network (DNN) with only the normalized pion transverse momentum spectra is also made. High-level features of pion spectra captured by a carefully-trained neural network were found to be able to distinguish the nature of the QCD transition even in a simulation scenario which is close to the experiments.
A generalized flux-vector splitting implicit Galerkin finite-element algorithm for the Euler equations in curvilinear coordinates for ideal and reacting gases is derived. For an arbitrary equation of state, the curvilinear-coordinate flux vector is split in kinematic and kinetic components, and the associated jacobian matrix eigenvalues explicitly depend on the metric data. After directional semidiscretization, the terminal ordinary differential-equation system is solved via a nonlinearly stable implicit Runge-Kutta scheme in concert with an accurate tensor matrix product factorization. The results for selected two-dimensional supersonic and axisymmetric hypersonic flows validate the algorithm and verify its robustness for curvilinear-coordinate computations. The evolution towards a steady state is achieved for large Courant numbers without indication of numerical instabilities.