Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Remarks on the MHD problem of generic magnetospheres and magnetotails

An understanding of the plasma physics problem of MHD flow around a magnetized sphere or cylinder in different Reynolds number regimes is shown to be crucial to magnetospheric and magnetotail studies. A significant amount of knowledge can be acquired from experimental and computational hydrodynamic precedents. It is believed that highly symmetric laminar models of the magnetosphere and magnetotail are unlikely in view of the high Reynolds numbers involved.

Montgomery, D.↗
NASA NTRS1990

One-dimensional cloud fluid model for propagating star formation

The aim of this project was to study the propagation of star formation (SF) with a self-consistent deterministic model for the interstellar gas. The questions of under what conditions does star formation propagate in this model and what are the mechanisms of the propagation are explored. Here, researchers used the deterministic Oort-type cloud fluid model of Scalo and Struck-Marcell (1984, also see the review of Struck-Marcell, Scalo and Appleton 1987). This cloud fluid approach includes simple models for the effects of cloud collisional coalescence or disruption, collisional energy dissipation, and cloud disruption and acceleration as the result of young star winds, HII regions and supernovae. An extensive one-zone parameter study is presented in Struck-Marcell and Scalo (1987). To answer the questions above, researchers carried out one-dimensional calculations for an annulus within a galactic disk, like the so-called solar neighborhood of the galactic chemical evolution. In the calculations the left-hand boundary is set equal to the right hand boundary. The calculation is obviously idealized; however, it is computationally convenient to study the first order effects of propagating star formation. The annulus was treated as if it were at rest, i.e., in the local rotating frame. This assumption may remove some interesting effects of a supersonic gas flow, but was necessary to maintain a numerical stability in the annulus. The results on the one-dimensional propagation of SF in the Oort cloud fluid model follow: (1) SF is propagated by means of hydrodynamic waves, which can be generated by external forces or by the pressure generated by local bursts. SF is not effectively propagated via diffusion or variation in cloud interaction rates without corresponding density and velocity changes. (2) The propagation and long-range effects of SF depend on how close the gas density is to the critical threshold value, i.e., on the susceptibility of the medium.

Titus, Timothy N.↗
NASA NTRS1991

Impact cratering calculations

In the course of carrying out finite difference calculations, it was discovered that for large craters, a previously unrecognized type of crater (diameter) growth occurred which was called lip wave propagation. This type of growth is illustrated for an impact of a 1000 km (2a) silicate bolide at 12 km/sec (U) onto a silicate half-space at earth gravity (1 g). The von Misses crustal strength is 2.4 kbar. The motion at the crater lip associated with this wave type phenomena is up, outward, and then down, similar to the particle motion of a surface wave. It is shown that the crater diameter has grown d/a of approximately 25 to d/a of approximately 4 via lip propagation from Ut/a = 5.56 to 17.0 during the time when rebound occurs. A new code is being used to study partitioning of energy and momentum and cratering efficiency with self gravity for finite-sized objects rather than the previously discussed planetary half-space problems. These are important and fundamental subjects which can be addressed with smoothed particle hydrodynamic (SPH) codes. The SPH method was used to model various problems in astrophysics and planetary physics. The initial work demonstrates that the energy budget for normal and oblique impacts are distinctly different than earlier calculations for silicate projectile impact on a silicate half space. Motivated by the first striking radar images of Venus obtained by Magellan, the effect of the atmosphere on impact cratering was studied. In order the further quantify the processes of meteor break-up and trajectory scattering upon break-up, the reentry physics of meteors striking Venus' atmosphere versus that of the Earth were studied.

Ahrens, Thomas J.↗
NASA NTRS1993

Testing the pressure-confined Ly alpha cloud model

The Ly-alpha absorption line forest, seen in quasar spectra, is generally interpreted as being due to cosmologically distributed 'clouds' of primordial gas. Analyses of the observations reveal that the number distribution can be described by power laws: dN/dz is proportional to (1 + z)(sup gamma) and dN/dN(sub HI) is proportional to N(sub HI)(sup -beyda), where N(sub HI) is the HI column density. The typical values for power law indices range between 2 is approximately less than gamma is approximately less than 2.6 and 1.7 is approximately less than gamma is approximately less than 1.9. One model postulates that the Ly-alpha clouds are optically thin entities, photoionized by the background UV flux, J(sub nu) is proportional to (1 + z)(sup j), and confined by an adiabatically evolving intercloud medium (ICM): P(z) is proportional to (1 + z)(sup 5). Analytic studies of this model suggest that the ensuing Ly-alpha line statistics can account for the observations (in particular, the dN/dz and the dN/dN(sub HI) distributions) if the cloud mass spectrum is a power law dN/dN is proportional to M(sup -delta), delta is approximately 1.9, and j is approximately 4. One of the simplifying assumptions incorporated into these studies is the existence of a large mass range for the clouds at all epochs, the validity of which is questionable. The pressure-confined model is investigated using a 1-D spherically symmetric hydrodynamical code to simulate cloud evolution over the epoch 1.8 less than z less than 6. This enables us to relax many of the assumptions incorporated in the analytic studies.

Babul, Arif↗
NASA NTRS2011

CFD Study of the J2-X LOX Two-Bladed Inducer in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of water flow experiments of the J-2X Oxidizer Turbopump (OTP) 2-bladed inducer with a circumferential groove that were conducted at Marshall Space Flight Center (MSFC). Sub-scale water flow testing results indicate that the circumferential groove greatly reduces synchronous cavitation and subsequent bearing loads at a minimal performance cost, but the energy reappears as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI-Stream CFD program were conducted in order to explore the root physical cause of the HOC. It was found that the axial recirculation pattern in the circumferential groove causes high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that the root cause of HOC is hydrodynamic. That is, the pearls are always present, even when cavitation is not. Comparisons to ongoing water flow experiments will be discussed, as well as predictions for the full-scale OTP inducer.

Coker, Robert↗
NASA NTRS2017

Aggregates: The Fundamental Building Blocks of Planetesimals?

The initial accretion of primitive asteroids (meteorite parent bodies) from freely-floating nebula particles remains problematic. Traditional growth-by-sticking models in turbulent nebulae encounter a formidable "meter-size barrier" due to both drift and destruction, or even a mmtocmsize "bouncing" barrier. Even if growth by sticking could somehow breach these barriers (perhaps if the actual sticking or strength is larger than current estimates, which are based on pure ice or pure silicate), turbulent nebulae present further obstacles through the 1-10km size range. On the other hand, nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids. Thus, the intensity of nebula turbulence (or "alpha ") is critical to the entire process. Theoretical understanding of nebula turbulence continues to evolve; while recent models of MRI (magneticallydriven) turbulence favor lowornoturbulence environments, purely hydrodynamic turbulence is making a comeback with three recently discovered mechanisms generating turbulence of moderate which do not rely on magnetic fields at all.

Cuzzi, J. N.↗
NASA NTRS2020

Intracranial Effects of Microgravity: A Prospective Longitudinal MRI Study

Purpose of Research: To investigate the intracranial effects of microgravity by measuring combined changes in intracranial volumetric parameters, pituitary morphologic structure, and aqueductal cerebrospinal fluid (CSF) hydrodynamics relative to spaceflight and to establish a comprehensive model of recovery after return to Earth.

Larry A Kramer↗
DOE OSTI2025

Sediment transport and morphodynamic changes due to upstream variation and sea level rise on the Western Mekong Delta, Vietnam

The hydrological and morphological dynamics across the Vietnam Mekong Delta (VMD) coastline are increasingly shaped by human interventions and climate change, affecting the sustainability of socio-economic environments. This study employs the TELEMAC-MASCARET modeling system (TELEMAC-2D, TOMAWAC, and SISYPHE modules) to simulate hydrodynamics, sediment transport, and morphological changes for a 10-year period under baseline (SC1) and future climate (SC2) scenarios, considering wave climatology of Northeast (NEM) and Southwest monsoons (SWM). The results of sedimentation and erosion were analyzed across five key coastal zones from Ca Mau Cape (KV1) to Ha Tien (KV5). Model results highlight distinct seasonal patterns of sediment redistribution shaped by monsoonal flows, wave action, and seabed morphology. During the NEM, sediment deposition is concentrated near estuarine zones, while the SWM favors offshore sediment accumulation due to perpendicular wave actions, and shallow, gently sloping bathymetry that dissipates wave energy. Coastal erosion predominantly occurs at convex shoreline features subjected to stronger coastal currents, particularly in fine-grained sediment areas. Under future sea level rise and land subsidence scenarios, reduced flow velocities lead to lower bottom shear stress, promoting deposition and mitigating erosion problems despite increased wave energy. Besides that, sediment transport analyses across 14 cross-sections (C-S) reveal a southward trend during the NEM and a northward reversal during the SWM, with upstream sediment reductions primarily influencing near-estuarine zones and exerting minimal impact beyond the mid-coastal regions. Furthermore, these findings advance the understanding of sediment redistribution processes under climate change and seasonal change, and provide a scientific basis for adaptive coastal management and sustainable development.

Climate change↗
DOE OSTI2023

Simple and Sophisticated Models of Taylor's Cylinder Impact

We present a study based on the simplest possible models for Taylor’s cylinder impact problem, in addition to examination of using convolutional neural networks (CNNs) to map cylinder profiles to strength calibrations. We find that the approximate treatments of Taylor and Hawkyard compare well with hydrodynamic simulations using an equivalent assumption of constant flow stress. The CNN models prove to be well suited to successfully infer parameterizations of the Preston-Tonks-Wallace model of plastic deformation based on the deformed profile of an impacted cylinder

36 MATERIALS SCIENCE↗
NASA NTRS1994

The structure of mass-loading shocks. 1: Comets

A new multifluid model to describe the solar wind-cometary ion plasma in the outer cometary coma is derived. This model is distinguished from previous multifluid models in that wave-particle effects are included explicitly. By considering hydrodynamic timescales, self-consistent dissipative terms (which correspond to the spatial diffusion of cometary ions) are derived, thereby allowing the structure of a cometary bow shock to be resolved. In the hypersonic limit, our diffusive multifluid model reveals an attractive connection to the one-fluid model formulated by Biermann et al. (1967), a connection which serves as an important guide when investigating our more complicated nonhypersonic diffusive model. The detailed model consists of cometary ions diffusing in a warm solar wind which is heated by the dissipation of Alfvenic turbulence (in situ and generated by the initial ring beam distribution). The quasi-parallel cometary shock is found to be at least an order of magnitude thicker than a quasi-perpendicular cometary shock. It is furthermore found that for undisturbed solar wind (sonic) Mach numbers in excess of approximately 6, the cometary bow shock (for both perpendicular and parallel cases) is smoothed completely by the cometary ions and no 'proton' or thermal subshock is necessary. Finally, a detailed comparison of the model with the plasma observations made at Halley by Giotto is undertaken. It is found that the observed shock structure, thickness, location, and plasma parameters compare well with the predictions of the theoretical model, particularly in the case of the quasi-parallel shock.

Zank, G. P.↗
NASA NTRS1972

Model development of supersonic trough wind with shocks

The time dependent one dimensional hydrodynamic equations describe the evolution of the thermal plasma flow along closed magnetic field lines outside of the plasmasphere. The convection of the supersonic polar wind onto a closed fieldline results in the assumed formation of collisionless plasma shocks. These shocks move earthward as the field line with its frozen-in plasma remains fixed or contracts with time to smaller L coordinates. The high equatorial plasma temperature (of the order of electron volts) produced by the shock process decreases with time if the flow is isothermal but it will increase if the contraction is under adiabatic conditions. Assuming adiabaticity a peak in the temperature forms at the equator in conjunction with a depression in the ion density. After an initial contraction, if the flux tube drifts to higher L coordinates the direction of the shock motion can be reversed so that the supersonic region will expand along the field line towards the state characterizing the supersonic polar wind. A rapid expansion will lower the equatorial density while the temperature decreases with time under adiabatic but not isothermal conditions.

Grebowsky, J. M.↗
NASA NTRS1978

Fluid mechanics of continuous flow electrophoresis

The following aspects of continuous flow electrophoresis were studied: (1) flow and temperature fields; (2) hydrodynamic stability; (3) separation efficiency, and (4) characteristics of wide gap chambers (the SPAR apparatus). Simplified mathematical models were developed so as to furnish a basis for understanding the phenomena and comparison of different chambers and operating conditions. Studies of the hydrodynamic stability disclosed that a wide gap chamber may be particularly sensitive to axial temperature variations which could be due to uneven heating or cooling. The mathematical model of the separation process includes effects due to the axial velocity, electro-osmotic cross flow and electrophoretic migration, all including the effects of temperature dependent properties.

Saville, D. A.↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗
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 OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2022

Validation of ablation model for polyethylene using pulsed x-ray and proton exposures

The surface erosion of polyethylene is interrogated using pulsed x rays at the Z Machine (Sandia National Laboratories) and with proton beams at the Gamble II generator (Naval Research Laboratory) to validate a coupled model for volumetric thermal ablation, photoionization, finite-rate decomposition, and molecular recombination of radicals. The intense radiation pulses (up to ∼1014W/m2 over tens of nanoseconds) are used to generate one-dimensional vapor flows with low ionization fractions and a simplified geometry compared to typical laser ablation, allowing for evaluation of the model under local thermal equilibrium conditions. Areal momentum carried by the ensuing uniaxial hydrodynamic shock is used to indicate the extent of ablation. The threshold fluence for ablation is found to be in close correspondence with the bulk melt transition, and reasonable agreement with the model is obtained for peak temperatures in polyethylene up to 5500 K and heating rates up to 1011K/s where thermal decomposition reactions are also active.

Moore, N. W. (ORCID:0000000272259674)↗
DOE OSTI2023

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗