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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Cooperative Education

Los Alamos National Laboratory (LANL) is a multidisciplinary national laboratory that conducts research and development in national security, engineering, materials science, computational modeling, and advanced manufacturing. The laboratory develops innovative technologies to address complex scientific and engineering challenges. This project focuses on the development and evaluation of high-performance absorbing structures through computational design, simulation, and engineering analysis. Absorbing structures are used in applications where damage mitigation, structural protection, and material efficiency are critical performance requirements. The increasing demand for lightweight, high-strength, and highly efficient structural systems has created a need for improved design methodologies capable of maximizing absorption while minimizing weight and material usage. The project utilizes advanced engineering software, including 3D CAD software and FEA, to generate and optimize structural concepts. Computational simulations are performed to evaluate structural behavior under loading conditions, while mathematical analyses are conducted using Python-based tools as well as established analytical equations from material and structural mechanics. The project benefits LANL by supporting the development of advanced design methodologies and improving the understanding of material and structural performance. During the internship term, a significant portion of the design development, simulation, and data analysis activities will be completed. Success of the project depends on collaboration among engineering mentors and technical staff members. Work will be conducted at Los Alamos National Laboratory using laboratory computing resources and engineering software.

42 ENGINEERING↗
DOE OSTI2022

Master State Distributed Estimator (masde)

The MSE technology provides a distributed state estimation scheme for power utilities by making it a real time dynamic system, versus the static snapshot method used today. This will allow utilities to verify power grid readings and identify false data on the communication network. The prevailing method of estimation takes all the data from the system into a single load flow equation at the utilities command center. This becomes a large algorithm that is very time consuming to solve, providing the engineers only snap shots of the system. The MFD distributes an algorithm to devises that are already installed at various locations on the power grid. This allows for much faster computer times because the algorithms contain magnitudes less data. When each distributed estimator completes its calculation, it sends the results back to the master device to check against all other distributed estimators.

Reen, DylanW.↗
DOE OSTI2025

Numerical schemes for 3-wave kinetic equations: A complete treatment of the collision operator

In our previous work Walton and Tran (2023), numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the long time dynamics of the equation given the multiple blow-up time phenomenon. In this second work in the series, we propose numerical treatments for the complete 3-wave kinetic equations, in which the complete, much more complicated collision operators are fully considered based on a novel conservative form of the equation. Here we then derive an implicit finite volume scheme to solve the equation. The new discretization uses an adaptive time-stepping method which allows for the simulations to be carried to very long times. Our computed solutions are compared with previously derived long-time asymptotic estimates for the decay rate of total energy of time-dependent solutions of 3-wave kinetic equations and found to be in excellent agreement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Effects of fragmentation on post-inflationary reheating

We consider the effects of fragmentation on the post-inflationary epoch of reheating. In simple single field models of inflation, an inflaton condensate undergoes an oscillatory phase once inflationary expansion ends. The equation of state of the condensate depends on the shape of the scalar potential, V(Φ), about its minimum. Assuming V(Φ) ~ Φ k , the equation of state parameter is given by w = P Φ /ρ Φ = (k - 2)/(k + 2). The evolution of condensate and the reheating process depend on k. For k ≥ 4, inflaton self-interactions may lead to the fragmentation of the condensate and alter the reheating process. Indeed, these self-interactions lead to the production of a massless gas of inflaton particles as w relaxes to 1/3. If reheating occurs before fragmentation, the effects of fragmentation are harmless. We find, however, that the effects of fragmentation depend sensitively to the specific reheating process. Reheating through the decays to fermions is largely excluded since perturbative couplings would imply that fragmentation occurs before reheating and in fact could prevent reheating from completion. Reheating through the decays to boson is relatively unaffected by fragmentation and reheating through scatterings results in a lower reheating temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Physical implications of the extrapolation and statistical bootstrap of nucleon structure function ratios F 2 n F 2 p for mirror nuclei He 3 and H 3

A nuclear physics example of statistical bootstrap is used on the MARATHON nucleon structure function ratio data in the quark momentum fraction regions x B → 0 and x B → 1. The extrapolated F 2 ratio as quark momentum fraction x B → 1 is $\frac{F^n_2}{F^p_2}$ → 0.4 ± 0.05 and this value is compared to theoretical predictions. The extrapolated ratio when x B → 0 favors the simple model of isospin symmetry with the complete dominance of sea quarks at low momentum fraction. At high- x B , the proton quark distribution function ratio d/u is derived from the F 2 ratio and found to be d/u → 1/6. Our extrapolated values for both the $\frac{F^n_2}{F^p_2}$ ratio and the d/u parton distribution function ratio are within uncertainties of perturbative QCD values from quark counting, helicity conservation arguments, and a Dyson-Schwinger equation with a contact interaction model. In addition, it is possible to match the statistical bootstrap value to theoretical predictions by allowing two compatible models to act simultaneously in the nucleon wave function. Finally, one such example is nucleon wave functions composed of a linear combination of a quark-diquark state and a three-valence quark correlated state with coefficients that combine to give the extrapolated F 2 ratio at x B = 1.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

DOE-ICoM/RIFT

Rapid Infrastructure Flood Tool (RIFT) is a two-dimensional hydrodynamic model based on the complete shallow water equations. RIFT has specifically been designed with rapid simulation in mind by utilizing commodity high performance computing technology and best-available nation-wide data. RIFT is used to predict the movement of water over land and resolve the spatial and temporal variability of flood depths, extent, and velocity. RIFT can be applied to many flood situations and has primarily been used to quantify flood extents from dam/levee failure or inland rainfall flooding.

Perkins, Bill [Pacific Northwest National Laborato↗
DOE OSTI2023

Boson star normal modes

Boson stars are gravitationally bound objects that arise in ultralight dark matter models and form in the centers of galactic halos or axion miniclusters. We systematically study the excitations of a boson star, taking into account the mixing between positive and negative frequencies introduced by gravity. We show that the spectrum contains zero-energy modes in the monopole and dipole sectors resulting from spontaneous symmetry breaking by the boson star background. We analyze the general properties of the eigenmodes and derive their orthogonality and completeness conditions which have non-standard form due to the positive-negative frequency mixing. The eigenvalue problem is solved numerically for the first few energy levels in different multipole sectors and the results are compared to the solutions of the Schrödinger equation in fixed boson star gravitational potential. The two solutions differ significantly for the lowest modes, but get close for higher levels. We further confirm the normal mode spectrum in 3D wave simulations where we inject perturbations with different multipoles. As an application of the normal mode solutions, we compute the matrix element entering the evaporation rate of a boson star immersed in a hot axion gas. The computation combines the use of exact wavefunctions for the low-lying bound states and of the Schrödinger approximation for the high-energy excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗
DOE OSTI2022

D 3 -creatine dilution for skeletal muscle mass measurement: historical development and current status

The French chemist Michel Eugène Chevreul discovered creatine in meat two centuries ago. Extensive biochemical and physiological studies of this organic molecule followed with confirmation that creatine is found within the cytoplasm and mitochondria of human skeletal muscles. Two groups of investigators exploited these relationships five decades ago by first estimating the creatine pool size in vivo with 14 C and 15 N labelled isotopes. Skeletal muscle mass (kg) was then calculated by dividing the creatine pool size (g) by muscle creatine concentration (g/kg) measured on a single muscle biopsy or estimated from the literature. This approach for quantifying skeletal muscle mass is generating renewed interest with the recent introduction of a practical stable isotope (creatine-(methyl-d 3 )) dilution method for estimating the creatine pool size across the full human lifespan. The need for a muscle biopsy has been eliminated by assuming a constant value for whole-body skeletal muscle creatine concentration of 4.3 g/kg wet weight. The current single compartment model of estimating creatine pool size and skeletal muscle mass rests on four main assumptions: tracer absorption is complete; tracer is all retained; tracer is distributed solely in skeletal muscle; and skeletal muscle creatine concentration is known and constant. Three of these assumptions are false to varying degrees. Not all tracer is retained with urinary isotope losses ranging from 0% to 9%; an empirical equation requiring further validation is used to correct for spillage. Not all tracer is distributed in skeletal muscle with non-muscle creatine sources ranging from 2% to 10% with a definitive value lacking. Lastly, skeletal muscle creatine concentration is not constant and varies between muscles (e.g. 3.89–4.62 g/kg), with diets (e.g. vegetarian and omnivore), across age groups (e.g. middle-age, ~4.5 g/kg; old-age, 4.0 g/kg), activity levels (e.g. athletes, ~5 g/kg) and in disease states (e.g. muscular dystrophies, <3 g/kg). Some of the variability in skeletal muscle creatine concentrations can be attributed to heterogeneity in the proportions of wet skeletal muscle as myofibres, connective tissues, and fat. These observations raise serious concerns regarding the accuracy of the deuterated-creatine dilution method for estimating total body skeletal muscle mass as now defined by cadaver analyses of whole wet tissues and in vivo approaches such as magnetic resonance imaging. A new framework is needed in thinking about how this potentially valuable method for measuring the creatine pool size in vivo can be used in the future to study skeletal muscle biology in health and disease.

body composition↗
DOE OSTI2025

EMPDF : inferring the Milky Way mass with data-driven distribution function in phase space

We introduce the emPDF (empirical distribution function), a novel dynamical modelling method that infers the gravitational potential from kinematic tracers with optimal statistical efficiency under the minimal assumption of steady state. emPDF determines the best-fitting potential by maximizing the similarity between instantaneous kinematics and the time-averaged phase-space distribution function (DF), which is empirically constructed from observation upon the theoretical foundation of oPDF (Han et al. 2016). This approach eliminates the need for presumed functional forms of DFs or orbit libraries required by conventional DF- or orbit-based methods. emPDF stands out for its flexibility, efficiency, and capability in handling observational effects, making it preferable to the popular Jeans equation or other minimal assumption methods, especially for the Milky Way (MW) outer halo where tracers often have limited sample size and poor data quality. We apply emPDF to infer the MW mass profile using Gaia DR3 data of satellite galaxies and globular clusters, obtaining enclosed masses of M (,r) = 26±8, 46±8, 90±13⁠, and 149±40 x 10 10 M ⊙ at r = 30, 50, 100⁠, and 200 kpc, respectively. These are consistent with the updated constraints from simulation-informed DF fitting (Li et al. 2020). While the simulation-informed DF offers superior precision owing to the additional information extracted from simulations, emPDF is independent of such supplementary knowledge and applicable to general tracer populations. emPDF is currently implemented for tracers with complete 6D kinematics within spherical potentials, but it can potentially be extended to address more general problems.

Astrophysics of Galaxies (astro-ph.GA)↗
DOE OSTI2026

Dark-matter-enhanced probe of relic neutrino clustering

We propose heavy decaying dark matter (DM) as a new probe of the cosmic neutrino background (C⁢𝜈⁢B). Heavy DM, with mass ≳ 10 9 GeV, decaying into neutrinos can be a new source of ultrahigh-energy (UHE) neutrinos. Including this contribution along with the measured astrophysical and predicted cosmogenic neutrino fluxes, we study the scattering of UHE neutrinos with the C⁢𝜈⁢B via standard weak interactions mediated by the 𝑍 boson. We solve the complete neutrino transport equation, taking into account both absorption and reinjection effects, to calculate the expected spectrum of UHE neutrino flux at future neutrino telescopes, such as the IceCube-Gen2 radio. We argue that such observations can be used to probe the C⁢𝜈⁢B properties and, in particular, local C⁢𝜈⁢B clustering. We find that, depending on the absolute neutrino mass and the DM mass and lifetime, a local C⁢𝜈⁢B overdensity ≳ 10 6 can be probed at the IceCube-Gen2 radio within ten years of data taking.

Dark matter↗
DOE OSTI2025

Global bases for nonplanar loop integrands, generalized unitarity, and the double copy to all loop orders

We introduce a constructive method for defining a global loop-integrand basis for scattering amplitudes, encompassing both planar and nonplanar contributions. Our approach utilizes a graph-based framework to establish a well-defined, non-redundant basis of integrands. This basis, constructed from a chosen set of non-redundant graphs together with a selection of irreducible scalar products, provides clear insights into various physical properties of scattering amplitudes and proves useful in multiple contexts, such as on-shell Ward identities and manifesting gauge-choice independence. A key advantage of our integrand basis is its ability to streamline the generalized unitarity method. Specifically, we can directly read off the coefficients of basis elements without resorting to ansätze or solving linear equations. This novel approach allows us to lift generalized unitarity cuts — expressed as products of tree amplitudes — to loop-level integrands, facilitating the use of the tree-level double copy to generate complete gravitational integrands at any loop order. This method circumvents the difficulties in identifying complete higher-loop-order gauge-theory integrands that adhere to the color-kinematics duality. Additionally, our cut-based organization is well-suited for expansion in hard or soft limits, aiding in the exploration of ultraviolet or classical limits of scattering amplitudes.

Effective Field Theories↗
DOE OSTI2021

Sublimation temperature sensor for temperature locale

The sublimation temperature sensor (or “sublime sensor”) provides a continuum of measurement locations in which certain maximum temperatures are achieved during a heat up/cool down cycle. A predetermined material is encapsulated within a vacuum-sealed, non-volatile long tube (i.e., both ends capped and L$\gg$ D). This assembly is then inserted and centered into a heated zone, such as a furnace, exhaust pipe, or reactor. As the temperature increases, the material will sublimate (i.e., a process of having both the solid and gaseous states of matter simultaneously present) and will begin to fill the void - moving outward in both directions toward the ends of the tube. Once beyond the elevated temperatures, the gas will de-sublimate (i.e., deposition) onto the inner wall of the tube. The desired result of the sensor is the ring of material that develops over a relatively short period of time. This material deposit can be equated with temperature at an exact location. There is no need to interpolate and/or extrapolate for the desired measurement. Accuracy has been recorded for temperature locations on the range of ±2 mm over a 1 m span. Likewise, the precision of the measurement is ±0.2% the overall sensor domain. Furthermore, individual tubes with unique materials and pressures can be bundled together to provide a complete temperature profile of the heated zone.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Complete definition of 𝑁 →Δ transition generalized parton distributions

We revisit the definition of the leading-twist chiral-even generalized parton distributions (GPDs) for 𝑁 →Δ baryon transitions. We identify and address deficiencies in previous definitions of the transition GPDs inspired by the transition form factors of the vector and axial vector currents. Through systematic analysis of all possible covariant structures, respecting discrete symmetries and the baryon spinor equations of motion, we derive complete sets of independent structures for the transition matrix elements of the vector and axial vector partonic operators. They contain additional structures proportional to the light-cone vector, corresponding to transition GPDs of vanishing first moment, which were not included in previous parametrizations. Their presence is confirmed independently by the light-front multipole expansion and the cross-channel SO(3) partial-wave analysis of the transition matrix elements. Our analysis provides a complete definition of the 𝑁 →Δ transition GPDs for use in theoretical and phenomenological studies.

baryons↗
DOE OSTI2024

Back-to-Back Inclusive Dijets in Deep Inelastic Scattering at Small x : Complete NLO Results and Predictions

We compute the back-to-back dijet cross section in deep inelastic scattering at small x to next-to-leading order (NLO) in the color glass condensate effective field theory. Our result can be factorized into a convolution of the Weizsäcker-Williams gluon transverse-momentum-dependent distribution function (WW gluon TMD) with a universal soft factor and an NLO coefficient function. The soft factor includes both double and single logarithms in the ratio of the relative transverse momentum P ⊥ of the dijet pair to the dijet momentum imbalance q ⊥ ; its renormalization group (RG) evolution is resummed into the Sudakov factor. Likewise, the WW TMD obeys a nonlinear RG equation in x that is kinematically constrained to satisfy both the lifetime and rapidity ordering of the projectile. Exact analytical expressions are obtained for the NLO coefficient function of transversely and longitudinally polarized photons. Our results allow for the first quantitative separation of the dynamics of Sudakov suppression from that of gluon saturation. They can be extended to other final states and provide a framework for precision tests of novel QCD many-body dynamics at the Electron-Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Chemically Recyclable Analogs of Styrene–Butadiene Copolymers Enabling Perfectly Linear Ethylene–Styrene Materials with Random Phenyl Distribution

Copolymerization of cyclopentene (CP) and 4-phenylcyclopentene (4PCP) at a full range of comonomer feed ratios is reported using Ru-based ring-opening metathesis polymerization (ROMP) yielding homogeneous copolymers analogous to poly(styrene-ran-1,4-butadiene) and poly(ethylene-ran-styrene) copolymers following hydrogenation under mild conditions. In all cases, total monomer conversions of 86%–92% yielded copolymers with compositions within 4% of monomer feeds. Analysis of equilibrium copolymerization thermodynamics, rarely performed on two cycloolefin monomers with low ring strain energies, provides rational design strategies for negotiating two monomers with different equilibrium monomer concentrations. Inverse-gated decoupled 13 C NMR analysis of dyad sequences on the resulting copolymer microstructures concludes a near-random distribution of comonomer units. The copolymers produced from ROMP have number-average molar masses up to 60 kg mol –1 , moderate dispersities (1.5 ± 0.1), and high trans olefin content (86% ± 2%) while glass-transitions temperatures follow the Fox equation and span the full range between homopolymer extremities of PCP (−96 °C) and P4PCP (17 °C). Unlike most prevulcanized elastomers, these materials undergo facile chemical recycling to monomer, producing complete ring-closing metathesis depolymerization (RCMD) of the polymer back to the CP comonomers. Quantitative olefin hydrogenation produced perfectly linear polyethylene with 4%–16% of the backbone carbons containing a phenyl pendant, analogous to ES copolymers with up to 71.5% w/w styrene units but with random distribution of the aromatic pendants. Thermal properties of these materials are discussed, which span from semicrystalline to amorphous, and with T g values notably less than the reported ES copolymer analogs at similar compositions.

animal feed↗
DOE OSTI2025

Electron and photon structure functions at two loops

We present a fully analytic computation of the complete electron and photon structure functions, or QED lepton parton distribution functions (PDFs) up to two-loop order. Our computation is performed using modern techniques of reduction to Master Integrals and solving them with the differential equation method. We obtain explicit expressions for the electron-in-electron, positron-in-electron, photon-in-electron, electron-in-photon, and photon-in-photon distributions at next-to-next-to-leading order (NNLO). Cross-checks against one-loop results, existing two-loop calculations, and a recent soft-collinear effective theory (SCET) analysis of the electron structure functions are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗