Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic profiling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

A Deep X-Ray Look at Abell 2142-Viscosity Constraints from Kelvin-Helmholtz Eddies, a Displaced Cool Peak that Makes a Warm Core, and A Possible Plasma Depletion Layer

cold fronts in detail. We find that the southern cold front exhibits well-developed Kelvin–Helmholtz (KH) eddies seen in the sky plane. Comparing their wavelength and amplitude with those in hydrodynamic simulations of cold fronts in viscous gas, and estimating the gas tangential velocity from centripetal acceleration, we constrain the effective viscosity to be at most 1/5 of Spitzer isotropic viscosity, but consistent with full Braginskii anisotropic viscosity for magnetized plasma. While the northwestern front does not show obvious eddies, its shape and the structure of its brightness profile suggest KH eddies seen in projection. The southern cold front continues in a spiral to the center of the cluster, ending with another cold front only 12 kpc from the gas density peak. The cool peak itself is displaced ∼30 kpc from the brightest cluster galaxy (BCG) (the biggest such offset among centrally peaked clusters), while the X-ray emission on a larger scale is still centered on the BCG, indicating that the BCG is at the center of the gravitational potential and the cool gas is sloshing in it. The specific entropy index of the gas in the peak (K≈49 keV sq.cm) makes A2142 a rare “warm core”; apparently the large displacement of the cool peak by sloshing is the reason. Finally, we find a subtle narrow, straight channel with a 10% drop in X-ray brightness, aligned with the southern cold front—possibly a plasma depletion layer in projection.

Wang, Qian H. S.↗
NASA NTRS2022

The Atmospheric Response to High Nonthermal Electron-beam Fluxes in Solar Flares. II. Hydrogen-broadening Predictions for Solar Flare Observations with the Daniel K. Inouye Solar Telescope

Redshifted components of chromospheric emission lines in the hard X-ray impulsive phase of solar flares have recently been studied through their 30 s evolution with the high resolution of the Interface Region Imaging Spectrograph. Radiative-hydrodynamic flare models show that these redshifts are generally reproduced by electron-beam-generated chromospheric condensations. The models produce large ambient electron densities, and the pressure broadening of the hydrogen Balmer series should be readily detected in observations. To accurately interpret the upcoming spectral data of flares with the DKIST, we incorporate nonideal, nonadiabatic line-broadening profiles of hydrogen into the RADYN code. These improvements allow time-dependent predictions for the extreme Balmer line wing enhancements in solar flares. We study two chromospheric condensation models, which cover a range of electron-beam fluxes (1 − 5 × 10^(11) erg s^(−1) cm^(−2)) and ambient electron densities (1 − 60 × 10^(13) cm^(−3)) in the flare chromosphere. Both models produce broadening and redshift variations within 10 s of the onset of beam heating. In the chromospheric condensations, there is enhanced spectral broadening due to large optical depths at Hα, Hβ, and Hγ, while the much lower optical depth of the Balmer series H12−H16 provides a translucent window into the smaller electron densities in the beam-heated layers below the condensation. The wavelength ranges of typical DKIST/ViSP spectra of solar flares will be sufficient to test the predictions of extreme hydrogen wing broadening and accurately constrain large densities in chromospheric condensations.

Adam F. Kowalski↗
NASA NTRS1983

Hydrodynamic models of the cometary atmosphere

Self-consistent transonic multifluid solutions for the dynamical and thermal structure of the two-phase dusty gas atmosphere of a comet containing a dirty H2O-clathrate ice nucleus have been obtained. This was done by solving a simultaneous set of differential equations representing conservation of number density, momentum, and energy, as well as the transfer of solar radiation responsible for the major photolytic processes and nuclear heating. The results presented include velocity and Mach number profiles of heavy species and the velocity of dust within the subsonic and supersonic regions, temperature profiles of the dust and gas in the supersonic region, and radial profiles of electron and ion densities in the collision-dominated region.

Mendis, D. A.↗
NASA NTRS1998

Hydrodynamic Analyses and Evaluation of New Fluid Film Bearing Concepts

Over the past several years, numerical and experimental investigations have been performed on a waved journal bearing. The research work was undertaken by Dr. Florin Dimofte, a Senior Research Associate in the Mechanical Engineering Department at the University of Toledo. Dr. Theo Keith, Distinguished University Professor in the Mechanical Engineering Department was the Technical Coordinator of the project. The wave journal bearing is a bearing with a slight but precise variation in its circular profile such that a waved profile is circumscribed on the inner bearing diameter. The profile has a wave amplitude that is equal to a fraction of the bearing clearance. Prior to this period of research on the wave bearing, computer codes were written and an experimental facility was established. During this period of research considerable effort was directed towards the study of the bearing's stability. The previously developed computer codes and the experimental facility were of critical importance in performing this stability research. A collection of papers and reports were written to describe the results of this work. The attached captures that effort and represents the research output during the grant period.

Keith, Theo G., Jr.↗
DOE OSTI2023

Halted-pendulum Relaxation: Application to White Dwarf Binary Initial Data

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Regulation of the central safety factor and normalized beta under low NBI torque in DIII-D

An algorithm has been designed to simultaneously control the central safety factor (q 0 ) and normalized beta (β N ) while ensuring near-zero torque from the neutral beam injection in DIII-D. Feedback control of q 0 and β N in tokamaks can be beneficial due to the close relationship that these variables have with plasma performance and magneto-hydrodynamic stability. In addition, low neutral-beam-torque conditions are of special interest in present devices because future burning-plasma tokamaks such as ITER will most likely operate at very low plasma rotation. The control synthesis of the algorithm presented in this work is based on a linearized, one-dimensional (1D) model of the current-profile dynamics coupled with a zero-dimensional (0D) plasma-energy balance. The actuators considered are neutral beam injection and electron-cyclotron heating and current drive, and discrete logic determines the neutral-beam injection powers that deliver near-zero torque. Here, the algorithm has been tested in nonlinear, 1D simulations using COTSIM (Control-Oriented Transport SIMulator) and in DIII-D experiments, demonstrating satisfactory performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Retention of Long-period Gas Giant Planets: Type II Migration Revisited

During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

Hydrodynamic air lubricated compliant surface bearing for an automotive gas turbine engine. 1: Journal bearing performance

A 38.1 mm (1.5 inch) diameter Hydresil Compliant Surface Air Lubricated Journal Bearing was designed and tested to obtain bearing performance characteristics at both room temperature and 315 C (600 F). Testing was performed at various speeds up to 60,000 rpm with varying loads. Rotating sensors provided an opportunity to examine the film characteristics of the compliant surface bearing. In addition to providing minimum film thickness values and profiles, many other insights into bearing operation were gained such as the influence of bearing fabrication accuracy and the influence of smooth foil deflection between the bumps.

Ruscitto, D.↗
DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of streamflow in a 30 km long reach spanning 5 years using OpenFOAM 5.x

Abstract. Developing accurate and efficient modeling techniques for streamflow at the tens-of-kilometers spatial scale and multi-year temporal scale is critical for evaluating and predicting the impact of climate- and human-induced discharge variations on river hydrodynamics. However, achieving such a goal is challenging because of limited surveys of streambed hydraulic roughness, uncertain boundary condition specifications, and high computational costs. We demonstrate that accurate and efficient three-dimensional (3-D) hydrodynamic modeling of natural rivers at 30 km and 5-year scales is feasible using the following three techniques within OpenFOAM, an open-source computational fluid dynamics platform: (1) generating a distributed hydraulic roughness field for the streambed by integrating water-stage observation data, a rough wall theory, and a local roughness optimization and adjustment strategy; (2) prescribing the boundary condition for the inflow and outflow by integrating precomputed results of a one-dimensional (1-D) hydraulic model with the 3-D model; and (3) reducing computational time using multiple parallel runs constrained by 1-D inflow and outflow boundary conditions. Streamflow modeling for a 30 km long reach in the Columbia River (CR) over 58 months can be achieved in less than 6 d using 1.1 million CPU hours. The mean error between the modeled and the observed water stages for our simulated CR reach ranges from −16 to 9 cm (equivalent to approximately ±7 % relative to the average water depth) at seven locations during most of the years between 2011 and 2019. We can reproduce the velocity distribution measured by the acoustic Doppler current profiler (ADCP). The correlation coefficients of the depth-averaged velocity between the model and ADCP measurements are in the range between 0.71 and 0.83 at 75 % of the survey cross sections. With the validated model, we further show that the relative importance of dynamic pressure versus hydrostatic pressure varies with discharge variations and topography heterogeneity. Given the model's high accuracy and computational efficiency, the model framework provides a generic approach to evaluate and predict the impacts of climate- and human-induced discharge variations on river hydrodynamics at tens-of-kilometers and decadal scales.

58 GEOSCIENCES↗
DOE OSTI2021

Synthetic Absorption Lines from Simulations of Multiphase Gas in Galactic Winds

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high-resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T$ ~ 10 4 K supernova-driven wind and a cool $T$ ~ 10 7 K cloud of interstellar material, we create mock observables, including the optical depth ($\tau$) and covering fraction ($C$ $f$ ) of six commonly observed ions (Si II , C II , Si IV , C IV , N V , and O VI ) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth and overestimate the "true" value of $\tau$ for ions with low optical depth relative to the simulated data. In conclusion, the empirically estimated covering fractions match our direct calculations very well for the low-ionization ions; for the high-ionization ions, the empirical covering fractions underestimate the directly computed values by up to ~40%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Three-dimensional Core-collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino radiation hydrodynamic simulations of the neutrino-driven phase performed with the Chimera code. We consider a 9.6 M ⊙ zero-metallicity progenitor starting from both 2D and 3D Chimera models and a 10 M ⊙ solar-metallicity progenitor starting from a 2D Chimera model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh–Taylor (R-T) "bullets" determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum 56 Ni velocities of ~1950 and ~1750 km s -1 at shock breakout from 2D and 3D 9.6 M ⊙ Chimera models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher-velocity structures is suppressed when the 2D Chimera model and 3D FLASH model meshes are aligned. The development of faster-growing spherical-bubble structures, as opposed to the slower-growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10 M ⊙ model, which achieves maximum 56 Ni velocities of ~2500 km s -1 at shock breakout.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Measuring the Mean Free Path of H i Ionizing Photons at 3.2 ≤ z ≤ 4.6 with DESI Y1 Quasars

The mean free path of ionizing photons in the intergalactic medium (IGM; λ $^{912}_{mpf}$ ) is a crucial quantity in modeling the ionization state of the IGM and the extragalactic ultraviolet background, and it is widely used in hydrodynamical simulations of galaxies and reionization. We construct the largest quasar spectrum data set to date—12,595 signal-to-noise ratio > 3 spectra—using the Y1 observations of the Dark Energy Spectroscopic Instrument, to make the most precise model-independent measurement of the mean free path at 3.2 ≤ z ≤ 4.6. By stacking the spectra in 17 redshift bins and modeling the Lyman-continuum profile, we get a redshift evolution λ $^{912}_{mpf}$ ∝ ( 1 + z ) − 4.27 at 2 ≤ z ≤ 5, which is much shallower than the previous estimate, λ $^{912}_{mpf}$ ∝ ( 1 + z ) − 5.4 . We then explore the sources of systematic bias, including the choice of the intrinsic quasar continuum, the consideration of Lyman-series opacity and Lyman-limit opacity evolution, and the definition of λ $^{912}_{mpf}$ . Combining our results with estimates of λ $^{912}_{mpf}$ at higher redshifts, we conclude at high confidence that the evolution in λ $^{912}_{mpf}$ steepens at z ≈ 5. We interpret this inflection as the transition from the end of H i reionization to a fully ionized plasma that characterizes the IGM of the past ∼10 billion yr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hydrodynamics of a multicomponent vesicle under strong confinement

Herein we numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Closed coronal structures. V - Gasdynamic models of flaring loops and comparison with SMM observations

A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.

Peres, G.↗
NASA NTRS1991

Hydrodynamical instabilities and mixing in SN 1987A - Two-dimensional simulations of the first 3 months

Results are presented from numerical simulations of the early stages of the explosion of SN 1987A. Using a two-dimensional cylindrical geometry version of a smooth particle hydrodynamics code, the explosion is followed for three months to investigate both the early hydrodynamical instabilities and the effect of the subsequent radioactive decay of Ni-56 and Co-56 with half-lives of 6.1 and 77.8 days, respectively. It is shown that the mixing induced by hydrodynamical instabilities occurring during the first few hours is substantially modified at later time by the radioactive decay of Ni-56 and Co-56. The inner cavity of the expanding supernova remnant fills up with nickel, its decay products thus forming a giant 'nickel bubble'. The peak velocity of the nickel increases by approximately 30 percent after the decays. While these results adequately model the core of the observed Fe line profiles, they fail to reproduce the high velocity wings of the spectra.

Herant, Marc↗
NASA NTRS2016

Simulations of the Mg II K and Ca II 8542 Lines From an Alfvén Wave-Heated Flare Chromosphere

We use radiation hydrodynamic simulations to examine two models of solar flare chromospheric heating: Alfven wave dissipation and electron beam collisional losses. Both mechanisms are capable of strong chromospheric heating, and we show that the distinctive atmospheric evolution in the mid-to-upper chromosphere results in Mg II k-line emission that should be observably different between wave-heated and beam-heated simulations. We also present Ca II 8542 A profiles that are formed slightly deeper in the chromosphere. The Mg II k-line profiles from our wave-heated simulation are quite different from those from a beam-heated model and are more consistent with Interface Region Imaging Spectrograph observations. The predicted differences between the Ca II 8542 A in the two models are small. We conclude that careful observational and theoretical study of lines formed in the mid-to-upper chromosphere holds genuine promise for distinguishing between competing models for chromospheric heating inflares.

Kerr, Graham S.↗
DOE OSTI2023

Impact of toroidal rotation on the resistive ballooning modes in ASDEX Upgrade tokamak

Abstract In this work, we investigate the behavior of instabilities appearing between type-I edge localized modes (ELMs), with increasing neutral beam injection (NBI) power concomitant increase in toroidal rotation, and compare it to the modeling result of the linear magneto-hydrodynamic (MHD) code CASTOR3D. An injection of one NBI beam, increasing toroidal rotation, results in the mode slowing down from 12 kHz to 7 kHz, and its associated radial displacement decreases from 5 mm to 3.5 mm. In addition, modes shift radially outwards towards higher q , decreasing their poloidal mode numbers. The mode velocity is measured to be close to the E × B velocity with significant uncertainties. Through a set of CASTOR3D simulations with varying profiles, resistivity has been identified as the primary contributor to the growth rates. Only a small stabilizing effect due to toroidal rotation has been observed. While experimental results show a decrease of mode frequency with rotation, the opposite trend is observed in modeling. Reasons for discrepancies between modeling and experiment are discussed. Nevertheless, a main contributor to the mode frequency has been identified to be rotation velocity. CASTOR3D classifies modes as resistive ballooning modes as they do not appear unstable in ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗