SEARCH · Search NASA

Results for “EQUATIONS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Astrophysical Equation-of-State Constraints on the Color-Superconducting Gap

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J⁢0348+0432, PSR J1624-2230, PSR J⁢0740+6620, GW170817). At baryon number chemical potential μ = 2.6 GeV we find a 95% upper limit on the color-flavor locked pairing gap Δ of 457 MeV using overly conservative assumptions and 216 MeV with more reasonable assumptions. Furthermore, this constraint may be strengthened by future astrophysical measurements as well as by future advances in high-density QCD calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

QCD Constraints on Isospin-Dense Matter and the Nuclear Equation of State

Understanding the behavior of dense hadronic matter is a central goal in nuclear physics as it governs the nature and dynamics of astrophysical objects such as supernovae and neutron stars. Because of the nonperturbative nature of quantum chromodynamics (QCD), little is known rigorously about hadronic matter in these extreme conditions. Here, lattice QCD calculations are used to compute thermodynamic quantities and the equation of state of QCD over a wide range of isospin chemical potentials with controlled systematic uncertainties. Agreement is seen with chiral perturbation theory when the chemical potential is small. Comparison to perturbative QCD at large chemical potential allows for an estimate of the gap in the superconducting phase, and this quantity is seen to agree with perturbative determinations. Since the partition function for an isospin chemical potential μ I bounds the partition function for a baryon chemical potential μ B = 3 μ I / 2 , these calculations also provide rigorous nonperturbative QCD bounds on the symmetric nuclear matter equation of state over a wide range of baryon densities for the first time. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Mean-field equation for phase-modulated optical parametric oscillator

The widely established techniques for the generation of ultrashort optical pulses rely on passive mode locking of lasers, with the output pulse duration and emission spectrum determined by the intrinsic lifetime of laser transition in the gain medium. Due to the instantaneous nature of nonlinear gain, optical parametric oscillators (OPOs) are capable of generating optical radiation in all timescales from continuous-wave (cw) to ultrashort femtosecond regime, if driven by laser pump sources in the corresponding time domain. In the ultrashort timescale, operation of OPOs conventionally relies on mode-locked pump lasers, with the concomitant disadvantages of large footprint and high cost. At the same time, the lack of gain storage mandates the use of synchronous pumping, resulting in increased complexity. In this paper, we present the concept of phase-modulated OPO driven by cw pump laser. The approach overcomes the traditional drawbacks of ultrafast OPOs, enabling femtosecond pulse generation without the need for synchronous pumping, resulting in a simplified, compact, and cost-effective architecture using cw input pump lasers. We derive a mean-field equation for a degenerate χ ( 2 ) OPO driven by a cw laser with intracavity electro-optic modulator (EOM), and also including dispersion compensation. The equation predicts the formation of stable femtosecond pulses ( < 200 fs ) , in both normal and anomalous dispersion regimes, with a controllable repetition rate determined by the frequency of the EOM. The remarkable functionality of the proposed scheme paves the way for the development of an alternative class of widely tunable coherent femtosecond light sources in both bulk and integrated format based on χ ( 2 ) OPOs using cw pump lasers. Published by the American Physical Society 2024

Sanchez, A. D. (ORCID:0000000327630339)↗

DOE OSTI2025

Conformal duality of the nonlinear Schrödinger equation: Theory and applications to parameter estimation

The nonlinear Schrödinger equation (NLSE) in one spatial dimension has stationary solutions similar to those of the linear Schrödinger equation (LSE) as well as more exotic solutions such as solitary waves and quantum droplets. Here, we present a newly discovered conformal duality which unifies the stationary and time-dependent traveling-wave solutions of the one-dimensional cubic-quintic NLSE, the cubic NLSE and LSE. Any two systems that are classified by the same single number called the cross ratio are related by this symmetry. Notably, the conformal duality can also be adapted in Newtonian mechanics and serves as a powerful tool for investigating physical systems that otherwise cannot be directly accessed in experiments. Further, we show that the conformal symmetry is a valuable resource to substantially improve NLSE parameter estimation from noisy empirical data by introducing an optimization afterburner. The new method therefore has far reaching practical applications for nonlinear physical systems. Published by the American Physical Society 2025

Reinhardt, David B. (ORCID:0009000409812838)↗

DOE OSTI2026

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗

DOE OSTI2024

On the Existence of Steady-State Solutions to the Equations Governing Fluid Flow in Networks

The steady-state solution of fluid flow in pipeline infrastructure networks driven by junction/node potentials is a crucial ingredient in various decision-support tools for system design and operation. While the nonlinear system is known to have a unique solution (when one exists), the absence of a definite result on the existence of solutions hobbles the development of computational algorithms, for it is not possible to distinguish between algorithm failure and non-existence of a solution. In this letter, we show that for any fluid whose equation of state is a scaled monomial, a unique solution exists for such nonlinear systems if the term solution is interpreted in terms of potentials and flows rather than pressures and flows. However, for gases following the CNGA equation of state, while the question of existence remains open, we construct an alternative system that always has a unique solution and show that the solution to this system is a good approximant of the true solution. Further, the existence result for flow of natural gas in networks also applies to other fluid flow networks such as water distribution networks or networks that transport carbon dioxide in carbon capture and sequestration. Most importantly, our result enables correct diagnosis of algorithmic failure, problem stiffness, and non-convergence in computational algorithms.

42 ENGINEERING↗

DOE OSTI2025

Hierarchical Network Partitioning for Solution of Potential-Driven, Steady-State Nonlinear Network Flow Equations

The solution of potential-driven steady-state flow in large networks is a task which manifests in various engineering applications, such as transport of natural gas or water through pipeline networks. The resultant system of nonlinear equations depends on the network topology, and in general, there is no numerical algorithm that offers guaranteed convergence to the solution (assuming a solution exists). Some methods offer guarantees in cases where the network topology satisfies certain assumptions, but these methods fail for larger networks. On the other hand, the Newton-Raphson algorithm offers a convergence guarantee if the starting point lies close to the (unknown) solution. It would be advantageous to compute the solution of the large nonlinear system through the solution of smaller nonlinear sub-systems wherein the solution algorithms (Newton-Raphson or otherwise) are more likely to succeed. Here, this letter proposes and describes such a procedure, a hierarchical network partitioning algorithm that enables the solution of large nonlinear systems corresponding to potential-driven steady-state network flow equations.

42 ENGINEERING↗

DOE OSTI2024

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum-inspired techniques with classical methods. It focuses on fixed-point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed-point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Influence of Equation-of-States on Supercritical CO2 Combustion Mixtures

Abstract Fossil fuel based direct-fired supercritical CO2 (sCO2) cycles are gaining the attention of industry, academia, and government due to their remarkable efficiency and carbon capture at high-source temperatures. Modeling plays an important role in the development of sCO2 combustors because experiments are very expensive at the designed operating conditions of these direct-fired cycles. Inaccurate density estimates are detrimental to the simulation output. Hence, this work focuses on comprehensive evaluation of the influence and applicability of various equation-of-states (EOS) which are being used in the supercritical combustion modeling literature. A state-of-the-art supercritical combustion modeling methodology is used to simulate counter-flow supercritical CO2 flames by using various equation-of-states. The results show that using the corresponding state principle to evaluate compressibility factor is not accurate. Also, van der Waal type EOSs predictions can be as accurate as complex Benedict–Webb–Rubin EOSs; hence van der Waal EOSs are more suitable to simulate sCO2 combustor simulations. Non-ideal effects are significant under the operating conditions considered in this work. The choice of EOS significantly influences the flame structure and heat release rate. Also, assuming the binary interaction parameter as zero is reasonable in sCO2 combustion simulations.

Energy & Fuels↗

DOE OSTI2021

Evaluating Diesel/Biofuel Blends Using Artificial Neural Networks and Linear/Nonlinear Equations

Abstract The use of biomass-derived additives in diesel fuel mixtures has the potential to increase the fuel’s efficiency, decrease the formation of particulate matter during its combustion, and retain the fuel’s behavior in cold weather. To this end, identifying compounds that enable these behaviors is paramount. The present work utilizes a series of linear and non-linear equations in series with artificial neural networks to predict the cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value of multi-component blends. Property values of pure components are predicted using artificial neural networks trained with existing experimental data, and these predictions and their expected errors are propagated through linear and non-linear equations to obtain property predictions for multi-component blends. Individual component property prediction errors, defined by blind prediction median absolute error, are 4.91 units, 7.84 units, 0.06 cSt, 4.00 °C, and 0.55 MJ/kg for cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value respectively. On average, property predictions for blends are shown to be accurate to within 6% of the blends’ experimental values. Further, a multitude of compounds expected to be produced from catalytically upgrading products of fast pyrolysis are evaluated with respect to their behavior in diesel fuel blends.

09 BIOMASS FUELS↗

DOE OSTI2022

Thermal equation of state for zoisite: Implications for the transportation of water into the upper mantle and the high-velocity anomaly in the Farallon plate

Seismic tomography studies have revealed a high-velocity anomaly at depths between 100 km and 300 km in the Farallon plate. However, the reasons for the high-velocity anomaly continue to be debated. An analysis of the mineral proportions of eclogites exhumed from the Farallon plate shows that the average amount of zoisite in eclogite is ∼16.0 vol%. Therefore, the presence of zoisite eclogite needs to be considered to explain the high-velocity anomaly of the Farallon plate. However, the thermal equation of state and stability of zoisite have not been fully investigated under high pressure−temperature (P-T) conditions. We investigated the high-pressure and high-temperature behavior of natural zoisite utilizing synchrotron single crystal-X-ray diffraction (XRD). The results indicate that zoisite is metastable up to 24.8 GPa and 700 K. We obtained the ambient unit-cell volume V0 = 901.26(3) Å3 by synchrotron single crystal-XRD measurement. We also fitted the pressure-volume-temperature data to a high-temperature Birch-Murnaghan equation of state and obtained the zero-pressure bulk modulus K0 = 134.7(8) GPa, the temperature derivative of the bulk modulus (∂K/∂T)P = −0.011(4) GPa/K, and the thermal expansion coefficients α0V = 1.9(7) × 10−5 K−1 and α1V = 3(2) × 10−8 K−2 when the pressure derivative of bulk modulus (K′0) is fixed at 4. By incorporating the results from previous studies, we calculated the density and bulk sound velocity profiles of zoisite eclogite along the Farallon plate geotherm. Finally, we infer that zoisite could carry water to depths of ∼300 km within cold subducting slabs and that the bulk sound velocity of typical zoisite eclogite with ∼61.0 vol% omphacite, ∼23.0 vol% garnet, and ∼16.0 vol% zoisite could cause the high-velocity anomalies at depths of 100−300 km in the Farallon plate.

Geology↗

DOE OSTI2020

Statistical Treatment of Inverse Problems Constrained by Differential Equations-based Models with Stochastic Terms

This paper introduces a statistical treatment of inverse problems constrained by models with stochastic terms. The solution of the forward problem is given by a distribution represented numerically by an ensemble of simulations. The goal is to formulate the inverse problem, in particular the objective function, to find the closest forward distribution (i.e., the output of the stochastic forward problem) that best explains the distribution of the observations in a certain metric. We use proper scoring rules, a concept employed in statistical forecast verification, namely energy, variogram, and hybrid (i.e., combination of the two) scores. We study the performance of the proposed formulation in the context of two applications: a coefficient field inversion for subsurface flow governed by an elliptic partial differential equation with a stochastic source and a parameter inversion for power grid governed by differential-algebraic equations. In both cases we show that the variogram and the hybrid scores produce better parameter inversion results than does the energy score, whereas the energy score leads to better probabilistic predictions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

KKT Preconditioners for PDE-Constrained Optimization with the Helmholtz Equation

This paper considers preconditioners for the linear systems that arise from optimal control and inverse problems involving the Helmholtz equation. Specifically, we explore an all-at-once approach. The main contribution centers on the analysis of two block preconditioners. Variations of these preconditioners have been proposed and analyzed in prior works for optimal control problems where the underlying partial differential equation is a Laplace-like operator. In this paper, we extend some of the prior convergence results to Helmholtz-based optimization applications. Our analysis examines situations where control variables and observations are restricted to subregions of the computational domain. We prove that solver convergence rates do not deteriorate as the mesh is refined or as the wavenumber increases. More specifically, for one of the preconditioners we prove accelerated convergence as the wavenumber increases. Additionally, in situations where the control and observation subregions are disjoint, we observe that solver convergence rates have a weak dependence on the regularization parameter. We give a partial analysis of this behavior. We illustrate the performance of the preconditioners on control problems motivated by acoustic testing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Neutron stars in $$f(\mathtt {R,L_m})$$ gravity with realistic equations of state: joint-constrains with GW170817, massive pulsars, and the PSR J0030+0451 mass-radius from NICER data

Abstract In this work, we investigate neutron stars (NS) in $$f(\mathtt {R,L_m})$$ f ( R , L m ) theory of gravity for the case $$f(\mathtt {R,L_m})= \mathtt {R}+ \mathtt {L_m}+ \sigma \mathtt {R}\mathtt {L_m}$$ f ( R , L m ) = R + L m + σ R L m , where $$\mathtt {R}$$ R is the Ricci scalar and $$\mathtt {L_m}$$ L m the Lagrangian matter density. In the term $$\sigma \mathtt {R}\mathtt {L_m}$$ σ R L m , $$\sigma $$ σ represents the coupling between the gravitational and particles fields. For the first time the hydrostatic equilibrium equations in the theory are solved considering realistic equations of state and NS masses and radii obtained are subject to joint constrains from massive pulsars, the gravitational wave event GW170817 and from the PSR J0030+0451 mass-radius from NASA’s Neutron Star Interior Composition Explorer ( NICER ) data. We show that in this theory of gravity, the mass-radius results can accommodate massive pulsars, while the general theory of relativity can hardly do it. The theory also can explain the observed NS within the radius region constrained by the GW170817 and PSR J0030+0451 observations for masses around $$1.4~M_{\odot }$$ 1.4 M ⊙ .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

The Bloch equation for spin dynamics in electron storage rings: Computational and theoretical aspects

In this paper, we describe our work on spin polarization in high-energy electron storage rings which we base on the Full Bloch equation (FBE) for the polarization density and which aims towards the [Formula: see text] option of the proposed Future Circular Collider (FCC-ee) and the proposed Circular Electron Positron Collider (CEPC). The FBE takes into account non spin-flip and spin-flip effects due to synchrotron radiation including the spin-diffusion effects and the Sokolov–Ternov effect with its Baier–Katkov generalization as well as the kinetic-polarization effect. This mathematical model is an alternative to the standard mathematical model based on the Derbenev–Kondratenko formulas. For our numerical and analytical studies of the FBE, we develop an approximation to the latter to obtain an effective FBE. This is accomplished by finding a third mathematical model based on a system of stochastic differential equations (SDEs) underlying the FBE and by approximating that system via the method of averaging from perturbative ODE theory. We also give an overview of our algorithm for numerically integrating the effective FBE. This discretizes the phase space using spectral methods and discretizes time via the additive Runge–Kutta (ARK) method which is a high-order semi-implicit method. We also discuss the relevance of the third mathematical model for spin tracking.

Physics↗

DOE OSTI2024

Numerical methods for fractional Fokker–Planck equation with multiplicative Marcus Lévy noises

The Fokker–Planck equation (FPE) is an important deterministic tool for investigating stochastic dynamical systems. In this paper, we consider the space-time fractional FPE driven by multiplicative Marcus Lévy noises. Efficient numerical schemes are presented to solve the equations. Stability and convergence of the methods are also discussed. We give some numerical experiments to validate our schemes, and examine the effects of parameters on solutions. Additionally, we analyze the maximal likely trajectories and the critical time for the change of the most probability location.

Mathematics↗

DOE OSTI2024

Simplified Universal Equations for Ionic Conductivity and Transference Number

Nernst-Einstein equation can provide a reasonable estimate of the ionic conductivity of dilute solutions. For concentrated solutions, alternate methods such as Green–Kubo relations and Einstein relations are more suitable to account for ion-ion interactions. Such computations can be expensive for multicomponent systems. Simplified mathematical expressions like the Nernst-Einstein equation do not exist for concentrated multicomponent mixtures. Newman’s treatment of multicomponent concentrated solutions yields a conductivity relation in terms of species concentration and Onsager phenomenological coefficients. However, the estimation of these phenomenological coefficients is not straightforward. Here, mathematical formulations that relate the phenomenological coefficients with the friction coefficients are developed, leading to simplified, ready-to-use expressions of conductivity and transference numbers that can be used for a wide range of ionic mixtures. This approach involves spectral decomposition of the matrix of Onsager phenomenological coefficients. The general analytical expressions for conductivity and transference number are simplified for binary electrolytes, and numerical solutions are provided for ternary and quaternary mixtures with ion dissociation.

Electrochemistry↗

DOE OSTI2019

H2OI95: A Stand-Alone Fortran Code for Evaluating the IAPWS-95 Equation-of-State Model for Water

H2OI95 is a stand-alone Fortran code for evaluating the IAPWS-95 equation-of-state model (Wagner and Pruss, 2002) for the thermodynamic properties of water. It further evaluates the corresponding thermochemical properties of water consistent with the CODATA recommendations (Cox et al., 1989). The IAPWS-95 model is based on a model equation for the dimensionless Helmholtz energy for which the primary variables are the inverse reduced temperature and reduced density. Here is the absolute temperature (K), is density (kg/m3), and the subscript refers to the critical point of water (647.096 K and 22.064 MPa pressure in this model, for which is 322 kg/m3). The code solves four basic types of problems, distinguished by the specified inputs: 1, Temperature (K) and density (kg/m3) or reduced density 2. Temperature (K) and pressure (MPa). 3. Temperature (K) on the saturation (liquid-vapor equilibrium) curve 4. Pressure (MPa) on the saturation curve Each type of problem is run using a corresponding input (text) file. All but the first type of problem require iteration. For example, to solve for desired temperature and pressure, the reduced density must be adjusted to give the desired pressure. Iteration is accomplished using the Newton-Raphson method, though the secant method is also used in solving the fourth type of problem. H2OI95 has been used to conduct numerical studies of convergence and the problem of multiple numerical solutions, only some of which correspond to valid results. Obtaining valid results depends mainly on appropriate choice of starting values. The default values in H2OI95 appear to consistently lead to generally desired results. With modification (not addressed here), H2OI95 can be used to support SUPCRT92 (Johnson et al., 1992) and similar codes that compute chemical thermodynamic properties of species and reactions over a wide range of temperature and pressure (273.16-1273K and 0-1000 MPa).

Wolery, ThomasJ.↗