Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete equation of state”

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 37 records · Page 2

Response analysis of flexible aircraft with active control

The small-perturbation equations of motion of a flexible aircraft with an active control technology (ACT) system are developed to evaluate the stability and performance of the controlled aircraft. The total aircraft system is formulated in state vector format and the system of equations is completed with fully unsteady and low-frequency aerodynamics for arbitrary, complex configurations based on a potential aerodynamic method. The ACT system equations are incorporated in the digital computer program FCAP (Flight Control Analysis Program) which can be used for the analysis of complete aircraft configurations, including control system, with either low-frequency or fully unsteady aerodynamics. The application of classical performance analyses including frequency response, poles and zeros, mean-square response, and time response in FCAP in state vector format is discussed.

Noll, R. B.↗

Equation-of-state spinning fluids in the Einstein-Cartan theory

The relativistic fluid equations may be completed in two physically distinct methods. One method assumes the mass, rho, (or particle number) is conserved, while the other method assumes an equation of state of the form P = P(rho). A variational principle for the mass conservation method both with and without an intrinsic spin for the fluid was constructed earlier (Ray and Smalley, 1982 and 1983). A variational principle for the fluid described by an equation of state both with and without spin is formulated. In all cases the variational principle is set in the Einstein-Cartan metric-torsion U4 geometry. The results for general relativity follow as a special case.

Ray, John R.↗

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗

Open‐Source Anaerobic Digestion Modeling Platform, Anaerobic Digestion Model No. 1 Fast (ADM1F)

An open‐source modeling platform, called Anaerobic Digestion Model No. 1 Fast (ADM1F), is introduced to achieve fast and numerically stable simulations of anaerobic digestion processes. ADM1F is compatible with an iPython interface to facilitate model configuration, simulation, data analysis, and visualization. Faster simulations and more stable results are accomplished by implementing an advanced open‐source library of numerical methods called Portable Extensive Toolkit for Scientific Computation (PETSc) to solve the ADM1 system of equations. Leveraging PETSc, ADM1F can consistently complete a steady‐state simulation under 0.2 s, over 99% faster than a benchmark ADM1 model implemented with MATLAB while achieving agreement of model outputs within 1% of those obtained with the benchmark model. For dynamic simulations, however, ADM1F has a computational speed advantage only when the influent characteristics update more frequently than every 4 h. The ability of ADM1F to be useful as a tool to study anaerobic digestion systems is demonstrated through two example implementations of ADM1F: (1) a two‐phase co‐digestion scenario evaluating the impact of the organic loading rate and the substrate composition on reactor performance and stability, and (2) a conventional digester scenario assessing the effectiveness of recovery strategies after disruptions that led to instability. These examples demonstrate how the high simulation speed and the convenience of the iPython interface allow ADM1F to complete complex analyses within minutes, much faster than computational strategies currently reported in the literature.

anaerobic co-digestion↗

Integrating out beyond tree level and relativistic superfluids

We revisit certain subtleties of renormalization that arise when one derives a low-energy effective action by integrating out the heavy fields of a more complete theory. Usually these subtleties are circumvented by matching some physical observables, such as scattering amplitudes, but a more involved procedure is required if one is interested in deriving the effective theory to all orders in the light fields (but still to fixed order in the derivative expansion). As a concrete example, we study the U(1) Goldstone low-energy effective theory that describes the spontaneously broken phase of a $\phi$ 4 theory for a complex scalar. Working to lowest order in the derivative expansion, but to all orders in the Goldstones, we integrate out the radial mode at one loop and express the low-energy effective action in terms of the renormalized couplings of the UV completion. This yields the one-loop equation of state for the superfluid phase of (complex) $\phi$ 4 . We perform the same analysis for a renormalizable scalar SO(N) theory at finite chemical potential, integrating out the gapped Goldstones as well, and confirm that the effective theory for the gapless Goldstone exhibits no obvious sign of the original SO(N) symmetry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Post-detonation fireball modeling: Validation of freeze out approximations

A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.

97 MATHEMATICS AND COMPUTING↗

Semidirect calculation of steady two- and three-dimensional flows

This paper describes a semidirect method for rapidly solving steady-state viscous flows described by the complete Navier-Stokes equations. The current results are for two-dimensional incompressible flows in general channels at arbitrary Reynolds numbers, but work in progress on compressible and three-dimensional flows is also described. The basic concept of semidirect methods is to use the recently developed fast (direct, or noniterative) linear solvers to solve linearized equations, which are then iterated to solve the nonlinearity. The method used here is an extension of the Split NOS method (Roache, 1975)

Roache, P. J.↗

Acoustic properties, elasticity, and equation of state of glycerol under pressure

Herein we employed high-pressure Brillouin scattering to study the pressure dependencies of acoustic modes of glycerol up to 14 GPa at 300 K. We observed longitudinal acoustic velocities and transverse acoustic velocities for the first time from 5 to 14 GPa. The results allow the determination of a complete set of elastic properties and an accurate determination of the pressure–volume (P–V) equation of state (EOS). EOS parameters, K 0 = 14.9 ± 1.8 GPa and K' 0 = 5.6 ± 0.5, were determined from fits to the data from ambient pressure to 14 GPa. Direct volume measurements of the P–V EOS are consistent with those determined by Brillouin scattering. A deviation from a Cauchy-like relationship for elastic properties was observed, and the pressure dependencies of the photoelastic constants and relaxation times were documented from 5 to 14 GPa. These results have broad implications for glass-forming liquids, viscoelastic theory, and mode coupling theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hypersonic shock structure with Burnett terms in the viscous stress and heat flux

The continuum Navier-Stokes and Burnett equations are solved for one-dimensional shock structure in various monatomic gases. A new numerical method is employed which utilizes the complete time-dependent continuum equations and obtains the steady-state shock structure by allowing the system to relax from arbitrary initial conditions. Included is discussion of numerical difficulties encountered when solving the Burnett equations. Continuum solutions are compared to those obtained utilizing the Direct Simulation Monte Carlo method. Shock solutions are obtained for a hard sphere gas and for argon from Mach 1.3 to Mach 50. Solutions for a Maxwellian gas are obtained from Mach 1.3 to Mach 3.8. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Shock density thickness, density asymmetry, and density-temperature separation are all more accurately predicted by the Burnett equations than by the Navier-Stokes equations.

Chapman, Dean R.↗

Stochastic time-optimal control problems

Two types of stochastic time-optimal controls in a one-dimensional setting are considered. Multidimensional problems, in the case of complete state information available and the system modeled by stochastic differential equations, are studied under the formulation of minimizing the expected transient-response time. The necessary condition of optimality is the satisfaction for the value function of a parabolic partial differential equation with boundary conditions. The sufficient condition of optimality is also provided, based on Dynkin's formula. Finally, three examples are given.

Zhang, W.↗

Comparison of Burnett, super-Burnett and Monte Carlo solutions for hypersonic shock structure

The continuum Navier-Stokes, Burnett, and super-Burnett equations are solved for one-dimensional shock structures in various monoatomic gases. Solutions for a hard sphere gas, argon, and a Maxwellian gas from Mach 1.3 to Mach 50 are obtained. A numerical method utilizing the complete time-dependent continuum equations and obtaining the steady-state shock structure by allowing the system to relax from arbitrary initial conditions is employed. Shock density, velocity, temperature, and entropy profiles are also obtained using the direct simulation Monte Carlo method, and these results are used as bases for comparison for continuum solution profiles. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Solutions to the super-Burnett equations with coefficients as presently derived, however, are considered to be inferior to those of the Burnett equations.

Fiscko, Kurt A.↗

Challenging Practices of Algebraic Battery Life Models through Statistical Validation and Model Identification via Machine-Learning

Various modeling techniques are used to predict the capacity fade of Li-ion batteries. Algebraic reduced-order models, which are inherently interpretable and computationally fast, are ideal for use in battery controllers, technoeconomic models, and multi-objective optimizations. For Li-ion batteries with graphite anodes, solid-electrolyte-interphase (SEI) growth on the graphite surface dominates fade. This fade is often modeled using physically informed equations, such as square-root of time for predicting solvent-diffusion limited SEI growth, and Arrhenius and Tafel-like equations predicting the temperature and state-of-charge rate dependencies. In some cases, completely empirical relationships are proposed. However, statistical validation is rarely conducted to evaluate model optimality, and only a handful of possible models are usually investigated. This article demonstrates a novel procedure for automatically identifying reduced-order degradation models from millions of algorithmically generated equations via bi-level optimization and symbolic regression. Identified models are statistically validated using cross-validation, sensitivity analysis, and uncertainty quantification via bootstrapping. On a LiFePO 4 /Graphite cell calendar aging data set, automatically identified models utilizing square-root, power law, stretched exponential, and sigmoidal functions result in greater accuracy and lower uncertainty than models identified by human experts, and demonstrate that previously known physical relationships can be empirically "rediscovered" using machine learning.

25 ENERGY STORAGE↗

Earth Filter Mathematica Toolbox V4 User’s Guide: Simplifying, Visualizing and Quick Modeling of All Fission Products from 60 Nuclear Sources

The Earth Filter Mathematica toolbox is a set of user-friendly routines that accurately calculate the ingrowth and decay of any of the 908 fission products for an infinite combination of sources over a large dynamic range and large time scales. This toolbox was developed by Dr. W. P. Brug at Los Alamos National Laboratory and eliminates the need for the user to define and solve differential equations by automatically scanning the ENDF-349 dataset and creating a complete set of differential equations (including delayed neutron branches to both the ground and metastable states) for every fission product. Additionally, this code allows the user to choose the data source for half-lives from several key references (i.e., TORI, IAEA, ENDF or NUDAT.) The complete set of ENDF-349 independent and cumulative fission yields for all 60 isotopic targets and 4 different neutron spectrums are included. This allows nearly any nuclear fission source to be accurately modeled as a linear combination of endpoint sources. The Earth Filter Toolbox has no limit as to the length of split list, and split lists of 5 or more endpoint sources can be quite common.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Transient and steady-state velocity of domain walls for a complete range of drive fields

Approximate analytic solutions for transient and steady-state 180 deg domain wall motion in bulk magnetic material are obtained from the dynamic torque equations with a Gilbert damping term. The results for the Walker region in which the transient solution approaches the familiar Walker steady-state solution are presented in a slightly new form for completeness. An analytic solution corresponding to larger drive fields predicts an oscillatory motion with an average value which decreases with drive field for reasonable values of the damping parameter. These results agree with those obtained by a computer solution of the torque equation and those obtained with the assumption of a very large anisotropy field.

Bourne, H. C., Jr.↗

Transient and steady-state velocity of domain walls for a complete range of drive fields

Approximate analytic solutions for transient and steady-state 180 deg domain-wall motion in bulk magnetic material are obtained from the dynamic torque equations with a Gilbert damping term. The results for the Walker region in which the transient solution approaches the familiar Walker steady-state solution are presented in a slightly new form for completeness. An analytic solution corresponding to larger drive fields predicts an oscillatory motion with an average value of the velocity which decreases with drive field for reasonable values of the damping parameter. These results agree with those obtained by others from a computer solution of the torque equation and those obtained by others with the assumption of a very large anisotropy field.

Bourne, H. C., Jr.↗

Fully Associative, Nonisothermal, Potential-Based Unified Viscoplastic Model for Titanium-Based Matrices

A number of titanium matrix composite (TMC) systems are currently being investigated for high-temperature air frame and propulsion system applications. As a result, numerous computational methodologies for predicting both deformation and life for this class of materials are under development. An integral part of these methodologies is an accurate and computationally efficient constitutive model for the metallic matrix constituent. Furthermore, because these systems are designed to operate at elevated temperatures, the required constitutive models must account for both time-dependent and time-independent deformations. To accomplish this, the NASA Lewis Research Center is employing a recently developed, complete, potential-based framework. This framework, which utilizes internal state variables, was put forth for the derivation of reversible and irreversible constitutive equations. The framework, and consequently the resulting constitutive model, is termed complete because the existence of the total (integrated) form of the Gibbs complementary free energy and complementary dissipation potentials are assumed a priori. The specific forms selected here for both the Gibbs and complementary dissipation potentials result in a fully associative, multiaxial, nonisothermal, unified viscoplastic model with nonlinear kinematic hardening. This model constitutes one of many models in the Generalized Viscoplasticity with Potential Structure (GVIPS) class of inelastic constitutive equations.

Source record↗

Relative Shock Effects in Mixed Powders of Calcite, Gypsum, and Quartz: A Calibration Scheme from Shock Experiments

The shock behavior of calcite and gypsum is important in understanding the Cretaceous/Tertiary event and other terrestrial impacts that contain evaporite sediments in their targets. Most interest focuses on issues of devolatilization to quantify the production of CO2 or SO2 to better understand their role in generating a temporary atmosphere and its effects on climate and biota [e.g., papers in 1,2,3,4]. Devolatilization of carbonate is also important because the dispersion and fragmentation of ejecta is strongly controlled by the expansion of large volumes of gas during the impact process as well [5,6]. Shock recovery experiments for calcite yield seemingly conflicting results: early experimental devolatilization studies [7,8,9] suggested that calcite was substantially outgassed at 30 GPa (> 50%). However, the recent petrographic work of [10,11,12] presented evidence that essentially intact calcite is recovered from 60 GPa experiments. [13] reported results of shock experiments on anhydrite, gypsum, and mixtures of those phases with silica. Their observations indicate little or no devolatilization of anhydrite shocked to 42 GPa and that the fraction of sulfur, by mass, that degassed is approx.10(exp -2) of theoretical prediction. In another (preliminary) report of shock experiments on calcite, anhydrite, and gypsum, [14] observe calcite recrystallization when shock loaded at 61 GPa, only intensive plastic deformation in anhydrite shock loaded at 63 GPa, and gypsum converted to anhydrite when shock loaded at 56 GPa. [15] shock loaded anhydrite and quartz to a peak pressure of 60 GPa. All of the quartz grains were trans-formed to glass and the platy anhydrite grains were completely pseudomorphed by small crystallized anhydrite grains. However, no evidence of interaction between the two phases could be observed and they suggest that recrystallization of anhydrite grains is the result of a solid state transformation. [16] reanalyzed the calcite and anhydrite shock wave experiments of [17] using improved equations of state of porous materials and vaporized products. They determined the pressures for incipient and complete vaporization to be 32.5 and 122 GPa for anhydrite and 17.8 and 54.1 GPa for calcite, respectively, a factor of 2 to 3 lower than reported earlier by [17].

Bell, Mary S.↗

Pressure-temperature equation of state of Al 2 ⁢O 3 up to 14 Mbar and 40 kK

Sapphire (Al 2 ⁢O 3 ), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al 2 ⁢ O 3 as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al 2 ⁢ O 3 dictating a comprehensive overhaul. We develop two advanced EoSs for Al 2 ⁢ O 3 the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al 2 ⁢O 3 up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Kalita, Patricia [Sandia National Laboratories (SN↗