Search NASA⌕ Search

SEARCH · Search NASA

Results for “velocity model”

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 361 records · Page 20

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

IGRINS RV: A Precision Radial Velocity Pipeline for IGRINS Using Modified Forward Modeling in the Near-infrared

Application of the radial velocity (RV) technique in the near-infrared is valuable because of the diminished impact of stellar activity at longer wavelengths, making it particularly advantageous for the study of late-type stars but also for solar-type objects. In this paper, we present the IGRINS RV open-source python pipeline for computing infrared RV measurements from reduced spectra taken with IGRINS, an R ≡ λ/Δλ ∼ 45,000 spectrograph with simultaneous coverage of the H band (1.49–1.80 μm) and K band (1.96–2.46 μm). Using a modified forward-modeling technique, we construct high-resolution telluric templates from A0 standard observations on a nightly basis to provide a source of common-path wavelength calibration while mitigating the need to mask or correct for telluric absorption. Telluric standard observations are also used to model the variations in instrumental resolution across the detector, including a yearlong period when the K band was defocused. Without any additional instrument hardware, such as a gas cell or laser frequency comb, we are able to achieve precisions of 26.8 m s{sup −1} in the K band and 31.1 m s{sup −1} in the H band for narrow-line hosts. These precisions are empirically determined by a monitoring campaign of two RV standard stars, as well as the successful retrieval of planet-induced RV signals for both HD 189733 and τ Boo A; furthermore, our results affirm the presence of the Rossiter–McLaughlin effect for HD 189733. The IGRINS RV pipeline extends another important science capability to IGRINS, with publicly available software designed for widespread use.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-dimensional wind-tunnel interference from measurements on two contours

This paper describes how wall-induced velocities near a model in a two-dimensional wind tunnel can be estimated from upwash distributions measured along two contours surrounding a model. The method is applicable to flows that can be represented by linear theory. It was derived by applying the Schwarz Integral Formula separately to the two contours and by exploiting the free-air relationship between upwashes along the contours. Advantages of the method are that only one flow quantity need by measured and no representation of the model is required. A weakness of the method is that it assumes streamwise interference velocity vanishes far upstream of the model. This method was applied to a simple theoretical model of flow in a solid-wall wind tunnel. The theoretical interference velocities and the velocities computed using the method were in excellent agreement. The method was then used to analyze experimental data acquired during adaptive-wall experiments at Ames Research Center. This analysis confirmed that the wall adjustments reduced wall-induced velocities near the model.

Schairer, E. T.↗

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗

Comparative study between experimental measurements and model predictions of bubble rise velocity in molten LiCl-KCl

This research investigated the shape and rise velocities of gas bubbles (helium, nitrogen, argon, and krypton) in molten LiCl-KCl at 500 °C using a specialized optical visualization apparatus. The apparatus features a custom quartz rectangular prism cell housed in a furnace, enabling the application of the shadowgraph measurement technique. Bubbles with diameters between 0.45 and 3.2 mm were produced via a capillary, and their equivalent diameters and rise velocities were analyzed using high-resolution imaging and centroid tracking algorithms. Important non-dimensional parameters were calculated to characterize behavior. Bubbles with 0.03 < Eo < 1.32, 17 < Re < 718, and Mo order of 10 -11 were generated in the molten LiCl-KCl, revealing that bubbles with diameters less than 1.46 mm remain spherical, whereas those with larger diameters are ellipsoidal. The transition and oscillation onset for bubble shape occurred at larger diameters than those previously observed in air–water systems. Comparative analysis of the gas-molten salt data against existing air–water bubble rise velocity correlations highlighted inaccuracies in predictions, particularly across the different bubble regimes. An existing correlation was calibrated to the molten salt data to improve the bubble velocity predictions. In conclusion, this study contributes essential experimental data for the design and safety evaluation of molten salt reactors and offers insights for the optimization of sparge systems for fission product removal.

Bubble rise velocity↗

3D Seismic Anatomy of a Watershed Reveals Climate–Topography Coupling That Drives Water Flowpaths and Bedrock Weathering

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east-west trending shale watershed in the Valley and Ridge province in the U.S.A. Using wave-equation travel-time tomography from a seismic array of >4000 geophones, we obtained a 3D P-wave velocity (Vp) model that resolves structures ~20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2700 m/s contour is greater under the N- as compared to S-facing hillslopes and under sub-planar as compared to concave-up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N- vs. S-facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P-wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. Here, the Vp result is best explained as the long-term integrated effect of groundwater flow-induced geochemical weathering of shale in response to climate-driven patterns of micro- and macro-topography.

54 ENVIRONMENTAL SCIENCES↗

Fluid-like shock compression of dilute polymer nanocomposites

The shock Hugoniot of heterogenous mixtures of discrete particles has been experimentally investigated for porous agglomerates and fluids, while the study of full-density solids has been primarily limited to compressed powders and high-dimensional composites. By dispersing ceria nanoparticles in a polyethylene matrix, we are able to examine the hydrodynamic behavior of a nonporous, heterogenous solid in thermal equilibrium during weak shock compression. Phase-driven discontinuities in the Hugoniot particle velocity–shock velocity (u−D) relationship of pure polyethylene are replicated in the nanocomposites but are shifted to lower velocity and to higher pressure with higher particle concentration. The results are explained using an isothermal, two-velocity fluid model under the hydrodynamic approximation. The model, which assumes a theoretical equation-of-state for ceria and either a low-order or high-order fit to the measured polyethylene Hugoniot, reasonably predicts the Hugoniot for two different polyethylene/ceria nanocomposites. Using the model, the mixture Hugoniot is shown to be insensitive to the Hugoniot of the stiffer constituent when the moduli are sufficiently disparate, while dependence on particle density and volume fraction is preserved through fluid-like motion.

Moore, Nathan W. [Sandia National Laboratories (SN↗

The DESI one-per cent survey: exploring the halo occupation distribution of luminous red galaxies and quasi-stellar objects with AbacusSummit

ABSTRACT We present the first comprehensive halo occupation distribution (HOD) analysis of the Dark Energy Spectroscopic Instrument (DESI) One-Percent Survey luminous red galaxy (LRG) and Quasi Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 h−1 Mpc in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of $f_\mathrm{sat} = 11\pm 1~{y{\ \mathrm{per\,cent}}}$, a mean halo mass of $\log _{10}\overline{M}_h/M_\odot =13.40^{+0.02}_{-0.02}$, and a linear bias of $b_\mathrm{lin} = 1.93_{-0.04}^{+0.06}$. For LRGs in 0.6 < z < 0.8, we find $f_\mathrm{sat}=14\pm 1~{{\ \mathrm{per\,cent}}}$, $\log _{10}\overline{M}_h/M_\odot =13.24^{+0.02}_{-0.02}$, and $b_\mathrm{lin}=2.08_{-0.03}^{+0.03}$. For QSOs, we infer $f_\mathrm{sat}=3^{+8}_{-2}\mathrm{per\,cent}$, $\log _{10}\overline{M}_h/M_\odot = 12.65^{+0.09}_{-0.04}$, and $b_\mathrm{lin} = 2.63_{-0.26}^{+0.37}$ in redshift range 0.8 < z < 2.1. Using these fits, we generate a large suite of high fidelity galaxy mocks, forming the basis of systematic tests for DESI Y1 cosmological analyses. We also study the redshift-evolution of the DESI LRG sample from z = 0.4 up to z = 1.1, revealling significant and interesting trends in mean halo mass, linear bias, and satellite fraction.

79 ASTRONOMY AND ASTROPHYSICS↗

Mapping the gas density with the kinematic Sunyaev-Zel’dovich and patchy screening effects: A self-consistent comparison

The secondary anisotropies of the cosmic microwave background (CMB) provide a wealth of astrophysical and cosmological information. Pairing measurements of the CMB temperature map obtained by DR5 of the Atacama Cosmology Telescope (ACT) with the imaging survey conducted by the Dark Energy Spectroscopic Instrument for the purposes of target selection, DECaLS DR9, we investigate two effects that are sensitive to the gas density 𝜏: kinematic Sunyaev-Zel’dovich (kSZ) and patchy screening or anisotropic screening (resulting from the Thomson scattering of CMB photons away from the line-of-sight by free electrons). In particular, we measure the stacked profiles of the gas density around luminous red galaxies (LRGs) at a mean redshift of 𝑧 ≈ 0.7. We detect the kSZ signal at 7.2⁢𝜎, and we find a signal at ∼ 4.1⁢𝜎 for the patchy screening estimator, which is in excess relative to the kSZ signal. We attribute this excess to contamination from CMB lensing. Here, we demonstrate the effect of lensing using 𝑁-body simulations, and we show that the screening signal is dominated by it. Accounting for lensing, our measurement places a 95% upper bound on the optical depth of the Extended DESI LRG sample of 𝜏 < 2.5 10 −4 for a mean value of the sample of 𝜏 ≈ 1.6 10 −4 . Furthermore, via hydro simulations, we show that the underlying optical depth signal measured by both effects (after removing the CMB lensing contribution) is in perfect agreement when adopting either a compensated aperture photometry (CAP) filter or a high-pass filter. Consistent with previous measurements, we see evidence for excess baryonic feedback around DESI LRGs in the patchy screening measurement. In the future, when both effects can be measured with high signal-to-noise, one can measure the amplitude ratio between them, which is proportional to the root-mean-square velocity of the host halo sample, and even place constraints on velocity-sensitive models such as modified gravity and phantom dark energy.

Astrophysical & cosmological simulations↗

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

A model to calculate fatigue damage caused by partial waking during wind farm optimization

Abstract. Wind turbines in wind farms often operate in waked or partially waked conditions, which can greatly increase the fatigue damage. Some fatigue considerations may be included, but currently a full fidelity analysis of the increased damage a turbine experiences in a wind farm is not considered in wind farm layout optimization because existing models are too computationally expensive. In this paper, we present a model to calculate fatigue damage caused by partial waking on a wind turbine that is computationally efficient and can be included in wind farm layout optimization. The model relies on analytic velocity, turbulence, and load models commonly used in farm research and design, and it captures some of the effects of turbulence on the fatigue loading. Compared to high-fidelity simulation data, our model accurately predicts the damage trends of various waking conditions. We also perform example wind farm layout optimizations with our presented model in which we maximize the annual energy production (AEP) of a wind farm while constraining the damage of the turbines in the farm. The results of our optimization show that the turbine damage can be significantly reduced, more than 10 %, with only a small sacrifice of around 0.07 % to the AEP, or the damage can be reduced by 20 % with an AEP sacrifice of 0.6 %.

17 WIND ENERGY↗

Back Trajectories during SAIL

This data set contains air mass back trajectories generated during the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign using the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT, Stein et al. 2015, Rolph et al. 2017). For each hour from January 1, 2021 to June 30, 2023, a 96-hour back trajectory was initiated at the SAIL sampling site, starting 100 m above ground level. The calculations used the GDAS meteorological data set with model vertical velocity. The files within this dataset show location and meteorological parameters of the back trajectory analysis. There are 21 columns within each data file. The first 2 columns pertain to back trajectory parameters: the trajectory number (which is always 1) and the meteorological grid number. The next 9 columns pertain to spatiotemporal information: the year, month, day, hour (0-23), and minute of the trajectory, the previous time of the trajectory (in negative hours), and the airmass location (latitude, longitude, and height in m AGL). The rest of the parameters within each file pertains to surface and meteorological characteristics of the airmass: pressure (in hPa), potential temperature (in K), temperature (in K), precipitation rate (in mm/hr), mixing depth (in m), relative humidity (in % relative to liquid water), specific humidity (in g/kg), H2O mixing ratio (in g/kg), underlying terrain height (in meters above sea level), and incoming radiation (in W/m^2).

back trajectory↗

Composition dependence of ion diffusion coefficients in gas mixtures at arbitrary field strengths

Expressions for the diffusion coefficient of ions in gas mixtures are obtained from momentum transfer theory, and are given in terms of the diffusion coefficients and drift velocities of the ions in the pure component gases. Blanc's law holds exactly at all field strengths if the mean free time between collisions is independent of velocity (Maxwell model), but otherwise there may be either positive or negative deviations from Blanc's law at high fields. Such deviations are of comparable magnitude for the diffusion coefficients and the mobility, but are not identical. Specific cases of inverse-power potentials are treated in further detail, and some numerical examples are given for rigid-sphere interactions.

Whealton, J. H.↗