Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRATIFIED FLOW”

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.

87 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

SST-TG-P7F4R3200: Decaying Stably-Stratified Turbulence (SST), Initialized Using Taylor-Green Vortices (TG) at Prandtl Number Pr=7, Froude Number Fr=4, Reynolds Number Re=3200

This dataset comprises direct numerical simulations (DNS) of decaying stably-stratified turbulence influenced by a linear background density gradient, initialized using an array of Taylor-Green vortices, extending the Pr=1 simulations performed in [Riley and de Bruyn Kops (2003)](https://doi.org/10.1063/1.1578077). The initial Prandtl, Froude, and Reynolds numbers are (Pr, Fr, Re) = (7, 4, 3200). A total of 15,250 snapshots are recorded at uniform time intervals, each with a spatial resolution of 1280x1280x640 grid points. Four flow variables are associated with each snapshot: the three velocity components (u,v,w) and the perturbed density field (rho) away from the background gradient. All fields are stored in binary format (32-bit little-endian), each with a size of 4 GB, yielding a total dataset size of 244 TB. Further details are referenced in the attached README file, and a current list of publications and associated analysis tools are provided at https://stratified-turbulence.github.io/web/.

42 ENGINEERING↗
DOE OSTI2022

A comprehensive model to capture electrical discharge and spark channel evolution during spark-ignition processes

Emissions reduction through engine efficiency improvements is a priority for automakers who have turned to unconventional engine operation such as highly dilute, boosted, and stratified charge. Given its importance to flame initiation and sustained turbulent flame propagation, reliable and accurate spark ignition models are necessary to design ignition systems that reduce cyclic indicated mean effective pressure (IMEP) variability and increase engine efficiency in these operation modes. In this paper, secondary electric circuit, short-circuit, blowout, and re-strike sub-models are added to the Lagrangian-Eulerian spark ignition (LESI) model to simulate electrical discharge and spark channel elongation in an inert cross-flow combustion vessel. First, the physics of the underlying sub-models are described and the governing equations discussed with the spark channel voltage expression playing a critical role. Then, the experimental and simulations setups are presented. The results section begins with the derivation of spark channel voltage from experimental results. Then, electrical discharge and spark channel elongation simulations in inert flows are carried out using LESI and compared against experimental results. Finally, the results validate the model's ability to accurately predict spark channel elongation, as well as the occurrence of short-circuits, blowouts, re-strikes, and end of discharge.

42 ENGINEERING↗
DOE OSTI2024

Response of salt intrusion in a tidal estuary to regional climatic forcing

Abstract Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Response Of Salt Intrusion In A Tidal Estuary To Regional Climatic Forcing (FVCOM Output for the baseline scenario)

Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels. These findings have been summarized into a journal manuscript. The data included in this information release request are FVCOM model output associated with the numerical simulations described in the journal manuscript.

Liu, Jinliang↗
DOE OSTI2023

Characterization of wind speed and directional shear at the AWAKEN field campaign site

The American wake experiment (AWAKEN) is taking place in northern Oklahoma, USA, close to the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) atmospheric observatory. The planning for the deployment of the instruments in this observational field campaign required an assessment of the wind characteristics of the site. This paper analyzes long-term data collected by instruments at the ARM SGP observatory to characterize the winds near the AWAKEN site. The analysis shows that this site experiences high wind shear and veer events with a large number of nocturnal low-level jets. A total of 7086 low-level jet wind profiles over 6 years are examined and found to be dominant from the south and southeast. Significant nocturnal wind veer is observed, which causes southerly wind near the surface to become westerly wind aloft. By identifying a strong relationship between atmospheric stability and wind shear, the wind shear at the site is predicted using the Monin–Obukhov similarity theory (MOST) and validated with the observational data collected by a scanning Doppler lidar. The results show that wind speed at a height of 91 m, a proxy hub height for wind turbines in this area, can be predicted from data collected at a height of 10 m with a bias of −0.35 and 0.65 m s−1 in unstable and stable atmospheric boundary layers, respectively. The bias of the predicted wind speed is mostly in the region of low wind speed, and wind speed above 5 m s−1 at a height of 91 m can be predicted with a bias of less than 0.2 m s−1, and the limitations of the MOST in predicting winds during the stably stratified boundary layer is well-observed.

17 WIND ENERGY↗
DOE OSTI2023

Propellant Mass Gauging in a Spherical Tank under Micro-Gravity Conditions Using Capacitance Plate Arrays and Machine Learning

Propellant mass gauging under micro-gravity conditions is a challenging task due to the unpredictable position and shape of the fuel body inside the tank. Micro-gravity conditions are common for orbiting satellites and rockets that operate on limited fuel supplies. Capacitance sensors have been investigated for this task in recent years; however, the effect of various positions and shapes of the fuel body is not analyzed in detail. In this paper, we investigate this with various fill types, such as annular, core-annular, and stratified fills at different positions. We compare the performance among several curve-fitting-based approaches and a machine-learning-based approach, the latter of which offers superior performance in estimating the fuel content.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Supernova Hydrodynamics: A Lab-scale Study of the Blast-driven Instability Using High-speed Diagnostics

A novel experimental approach to study the blast-driven instability at a nondiffuse, gaseous interface with a density gradient is presented. Under Euler similarity, this approach enables study of dissipative-scale hydrodynamics relevant to many astrophysical and laboratory high energy density phenomena in a well-resolved manner. The instability is initiated by passing a Taylor–Sedov blast wave originating from a controlled detonation through a perturbed and stably stratified interface between two gases. The facility and driving blast wave are characterized to obtain repeatable conditions and capture large ensembles of time-resolved Mie scattering imaging that show consistent hydrodynamic development. We analyze the instability evolution between different gas pairs to demonstrate the wide range of development and turbulent behavior that may occur between different supernova layers. The mean evolution of the hydrodynamic instability is compared to a buoyancy–drag model that is frequently used to estimate perturbation growth in supernova mixing research. We propose a time delay to this model in order to reproduce the measured instability behavior and demonstrate model robustness in handling flows driven by a time-varying acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Time-resolved particle image velocimetry measurements of the turbulent Richtmyer–Meshkov instability

Experiments are presented on the Richtmyer–Meshkov instability (RMI) with a three-dimensional, multi-mode initial perturbation. The experiments use a vertical shock tube, where a stably stratified interface is formed between air and sulphur hexafluoride (SF 6 ) via counterflow. A perturbation is imposed at the interface by vertical oscillation of the gas column, forming Faraday waves. The interface is accelerated by a Mach 1.17 (in air) shock wave, and the development of the mixing region between the gases is investigated using particle image velocimetry. Following shock acceleration, a reflected shock wave from the bottom of the shock tube interacts with the mixing layer a second time (reshock). The experiment is initialized with both high and low amplitude perturbations to examine the effect of the perturbation amplitude on measured quantities. The instability growth exponent (θ) is determined from the kinetic energy field using the width of the mixing layer and the decay of kinetic energy, which are found to be in agreement when the flow is most strongly excited. A growth exponent of θ ≈ 0.5 is found for all cases except the high-amplitude reshocked regime (where θ ≈ 0.33). High-amplitude experiments exhibit the transitional outer Reynolds number (Re ≡ hh˙/ν >10 4 ) required for mixing transition following the incident shock, and both experiments are elevated well above this threshold following reshock. However, neither set of experiments meet the more stringent requirements proposed by Zhou et al. (Phys. Rev.E, vol. 67, issue 5, 2003) which include the time dependent aspect of the RMI, an observation which is also made when examining the spectra.

42 ENGINEERING↗
DOE OSTI2020

Spark-ignited direct-injection engine combustion systems

A direct-injection stratified charge internal combustion engine includes a combustion cylinder to receive an air-fuel mixture, and an air intake port to inlet air into the combustion cylinder. The direct-injection engine also includes a fuel injector configured to deliver fuel within the cylinder in a spray pattern substantially aligned to a cylinder central axis to create the air-fuel mixture. A spark igniter is located within a path of the spray pattern to ignite combustion of the air-fuel mixture. The direct-injection engine further includes a movable piston defining a lower boundary of the combustion cylinder to contain the combustion of the air-fuel mixture. The piston is configured to include a bowl portion having local geometric features located on an intake port side of the combustion cylinder to redirect fluid flow towards a vortex in fluid communication with a combustion location near the cylinder central axis.

42 ENGINEERING↗
DOE OSTI2022

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗
DOE OSTI2020

Potential Vorticity Mixing in a Tangled Magnetic Field

A theory of potential vorticity (PV) mixing in a disordered (tangled) magnetic field is introduced. The analysis is in the context of β -plane MHD, with a special focus on the physics of momentum transport in the stably stratified, quasi-2D solar tachocline. A physical picture of mean PV evolution by vorticity advection and tilting of magnetic fields is proposed. In the case of weak field perturbations, quasi-linear theory predicts that the Reynolds and magnetic stresses balance as turbulence Alfvénizes for a larger mean magnetic field. Jet formation is explored quantitatively in the mean field–resistivity parameter space. However, since even a modest mean magnetic field leads to large magnetic perturbations for large magnetic Reynolds number, the physically relevant case is that of a strong but disordered field. We show that numerical calculations indicate that the Reynolds stress is modified well before Alfvénization—i.e., before fluid and magnetic energies balance. To understand these trends, a double-average model of PV mixing in a stochastic magnetic field is developed. Calculations indicate that mean-square fields strongly modify Reynolds stress phase coherence and also induce a magnetic drag on zonal flows. The physics of transport reduction by tangled fields is elucidated and linked to the related quench of turbulent resistivity. We propose here a physical picture of the system as a resisto-elastic medium threaded by a tangled magnetic network. Applications of the theory to momentum transport in the tachocline and other systems are discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Blockage and speedup in the proximity of an onshore wind farm: A scanning wind LiDAR experiment

To maximize the profitability of wind power plants, wind farms are often characterized by high wind turbine density leading to operations with reduced turbine spacing. As a consequence, the overall wind farm power capture is hindered by complex flow features associated with flow modifications induced by the various wind turbine rotors. In addition to the generation of wakes, the velocity of the incoming wind field can reduce due to the increased pressure in the proximity of a single turbine rotor (named induction); a similar effect occurs at the wind-farm level (global blockage), which can have a noticeable impact on power production. On the other hand, intra-wind-farm regions featuring increased velocity compared to the freestream (speedups) have also been observed, which can be a source for a potential power boost. To quantify these rotor-induced effects on the incoming wind velocity field, three profiling LiDARs and one scanning wind LiDAR were deployed both before and after the construction of an onshore wind turbine array. The different wind conditions are classified according to the ambient turbulence intensity and streamwise/spanwise spacing among wind turbines. The analysis of the mean velocity field reveals enhanced induction and speedup under stably stratified atmospheric conditions. Additionally, a reduced horizontal area between adjacent turbines has a small impact on the induction zone but increases significantly the speedup between adjacent rotors.

17 WIND ENERGY↗
DOE OSTI2025

Stratified wind from a super-Eddington X-ray binary is slower than expected

Accretion disks in strong gravity ubiquitously produce winds, seen as blueshifted absorption lines in the X-ray band of both stellar mass X-ray binaries (black holes and neutron stars) and supermassive black holes. Some of the most powerful winds (termed Eddington winds) are expected to arise from systems in which radiation pressure is sufficient to unbind material from the inner disk (L ≳ L Edd ). These winds should be extremely fast and carry a large amount of kinetic power, which, when associated with supermassive black holes, would make them a prime contender for the feedback mechanism linking the growth of those black holes with their host galaxies. Here we show the XRISM Resolve spectrum of the galactic neutron star X-ray binary, GX 13+1, which reveals one of the densest winds ever seen in absorption lines. This Compton-thick wind significantly attenuates the flux, making it appear faint, although it is intrinsically more luminous than usual (L ≳ L Edd ). However, the wind is extremely slow, more consistent with the predictions of thermal-radiative winds launched by X-ray irradiation of the outer disk than with the expected Eddington wind driven by radiation pressure from the inner disk. This puts new constraints on the origin of winds from bright accretion flows in binaries, but also highlights the very different origin required for the ultrafast (v ~ 0.3c) winds seen in recent Resolve observations of a supermassive black hole at a similarly high Eddington ratio.

Audard, Marc [Univ. of Geneva, Versoix (Switzerlan↗
DOE OSTI2024

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Does the rotational direction of a wind turbine impact the wake in a stably stratified atmospheric boundary layer?

Stably stratified atmospheric boundary layers are often characterized by a veering wind profile, in which the wind direction changes clockwise with height in the Northern Hemisphere. Wind-turbine wakes respond to this veer in the incoming wind by stretching from a circular shape into an ellipsoid. We investigate the relationship between this stretching and the direction of the turbine rotation by means of large-eddy simulations. Clockwise rotating, counterclockwise rotating, and non-rotating actuator disc turbines are embedded in wind fields of a precursor simulation with no wind veer and in wind fields with a Northern Hemispheric Ekman spiral, resulting in six combinations of rotor rotation and inflow wind condition. The wake strength, extension, width, and deflection depend on the interaction of the meridional component of Ekman spiral with the rotational direction of the actuator disc, whereas the direction of the disc rotation only marginally modifies the wake if no veer is present. The differences result from the amplification or weakening/reversion of the spanwise and the vertical wind components due to the effect of the superposed disc rotation. They are also present in the streamwise wind component of the wake and in the total turbulence intensity. In the case of an counterclockwise rotating actuator disc, the spanwise and vertical wind components increase directly behind the rotor, resulting in the same rotational direction in the whole wake while its strength decreases downwind. In the case of a clockwise rotating actuator disc, however, the spanwise and vertical wind components of the near wake are weakened or even reversed in comparison to the inflow. This weakening/reversion results in a downwind increase in the strength of the flow rotation in the wake or even a different rotational direction in the near wake in comparison to the far wake. The physical mechanism responsible for this difference can be explained by a simple linear superposition of a veering inflow with a Rankine vortex.

17 WIND ENERGY↗