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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Simplified Computation Of Blade-Surface Coordinates

Coordinates of points on complicated three-dimensional surface of turbo-machine blade computed by method partly simplifying three-dimensional problem to two-dimensional problem. Method implemented in computer program enabling rapid hydrodynamical evaluation and optimization of blade surface prior to definition of final design of, and analysis of stresses in blade. By incorporating program into procedure for designing turbomachine, one reduces time spent on design iterations.

Chen, Wei-Chung↗
NASA NTRS1994

X-ray clusters from a high-resolution hydrodynamic PPM simulation of the cold dark matter universe

A new three-dimensional hydrodynamic code based on the piecewise parabolic method (PPM) is utilized to compute the distribution of hot gas in the standard Cosmic Background Explorer (COBE)-normalized cold dark matter (CDM) universe. Utilizing periodic boundary conditions, a box with size 85 h(exp-1) Mpc, having cell size 0.31 h(exp-1) Mpc, is followed in a simulation with 270(exp 3)=10(exp 7.3) cells. Adopting standard parameters determined from COBE and light-element nucleosynthesis, Sigma(sub 8)=1.05, Omega(sub b)=0.06, we find the X-ray-emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution, and estimated sizes, as well as the evolution of these quantities with redshift. The results, which are compared with those obtained in the preceding paper (Kang et al. 1994a), may be used in conjuction with ROSAT and other observational data sets. Overall, the results of the two computations are qualitatively very similar with regard to the trends of cluster properties, i.e., how the number density, radius, and temeprature depend on luminosity and redshift. The total luminosity from clusters is approximately a factor of 2 higher using the PPM code (as compared to the 'total variation diminishing' (TVD) code used in the previous paper) with the number of bright clusters higher by a similar factor. The primary conclusions of the prior paper, with regard to the power spectrum of the primeval density perturbations, are strengthened: the standard CDM model, normalized to the COBE microwave detection, predicts too many bright X-ray emitting clusters, by a factor probably in excess of 5. The comparison between observations and theoretical predictions for the evolution of cluster properties, luminosity functions, and size and temperature distributions should provide an important discriminator among competing scenarios for the development of structure in the universe.

Bryan, Greg L.↗
NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗
NASA NTRS1996

Numerical Investigation of Two-Phase Flows With Charged Droplets in Electrostatic Field

A numerical method to solve two-phase turbulent flows with charged droplets in an electrostatic field is presented. The ensemble-averaged Navier-Stokes equations and the electrostatic potential equation are solved using a finite volume method. The transitional turbulence field is described using multiple-time-scale turbulence equations. The equations of motion of droplets are solved using a Lagrangian particle tracking scheme, and the inter-phase momentum exchange is described by the Particle-In-Cell scheme. The electrostatic force caused by an applied electrical potential is calculated using the electrostatic field obtained by solving a Laplacian equation and the force exerted by charged droplets is calculated using the Coulombic force equation. The method is applied to solve electro-hydrodynamic sprays. The calculated droplet velocity distributions for droplet dispersions occurring in a stagnant surrounding are in good agreement with the measured data. For droplet dispersions occurring in a two-phase flow, the droplet trajectories are influenced by aerodynamic forces, the Coulombic force, and the applied electrostatic potential field.

Kim, Sang-Wook↗
NASA NTRS2000

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

3D Cell-centered hydrodynamics with subscale closure model and multi-material remap

We extend a higher-order finite volume cell-centered hydrodynamic (CCH) formulation to include an interface-aware subscale closure model and a multi-material remap for simulating 3D compressible hydrodynamic problems within an arbitrary Lagrangian-Eulerian (ALE) framework. This CCH formulation involves a multidirectional approximate Riemann solution using quadratic polynomial reconstructions of the stress tensor and the velocity. At the subscale level, we determine pair-wise material interactions by solving a distinct approximate Riemann problem at the common interface, using the volume of fluids (VOF) method to find the interface. Material interactions are constrained to ensure smooth pressure equilibration among materials. The accuracy and robustness of the ALE method is demonstrated by simulating a suite of 3D Cartesian multi-material problems covering both gas and solid dynamics, where each test case has two or more materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1970

Lagrangian description of warm plasmas

Efforts are described to extend the averaged Lagrangian method of describing small signal wave propagation and nonlinear wave interaction, developed by earlier workers for cold plasmas, to the more general conditions of warm collisionless plasmas, and to demonstrate particularly the effectiveness of the method in analyzing wave-wave interactions. The theory is developed for both the microscopic description and the hydrodynamic approximation to plasma behavior. First, a microscopic Lagrangian is formulated rigorously, and expanded in terms of perturbations about equilibrium. Two methods are then described for deriving a hydrodynamic Lagrangian. In the first of these, the Lagrangian is obtained by velocity integration of the exact microscopic Lagrangian. In the second, the expanded hydrodynamic Lagrangian is obtained directly from the expanded microscopic Lagrangian. As applications of the microscopic Lagrangian, the small-signal dispersion relations and the coupled mode equations are derived for all possible waves in a warm infinite, weakly inhomogeneous magnetoplasma, and their interactions are examined.

Kim, H.↗
NASA NTRS1992

A new scheme for multidimensional line transfer. II - ETLA method in one dimension with application to iron K-alpha lines

The 1D version of the radiation hydrodynamics computer code ALTAIR is discribed. It is an implementation of the equivalent two-level atom method, and also incorporates a variety of other iterative techniques to achieve a fast, accurate solution of the problem of coupled atomic kinetics and radiative transfer. All the techniques are described in sufficient detail to enable their use by others. The application of the program is illustrated by a set of calculations of the iron X-ray emission spectrum from intercloud hot gas that may exist in the broad-line region of active galactic nuclei. These calculations demonstrate that such spectra may be a rich load to mine for information about the hot gas in AGNs, and also that the concept of 'the iron K-alpha line' inadequately depicts the spectrum.

Castor, John I.↗
NASA NTRS1984

Real Surface Effects in Elastohydrodynamic Lubrication

Following a general discussion of real effects in EHL, consideration is given in detail to the role played by stochastic roughness superimposed on the nominal solid boundaries. Particular attention is given to the full-film EHL regime where incipient asperity contact bears a negligible fraction of the load. Flow, from which an averaged Reynolds equation can be formed, is nonetheless modified by the amplitude and texture of the roughness patterns. In describing these effects by means of certain well-defined flow factors, it is found that only two extra parameters are needed - the rms surface height and the ratio of the correlation lengths in the two principal roughness directions. In a lowest order perturbation expansion of these flow factors in powers of the ratio of rms roughness to nominal film thickness, no other properties of the roughness appear. In most cases of practical interest, the factors describing Poiseuille flow are separable into the sum of two single-surface flow factors which means that a combination of a single equivalent rough surface versus an ideal smooth surface can always be found. For flow entrained by slip velocity the factors separate instead into a difference of single-surface factors and it becomes significant which of the two surfaces carries the equivalent roughness. Results are discussed of some applications of the averaged Reynolds equation based on the flow factor method to the EHL line contact problem. Finally, the partial EHL regime is considered where comparable load fractions are carried by the hydrodynamic film and by incipient mechanical contact. An extension of the method into this regime by combining it with asperity contact models appears most encouraging.

John H Tripp↗
NASA NTRS1975

Evolution of rotating interstellar clouds. I - Numerical techniques

A method is described for the numerical calculation of the hydrodynamic evolution of a self-gravitating configuration in two space dimensions with assumed axial symmetry. The calculation is formulated in cylindrical coordinates with respect to a moving Eulerian grid and is solved using explicit hydrodynamics combined with implicit radiative transfer. The physics included is appropriate for calculation of the collapse of a rotating protostellar cloud. The gravitational field is obtained by means of an alternating-direction iterative technique. Numerical tests to demonstrate the correctness of the method are presented for special cases.

Black, D. C.↗
DOE OSTI2026

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗
DOE OSTI2021

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An efficient hybrid downscaling framework to estimate high-resolution river hydrodynamics

Flow depth and velocity are the most important hydrodynamic variables that govern various river functions, including water resources, navigation, sediment transport, and biogeochemical cycling. Existing high-resolution flow depth simulations rely on either computationally expensive river hydrodynamic models (RHMs) or data-driven models with formidable training costs, whereas data-driven modeling of flow velocity has rarely been explored. Here, using the hybrid Low-fidelity, Spatial analysis, and Gaussian process learning (LSG) model, we developed a downscaling approach to construct high-resolution flow depth and velocity from a two-dimensional (2-D) RHM simulation at coarse resolution. The LSG models were trained and tested in an urban watershed in Houston using two different hurricane-driven flood events. The high-resolution (as fine as 30 m resolution) and low-resolution (mostly 1000 m resolution) meshes include 664 724 and 14 536 grid cells, respectively. The results showed that through downscaling, the simulation errors were reduced to less than one-fourth and one-third of the errors of the low-resolution 2-D RHM for flow depth and velocity, respectively. Our analysis further revealed that the dominant uncertainty sources of the downscaled hydrodynamics are different, with flow velocity dominated by the dimensionality reduction error, which we reduced by using a regionalized training procedure. The downscaling approach achieves an 84-fold acceleration in computational time compared to the high-resolution 2-D RHM, making high-fidelity ensemble flood modeling feasible. More importantly, the developed method provides an opportunity to couple large-scale hydrodynamical processes with local physical, chemical, and biological processes in river models.

Tan, Zeli [Pacific Northwest National Laboratory (↗
NASA NTRS1985

Flight vehicle propellant surface motion by numerical method

The marker and cell technique is applied to three idealized hydrodynamic situations of interest to the Space Shuttle External Tank (ET). These applications illustrate the versatility of this technique as well as its usefulness in solving similar problems in the aerospace industry. The codes used are the simplified marker and cell (SMAC) code and the variable grid simplified marker and cell code named ERIE. The examples focus on the dropout height in the LH2 tank, the slosh behavior of propellant in the bottom dome of the ET during thrust termination, both of which are required to accurately determine the amount of usable propellants, and the effect on liquid surface due to LO2 tank wall motions occurring during Space Shuttle lift-off. Comparisons with experimental results and linear theory are also presented.

Mehta, G.↗
DOE OSTI2023

Review of TEAMER Awards for WEC-Sim Support: Preprint

Testing Expertise and Access for Marine Energy Research (TEAMER) is a U.S. Department of Energy Water Power Technologies Office sponsored program, overseen by the Pacific Ocean Energy Trust, which aims to advance the state of marine energy technologies. The program connects technology developers with experts at U.S. facilities, including numerical modeling and analysis facilities, like WEC-Sim. The WEC-Sim facility is supported by the WEC-Sim development team at Sandia National Laboratories and the National Renewable Energy Laboratory. WEC-Sim (Wave Energy Converter SIMulator) is an open-source software for simulating wave energy converters. WEC-Sim can model the multi-body dynamics of devices comprised of bodies, joints, power take-off systems, and mooring systems. Since TEAMER's first round of support in 2020, there have been eighteen TEAMER awards focused on numerical model development in WEC-Sim. TEAMER awards for WEC-Sim support have modeled a wide range of wave energy converter archetypes, including point absorbers, attenuators, oscillating water columns, and many other novel architectures. A wide variety of studies have been conducted, leading to important insights for TEAMER partners and software improvements for WECSim. This article highlights several successful WEC-Sim TEAMER awards. The awards described herein include TEAMER recipients Ocean Motion Technologies, AquaHarmonics, iProTech, East Carolina University, Virginia Tech, Maiden Wave Energy, and the University of Massachusetts Dartmouth. The awards of these seven partners contain a wide range of investigations and cover the creation of baseline hydrodynamic models, PTO modeling, geometry optimization in both boundary element methods and WECSim, and model tuning and validation.

industry support↗
DOE OSTI2021

Mitigating baryonic effects with a theoretical error covariance

ABSTRACT One of the primary sources of uncertainties in modelling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for stage-IV surveys requires knowledge of the matter power spectrum deep in the non-linear regime at the per cent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an augmented covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a ΛCDM model using the extended halo model code HMCode to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on Ωm and σ8 of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

79 ASTRONOMY AND ASTROPHYSICS↗