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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SANTA BARBARA CLUSTER COMPARISON TEST WITH DISPH

The Santa Barbara cluster comparison project revealed that there is a systematic difference between entropy profiles of clusters of galaxies obtained by Eulerian mesh and Lagrangian smoothed particle hydrodynamics (SPH) codes: mesh codes gave a core with a constant entropy, whereas SPH codes did not. One possible reason for this difference is that mesh codes are not Galilean invariant. Another possible reason is the problem of the SPH method, which might give too much “protection” to cold clumps because of the unphysical surface tension induced at contact discontinuities. In this paper, we apply the density-independent formulation of SPH (DISPH), which can handle contact discontinuities accurately, to simulations of a cluster of galaxies and compare the results with those with the standard SPH. We obtained the entropy core when we adopt DISPH. The size of the core is, however, significantly smaller than those obtained with mesh simulations and is comparable to those obtained with quasi-Lagrangian schemes such as “moving mesh” and “mesh free” schemes. We conclude that both the standard SPH without artificial conductivity and Eulerian mesh codes have serious problems even with such an idealized simulation, while DISPH, SPH with artificial conductivity, and quasi-Lagrangian schemes have sufficient capability to deal with it.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Novel volumetric method for highly repeatable injection in microchip electrophoresis

A novel injector for microchip electrophoresis (MCE) has been designed and evaluated that achieves very high repeatability of injection volume suitable for quantitative analysis. It eliminates the injection biases in electrokinetic injection and the dependence on pressure and sample properties in hydrodynamic injection. The microfluidic injector, made of poly(dimethylsiloxane) (PDMS), operates similarly to an HPLC injection valve. It contains a channel segment (chamber) with a well-defined volume that serves as an “injection loop”. Using on-chip microvalves, the chamber can be connected to the sample source during the “loading” step, and to the CE separation channel during the “injection” step. Once the valves are opened in the second state, electrophoretic potential is applied to separate the sample. For evaluation and demonstration purposes, the microinjector was connected to a 75 μm ID capillary and UV absorbance detector. For single compounds, a relative standard deviation (RSD) of peak area as low as 1.04% (n = 11) was obtained, and for compound mixtures, RSD as low as 0.40% (n = 4) was observed. Using the same microchip, the performance of this new injection technique was compared to hydrodynamic injection and found to have improved repeatability and less dependence on sample viscosity. Furthermore, a non-radioactive version of the positron-emission tomography (PET) imaging probe, FLT, was successfully separated from its known 3 structurally-similar byproducts with baseline resolution, demonstrating the potential for rapid, quantitative analysis of impurities to ensure the safety of batches of short-lived radiotracers. Both the separation efficiency and injection repeatability were found to be substantially higher when using the novel volumetric injection approach compared to electrokinetic injection (performed in the same chip). This novel microinjector provides a straightforward way to improve the performance of hydrodynamic injection and enables extremely repeatable sample volume injection in MCE. It could be used in any MCE application where volume repeatability is needed, including the quantitation of impurities in pharmaceutical or radiopharmaceutical samples.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Revisiting non-relativistic limits

Here, we show that the full spurionic symmetry of Galilean-invariant field theories can be deduced when those theories are the limits of relativistic parents. Under the limit, the non-relativistic daughter couples to Newton-Cartan geometry together with all of the symmetries advocated in previous work, including the recently revived Milne boosts. Our limit is a covariant version of the usual one, where we start with a gapped relativistic theory with a conserved charge, turn on a chemical potential equal to the rest mass of the lightest charged state, and then zoom in to the low energy sector. This procedure gives a simple physical interpretation for the Milne boosts. Our methods even apply when there is a magnetic moment, which is known to modify the non-relativistic symmetry transformations. We focus on two examples. Free scalars are used to demonstrate the basic procedure, whereas hydrodynamics is used in order to exhibit the power of this approach in a fully dynamical setting, correcting several inaccuracies in the existing literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

A simple method to verify the opacity and equation of state of high-Z plasmas

Accurate opacity and equation of state (EOS) data are important parameters for reliable modeling and understanding many phenomena of high energy density physics. However, they are difficult to obtain either theoretically or experimentally especially for high-Z plasmas. In this paper, we proposed a simple method to verify the opacity and EOS of high-Z plasmas in rarefactive states, which is build on the two key aspects of ablation processes, i.e., the propagation of heatwave into materials and the hydrodynamic expansion of the blow-off plasmas. Based on the integrated analysis of the Au foil burnthrough and foam-confined plasma expansion experiments performed on the Shenguang II and Shenguang III prototype laser facilities, we suggested two semi-empirical multipliers to adjust the opacity and EOS of Au plasma. These adjusted opacity and EOS may then provide a foundation for other radiation hydrodynamic simulations under similar drive conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Weighing galaxy clusters with gas. II. On the origin of hydrostatic mass bias in ΛCDM galaxy clusters

The use of galaxy clusters as cosmological probes hinges on our ability to measure their masses accurately and with high precision. Hydrostatic mass is one of the most common methods for estimating the masses of individual galaxy clusters, which suffer from biases due to departures from hydrostatic equilibrium. Using a large, mass-limited sample of massive galaxy clusters from a high-resolution hydrodynamical cosmological simulation, in this work we show that in addition to turbulent and bulk gas velocities, acceleration of gas introduces biases in the hydrostatic mass estimate of galaxy clusters. In unrelaxed clusters, the acceleration bias is comparable to the bias due to non-thermal pressure associated with merger-induced turbulent and bulk gas motions. In relaxed clusters, the mean mass bias due to acceleration is small (≲ 3%), but the scatter in the mass bias can be reduced by accounting for gas acceleration. Additionally, this acceleration bias is greater in the outskirts of higher redshift clusters where mergers are more frequent and clusters are accreting more rapidly. Since gas acceleration cannot be observed directly, it introduces an irreducible bias for hydrostatic mass estimates. This acceleration bias places limits on how well we can recover cluster masses from future X-ray and microwave observations. We discuss implications for cluster mass estimates based on X-ray, Sunyaev-Zel'dovich effect, and gravitational lensing observations and their impact on cluster cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Event-shape fluctuations and flow correlations in ultra-relativistic heavy-ion collisions

I review recent measurements of a large set of flow observables associated with event-shape fluctuations and collective expansion in heavy ion collisions. First, these flow observables are classified and experiment methods are introduced. The experimental results for each type of observables are then presented and compared to theoretical calculations. A coherent picture of initial condition and collective flow based on linear and non-linear hydrodynamic responses is derived, which qualitatively describe most experimental results. I discuss new types of fluctuation measurements that can further our understanding of the event-shape fluctuations and collective expansion dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Range wide genetic differentiation in the bull kelp Nereocystis luetkeana with a seascape genetic focus on the Salish Sea

Introduction: In temperate regions, one of the most critical determinants of present range-wide genetic diversity was the Pleistocene climate oscillations, the most recent one created by the last glacial maximum (LGM). This study aimed to describe N. luetkeana genetic structure across its entire range (Alaska to California) and test different models of population connectivity within the Salish Sea. This region was colonized after the LGM and has been under increased disturbance in recent decades. Methods: We utilized microsatellite markers to study N. luetkeana genetic diversity at 53 sites across its range. Using higher sampling density in the Salish Sea, we employed a seascape genetics approach and tested isolation by hydrodynamic transport and environment models. Results: At the species distribution scale, we found four main groups of genetic co-ancestry, Alaska; Washington with Vancouver Island’s outer coast and Juan de Fuca Strait; Washington’s inner Salish Sea; and Oregon with California. The highest allelic richness (AR) levels were found in California, near the trailing range edge, although AR was also high in Alaska. The inner Salish Sea region had the poorest diversity across the species distribution. Nevertheless, a pattern of isolation by hydrodynamic transport and environment was supported in this region. Discussion: The levels of allelic, private allele richness and genetic differentiation suggest that during the LGM, bull kelp had both northern and southern glacial refugia in the Prince of Wales Island-Haida Gwaii region and Central California, respectively. Genetic diversity in Northern California sites seems resilient to recent disturbances, whereas the low levels of genetic diversity in the inner Salish Sea are concerning.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1983

Flare loop radiative hydrodynamics. III - Nonlocal radiative transfer effects

The study has three goals. The first is to demonstrate that processes exist whose intrinsic nonlocal nature cannot be represented by local approximations. The second is to elucidate the physical nature and origins of these nonlocal processes. The third is to suggest that the methods and results described here may prove useful in constructing semiempirical models of the chromosphere by means more efficient than trial and error. Matrices are computed that describe the effect of a temperature perturbation at an arbitrary point in the loop on density, hydrogen ionized fraction, total radiative loss rate, and radiative loss rate of selected hydrogen lines and continua at all other points. It is found that the dominant nonlocal radiative transfer effects can be separated into flux divergence coefficient effects and upper level population effects. The former are most important when the perturbation takes place in a region of significant opacity. Upper level population effects arise in both optically thick and thin regions in response to nonlocal density, ionization, and interlocking effects.

Canfield, R. C.↗
NASA NTRS1990

On the numerical solution of the dynamically loaded hydrodynamic lubrication of the point contact problem

The transient analysis of hydrodynamic lubrication of a point-contact is presented. A body-fitted coordinate system is introduced to transform the physical domain to a rectangular computational domain, enabling the use of the Newton-Raphson method for determining pressures and locating the cavitation boundary, where the Reynolds boundary condition is specified. In order to obtain the transient solution, an explicit Euler method is used to effect a time march. The transient dynamic load is a sinusoidal function of time with frequency, fractional loading, and mean load as parameters. Results include the variation of the minimum film thickness and phase-lag with time as functions of excitation frequency. The results are compared with the analytic solution to the transient step bearing problem with the same dynamic loading function. The similarities of the results suggest an approximate model of the point contact minimum film thickness solution.

Lim, Sang G.↗
DOE OSTI2014

Time-resolved characterization and energy balance analysis of implosion core in shock-ignition experiments at OMEGA

Time-resolved temperature and density conditions in the core of shock-ignition implosions have been determined for the first time. The diagnostic method relies on the observation, with a streaked crystal spectrometer, of the signature of an Ar tracer added to the deuterium gas fill. The data analysis confirms the importance of the shell attenuation effect previously noted on time-integrated spectroscopic measurements of thick-wall targets [R. Florido et al., Phys. Rev. E 83, 066408 (2011)]. This effect must be taken into account in order to obtain reliable results. The extracted temperature and density time-histories are representative of the state of the core during the implosion deceleration and burning phases. As a consequence of the ignitor shock launched by the sharp intensity spike at the end of the laser pulse, observed average core electron temperature and mass density reach T ~ 1100 eV and ρ ~ 2 g/cm 3 ; then temperature drops to T ~ 920 eV while density rises to ρ ~ 3.4 g/cm 3 about the time of peak compression. Compared to 1D hydrodynamic simulations, the experiment shows similar maximum temperatures and smaller densities. Simulations do not reproduce all observations. Differences are noted in the heating dynamics driven by the ignitor shock and the optical depth time-history of the compressed shell. Time-histories of core conditions extracted from spectroscopy show that the implosion can be interpreted as a two-stage polytropic process. Furthermore, an energy balance analysis of implosion core suggests an increase in total energy greater than what 1D hydrodynamic simulations predict. This new methodology can be implemented in other ICF experiments to look into implosion dynamics and help to understand the underlying physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

EXPLICIT-IMPLICIT SCHEME FOR RELATIVISTIC RADIATION HYDRODYNAMICS

We propose an explicit-implicit scheme for numerically solving special relativistic radiation hydrodynamic equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, zeroth and first moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion approximation. For an hyperbolic term, of which the timescale is the light crossing time when the flow velocity is comparable to the speed of light, is explicitly solved using an approximate Riemann solver. Source terms describing an exchange of energy and momentum between the matter and the radiation via the gas-radiation interaction are implicitly integrated using an iteration method. The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in optically thick media, where heating/cooling and scattering timescales could be much shorter than the dynamical timescale. We show that our numerical code can pass test problems of one- and two-dimensional radiation energy transport, and one-dimensional radiation hydrodynamics. Our newly developed scheme could be useful for a number of relativistic astrophysical problems. We also discuss how to extend our explicit-implicit scheme to the relativistic radiation magnetohydrodynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SIMULATING THE FORMATION OF MASSIVE PROTOSTARS. I. RADIATIVE FEEDBACK AND ACCRETION DISKS

We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M {sub ⊙}. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2–10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M {sub ⊙} simulations, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M {sub ⊙} simulation shows a star with a mass of 5.48 M {sub ⊙} and a disk of mass 3.3 M {sub ⊙}, while our 100 M {sub ⊙} simulation forms a 28.8 M {sub ⊙} mass star with a 15.8 M {sub ⊙} disk over the course of 41.6 kyr, and our 200 M {sub ⊙} simulation forms a 43.7 M {sub ⊙} star with an 18 M {sub ⊙} disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear to be limited by their own luminosities.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2011

Improved Calibration of Modeled Discharge and Storage Change in the Atchafalaya Floodplain Using SAR Interferometry

This study focuses on the feasibility of using SAR interferometry to support 2D hydrodynamic model calibration and provide water storage change in the floodplain. Two-dimensional (2D) flood inundation modeling has been widely studied using storage cell approaches with the availability of high resolution, remotely sensed floodplain topography. The development of coupled 1D/2D flood modeling has shown improved calculation of 2D floodplain inundation as well as channel water elevation. Most floodplain model results have been validated using remote sensing methods for inundation extent. However, few studies show the quantitative validation of spatial variations in floodplain water elevations in the 2D modeling since most of the gauges are located along main river channels and traditional single track satellite altimetry over the floodplain are limited. Synthetic Aperture Radar (SAR) interferometry recently has been proven to be useful for measuring centimeter-scale water elevation changes over the floodplain. In the current study, we apply the LISFLOOD hydrodynamic model to the central Atchafalaya River Basin, Louisiana, during a 62 day period from 1 April to 1 June 2008 using two different calibration schemes for Manning's n. First, the model is calibrated in terms of water elevations from a single in situ gauge that represents a more traditional approach. Due to the gauge location in the channel, the calibration shows more sensitivity to channel roughness relative to floodplain roughness. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. Since SAR interferometry receives strongly scatters in floodplain due to double bounce effect as compared to specular scattering of open water, the calibration shows more dependency to floodplain roughness. An iterative approach is used to determine the best-fit Manning's n for the two different calibration approaches. Results suggest similar floodplain roughness but slightly different channel roughness. However, application of SAR interferometry provides a unique view of the floodplain flow gradients, not possible with a single gauge calibration. These gradients, allow improved computation of water storage change over the 46-day simulation period. Overall, the results suggest that the use of 2D SAR water elevation changes in the Atchafalaya basin offers improved understanding and modeling of floodplain hydrodynamics.

Jung, Hahn Chul↗
DOE OSTI2025

Rippled shock propagation in a laser-driven target at multimegabar pressures

The evolution of non-uniform shocks produced by modulated laser irradiation or surface perturbations is relevant to studies of inertial confinement fusion and material properties at high-energy-density conditions. We present results from an experiment conducted at the OMEGA EP laser facility, where a 300 GPa shock was driven into a fused silica sample with pre-fabricated single-mode surface modulations. Using time-resolved optical velocimetry, we captured the continuous evolution of rippled shock motion, enabling a comprehensive mapping of the spatial amplitude history from formation to phase reversal in a single experiment. Initially, the ablation-driven shock inherits a fraction of the surface modulation amplitude from the sample, which subsequently grows before decaying, ultimately leading to the flattening of the rippled shock and a phase reversal. We find that two-dimensional inviscid hydrodynamic simulation of the experiment is able to qualitatively capture many aspects of the rippled shock evolution but over-predicts the initial amplitude growth. This experimental platform, capable of accommodating varying ripple wavelengths, lays the groundwork for a potential viscometry method at extreme pressures, where viscous effects manifest as differences in shock flattening times between rippled shocks of two distinct wavelengths propagating through the sample.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Two-dimensional simulations of thermonuclear burn in ignition-scale inertial confinement fusion targets under compressed axial magnetic fields

We report for the first time on full 2-D radiation-hydrodynamic implosion simulations that explore the impact of highly compressed imposed magnetic fields on the ignition and burn of perturbed spherical implosions of ignition-scale cryogenic capsules. Using perturbations that highly convolute the cold fuel boundary of the hotspot and prevent ignition without applied fields, we impose initial axial seed fields of 20–100 T (potentially attainable using present experimental methods) that compress to greater than 4 × 10{sup 4} T (400 MG) under implosion, thereby relaxing hotspot areal densities and pressures required for ignition and propagating burn by ∼50%. The compressed field is high enough to suppress transverse electron heat conduction, and to allow alphas to couple energy into the hotspot even when highly deformed by large low-mode amplitudes. This might permit the recovery of ignition, or at least significant alpha particle heating, in submarginal capsules that would otherwise fail because of adverse hydrodynamic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Extension of a coarse grained particle method to simulate heat transfer in fluidized beds

The heat transfer in a gas-solids fluidized bed is simulated with computational fluid dynamic-discrete element method (CFD-DEM) and coarse grained particle method (CGPM). In CGPM fewer numerical particles and their collisions are tracked by lumping several real particles into a computational parcel. Here, the assumption is that the real particles inside a coarse grained particle (CGP) are made from same species and share identical physical properties including density, diameter and temperature. The parcel-fluid convection term in CGPM is calculated using the same method as in DEM. For all other heat transfer mechanisms, we derive in this study mathematical expressions that relate the new heat transfer terms for CGPM to those traditionally derived in DEM. This newly derived CGPM model is verified and validated by comparing the results with CFD-DEM simulation results and experiment data. The numerical results compare well with experimental data for both hydrodynamics and temperature profiles. Finally, the proposed CGPM model can be used for fast and accurate simulations of heat transfer in large scale gas-solids fluidized beds.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

Coupled Geomechanical Modeling to Assess Cap Rock Integrity and Mechanical Fault Stability: Application to Farnsworth Field Unit Project

Abstract The SWP project is located in a mature waterflood undergoing conversion to CO2-WAG operations in the Farnsworth, Texas, USA. Anthropogenic CO2 is sourced from a fertilizer and an ethanol plant. This work utilizes Farnsworth’s full-field, history-matched, compositional hydrodynamic coupled geomechanical model for assessing the impact of stress changes observed through the history matched field life. Production and injection induced stress changes, fault stability and caprock integrity investigations are performed to project the potential for fault reactivation and the loss of caprock integrity under shear failure. A static mechanical earth model (MEM) was constructed for use in transient coupled geomechanical model based on the existing Southwest Regional Partnership (SWP) geological model. The static MEM inherits the stratigraphic and structural features of the geologic model and incorporates additional overburden, underburden, and sideburden formations required to impose mechanical boundary conditions. Mechanical properties were distributed in the 3D MEM through integration of geophysical logs and 3-dimensional seismic elastic inversion properties using a combination of Bayesian and stochastic interpolation methods. These data are further enhanced by lab derived strength and failure criteria for the caprock interval. Additionally, interpreted faults and other geological features were included as part of the static structural framework to fully represent subsurface hydraulic and mechanical systems and appropriately integrate heterogeneity. Two-way coupling of hydrodynamic flow and geomechanical simulations incorporates the Kozeny-Carman relationship for updating permeability and is history matched through primary, secondary (waterflooding) and tertiary (CO2 WAG) recovery periods before performing a twenty year forecast. Two-way simulations are performed to understand the effective stress perturbations imposed by field operations: water injection and the more recent implementation of WAG. Evaluation of Mohr circles with liquid production and field pressure charts as well as slip tendency and distance to failure metrics indicate that neither our faults or cap rock are critically stressed. This paper presents the results of the Farnsworth initial attempts to integrate the seismic driven 3D MEM into coupled hydrodynamic geomechanical history match simulation workflow.

02 PETROLEUM↗