Search NASA⌕ Search

SEARCH · Search NASA

Results for “sound speed”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

How slowly can the early Universe expand?

When the expansion of the Universe is dominated by a perfect fluid with equation of state parameter w and a sound speed c s satisfying w = $c^{2}_{s}$ ≤ 1 , the Hubble parameter H and time t satisfy the bound Ht ≥ 1/3 . There has been recent interest in “ultraslow” expansion laws with Ht <1/3 (sometimes described as “fast expanding” models). Here, we examine various models that can produce ultraslow expansion: scalar fields with negative potentials, barotropic fluids, braneworld models, or a loitering phase in the early Universe. Scalar field models and barotropic models for ultraslow expansion are unstable to evolution toward w = 1 or w → ∞ in the former case and w → ∞ in the latter case. Braneworld models can yield ultraslow expansion but require an expansion law beyond the standard Friedman equation. Loitering early universe models can produce a quasistatic expansion phase in the early Universe but require an exotic negative-density component. These results suggest that appeals to an ultraslow expansion phase in the early Universe should be approached with some caution, although the loitering early universe may be worthy of further investigation. These results do not apply to ultraslow contracting models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Framework for phase transitions between the Maxwell and Gibbs constructions

By taking the nucleon-to-quark phase transition within a neutron star as an example, we present a thermodynamically consistent method to calculate the equation of state of ambient matter so that transitions that are intermediate to those of the familiar Maxwell and Gibbs constructions can be described. This method does not address the poorly known surface tension between the two phases microscopically (as, for example, in the calculation of the core pasta phases via the Wigner-Seitz approximation) but instead combines the local and global charge neutrality conditions characteristic of the Maxwell and Gibbs constructions, respectively. Overall charge neutrality is achieved by dividing the leptons to those that obey local charge neutrality (Maxwell) and those that maintain global charge neutrality (Gibbs). The equation of state is obtained by using equilibrium constraints derived from minimizing the total energy density. The results of this minimization are then used to calculate neutron star mass-radius curves, tidal deformabilities, equilibrium and adiabatic sound speeds, and nonradial g-mode oscillation frequencies for several intermediate constructions. Various quantities of interest transform smoothly from their Gibbs structures to those of Maxwell as the local-to-total electron ratio η, introduced to mimic the hadron-to-quark interface tension from 0 (Gibbs) to ∞ (Maxwell), is raised from 0 to 1. As a result, a notable exception is the g-mode frequency for the specific case of η = 1 for which a gap appears between the quark and hadronic branches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Equation of state at neutron-star densities and beyond from perturbative QCD

We explore the consequences of imposing robust thermodynamic constraints arising from perturbative quantum chromodynamics (QCD) when inferring the dense-matter equation-of-state (EOS). We find that the termination density, up to which the EOS modeling is performed in an inference setup, strongly affects the constraining power of the QCD input. This sensitivity in the constraining power arises from EOSs that have a specific form, with drastic softening immediately above the termination density followed by a strong stiffening. Here we also perform explicit modeling of the EOS down from perturbative-QCD densities to construct a new QCD likelihood function that incorporates additional perturbative-QCD calculations of the sound speed and is insensitive to the termination density, which we make publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗
DOE OSTI2022

Ablating Ion Velocity Distributions in Short-Pulse-Heated Solids via X-Ray Doppler Shifts

Solids ablate under laser irradiation, but experiments have not previously characterized the initiation of this process at ultrarelativistic laser intensities. Here, we present first measurements of bulk ion velocity distributions as ablation begins, captured as a function of depth via Doppler-shifted x-ray line emission from two viewing angles. Bayesian analysis indicates that bulk ions are either nearly stationary or flowing outward at the plasma sound speed. The measurements quantitatively constrain the laser-plasma ablation mechanism, suggesting that a steplike electrostatic potential structure drives solid disassembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Towards performing high-resolution inelastic X-ray scattering measurements at hard X-ray free-electron lasers coupled with energetic laser drivers

High-resolution inelastic X-ray scattering is an established technique in the synchrotron community, used to investigate collective low-frequency responses of materials. When fielded at hard X-ray free-electron lasers (XFELs) and combined with high-intensity laser drivers, it becomes a promising technique for investigating matter at high temperatures and high pressures. This technique gives access to important thermodynamic properties of matter at extreme conditions, such as temperature, material sound speed, and viscosity. The successful realization of this method requires the acquisition of many identical laser-pump/X-ray-probe shots, allowing the collection of a sufficient number of photons necessary to perform quantitative analyses. Here, a 2.5-fold improvement in the energy resolution of the instrument relative to previous works at the Matter in Extreme Conditions (MEC) endstation, Linac Coherent Light Source (LCLS), and the High Energy Density (HED) instrument, European XFEL, is presented. Some aspects of the experimental design that are essential for improving the number of photons detected in each X-ray shot, making such measurements feasible, are discussed. A careful choice of the energy resolution, the X-ray beam mode provided by the XFEL, and the position of the analysers used in such experiments can provide a more than ten-fold improvement in the photometrics. The discussion is supported by experimental data on 10 µm-thick iron and 50 nm-thick gold samples collected at the MEC endstation at the LCLS, and by complementary ray-tracing simulations coupled with thermal diffuse scattering calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multidimensional Tests of a Finite-Volume Solver for MHD With a Real-Gas Equation of State

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

stochprop Documentation (Release 1.0)

Simulations of infrasonic propagation in the atmosphere typically utilize a single atmospheric specification describing the acoustic sound speed, ambient winds, and density as a function of altitude. Due to the dynamic and sparsely sampled nature of the atmosphere, there is a notable amount of uncertainty in the atmospheric state at a given location and time so that a more robust analysis of infrasonic propagation requires inclusion of this uncertainty. This Python library, stochprop, has been implemented using methods developed jointly by infrasound scientists at Los Alamos National Laboratory (LANL) and the University of Mississippi’s National Center for Physical Acoustics (NCPA). This software library includes methods to quantify variability in the atmospheric state, identify typical seasonal variability in the atmospheric state and generate suites of representative atmospheric states during a given season, as well as perform uncertainty analysis on a specified atmospheric state given some level of uncertainty. These methods have been designed to interface between propagation modeling capabilities such as InfraGA/GeoAc and NCPAprop and signal analysis methods in the LANL InfraPy tool.

58 GEOSCIENCES↗
DOE OSTI2022

Sesame Equation of State for Kel F-800

I describe the construction of new equations of state for Kel F-800, as both linear and quadratic fits to shock data and Sesame 97200. This is the first SESAME table for Kel F-800 of which I’m aware. The SESAME table recovers shock data by construction and does reasonably well with thermal data, but differs with Brillouin scattering results for sound speed and density under static compression.

42 ENGINEERING↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
DOE OSTI2023

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

High explosives are almost always evolving in some form of metal containment vessel. This fact requires that a materials model for the evolution of the metal containment vessel be part of any simulation of the HE. Since the actual form of the strength model to be used is an open question, we incorporate the material strength properties in a relatively agnostic fashion, which will accommodate many strength models. The only restriction is that the stress components be a function of density, specific internal energy and velocities. In addition, the timescales and rates of the chemical reactions in the HE vary between thousands of seconds and nanoseconds. A hydrodynamics capable of operating in anticipation of the eventual violent release of energy via these chemical reactions is a necessity. The central issue in creating a hydrodynamics capable of spanning these timescales is to forgo the use of a Courant time step control necessary in an explicit hydrodynamics. The Courant time-step arises because of characteristic velocities associated with the material, such as the sound speed or, in materials with strength, characteristic longitudinal and transverse speeds associated with compression and shear. The method used to circumvent the need for a Courant time-step limit is to develop an implicit calculation of the advanced particle pressure and shear modulus derived from the conservation laws of mass, momentum and specific internal energy. This process creates the implicit hydrodynamics needed for the HE calculations of interest to this study.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Aging of Limestones and Silane–Siloxane-Based Protective Hydrophobics: The Impact of Heating–Cooling and Freeze–Thaw Cycles

Stones are traditionally used in construction and architectural applications as building elements due to their aesthetic and technical/structural performance. Like other environmental factors (rain, humidity, moisture, salt presence, biological activity, etc.), heating–cooling and freeze–thaw cycles significantly threaten the longevity of stone materials. Hence, considering the socio-economic and cultural value of stones, preventive actions such as hydrophobic coatings are applied to prevent or mitigate damage. The scope of this study is the performance assessment of limestones with different characteristics and the efficiency of various commercial silane/siloxane-based hydrophobic coatings when exposed to thermal variation and freeze–thaw. For that purpose, the standards EN 14066:2013 (determination of resistance to aging by thermal shock) and EN 12371:2010 (determination of frost resistance) were followed. Open porosity and static contact angles were estimated to assess the stone durability and water protection capabilities of the hydrophobics. Additionally, sound speed propagation velocity, quality of building material index, elastic modulus and flexural strength were measured to evaluate the variation of mechanical properties. Static contact angle revealed that the coatings maintained an efficient level of hydrophobicity even after thermal-shock and freeze–thaw weathering tests. The study also revealed a critical interaction between freeze–thaw cycles, hydrophobic coatings and structural integrity of the stones, mostly on more porous ones. When they are subjected to harsh environmental conditions, untreated porous limestones keep structural cohesion, allowing for the natural absorption and release of water during freezing and thawing. On the contrary, when limestones are treated, the hydrophobic coatings can moderately obstruct the water release due to the partial saturation of the porous framework by the products. It also probably resulted from the different mechanical behavior between the inner matrix and layer of stone coated, resulting in a premature breakout and mechanical damage of the stone.

Lisci, Carla (ORCID:000000032179130X)↗
DOE OSTI2022

Shock Hugoniot of Forged and Additively Manufactured 304L Stainless Steel

The purpose of this research was to measure the equation of state for additively manufactured (AM) and forged 304L stainless steel using a novel experimental technique. An understanding of the dynamic behavior of AM metals is integral to their timely adoption into various applications. The Hugoniot of the AM 304L was compared to that of the forged 304L at particle velocities where the material retains a two-wave structure. This comparison enabled us to determine the sensitivity of the equation of state to microstructure as varied due to processing. Our results showed that there was a measurable difference in the measured shock velocity between the AM and forged 304L. The shock wave velocities for the AM 304L were found to be ~3% slower than those for the forged 304L at similar particle velocities. To understand these differences, properties such as densities, sound speeds, and texture were measured and compared between the forged and AM materials. Our results showed that no measurable difference was found in these properties. Additionally, it is possible that differing elastic wave amplitudes may influence shock velocity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Data-Driven Tailoring Optimization of Thermoset Polymers Using Ultrasonics and Machine Learning

Thermoset polymers are highly demanded for their structural robustness, thermal stability, and chemical resistance. Tailoring the properties of these polymers for high-performance applications is often preferred to designing brand-new polymers. However, the traditional destructive techniques used to characterize their properties as a function of manufacturing parameters are expensive and time-consuming. A novel non-destructive, data-driven method leveraging ultrasonics and machine learning techniques to tailor the properties of thermosets as a function of the manufacturing parameters is demonstrated. Thermoset epoxy samples with varying curing temperatures (15–40 °C) and curing agent amounts (±40%) were manufactured and tested. Their curing kinetics were monitored by determining the sound speed in the material in real time, while the longitudinal modulus of the samples was determined post-cure. Machine learning models were developed using a k-nearest neighbors algorithm. These models were implemented to predict the curing and final elastic properties using the manufacturing parameters, i.e., stoichiometry and curing temperature, and vice versa. Understanding and modeling how these parameters affect the cure kinetics and final properties will allow for efficient and reliable optimization of thermoset tailoring and manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the Generalizability of Time-of-Flight Convolutional Neural Networks for Noninvasive Acoustic Measurements

Bulk wave acoustic time-of-flight (ToF) measurements in pipes and closed containers can be hindered by guided waves with similar arrival times propagating in the container wall, especially when a low excitation frequency is used to mitigate sound attenuation from the material. Convolutional neural networks (CNNs) have emerged as a new paradigm for obtaining accurate ToF in non-destructive evaluation (NDE) and have been demonstrated for such complicated conditions. However, the generalizability of ToF-CNNs has not been investigated. In this work, we analyze the generalizability of the ToF-CNN for broader applications, given limited training data. We first investigate the CNN performance with respect to training dataset size and different training data and test data parameters (container dimensions and material properties). Furthermore, we perform a series of tests to understand the distribution of data parameters that need to be incorporated in training for enhanced model generalizability. This is investigated by training the model on a set of small- and large-container datasets regardless of the test data. We observe that the quantity of data partitioned for training must be of a good representation of the entire sets and sufficient to span through the input space. The result of the network also shows that the learning model with the training data on small containers delivers a sufficiently stable result on different feature interactions compared to the learning model with the training data on large containers. To check the robustness of the model, we tested the trained model to predict the ToF of different sound speed mediums, which shows excellent accuracy. Furthermore, to mimic real experimental scenarios, data are augmented by adding noise. We envision that the proposed approach will extend the applications of CNNs for ToF prediction in a broader range.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Correlations between the Neutron Star Mass–Radius Relation and the Equation of State of Dense Matter

We develop an analytic method of inverting the Tolman–Oppenheimer–Volkoff relations to high accuracy. In principle, a specified energy density–pressure relation gives a unique mass–radius (M–R) relation and vice versa. Our method is developed from the strong correlations that are shown to exist between the neutron star mass–radius curve and the equation of state (EOS) or pressure–energy density relation. Selecting points that have masses equal to fixed fractions of the maximum mass, we find a semi-universal power-law relation between the central energy densities, pressures, sound speeds, chemical potentials, and number densities of those stars, with the maximum mass and the radii of one or more fractional maximum mass points. Rms fitting accuracies, for EOSs without large first-order phase transitions, are typically 0.5% for all quantities at all mass points. The method also works well, although less accurately, in reconstructing the EOS of hybrid stars with first-order phase transitions. These results permit, in effect, an analytic method of inverting an arbitrary M–R curve to yield its underlying EOS. We discuss applications of this inversion technique to the inference of the dense matter EOS from measurements of neutron star masses and radii as a possible alternative to traditional Bayesian approaches.

Bayesian statistics↗