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 55 records · Page 3

Enhanced Tensor Completion Based Approaches for State Estimation in Distribution Systems

Grid state estimation is essential for effective control and management of distribution systems. While weighted least squares has been the conventional method for state estimation, sparsity-aware methods have become popular due to their superior performance with limited data. Matrix completion and compressed sensing-based state estimation approaches exploit the underlying smoothness in the state variables. However, classic matrix completion methods do not take into account the temporal correlation of system states. Compressed sensing methods, on the other hand, require an appropriate choice of sparsifying basis that may not be easy to identify. This paper proposes a blocktensor completion based framework which uses an alternative approach to estimate voltage phasor, power injections and branch currents. This approach utilizes the temporal correlation of the system states in a tensor trace-norm minimization formulation with power flow equations as constraints. Herein, feature scaling is introduced in the problem formulation to benefit from the improved sensitivity of the tensor trace norm to the matrix columns in the scaled unfoldings of the tensor. Weighted tensor norm is utilized to exploit the structures of the different unfoldings of the state measurement tensor to improve the voltage estimation. The estimation accuracy is further improved by alternatively estimating the tensor columns and increasing the available data at each stage in the tensor completion process. The proposed methods are evaluated on the IEEE-33, 37 test systems and a 100- node test system. The proposed methods are shown to provide significant performance gains relative to the classic matrix and tensor completion based approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗

Bilocal Field Theory for Composite Scalar Bosons

We give a bilocal field theory description of a composite scalar with an extended binding potential that reduces to the Nambu–Jona-Lasinio (NJL) model in the pointlike limit. This provides a description of the internal dynamics of the bound state and features a static internal wave function, $\phi(\vec{r})$, in the center-of-mass frame that satisfies a Schrödinger–Klein–Gordon equation with eigenvalues m 2 . We analyze the “coloron” model (single perturbative massive gluon exchange) which yields a UV completion of the NJL model. This has a BCS-like enhancement of its interaction, ∝N c the number of colors, and is classically critical with g critical remarkably close to the NJL quantum critical coupling. Negative eigenvalues for m 2 lead to spontaneous symmetry breaking, and the Yukawa coupling of the bound state to constituent fermions is emergent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The quantum p -spin glass model: a user manual for holographers

We study a large-N bosonic quantum mechanical sigma-model with a spherical target space subject to disordered interactions, more colloquially known as the p-spin spherical model. Replica symmetry is broken at low temperatures and for sufficiently weak quantum fluctuations, which drives the system into a spin glass (SG) phase. The first half of this paper is dedicated to a discussion of this model's thermodynamics, with particular emphasis on the marginally stable SG. This phase exhibits an emergent conformal symmetry in the strong coupling regime, which dictates its thermodynamic properties. It is associated with an extensive number of nearby states in the free energy landscape. We discuss in detail an elegant approximate solution to the SG equations, which interpolates between the conformal regime and an ultraviolet-complete short distance solution. In the second half of this paper we explore the real-time dynamics of the model and uncover quantum chaos as measured by out-of-time-order four-point functions, both numerically and analytically. We find exponential Lyapunov growth, which intricately depends on the model's couplings and becomes strongest in the quantum critical regime. We emphasize that the SG phase also exhibits quantum chaos, albeit with parametrically smaller Lyapunov exponent than in the replica symmetric phase. An analytical calculation in the marginal SG phase suggests that this Lyapunov exponent vanishes in a particular infinite coupling limit. Finally, we comment on the potential meaning of these observations from the perspective of holography.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nuclear Charge Radii of Silicon Isotopes

The nuclear charge radius of 32 Si was determined using collinear laser spectroscopy. The experimental result was confronted with ab initio nuclear lattice effective field theory, valence-space in-medium similarity renormalization group, and mean field calculations, highlighting important achievements and challenges of modern many-body methods. The charge radius of 32 Si completes the radii of the mirror pair 32 Ar – 32 Si, whose difference was correlated to the slope L of the symmetry energy in the nuclear equation of state. Furthermore, our result suggests L ≤ 60 MeV, which agrees with complementary observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gluon saturation in proton and its contribution to single inclusive soft gluon production in high energy proton-nucleus collisions

The leading order single inclusive soft gluon production in high energy proton-nucleus (pA) collisions has been studied by various approaches for more than two decades. The first correction due to the gluon saturation in proton was analytically attempted recently through a diagrammatic approach in which only partial results were obtained. An important feature of the first saturation correction is that it includes both initial state and final state interactions. In this paper, we present the complete result derived from the Color Glass Condensate framework. Our approach is to analytically solve the classical Yang-Mills equations in the dilute-dense regime and then use the Lehmann-Symanzik-Zimmermann (LSZ) reduction formula to obtain gluon production from classical gluon fields.

Li, Ming↗

Samarium: from a distorted-fcc phase to melting under dynamic compression using in-situ x-ray diffraction

Abstract Lattice and electronic structure interactions for f -electrons are fundamental challenges for lanthanide equation of state development. Difficulties in first-principles calculations, such as density functional theory (DFT), emphasize the need for well-characterized experimental data. Here, we measure in-situ x-ray diffraction of shocked samarium (Sm) and temperature along the Hugoniot for the first time, providing direct evidence for phase transitions. We report direct evidence of a distorted fcc (dfcc) phase at 23 GPa. Shocked samarium melts from the dfcc phase starting at 33 GPa (1333 K), with complete melt at 40 GPa (1468 K). Previous work indicated shock melt at 27 GPa (1200 K), underscoring the significance of x-ray measurements for detecting phase transitions. Interestingly, our observed melting is in sharp contrast with the melting reported by a diamond anvil cell study. These experimental data can tightly constrain first principles calculations and serve as key touchstones for equation of state modeling.

36 MATERIALS SCIENCE↗

Equation of state of boron carbide B 4 ⁢C

We present the results of recent experiments conducted on the Sandia Z machine and a new tabular equation of state for B 4 ⁢C. The equation of state was calibrated to a combination of density functional calculations reported here and fits to preexisting data. It was constructed partly to recover the effects of a shock-driven, polymorphic phase transition of unknown character beginning at particle velocities of just under 3 km/s (shock pressures of 95 GPa). Some of the Z experiments included sound speeds determined by the overtaking rarefaction method, from which we calculate the Grüneisen parameter and compare with previous experiments conducted at the OMEGA laser [Fratanduono et al ., Phys. Rev. B 94 , 184107 (2016)], our own first principles calculations, and another recent tabular equation of state [Zhang et al ., Phys. Rev. E 102 , 053203 (2020)]. We also compare our results with previous static compression, thermophysical, and melt studies, finding mixed consistency. We predict the onset and completion of shock melting at 225 and 265 GPa, respectively, and predict a melt curve that is largely flat to pressures of several hundred GPa.

36 MATERIALS SCIENCE↗

The extended Vinet analytical equations of state

When testing new materials models, it is often not possible to use tabulated Equations of States (EOS) since the model might only be available in a smaller specialized research code without capability to read in tables. In such cases a versatile and easy to implement analytical EOS is needed. In this report I derive, give, and discuss all formulas necessary for implementing the thermodynamically complete extended Vinet EOS 3 into codes. While similar to a Mie-Grüneisen (MG) EOS, this Vinet EOS is based on an isotherm and substantially easier to correctly implement than a MG EOS. In the extended Vinet the isotherm’s density dependence is given as a power series in density to accommodate for higher density behavior which normally requires the Hugoniot based MG.

36 MATERIALS SCIENCE↗

DESI2024: Dark Energy Spectroscopic Instrument (DESI) survey, presentation of Year 1 results

The series of papers included in this special issue represent the first cosmological measurements from the Dark Energy Spectroscopic Instrument (DESI). These papers present the validation of the DESI clustering catalogs, the methodology for the cosmology measurements, and the full cosmology implications. These measurements were derived from a catalog of approximately 6.4 million extragalactic sources in a redshift range of 0.1 < z < 4.2. Owing to advances in fiber positioner technology and to a larger telescope, this sample was acquired in only one year of observation, yet is more than double the size of the galaxy and quasar samples from the first four phases of the Sloan Digital Sky Survey (SDSS) that spanned 20 years. With these new galaxy and quasar catalogs, the DESI collaboration completed percent-level constraints on the cosmic expansion history and the growth of structure. In turn, these new measurements provided a new platform from which to test the equation of state of dark energy, the summed mass of the three neutrino mass eigenstates, global curvature, and General Relativity. By surpassing the precision of the previous generation spectroscopic survey, the results in this special issue mark the transition in the field of cosmology to Stage-IV dark energy science.

baryon acoustic oscillations↗

Precise measurement of nuclear interaction cross sections towards neutron-skin determination with R 3 B

The R 3 B (Reactions with Relativistic Radioactive Beams) experiment as a major instrument of the NUSTAR collaboration for the research facility FAIR in Darmstadt is designed for kinematically complete studies of reactions with high-energy radioactive beams. Part of the broad physics program of R 3 B is to constrain the asymmetry term in the nuclear equation-of-state and hence improve the description of highly asymmetric nuclear matter (e.g., in neutron stars). For a precise determination of the neutron-skin thickness – an observable which is directly correlated with the symmetry energy in theoretical calculations – by measuring absolute fragmentation cross sections, it is essential to quantify the uncertainty and challenge the reaction model under stable conditions. During the successful FAIR Phase-0 campaign of R 3 B, we precisely measured the energy dependence of total interaction cross sections in 12 C + 12 C collisions, for a direct comparison with calculations based on the eikonal reaction theory.

Ponnath, L. (ORCID:0000000286742624)↗

Polyimide dynamically compressed to decomposition pressures: Two-wave structures captured by velocimetry and modeling

We performed a series of six plate impact experiments on polyimide and modeled them using new reactant and products equations of state combined with an Arrhenius rate model. The first experiment was diagnosed with embedded electromagnetic velocity gauges through which we directly observed attenuation of the lead shock to an approximately constant state over a propagation distance of roughly 4 mm. Simulated gauge profiles were in excellent qualitative agreement with experiment and suggested a sluggish chemical reaction that did not proceed to completion. The remaining five experiments were conducted in a transmission geometry and diagnosed velocimetrically at the sample/window interface. All five of these yielded profiles with a sharp shock followed by a more gradual approach to maximum interface velocity that was “rounded” to varying degree. These profiles proved difficult to interpret unambiguously due to the convolution of the reactive wave upon first shock with reflection of the lead wave and reshock or release by the window. Comparison with thermochemical calculations strongly suggests that the point of maximum interface velocity corresponds to the equilibrium reshock or release locus. We discuss the implications of this point for the practice of impedance matching based on the reflected Hugoniot of reactive materials such as polymers. The reactant and thermochemical products equations of state are developmental SESAME tables 97710 and 97720, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo

Abstract We present a discussion of recent progress in excited‐state‐specific quantum chemistry and quantum Monte Carlo alongside a demonstration of how a combination of methods from these two fields can offer reliably accurate excited state predictions across singly excited, doubly excited, and charge transfer states. Both of these fields have seen important advances supporting excited state simulation in recent years, including the introduction of more effective excited‐state‐specific optimization methods, improved handling of complicated wave function forms, and ways of explicitly balancing the quality of wave functions for ground and excited states. To emphasize the promise that exists at this intersection, we provide demonstrations using a combination of excited‐state‐specific complete active space self‐consistent field theory, selected configuration interaction, and state‐specific variance minimization. These demonstrations show that combining excited‐state‐specific quantum chemistry and variational Monte Carlo can be more reliably accurate than either equation of motion coupled cluster theory or multi‐reference perturbation theory, and that it can offer new clarity in cases where existing high‐level methods do not agree. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry

Otis, Leon↗

Automated Detection of Instability-Inducing Channel Geometry Transitions in Saint-Venant Simulation of Large-Scale River Networks

A new sweep-search algorithm (SSA) is developed and tested to identify the channel geometry transitions responsible for numerical convergence failure in a Saint-Venant equation (SVE) simulation of a large-scale open-channel network. Numerical instabilities are known to occur at “sharp” transitions in discrete geometry, but the identification of problem locations has been a matter of modeler’s art and a roadblock to implementing large-scale SVE simulations. The new method implements techniques from graph theory applied to a steady-state 1D shallow-water equation solver to recursively examine the numerical stability of each flowpath through the channel network. The SSA is validated with a short river reach and tested by the simulation of ten complete river systems of the Texas–Gulf Coast region by using the extreme hydrological conditions recorded during hurricane Harvey. The SSA successfully identified the problematic channel sections in all tested river systems. Subsequent modification of the problem sections allowed stable solution by an unsteady SVE numerical solver. The new SSA approach permits automated and consistent identification of problem channel geometry in large open-channel network data sets, which is necessary to effectively apply the fully dynamic Saint-Venant equations to large-scale river networks or for city-wide stormwater networks.

54 ENVIRONMENTAL SCIENCES↗

Characterizing Topological Order with Matrix Product Operators

One of the most striking features of gapped quantum phases that exhibit topological order is the presence of long-range entanglement that cannot be detected by any local order parameter. The formalism of projected entangled-pair states is a natural framework for the parameterization of gapped ground state wavefunctions which allows one to characterize topological order in terms of the virtual symmetries of the local tensors that encode the wavefunction. In their most general form, these symmetries are represented by matrix product operators acting on the virtual level, which leads to a set of algebraic rules characterizing states with topological quantum order. This construction generalizes the concepts of $\mathsf G$- and twisted injectivity; the corresponding matrix product operators encode all topological features of the theory and provide a complete picture of the ground state manifold on the torus. We show how the string-net models of Levin and Wen fit within this formalism and in doing so provide a particularly intuitive interpretation of the pentagon equation for F-symbols as the pulling of matrix product operators through the string-net tensor network. Our approach paves the way to finding novel topological phases beyond string nets and elucidates the description of topological phases in terms of entanglement Hamiltonians and edge theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accurate effective fluid approximation for ultralight axions

Ultralight axions are theoretically interesting and phenomenologically rich dark sector candidates, but they are difficult to track across cosmological timescales because of their fast oscillations. We resolve this problem by developing a novel method to evolve them efficiently and accurately. We first construct an exact effective fluid which at late times matches the axion but which evolves in a simple way. We then approximate this evolution with a carefully chosen equation of state and sound speed. With our scheme we find that we can obtain subpercent accuracy for the linear theory suppression of axion density fluctuations relative to that of cold dark matter without tracking even a single complete oscillation of the axion field. Furthermore, we use our technique to test other approximation schemes and to provide a fitting formula for the transfer function for the matter power spectrum in linear theory in axion models. Implementing our approach in existing cosmological axion codes is straightforward and will help unleash the potential of high-precision next-generation experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Block smoothers and generalized ideal interpolation in AMG (Final Report)

The Pennsylvania State University (“Subcontractor”) worked on developing new parallel algebraic multilevel methods suitable for solving PDEs. Specifically, work on the design of multigrid solvers for coupled systems of partial differential equations arising in numerical modeling of various applications was completed. A main emphasis was on the design of new ideal algebraic multigrid interpolation for problems such as Maxwell’s equations where block smoothers are needed and the standard form of ideal interpolation is not an effective choice.

97 MATHEMATICS AND COMPUTING↗

Near-optimal combination of high performance and insensitivity in a shock compressed high explosive single crystal

Achieving the desired combination of superior detonation performance and insensitivity to shock initiation has been a long-standing goal in high explosive (HE) science and technology. Having previously established the shock insensitivity of 1,1-diamino-2,2-dinitroethene (also known as DADNE or FOX-7) single crystals to 20 GPa (extended to 25 GPa in this work), the FOX-7 detonation response was determined through wave profile measurements in ~250 μm thick single crystals shock compressed to 64 GPa. Quite unexpectedly, FOX-7 demonstrated the classic Chapman–Jouguet (C–J) detonation response—reaction completion in the detonation front (<0.7 ns) at pressures of 44 GPa and higher—not observed in other insensitive high explosives. Here, the experimentally determined C–J pressure (35 GPa), detonation wave velocities and the detonation products equation of state—together with shock insensitivity to 25 GPa—demonstrate that FOX-7 single crystals display a near-optimal combination of high performance and shock insensitivity, not observed in another HE crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗