SEARCH · Search NASA

Results for “hydrodynamics methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A Novel Framework for Parametric Analysis of Coastal Transition Zone Modeling

Abstract Vulnerability of coastal regions to extreme events motivates an operational coupled inland‐coastal modeling strategy focusing on the coastal transition zone (CTZ), an area between the coast and upland river. To tackle this challenge, we propose a top‐down framework for investigating the contribution of different processes to the hydrodynamics of CTZs with various geometrical shapes, different physical properties, and under several forcing conditions. We further propose a novel method, called tidal vanishing point (TVP), for delineating the extent of CTZs through the upland. We demonstrate the applicability of our framework over the United States East and Gulf coasts. We categorize CTZs in the region into three classes, namely, without estuary (direct river–coast connection), triangular‐, and trapezoidal‐shaped estuary. The results show that although semidiurnal tidal constituents are dominant in most cases, diurnal tidal constituents become more prevalent in the river segment as the discharge increases. Also, decreasing the bed roughness value promotes more significant changes in the results than increasing it by the same value. Additionally, the estuary promotes tidal energy attenuation and consequently decreases the reach of tidal signals through the upland. The proposed framework is generic and extensible to any coastal region.

Chegini, Taher↗

NASA NTRS1974

Approach to the origin of turbulence on the basis of two-point kinetic theory

Equations for the fluctuation correlation in an incompressible shear flow are derived on the basis of kinetic theory, utilizing the two-point distribution function which obeys the BBGKY hierarchy equation truncated with the hypothesis of 'ternary' molecular chaos. The step from the molecular to the hydrodynamic description is accomplished by a moment expansion which is a two-point version of the thirteen-moment method, and which leads to a series of correlation equations, viz., the two-point counterparts of the continuity equation, the Navier-Stokes equation, etc. For almost parallel shearing flows the two-point equation is separable and reduces to two Orr-Sommerfeld equations with different physical implications.

Tsuge, S.↗

NASA NTRS1996

Interplanetary and Interstellar Dust Particles: Reentry Heating and Capture in Aerogel

Studies of interplanetary and interstellar dust can provide significant information on the evolution of the solar system or stars, respectively. However, for reliable analysis it is crucial to know how the particles have been modified during reentry (in the case of interplanetary dust particles, or [DP's) and impact into the capture medium. In the case of stratospheric capture, particles will be heated by atmospheric drag. Subsequent capture of the particles will result in heating, ablation, accretion of the capture medium and possible fragmentation. Modeling of these processes is a useful way of assessing their effects on the interpretation of the compositional data for these particles. Previous work on reentry heating has shown that heat diffusion alone cannot adequately account for temperature gradients observed in IDP's. In fact, for any reasonable thermal parameters, calculations show the particles to be nearly isothermal. Here we extend those studies to include phase transitions. These preliminary results are promising and show a significant temperature gradient for a 100 micron (diameter) particle. The actual capture of the particles in silica aerogel is being modeled using a comprehensive shock hydrodynamic code (called CTH). Various options of this code were investigated to attempt to make the most appropriate choice of methods of impact, equations of state, and processes of energy transfer from capture material to particle. The initial calculations with the code used only 'reasonable' estimates for the physical parameters of silica aerogel. Through the literature searches and personal contacts with the knowledgeable scientists, the best possible mechanical and thermal data have been made available for these computations.

Szydlik, Paul P.↗

NASA NTRS1980

Clear air turbulence forecasting techniques

A method to improve clear air turbulence (CAT) forecasting by more effectively using the currently operational rawinsonde (RW) system is discussed. The method is called the Diagnostic Richardson Number Tendency (DRT) technique. The technique does not attempt to use the RW as a direct detector of the turbulent motion or even of the CAT mechanism structure but rather senses the synoptic scale centers of action which provide the energy to the CAT mechanism at the mesoscale level. The DRT algorithm is deterministic rather than statistical in nature, using the hydrodynamic equations (equations of motion) relevant to the synoptic scale. However, interpretation, by necessity, is probabilistic. What is most important with respect to its operational implementation is that this method uses the same input data as currently used by the operational National Meteorological Center prognostic models.

Keller, J. L.↗

DOE OSTI2020

Cross-code comparison of the impact of the fill tube on high yield implosions on the National Ignition Facility

Fill tubes are used to inject deuterium and tritium fuel into inertial confinement fusion capsules fielded on the National Ignition Facility. These fill tubes have been shown to have a detrimental effect on capsule performance, primarily by introducing a low-density pathway into the central fuel region that enables the jetting of ablation material into the hot spot. Due to the complexity of the highly nonlinear flow associated with the fill tube and the challenge of diagnosing the evolution of the fill tube jet late in the implosion experiments, the uncertainty in how this perturbation source evolves is great. In this work, we report on the results of a detailed code comparison performed to understand uncertainties in computational modeling of the impact of fill tubes on implosion performance. The study employed two radiation-hydrodynamics codes, HYDRA and xRAGE, which employ very different meshing strategies and hydrodynamics solvers, as well as two radiation transport methodologies, discrete ordinates and multi-group diffusion. Our results demonstrate generally good agreement between codes through most of the implosion although they indicate sensitivity to opacity averaging methods. Late in the implosion, differences arise in the distribution and amount of contaminant although these differences have a remarkably small impact on the amount of yield reduction due to the fill tube. While these results demonstrate sensitivity in fill tube modeling to algorithmic choices, the observed differences between codes are small relative to known sensitivities due to expected variations in the fill tube geometry. Finally, we have developed a methodology for performing multi-group diffusion simulations that show good agreement with the more accurate discrete ordinates method

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Machine Learning-Driven Conservative-to-Primitive Conversion in Hybrid Piecewise Polytropic and Tabulated Equations of State

We present a novel machine learning (ML)-based method to accelerate conservative-to-primitive inversion, focusing on hybrid piecewise polytropic and tabulated equations of state. Traditional root-finding techniques are computationally expensive, particularly for large-scale relativistic hydrodynamics simulations. To address this, we employ feedforward neural networks (NNC2PS and NNC2PL), trained in PyTorch (2.0+) and optimized for GPU inference using NVIDIA TensorRT (8.4.1), achieving significant speedups with minimal accuracy loss. The NNC2PS model achieves 𝐿 1 and 𝐿 ∞ errors of 4.54 × 10 −7 and 3.44 × 10−6, respectively, while the NNC2PL model exhibits even lower error values. TensorRT optimization with mixed-precision deployment substantially accelerates performance compared to traditional root-finding methods. Specifically, the mixed-precision TensorRT engine for NNC2PS achieves inference speeds approximately 400 times faster than a traditional single-threaded CPU implementation for a dataset size of 1,000,000 points. Ideal parallelization across an entire compute node in the Delta supercomputer (dual AMD 64-core 2.45 GHz Milan processors and 8 NVIDIA A100 GPUs with 40 GB HBM2 RAM and NVLink) predicts a 25-fold speedup for TensorRT over an optimally parallelized numerical method when processing 8 million data points. Moreover, the ML method exhibits sub-linear scaling with increasing dataset sizes. We release the scientific software developed, enabling further validation and extension of our findings. By exploiting the underlying symmetries within the equation of state, these findings highlight the potential of ML, combined with GPU optimization and model quantization, to accelerate conservative-to-primitive inversion in relativistic hydrodynamics simulations.

conservative-to-primitive conversion↗

DOE OSTI2022

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Simulation of non-linear bearing forces for post-stability investigation

Different types of bearing designs were developed to improve dynamic properties of rotor-bearing systems. Elliptical bearings, multisleeve bearings, tilting pad and other designs such as herringbone groove were utilized to increase resistance to the onset of self excited vibrations. Experimental trials are costly, two alternative methods are used to gain a qualitative insight. The first one creates mathematical model and applies both a digital or an analog computer simulation. The second one investigates phenomena occurring on the laboratory rig with the bearing replaced by an electronic simulating device, working in a feedback loop, which produces forces,which are functions of journal displacement and velocity. The simulated hydrodynamic forces are produced according to assumed characteristics matched to the bearing type. The principal benefit of the analog simulation is that nonlinear characteristics of a subsystem are precisely identified and mathematical methods applied for a wide class of problems are checked on the experimental installation.

Parszewski, Z. A.↗

DOE OSTI2020

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Verification and Validation of Flow Confinement Modeling using OpenFAST Source Panel Method

Blockage or flow confinement effects in scale wind and water turbine testing in experimental facilities, including wind and water tunnels and tow tanks, can significantly alter incident aerodynamic and hydrodynamic loads, turbine performance (thrust, torque, power output), and wake dynamics compared to unconfined conditions. As scale model tests are key sources of OpenFAST model validation measurements, directly modeling flow confinement and its effects in OpenFAST enables improved model validat

Kim, Dongyoung [Sandia National Laboratories (SNL-↗

NASA NTRS1982

PAN AIR summary document (version 1.0)

The capabilities and limitations of the panel aerodynamics (PAN AIR) computer program system are summarized. This program uses a higher order panel method to solve boundary value problems involving the Prandtl-Glauert equation for subsonic and supersonic potential flows. Both aerodynamic and hydrodynamic problems can be solved using this modular software which is written for the CDC 6600 and 7600, and the CYBER 170 series computers.

Derbyshire, T.↗

NASA NTRS1959

A Generalized Hydrodynamic-Impact Theory for the Loads and Motions of Deeply Immersed Prismatic Bodies

A theory is derived for determining the loads and motions of a deeply immersed prismatic body. The method makes use of a two-dimensional water-mass variation and an aspect-ratio correction for three-dimensional flow. The equations of motion are generalized by using a mean value of the aspect-ratio correction and by assuming a variation of the two-dimensional water mass for the deeply immersed body. These equations lead to impact coefficients that depend on an approach parameter which, in turn, depends upon the initial trim and flight-path angles. Comparison of experiment with theory is shown at maximum load and maximum penetration for the flat-bottom (0 deg dead-rise angle) model with bean-loading coefficients from 36.5 to 133.7 over a wide range of initial conditions. A dead-rise angle correction is applied and maximum-load data are compared with theory for the case of a model with 300 dead-rise angle and beam-loading coefficients from 208 to 530.

Markey, Melvin F.↗

NASA NTRS2016

Patterns of Activity in A Global Model of A Solar Active Region

In this work we investigate the global activity patterns predicted from a model active region heated by distributions of nanoflares that have a range of frequencies. What differs is the average frequency of the distributions. The activity patterns are manifested in time lag maps of narrow-band instrument channel pairs. We combine hydrodynamic and forward modeling codes with a magnetic field extrapolation to create a model active region and apply the time lag method to synthetic observations. Our aim is not to reproduce a particular set of observations in detail, but to recover some typical properties and patterns observed in active regions. Our key findings are the following. (1) Cooling dominates the time lag signature and the time lags between the channel pairs are generally consistent with observed values. (2) Shorter coronal loops in the core cool more quickly than longer loops at the periphery. (3) All channel pairs show zero time lag when the line of sight passes through coronal loop footpoints. (4) There is strong evidence that plasma must be re-energized on a timescale comparable to the cooling timescale to reproduce the observed coronal activity, but it is likely that a relatively broad spectrum of heating frequencies are operating across active regions. (5) Due to their highly dynamic nature, we find nanoflare trains produce zero time lags along entire flux tubes in our model active region that are seen between the same channel pairs in observed active regions.

Bradshaw, S. J.↗

DOE OSTI2024

, Experimental Investigation of Electrically Charged Water Droplets Two-Phase Cross-Flow Interactions with Humid

Electrically charged water droplets can capture water vapor molecules in surrounding moist air and promote vapor condensation. This concept unfolds an alternative to mechanical cooling-based air dehumidification systems that separate sensible and latent cooling loads in HVAC applications. The effect of dielectrophoresis, electro-hydrodynamic, and micro-diffusive flows due to highly electrically charged water droplets traveling through moist air were investigated as new methods to reduce the air’s water vapor content. This paper presents new experimental data on humidification and power consumption for a system consisting of highly energized electrosprays in combination with mist eliminators. Electrosprays produced highly charged water spray cones crossing the air perpendicularly and upwardly. The moist air stream was circulated in a horizontal duct at 22C, 80% relative humidity, and 0.24 ms-1 airspeed. Mist eliminators made of stainless wire clothes of 36 and 25 μm pore diameter were used to screen the droplets growing during their flight trajectory. The mist eliminators blocked some of the largest diameter droplets. Still, they could not prevent re-evaporation of the droplets into the bulk air from the droplets deposited on their wire clothes. The equilibrium between the dielectrophoresis condensation from the electrosprays and the droplets’ re-evaporation phenomena resulted in reduced humidification when compared to conventional-type spray evaporative cooling systems. However, to dehumidify the air, the droplets injected by the electrosprays must be separated entirely and promptly removed from the bulk airflow at the end of their flight trajectory.

Electrostatic droplets, Electrospray, humidificati↗

NASA NTRS2003

Scaling of Two-Phase Flows to Partial-Earth Gravity

A report presents a method of scaling, to partial-Earth gravity, of parameters that describe pressure drops and other characteristics of two-phase (liquid/ vapor) flows. The development of the method was prompted by the need for a means of designing two-phase flow systems to operate on the Moon and on Mars, using fluid-properties and flow data from terrestrial two-phase-flow experiments, thus eliminating the need for partial-gravity testing. The report presents an explicit procedure for designing an Earth-based test bed that can provide hydrodynamic similarity with two-phase fluids flowing in partial-gravity systems. The procedure does not require prior knowledge of the flow regime (i.e., the spatial orientation of the phases). The method also provides for determination of pressure drops in two-phase partial-gravity flows by use of a generalization of the classical Moody chart (previously applicable to single-phase flow only). The report presents experimental data from Mars- and Moon-activity experiments that appear to demonstrate the validity of this method.

Hurlbert, Kathryn M.↗

DOE OSTI2023

Quantification of hydraulic trait control on plant hydrodynamics and risk of hydraulic failure within a demographic structured vegetation model in a tropical forest (FATES–HYDRO V1.0)

Abstract. Vegetation plays a key role in the global carbon cycle and thus is an important component within Earth system models (ESMs) that project future climate. Many ESMs are adopting methods to resolve plant size and ecosystem disturbance history, using vegetation demographic models. These models make it feasible to conduct more realistic simulation of processes that control vegetation dynamics. Meanwhile, increasing understanding of the processes governing plant water use, and ecosystem responses to drought in particular, has led to the adoption of dynamic plant water transport (i.e., hydrodynamic) schemes within ESMs. However, the extent to which variations in plant hydraulic traits affect both plant water stress and the risk of mortality in trait-diverse tropical forests is understudied. In this study, we report on a sensitivity analysis of an existing hydrodynamic scheme (HYDRO) model that is updated and incorporated into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (FATES–HYDRO V1.0). The size- and canopy-structured representation within FATES is able to simulate how plant size and hydraulic traits affect vegetation dynamics and carbon–water fluxes. To better understand this new model system, and its functionality in tropical forest systems in particular, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We assembled 942 observations of plant hydraulic traits on 306 tropical plant species for stomata, leaves, stems, and roots and determined the best-fit statistical distribution for each trait, which was used in model parameter sampling to assess the parametric sensitivity. We showed that, for simulated leaf water potential and loss of hydraulic conductivity across different plant organs, the four most important traits were associated with xylem conduit taper (buffers increasing hydraulic resistance with tree height), stomatal sensitivity to leaf water potential, maximum stem hydraulic conductivity, and the partitioning of total hydraulic resistance above vs. belowground. Our analysis of individual ensemble members revealed that trees at a high risk of hydraulic failure and potential tree mortality generally have a lower conduit taper, lower maximum xylem conductivity, lower stomatal sensitivity to leaf water potential, and lower resistance to xylem embolism for stem and transporting roots. We expect that our results will provide guidance on future modeling studies using plant hydrodynamic models to predict the forest responses to droughts and future field campaigns that aim to better parameterize plant hydrodynamic models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A Gigaparsec-scale Hydrodynamic Volume Reconstructed with Deep Learning

The next generation of spectroscopic surveys will map the large-scale structure of the Universe at high redshifts (2 ≤ z ≤ 5) using millions of quasar spectra, enabling major advances in constraining both the standard cosmological model and its extensions. Robust cosmological analyses of these data sets require numerical simulations that both cover gigaparsec volumes and resolve features on ∼10 kpc scales and smaller. However, running such large-volume, high-resolution hydrodynamic simulations is computationally prohibitive. We present a generative deep learning model that enhances a low-resolution, gigaparsec-scale (960 h −1 Mpc) hydrodynamic simulation using a smaller (80 h −1 Mpc) high-resolution input hydrodynamic simulation as training data. The resulting enhanced simulation reproduces the line-of-sight power spectrum to within ∼10% and the three-dimensional power spectrum at the ∼20% level at intermediate to small scales (k ≲ 2 h Mpc −1 ). Our method shows strong promise for producing realistic simulations for cosmological analyses with current surveys such as the Dark Energy Spectroscopic Instrument and upcoming next-generation experiments, but further improvements are needed to accurately recover the large-scale modes. We publicly release the enhanced hydrodynamic simulation, along with a halo catalog from a companion N-body dark matter simulation to support the calibration of data analysis pipelines for these large-scale surveys.

Convolutional neural networks↗

DOE OSTI2020

Drive-pressure optimization in ramp-wave compression experiments through differential evolution

Ramp-wave dynamic-compression experiments are used to examine quasi-isentropic loading paths in materials. The gradual and continuous increase in pressure created by ramp waves make these types of experiments ideal for studying nonequilibrium material behavior, such as solidification kinetics. In ramp-wave compression experiments, the input drive pressure to the experimental setup may be exerted through one of a number of different mechanisms (e.g., magnetic fields, gas-gun-driven impactors, or high-energy lasers) and is generally required for simulating such experiments. Yet, regardless of the specific mechanism, this drive pressure cannot be measured directly (measurements are generally taken at a location near the back of the experimental setup through a transparent window), leading to an inverse problem where one must determine the drive pressure at the front of the experimental setup (i.e., the input) that corresponds to the particle velocity (the output) measured near the back of the experimental setup. Furthermore, we solve this inverse problem using a heuristic optimization algorithm, known as differential evolution, coupled with a multiphysics, hydrodynamics code that simulates the compression of the experimental setup. By running many rounds of forward simulations of the experimental setup, our optimization process iteratively searches for a drive pressure that is optimized to closely reproduce the experimentally measured particle velocity near the back of the experimental setup. While our optimization methodology requires a significant number of hydrodynamics simulations to be conducted, many of these can be performed in parallel, which greatly reduces the time cost of our methodology. One novel aspect of our method for determining the drive pressure is that it does not require physical modeling of the drive mechanism and can thus be broadly applied to many types of ramp-compression experiments, regardless of the drive mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗