SEARCH · Search NASA
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.
This report details the reduced NB5 velocity selector shielding box that has been isolated from the complete HFIR CGH model in order to speed up detailed calculations of the velocity selector shielding. This shielding model is intended for outsourcing shielding optimization calculations.
Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.
Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.
Abstract Modeling the Milky Way stellar halo requires well-determined density and velocity anisotropy profiles. However, it has been challenging to gather a large sample of stars with six-dimensional data that extend beyond 40 kpc to map the outer halo. Our work investigates the velocity anisotropy in the Milky Way stellar halo with more than 10,000 blue horizontal-branch stars, combining Gaia astrometric data and spectroscopic data from SEGUE, DESI, and LAMOST. This large sample allows us to obtain a detailed profile of up to ∼50 kpc. Radial velocities are predominant in the inner halo (<30 kpc), and the anisotropy presents a smooth decrease before rapidly dropping to negative values, becoming dominated by tangential dispersion velocities. Removing the main known accreted structures of the Milky Way makes the anisotropy profile radially dominated at all radii. Our profile clearly shows an increase in the anisotropy from the center of the Galaxy, in accordance to the simulations. We also investigate the correlation of anisotropy with metallicity and with color. The lack of correlation between metallicity and anisotropy in our clean sample reinforces that this relation is driven by merger events. The initial exploration with color indicates a relation between kinematics and age, showing that older stars are dynamically colder and present less radial orbits than younger stars in the inner halo.
Compact heat exchangers (HXs) have gained attention in recent years in various fields such as solar and nuclear power generation, oil and gas, and refrigeration due to their low cost, high power density, and robustness in high-pressure and/or high-temperature environments. However, the large difference between a compact HX's overall dimensions (~m) and the much smaller scale of its channels (~mm) makes it challenging to model the entire HX at once, due to computational limitations. In this work, we treat the channeled region of a compact HX as a porous medium (PM) to circumvent the need to model/mesh each individual channel. This allows us to simulate the entire HX, including both the header and channeled regions while maintaining the computational cost at a practical level. Although the porous medium approach has been used to model heat exchangers, its validity is still questionable because (1) the resistance coefficients are heavily data-based and thus difficult to be applied to new heat exchangers and (2) the validation has been focused on matching the overall pressure drop in the channel region, which does not address whether such model can predict detailed pressure and velocity field. For the first time, this work addresses under what circumstances and with what uncertainty does the PM approach work for hydrodynamics modeling in compact HXs. By answering these questions, we introduce the PM approach as a powerful tool for HX hydrodynamics modeling that can predict not only the overall pressure drop but also the detailed pressure and velocity distributions.
Nitric oxide molecular tagging velocimetry (NO MTV) is used to characterize the wake behind a free-flight spherical model in the test section of a free-piston reflected shock tunnel using a burst-mode laser operated at 100 kHz. A novel (to our knowledge) multi-delay timing scheme was implemented to measure velocity in a varying collisional environment in the wake of the free-flight model. Four simultaneous velocity profiles were measured in the wake of the model from -600 to 3600 m/s for flow enthalpies of 10.3 and 12.0 MJ/kg. Finally, the measured velocity distributions show good agreement when compared to computational fluid dynamics (CFD) modeling.
The ocean circulation around and over the Seychelles Plateau is characterized using 35 months of temperature and velocity measurements along with a numerical model. The results here provide the first documented description of the ocean circulation atop the Seychelles Plateau. The Seychelles Plateau is an unusually broad (~200 km), shallow (~50 m) plateau, dropping off steeply to the abyss. It is situated in a dynamic location (3.5-5.5S, 54-57E) in the south-western tropical Indian Ocean where northwesterly winds are present during austral summer and become southeasterly in austral winter, following the reversal of the western Indian Ocean monsoon winds. Measurements around the Inner Islands, on the Seychelles Plateau, have been carried out since 2015. Velocity measurements show that most of the depth-averaged current variance on the Seychelles Plateau arises from near-inertial oscillations and lower-frequency variability. Lower-frequency variability encompasses seasonal and intraseasonal variability, the latter of which includes the effects of mixed Rossby-gravity waves and mesoscale eddies. Here, a global 0.1-degree numerical ocean simulation is used in conjunction with these observations to describe the regional circulation around the Seychelles Plateau. Atop the Seychelles Plateau, circulation is dominated by ageostrophic processes consistent with Ekman dynamics, while around the Seychelles Plateau, both geostrophic and ageostrophic processes are important and vary seasonally. Stratification responds to the sea surface height semi-annual signal which is due to Ekman pumping-driven upwelling (related to the Seychelles-Chagos Thermocline Ridge) and the arrival of an annual downwelling Rossby wave.
In this work, we present a detailed implementation and validation of the droplet modeling framework proposed by Dahms and Oefelein (2016) into the engine commercial CFD software CONVERGE using the User Defined Function (UDF) interface. The model accounts for the nonlinear deformation and oscillation experienced by liquid spray droplet injected into high pressure and temperature. Lagrangian spray simulations of Engine Combustion Network (ECN) Spray A are performed. Model validation against standard experimental measurements of liquid velocity, vapor mixture fraction is conducted. To perform more rigorous model validation, new experimental measurements based on Diffused Back Illumination (DBI) are introduced. The new measurements are processed for Projected Liquid Volume (PLV), which offers as close as possible one-to-one model validation for liquid penetration while offering new insights into the spray physics. Comparison with a One-D model based on adiabatic mixing theory by Siebers (1999) and Desantes et al. (2007) are also conducted. Through these model validation exercises, it is shown that the new framework improves liquid-phase penetration predictions, following a tendency for enhanced evaporation, compared to the standard approach for both Reynolds Average Navier Stokes (RANS) and Large Eddy Simulation (LES). Additionally, at the liquid length, maximum mixture fraction values predicted by the new approach are in good agreement those of an adiabatic mixing model. Qualitative analysis of the spray behaviors during the early stage of the injection process reveals that the proposed framework predicts significant increase in droplet evaporation rate with lower drop drag compared to the current standard approach.
We extend the scale-dependent Gaussian Streaming Model (GSM) to produce analytical predictions for the anisotropic redshift-space correlation function for biased tracers in modified gravity models. Employing the Convolution Lagrangian Perturbation Theory (CLPT) re-summation scheme, with a local Lagrangian bias schema provided by the peak-background split formalism, we predict the necessary ingredients that enter the GSM, the real-space halo pairwise velocity and the pairwise velocity dispersion. We further consider effective field theory contributions to the pairwise velocity dispersion in order to model correctly its large scale behavior. We apply our method on two widely-considered modified gravity models, the chameleon-screened ƒ(R) Hu-Sawicki model and the nDGP Vainshtein model and compare our predictions against state-of-the-art N-body simulations for these models. We demonstrate that the GSM approach to predict the monopole and the quadrupole of the redshift-space correlation function for halos, gives very good agreement with the simulation data, for a wide range of screening mechanisms, levels of screening and halo masses at z = 0.5 and z = 1. Our work shows the applicability of the GSM, for cosmologies beyond GR, demonstrating that it can serve as a powerful predictive tool for the next stage of cosmological surveys like DESI, Euclid, LSST and WFIRST.
Continuous measurement and monitoring of surface water coverage in non-perennial streams are essential for understanding the exchange fluxes between surface and subsurface waters under both inundated and non-inundated conditions. In this study, a wildlife camera photo-based framework was developed to monitor small stream water inundation, depth, discharge, and velocity. Two advanced machine learning models, YOLOv8 and Mask2Former, were utilized to efficiently analyze images captured by wildlife cameras. The accuracy of the framework was validated against on-site depth measurements at six sites in the Yakima River Basin, along with the gage height, discharge, and velocity data from four USGS sites. This approach facilitates long-term, continuous monitoring and quantification of river intermittency and water availability with high precision and low cost, thereby advancing river ecosystem research and management.
Kilonovae are the radioactively powered optical/IR transients following a binary neutron star or neutron star black hole merger. A broad grid (varying ejecta mass, velocity, composition) of 900 2D kilonova models were simulated with Monte Carlo on Grizzly, using MPI+OpenMP. The set of models are being compared to observation AT2017gfo (the electromagnetic counterpart of gravitational wave GW170817), and hint at potential fidelity improvements.
Abstract In droplet-on-demand liquid metal jetting (DoD-LMJ) additive manufacturing, complex physical interactions govern the droplet characteristics, such as size, velocity, and shape. These droplet characteristics, in turn, determine the functional quality of the printed parts. Hence, to ensure repeatable and reliable part quality it is necessary to monitor and control the droplet characteristics. Existing approaches for in-situ monitoring of droplet behavior in DoD-LMJ rely on high-speed imaging sensors. The resulting high volume of droplet images acquired is computationally demanding to analyze and hinders real-time control of the process. To overcome this challenge, the objective of this work is to use time series data acquired from an in-process millimeter-wave sensor for predicting the size, velocity, and shape characteristics of droplets in DoD-LMJ process. As opposed to high-speed imaging, this sensor produces data-efficient time series signatures that allows rapid, real-time process monitoring. We devise machine learning models that use the millimeter-wave sensor data to predict the droplet characteristics. Specifically, we developed multilayer perceptron-based non-linear autoregressive models to predict the size and velocity of droplets. Likewise, a supervised machine learning model was trained to classify the droplet shape using the frequency spectrum information contained in the millimeter-wave sensor signatures. High-speed imaging data served as ground truth for model training and validation. These models captured the droplet characteristics with a statistical fidelity exceeding 90%, and vastly outperformed conventional statistical modeling approaches. Thus, this work achieves a practically viable sensing approach for real-time quality monitoring of the DoD-LMJ process, in lieu of the existing data-intensive image-based techniques.
Cosmic strings are predicted by many Standard Model extensions involving the cosmological breaking of a symmetry with nontrivial first homotopy group and represent a potential source of primordial gravitational waves (GWs). Present efforts to model the GW signal from cosmic strings are often based on minimal models, such as, e.g., the Nambu-Goto action that describes cosmic strings as exactly one-dimensional objects without any internal structure. In order to arrive at more realistic predictions, it is therefore necessary to consider nonminimal models that make an attempt at accounting for the microscopic properties of cosmic strings. With this goal in mind, we derive in this paper the GW spectrum emitted by current-carrying cosmic strings (CCCSs), which may form in a variety of cosmological scenarios. Our analysis is based on a generalized version of the velocity-dependent one-scale (VOS) model, which, in addition to the mean velocity and correlation length of the string network, also describes the evolution of a chiral (light-like) current. Here we are able to show, the solutions of the VOS equations imply a temporarily growing fractional cosmic-string energy density, Ω cs . This results in an enhanced GW signal across a broad frequency interval, whose boundaries are determined by the times of generation and decay of cosmic-string currents. Our findings have important implications for GW experiments in the Hz to MHz band and motivate the construction of realistic particle physics models that give rise to large currents on cosmic strings.
The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.
The additive manufacturing (AM) of metals is becoming an increasingly important production process with the potential to replace traditional techniques such as casting. Laser Powder Bed Fusion (LPBF) is used in many applications to print metal parts from powder. The metal powder is heated locally with sufficient laser radiation that the liquid melt easily reaches its boiling temperature, which leads to a metallic vapour jet that can entrain both powder bed particles and molten droplets. The small size of laser-matter interaction site makes a detailed experimental analysis of the process challenging. Synchrotron X-ray imaging experiments are one of the few methods which can capture the dynamic melting and solidification processes. Comparing such experiments with computer simulations of the process is an important approach in order to better understand the manufacturing process and to analyse the influence of process parameters on the evaporated gas jet and the subsequent impact on particle ejection, leading to potentially reduced AM component quality. The melting and solidification of the metal powder is simulated using an Eulerian multiphase approach based on a control volume discretization of powder bed and substrate and a volume of liquid separation from melt and gas phase. The gas phase modelled as an ideal gas reaches velocities up to 100 m/s. Lagrangian particle tracking in the simulation demonstrates that the velocity fields calculated by the Eulerian multi-phase approach in combination with a standard drag-force model lead to particle accelerations in good agreement with those measured experimentally. In order to avoid numerical laborious Lagrangian calculations, a direct method to compare an Eulerian multiphase simulation with synchrotron X-ray experiments was introduced and validated. This approach is used to analyse the influence of process parameters including laser power and laser speed on the maximal acceleration of particles from the melt pool area. While the particle acceleration increases linearly with line energy in the conduction mode, a linear decrease of the acceleration with increasing line energy can be found in the transition mode before the accel-eration increases again with line energy in the keyhole mode.
Dense, granular flows were examined to effectively capture and model bulk viscous properties in thin packed beds. A modified Couette cell with particle image velocimetry was used to experimentally determine pseudo-viscosity properties of four particulate media with varying morphologies: (1) iron oxide-coated SiO 2 particles, (2) CARBOBEAD CP30-60 particles, (3) CARBOBEAD CP40-100 particles, and (4) Al 2 O 3 beads. The pseudo-viscosity functions were fitted using a power law to correlate the measured shear stress as a function of measured shear rate. The pseudo-viscous functions were used as inputs to computation fluid dynamics models for a single-phase viscous fluid to predict granular flow profiles. Steady-state free surface velocity profiles at angular velocities <7 rad/s predicted by the model were in good agreement with the experimental particle image velocimetry measurements, resulting in Pearson correlation coefficients of 0.97 for iron-oxide coated SiO 2 particles and 0.95 for CP30-60 particles. As a result, this alternative approach to measuring pseudo-viscous properties under shearing and modeling bulk transport behavior of granular flow using computation fluid dynamics model offered significant reduction in computational load compared to discrete element methods.
ABSTRACT We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\, h^{-1}{\rm Mpc}$ with a model based on the aemulus cosmological emulator to measure the growth rate of cosmic structure, parametrized by fσ8. We obtain a measurement of fσ8(z = 0.737) = 0.408 ± 0.038, which is 1.4σ lower than the value expected from 2018 Planck data for a flat ΛCDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1\text{--}60\, h^{-1}{\rm Mpc}$ modelled by the aemulus cosmological emulator and find a 4.5σ tension in the amplitude of the halo velocity field with the Planck + ΛCDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used in the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales.
Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.