Search NASA⌕ Search

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.

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

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1947

On the Theory of the Unsteady Motion of an Airfoil

The paper presents a systematical analysis of the problem of the determination of the unsteady motion about an airfoil moving in an infinite fluid that contains a system of vortices and the determination of the hydrodynamical forces acting on the airfoil. The hydrodynamical problem is reduced to the determination of the function f (xi) which transforms conformally the external region of the airfoil into the interior of a circle. The proposed methods of determining the irrotational motion of a fluid that is produced by any motion of the airfoil are especially simple and effective if the function f (xi) is rational. As an example the flow is determined for the case of an arbitrary motion of an airfoil of the Joukowsky type. The formulas obtained for the determination of the hydrodynamical forces by means of contour integration are similar to those given by S. Chaplygin. These formulas are used to determine the force acting on the airfoil in the cases where the unsteady motion is potential throughout and the circulation about the airfoil is constant and also when the fluid contains a system of vortices. A full discussion is given of the concept of virtual masses together with practical formulas for computing the virtual mass coefficients.

Sedov, L. I.↗
DOE OSTI2016

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.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Shock structure in jets of arbitrary exit geometry

The shock structure and hydrodynamic stability of single and twin jets with arbitrary exit geometry are studied. Both vortex sheets and realistic mean profiles are considered for the jet shear layer. The boundary element method is used to predict the shock spacing and screech tones in a vortex sheet model of a single jet. Results are given for circular, elliptic and rectangular jets and the results compared with analysis and experiment. A finite difference technique is used to examine the effect of finite mixing layer thickness for a single jet. The hydrodynamic stability of twin circular jets is also considered. The variation of the stability characteristics are presented for the vortex sheet case for several separation distances. The formulation of the problem to account for finite thickness effects is also given.

Morris, P. J.↗
DOE OSTI2023

Reconstructing Ly$α$ Fields from Low-resolution Hydrodynamical Simulations with Deep Learning

Hydrodynamical cosmological simulations are a powerful tool for accurately predicting the properties of the intergalactic medium (IGM) and for producing mock skies that can be compared against observational data. However, the need to resolve density fluctuation in the IGM puts a stringent requirement on the resolution of such simulations, which in turn limits the volumes that can be modeled, even on the most powerful supercomputers. In this work, we present a novel modeling method that combines physics-driven simulations with data-driven generative neural networks to produce outputs that are qualitatively and statistically close to the outputs of hydrodynamical simulations employing eight times higher resolution. We show that the Ly$α$ flux field, as well as the underlying hydrodynamic fields, have greatly improved statistical fidelity over a low-resolution simulation. Importantly, the design of our neural network allows for sampling multiple realizations from a given input, enabling us to quantify the model uncertainty. Using test data, we demonstrate that this model uncertainty correlates well with the true error of the Ly$α$ flux prediction. Ultimately, our approach allows for training on small simulation volumes and applying it to much larger ones, opening the door to producing accurate Ly$α$ mock skies in volumes of Hubble size, as will be probed with DESI and future spectroscopic sky surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Mature landfill leachate utilization using a cost-effective hybrid method

Highlights: • Co-digestion of sewage sludge and cavitated mature landfill leachate is studied. • The use of a small CL dose and optimal co-digestion HRT of 20 days is suggested. • Biogas production rate increases when optimal HRT is ensured. • FT-IR/PAS analysis seems to indicate that no new toxic compounds are formed. • The method is recommended as cost-effective for utilizing mature leachate. - Abstract: The main goal of this study was to investigate the co-digestion of sewage sludge and mature landfill leachate pretreated through hydrodynamic cavitation. The process efficiency was analyzed from the aspects of organics removal, biogas production, kinetics and digestate quality. Energy efficiency of the overall treatment was evaluated as well. A dose of hydrodynamically cavitated leachate of 5% v/v was used as a substrate for co-digestion with sewage sludge at a hydraulic retention time of 20 days. This improved the kinetics, biogas production and corresponding yields, as well as digestate quality. Fourier transform infrared photoacoustic spectroscopy analysis seemed to indicate that no new toxic compounds were formed with co-digestion. The investigated two-stage treatment is recommended as a safe and cost-effective method of utilizing mature leachate.

09 BIOMASS FUELS↗
DOE OSTI2017

Modeling electrokinetic flows by consistent implicit incompressible smoothed particle hydrodynamics

We present an efficient implicit incompressible smoothed particle hydrodynamics (I2SPH) discretization of Navier-Stokes, Poisson-Boltzmann, and advection-diffusion equations subject to Dirichlet or Robin boundary conditions. It is applied to model various two and three dimensional electrokinetic flows in simple or complex geometries. The I2SPH's accuracy and convergence are examined via comparison with analytical solutions, grid-based numerical solutions, or empirical models. The new method provides a framework to explore broader applications of SPH in microfluidics and complex fluids with charged objects, such as colloids and biomolecules, in arbitrary complex geometries.

Pan, Wenxiao↗
DOE OSTI2023

Parallel Implicit Hydrodynamics with Material Strength for High Explosive Burn Calculations

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

36 MATERIALS SCIENCE↗
DOE OSTI2018

Fluctuating hydrodynamics of reactive liquid mixtures

Fluctuating hydrodynamics (FHD) provides a framework for modeling microscopic fluctuations in a manner consistent with statistical mechanics and nonequilibrium thermodynamics. This paper presents an FHD formulation for isothermal reactive incompressible liquid mixtures with stochastic chemistry. Fluctuating multispecies mass diffusion is formulated using a Maxwell-Stefan description without assuming a dilute solution, and momentum dynamics is described by a stochastic Navier-Stokes equation for the fluid velocity. We consider a thermodynamically consistent generalization for the law of mass action for non-dilute mixtures and use it in the chemical master equation (CME) to model reactions as a Poisson process. The FHD approach provides remarkable computational efficiency over traditional reaction-diffusion master equation methods when the number of reactive molecules is large, while also retaining accuracy even when there are as few as ten reactive molecules per hydrodynamic cell. We present a numerical algorithm to solve the coupled FHD and CME equations and validate it on both equilibrium and nonequilibrium problems. We simulate a diffusively driven gravitational instability in the presence of an acid-base neutralization reaction, starting from a perfectly flat interface. We demonstrate that the coupling between velocity and concentration fluctuations dominates the initial growth of the instability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗
DOE OSTI2023

Alpha heating of indirect-drive layered implosions on the National Ignition Facility

We report in order to understand how close current layered implosions in indirect-drive inertial confinement fusion are to ignition, it is necessary to measure the level of alpha heating present. To this end, pairs of experiments were performed that consisted of a low-yield tritium–hydrogen–deuterium (THD) layered implosion and a high-yield deuterium–tritium (DT) layered implosion to validate experimentally current simulation-based methods of determining yield amplification. The THD capsules were designed to reduce simultaneously DT neutron yield (alpha heating) and maintain hydrodynamic similarity with the higher yield DT capsules. The ratio of the yields measured in these experiments then allowed the alpha heating level of the DT layered implosions to be determined. The level of alpha heating inferred is consistent with fits to simulations expressed in terms of experimentally measurable quantities and enables us to infer the level of alpha heating in recent high-performing implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2015

Total Dissolved Gas Prediction and Optimization in RIVERWARE

Management and operation of dams within the Columbia River Basin (CRB) provides the region with irrigation, hydropower production, flood control, navigation, and fish passage. These various system-wide demands can require unique dam operations that may result in both voluntary and involuntary spill, thereby increasing tailrace levels of total dissolved gas (TDG) which can be fatal to fish. Appropriately managing TDG levels within the context of the systematic demands requires a predictive framework robust enough to capture the operationally related effects on TDG levels. Development of the TDG predictive methodology herein attempts to capture the different modes of hydro operation, thereby making it a viable tool to be used in conjunction with a real-time scheduling model such as RiverWare. The end result of the effort will allow hydro operators to minimize system-wide TDG while meeting hydropower operational targets and constraints. The physical parameters such as spill and hydropower flow proportions, accompanied by the characteristics of the dam such as plant head levels and tailrace depths, are used to develop the empirically-based prediction model. In the broader study, two different models are developed a simplified and comprehensive model. The latter model incorporates more specific bubble physics parameters for the prediction of tailrace TDG levels. The former model is presented herein and utilizes an empirically based approach to predict downstream TDG levels based on local saturation depth, spillway and powerhouse flow proportions, and entrainment effects. Representative data collected from each of the hydro projects is used to calibrate and validate model performance and the accuracy of predicted TDG uptake. ORNL, in conjunction with IIHR - Hydroscience & Engineering, The University of Iowa, carried out model adjustments to adequately capture TDG levels with respect to each plant while maintaining a generalized model configuration. Validation results indicate excellent model performance with coefficient of determination values exceeding 92% for all sites. This approach enables model extension to an increasingly wider array of hydropower plants, i.e., with the proper data input, TDG uptake can be calculated independent of actual physical component design. The TDG model is used as a module in the systematic optimization framework of RiverWare, a river and reservoir modeling tool used by federal agencies, public utility districts, and other dam owners and operators to forecast, schedule, and manage hydropower assets. The integration and testing of the TDG module within RiverWare, led by University of Colorado's Center for Advanced Decision Support for Water and Environmental Systems (CADSWES), will allow users to generate optimum system schedules based on the minimization of TDG. Optimization analysis and added value will be quantified as system wide reductions in TDG achieved while meeting existing hydropower constraints. Future work includes the development of a method to predict downstream reservoir forebay TDG levels as a function of upstream reservoir tailrace TDG values based on river hydrodynamics, hydro operations, and reservoir characteristics. Once implemented, a holistic model that predicts both TDG uptake and transport will give hydropower operators valuable insight into how system-wide environmental effects can be mitigated while simultaneously balancing stakeholder interests.

13 HYDRO ENERGY↗
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 OSTI2014

A coarse-graining approach for molecular simulation that retains the dynamics of the all-atom reference system by implementing hydrodynamic interactions

We report on a new approach for deriving coarse-grained intermolecular forces that retains the frictional contribution that is often discarded by conventional coarse-graining methods. The approach is tested for water and an aqueous glucose solution, and the results from the new implementation for coarse-grained molecular dynamics simulation show remarkable agreement with the dynamics obtained from reference all-atom simulations. The agreement between the structural properties observed in the coarse-grained and all-atom simulations is also preserved. We discuss how this approach may be applied broadly to any existing coarse-graining method where the coarse-grained models are rigorously derived from all-atom reference systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Modeling UV Radiation Feedback from Massive Stars. I. Implementation of Adaptive Ray-tracing Method and Tests [Modeling UV Radiation Feedback from Massive Stars: I. Implementation of Adaptive Ray Tracing Algorithm]

Here, we present an implementation of an adaptive ray-tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses nonblocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems, including the propagation of radiation in vacuum and the expansions of various types of H II regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $\sim {10}^{3}$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is 12% when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with those from the M 1 closure relation. Although the ART and M 1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.

79 ASTRONOMY AND ASTROPHYSICS↗