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 127 records · Page 7

Online and noninvasive monitoring of battery health at negative-half cell in all-vanadium redox flow batteries using ultrasound

Hydrogen evolution is one of the major side reactions that is detrimental to the health of all-vanadium redox flow batteries, especially for long-term cycling. Effective, low-cost, and accurate online prediction and detection methods for hydrogen generation are not yet available. In this work, we designed an online, noninvasive ultrasonic probing approach for monitoring the state of charge (SoC), predicting the hydrogen generation, and detecting hydrogen gas bubbles in anolyte solutions. The technique employs a pulse-echo method to measure the sound speed and the acoustic attenuation coefficient of the anolyte solution. Through static offline experiments and online in operando experiments, we have demonstrated that when hydrogen gas is generated in anolyte solutions, large variations are observed in both sound speed and acoustic attenuation coefficient measurements. We found that the variations of acoustic attenuation coefficient are highly correlated (correlation coefficients >0.9) with the gas flow rate. In conclusion, the designed acoustic method can monitor the SoC of anolyte, predict the hydrogen generation, and detect the presence of gas bubbles in an anolyte solution and, thus, provide information about the state of health for operation and management of flow battery systems.

25 ENERGY STORAGE↗

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(exp -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗

Protostellar Collapse with a Shock

We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(sup -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.

Tsai, John C.↗

Characterization of shallow ocean sediments using the airborne electromagnetic method

Experimental airborne electromagnetic (AEM) survey data collected in Cape Cod Bay are used to derive continuous profiles of water depth, electrical depth, water conductivity, and bottom sediment conductivity. Through a few well-known empirical relationships, the conductivities are used, in turn, to derive density, porosity, sound speed, and acoustic reflectivity of the ocean bottom. A commercially available Dighem III AEM system was used for the survey without any significant modification. The helicopter-borne system operated at 385 and 7200 Hz; both were in a horizontal coplanar configuration. The interpreted profiles show good agreement with available ground truth data. Where no such data are available, the results appear to be very reasonable. Compared with the shipborne electrode array method, the AEM method can determine the necessary parameters at a much higher speed with a better lateral resolution over a wide range of water depths from 0 to perhaps 100 m. The bottom sediment conductivity that can be measured by the AEM method is closely related to physical properties of sediments, such as porosity, density, sound speed, and, indirectly, sediment types that might carry broad implications for various offshore activities.

Won, I. J.↗

Acoustic propagation in a thermally stratified atmosphere

Acoustic propagation in an atmosphere with a specific form of temperature profile has been investigated by analytical means. The temperature profile used is representative of an actual atmospheric profile and contains three free parameters. Both lapse and inversion cases have been considered. Although ray solution have been considered the primary emphasis has been on solutions of the acoustic wave equation with point force where the sound speed varies with height above the ground corresponding to the assumed temperature profile. The method used to obtain the solution of the wave equation is based on Hankel transformation of the wave equation, approximate solution of the transformed equation for wavelength small compared to the scale of the temperature (or sound speed) profile, and approximate or numerical inversion of the Hankel transformed solution. The solution displays the characteristics found in experimental data but extensive comparison between the models and experimental data has not been carried out.

Vanmoorhem, W. K.↗

Acoustic propagation in a thermally stratified atmosphere

Acoustic propagation in an atmosphere with a specific form of a temperature profile has been investigated by analytical means. The temperature profile used is representative of an actual atmospheric profile and contains three free parameters. Both lapse and inversion cases have been considered. Although ray solutions have been considered, the primary emphasis has been on solutions of the acoustic wave equation with point source where the sound speed varies with height above the ground corresponding to the assumed temperature profile. The method used to obtain the solution of the wave equation is based on Hankel transformation of the wave equation, approximate solution of the transformed equation for wavelength small compared to the scale of the temperature (or sound speed) profile, and approximate or numerical inversion of the Hankel transformed solution. The solution displays the characteristics found in experimental data but extensive comparison between the models and experimental data has not been carried out.

Vanmoorhem, W. K.↗

Data-driven acoustic measurement of moisture content in flowing biomass

Measuring the moisture content in flowing biomass is critical to processes such as liquid biofuel conversion, such as biogasoline, biodiesel, bio jet kerosene, etc. However, biomass tends to flow in aggregates, which results in significant inhomogeneities in the amount of biomass flowing in front of a sensor at a given time, and there can be significant overlap in the material properties of dry vs wet biomass, leading to poor signal-to-noise ratio. We present a technique for identifying biomass moisture content using a series of acoustic pitch-catch measurements to quantify the sound speed and acoustic amplitude through the biomass, in conjunction with classical machine learning techniques, including Naive Bayes, Random Forest, and K-Nearest Neighbors classification. We amplify the differences between the acoustic measurements in different moisture levels by collecting a series of pulse-echo measurements, which we sort in order of ascending sound speed. We test the accuracy of the technique on experimentally-prepared batches of corn stover biomass with specified moisture levels and measure the average error in the estimated moisture level as a function of the number of pitch-catch measurements used. We observe average estimation errors as low as 6.7% by increasing the number of measurements and optimizing the hyperparameters. This work presents a novel method determining moisture content in flowing biomass with inhomogeneous flow. Additionally, this technique has application in optimizing biomass conversion processes, as well as other fields including, paper production, natural fiber processing, and mineral extraction.

09 BIOMASS FUELS↗

Aircraft laser sensing of sound velocity in water - Brillouin scattering

A real-time data source for sound speed in the upper 100 m has been proposed for exploratory development. This data source is planned to be generated via a ship- or aircraft-mounted optical pulsed laser using the spontaneous Brillouin scattering technique. The system should be capable (from a single 10 ns 500 mJ pulse) of yielding range resolved sound speed profiles in water to depths of 75-100 m to an accuracy of 1 m/s. The 100 m profiles will provide the capability of rapidly monitoring the upper-ocean vertical structure. They will also provide an extensive, subsurface-data source for existing real-time, operational ocean nowcast/forecast systems.

Hickman, G. D.↗

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ultrasound Characterization of Plastic Bonded Explosives With Varied Binder Content and Density

Advanced explosive formulations and methods of production require rapid characterization for both processes and materials. Here, we present results for ultrasound characterization of a plastic bonded explosive (PBX) formulation as a function of binder weight percent and pressing density using through-transmission sound speed measurements. The explosive in this study was composed of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) combined with the fluoropolymer binder Kynar Ultraflex B. Binder content ranged from 2.6 wt.% to 11.5 wt.% with associated explosive formulation densities from 1.69 g/cc to 1.89 g/cc. We find linear and nonlinear changes in longitudinal and shear sound speeds, peak frequencies, and measured elastic properties for the different PBX formulations. These measured changes can differentiate the tested formulations, highlighting the potential for ultrasound as a nondestructive characterization method for newly-formulated PBXs. We also demonstrate through-transmission measurements using different coupling media and present results from an ultrasound technique called resonant ultrasound spectroscopy (RUS). Both transmission and RUS approaches may allow for rapid and even automated characterization of PBX components suitable for high-throughput manufacturing operations.

Engineering↗

Is there a seismic signature of the sun's magnetic field?

Both weak magnetic fields and latitudinally dependent acoustic perturbations remove the degeneracy of the azimuthal quantum number, m, of acoustic modes of otherwise spherically symmetrical solar model. In the case of acoustic perturbations, the degeneracy is removed because the range of latitude in which a mode propagates depends on m, and therefore modes of like principle order n and degree l sample the aspherical scalar sound speed distribution differently. In the magnetic case, the removal of the degeneracy is caused by the same geometrical effect, and is influenced by the anisotropy of the Lorentz forces. Asymptotic analysis is used to show that the frequency splittings cannot be unambiguously attributed to the direct effect of a magnetic field, and that the effect of such a field on frequency splittings can be reproduced by a perturbation to the sound speed.

Zweibel, Ellen G.↗

3D Relativistic Magnetohydrodynamic Simulations of Magnetized Spine-Sheath Relativistic Jets

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗

The antesonic condition for the explosion of core-collapse supernovae – II. Rotation and turbulence

In the problem of steady free fall on to a standing shockwave around a central mass, the ‘antesonic’ condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\le 3/16$, where cT is the isothermal sound speed in the subsonic post-shock flow, and vesc is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and the shock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models of accretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a ~5 percent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by ~20 to 40 percent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multidimensional simulations explode more easily than their one-dimensional counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG↗

Effects of heavy-element settling on solar neutrino fluxes and interior structure

We consider the effects of gravitational settling of both He and heavier elements on the predicted solar neutrino fluxes and interior sound speed and density profiles. We find that while the structural changes that result from the inclusion of both He and heavy-element settling are only slightly larger than the changes resulting from the inclusion of He settling alone, the additional increases in expected neutrino fluxes are of comparable size. Our preferred model with both He and heavy-element settling has neutrino count rates of 9.0 SNU for Cl-37 detectors and 137 SNU for Ga-71 detectors, as compared to 7.1 and 127 SNU for a comparable model without any diffusive separation, or 8.0 and 132 SNU for a model that includes He settling alone. We suggest that the correction factors by which the predicted neutrino fluxes of solar models calculated without including the effects of diffusion should be multiplied are 1.25 +/- 0.08 for Cl detectors, 1.07 +/- 0.02 for Ga detectors, and 1.28 +/- 0.09 for the B-8 flux (1 sigma errors). Comparison of internal sound speed and density profiles strongly suggests that the additional changes in calculated p-mode oscillation frequencies due to the inclusion of heavy-element settling will be small compared to the changes that result from He settling alone, especially for the higher degree modes. All models with diffusive separation give much better agreement with the observed depth of the convection zone than do nondiffusive models. The model that includes both He and heavy-element settling requires an initial He mass fraction Y = 0.280 and has a surface He abundance of Y = 0.251 at the solar age.

Proffitt, Charles R.↗

Volatile loss from accreting icy protoplanets

A large self-gravitating body does not easily lose significant mass because the escape velocity is much larger than the sound speed of atmosphere-forming species under ambient thermal conditions. The most significant exceptions to this are giant impacts or impact jetting by fast-moving projectiles. A very small object (e.g. a comet) also does not easily lose significant volatile mass upon formation because the energy release associated with its accretion is so small. (It can however lose a great deal of mass if it is subsequently moved closer to the Sun.) I argue that there is an intermediate mass range (corresponding to bodies with radii of approximately 300-800 km) for which the ambient steady-state mass loss is a maximum. By ambient, I mean those conditions pertaining to the formation region of the body. By steady state, I mean to exclude infrequent traumas (giant impacts). The existence of a preferred intermediate mass arises through the competition of growing gravitational containment and growing energy release by accretion; it corresponds typically to GM/(Rc(sub s)(exp 2)) approximately equals 2 to 4, where M is the protoplanet mass of radius R, and c(sub s) is the sound speed. Several factors determine the amount of volatile loss is this vulnerable zone during accretion but in general the loss is a substantial fraction of the volatiles, sometimes approaching 100 percent. The principal implication is that bodies larger than a few hundred kilometers in radius will not have a 'primitive' (i.e. cometary) composition. This is relevant for understanding Triton, Pluto, Charon, and perhaps Chiron.

Stevenson, D. J.↗

Progress in modeling atmospheric propagation of sonic booms

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by humidity and turbulence. A realistic atmosphere is invariably somewhat turbulent, and may be characterized by an ambient fluid velocity v and sound speed c that vary from point to point. The absolute humidity will also vary from point to point, although possibly not as irregularly. What is ideally desired is a relatively simple scheme for predicting the probable spreads in key sonic boom signature parameters. Such parameters could be peak amplitudes, rise times, or gross quantities obtainable by signal processing that correlate well with annoyance or damage potential. The practical desire for the prediction scheme is that it require a relatively small amount of knowledge, possibly of a statistical nature, concerning the atmosphere along, the propagation path from the aircraft to the ground. The impact of such a scheme, if developed, implemented, and verified, would be that it would give the persons who make planning decisions a tool for assessing the magnitude of environmental problems that might result from any given overflight or sequence of overflights. The technical approach that has been followed by the author and some of his colleagues is to formulate a hierarchy of simple approximate models based on fundamental physical principles and then to test these models against existing data. For propagation of sonic booms and of other types of acoustic pulses in nonturbulent model atmospheres, there exists a basic overall theoretical model that has evolved as an outgrowth of geometrical acoustics. This theoretical model depicts the sound as propagating within ray tubes in a manner analogous to sound in a waveguide of slowly varying cross-section. Propagation along the ray tube is quasi-one-dimensional, and a wave equation for unidirectional wave propagation is used. A nonlinear term is added to this equation to account for nonlinear steepening, and the formulation has been carried through to allow for spatially varying sound speed, ambient density, and ambient wind velocities. The model intrinsically neglects diffraction, so it cannot take into account what has previously been mentioned in the literature as possibly important mechanisms for turbulence-related distortion. The model as originally developed could predict an idealized N-waveform which often agrees with data in terms of peak amplitude and overall positive phase duration. It is possible, moreover, to develop simple methods based on the physics of relaxation processes for incorporating molecular relaxation into the quasi-one-dimensional model of nonlinear propagation along ray tubes.

Pierce, Allan D.↗

Universal relations for neutron star 𝑓-mode and 𝑔-mode oscillations

Among the various oscillation modes of neutron stars, 𝑓- and 𝑔- modes are the most likely to be ultimately observed in binary neutron star mergers due to their relatively large coupling and shared frequencies with tidal excitations. The 𝑓-mode frequency and damping time are known to correlate in normal neutron stars with their compactness, and previous fits to hadronic stars are extended and shown to be valid for an extremely broad sampling of equations of state using a piecewise polytropic parametrization scheme for hadrons and a constant sound-speed parametrization for quark matter. Separate fits applicable to quark (self-bound) stars are improved. Much more significant correlations exist with tidal deformability, and therefore with moment of inertia and quadrupole moment. It is conclusively demonstrated that these correlations are the same for all types of stars, whether hadronic, hybrid, or pure quark, and its accuracy is quantified. A novel 1-node branch of the 𝑓-mode that occurs in low-mass hybrid stars in a narrow mass range just beyond the critical mass necessary for a phase transition to appear is identified. This 1-node branch shows the largest, but still small, deviations from the universal correlation for any configuration. It is characterized by a nonmonotonic relation between neutron star mass and 𝑓-mode frequency, in contrast to the behavior otherwise observed in normal, quark and hybrid stars. The 𝑔-mode only exists in matter with a nonbarotropic equation of state involving temperature, chemical potential or composition (such as being out of beta equilibrium), or a phase transition in barotropic matter. Here, the 𝑔-mode therefore could serve as a probe for studying phase transitions in hybrid stars. In contrast with the 𝑓-mode, 𝑔-mode frequencies do not correlate well with tidal deformability, but depend strongly on properties of the transition (the density and the magnitude of the discontinuity) at the transition. Imposing causality and maximum mass constraints, a fit involving neutron star and phase transition properties is found and the 𝑔-mode frequency is determined to have an upper bound of about 1.25 kHz. However, if the sound speed 𝑐𝑠 in the inner core at densities above the phase transition density is restricted to 𝑐$^{2}_{𝑠}$ ≤1/3, 𝑔-mode frequencies can only reach about 0.8 kHz, which are significantly lower than 𝑓-mode frequencies (1.3–2.8 kHz). 𝑔-mode gravitational wave damping times are found to be extremely long, >10 4 s (102 s) in the inner core with 𝑐$^{2}_{𝑠}$ ≤1/3 (1), in comparison with 𝑓-mode damping times (0.1–1 s).

Composition of astronomical objects↗