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.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

First Cloud-to-Ground Lightning Timing Study

NASA's LSP, GSDO and other programs use the probability of cloud-to-ground (CG) lightning occurrence issued by the 45th Weather Squadron (45 WS) in their daily and weekly lightning probability forecasts. These organizations use this information when planning potentially hazardous outdoor activities, such as working with fuels, or rolling a vehicle to a launch pad, or whenever personnel will work outside and would be at-risk from lightning. These organizations would benefit greatly if the 45 WS could provide more accurate timing of the first CG lightning strike of the day. The Applied Meteorology Unit (AMU) has made significant improvements in forecasting the probability of lightning for the day, but forecasting the time of the first CG lightning with confidence has remained a challenge. To address this issue, the 45 WS requested the AMU to determine if flow regimes, wind speed categories, or a combination of the two could be used to forecast the timing of the first strike of the day in the Kennedy Space Center (KSC)/Cape Canaveral Air Force Station (CCAFS) lightning warning circles. The data was stratified by various sea breeze flow regimes and speed categories in the surface to 5,000-ft layer. The surface to 5,000-ft layer was selected since that is the layer the 45 WS uses to predict the behavior of sea breeze fronts, which are the dominant influence on the occurrence of first lightning in Florida during the warm season. Due to small data sample sizes after stratification, the AMU could not determine a statistical relationship between flow regimes or speed categories and the time of the first CG strike.. As expected, although the amount and timing of lightning activity varies by time of day based on the flow regimes and speed categories, there are extended tails of low lightning activity making it difficult to specify times when the threat of the first lightning flash can be avoided. However, the AMU developed a graphical user interface with input from the 45 WS that allows forecasters to visualize the climatological frequencies of the timing of the first lightning strike. This tool should contribute directly to the 45 WS goal of improving lightning timing capability for its NASA, US Air Force and commercial customers.

statistical analysis↗
NASA NTRS1991

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗
NASA NTRS1989

Fuel-air mixing and distribution in a direct-injection stratified-charge rotary engine

A three-dimensional model for flows and combustion in reciprocating and rotary engines is applied to a direct-injection stratified-charge rotary engine to identify the main parameters that control its burning rate. It is concluded that the orientation of the six sprays of the main injector with respect to the air stream is important to enhance vaporization and the production of flammable mixture. In particular, no spray should be in the wake of any other spray. It was predicted that if such a condition is respected, the indicated efficiency would increase by some 6 percent at higher loads and 2 percent at lower loads. The computations led to the design of a new injector tip that has since yielded slightly better efficiency gains than predicted.

Abraham, J.↗
NASA NTRS2006

Using Flow Regime Lightning and Sounding Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourne, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (http://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in East Central Florida, especially during the warm season months of May September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Currently, the forecasters create each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent is to improve consistency between forecasters while allowing them to focus on the mesoscale detail of the forecast, ultimately benefiting the end-users of the product. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) in which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities, or number of strikes per specified area. The soundings used to determine the flow regimes were taken at Miami (MIA), Tampa (TBW), and Jacksonville (JAX), FL, and the lightning data for the strike densities came from the National Lightning Detection Network (NLDN). The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at FSU and NWS TAE provided this data and supporting software for the work performed by the AMU.

Lambert, Winifred↗
NASA NTRS2007

Using Flow Regime Lightning and Sounding Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourne, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (httl://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in East Central Florida, especially during the warm season months of May September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Until recently, the forecasters created each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent was to improve consistency between forecasters while allowing them to focus on the mesoscale detail of the forecast. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) in which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities, or number of strikes per specified area. The soundings used to determine the flow regimes were taken at Miami (MIA), Tampa (TBW), and Jacksonville (JAX), FL, and the lightning data for the strike densities came from the National Lightning Detection Network (NLDN). The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at FSU and NWS TAE provided this data and supporting software for the work performed by the AMU.

Lambert, Winifred↗
NASA NTRS1989

On the role of ozone in the stability of Rossby normal modes

The role of ozone in the linear stability of Rossby normal modes is examined in a continuously stratified, extratropical baroclinic atmosphere. The flow is described by coupled equations for the quasi-geostraphic potential vorticity and ozone volume mixing ratio. A perturbation analysis is carried out under the assumption of weak diabatic heating, which is generated by Newtonian cooling and dynamics-ozone interaction. An expression for the propagation and growth characteristics is obtained analytically in terms of the vertically averaged wave activity, which depends explicitly on the wave spatial structure, photochemistry, and basic state distributions of wind, temperature, and ozone mixing ratio. Calculations show that stationary internal modes, whose amplitudes are largest in the stratosphre, are destabilized by dynamics-ozone interaction and Newtonian cooling, with e-folding times on the order of 20-40 days.

Nathan, Terrence R.↗
NASA NTRS2001

Hysteresis and Wavenumber Vacillation in Unstable Baroclinic Flows

Hysteresis and wavenumber vacillation are studied numerically in a weakly stratified quasigeostrophic model. In general, the amplitude of the most unstable wave increases, as the flow becomes more unstable. When the wave becomes saturated, the next longer wave will grow at the expanse of the most unstable wave and becomes the dominant wave. However, once the longwave state is established, it may remain in that regime even as the instability is decreased beyond the threshold where it first developed, thus constituting a hysteresis loop. In a highly unstable case, the flow may not show a preference for any single wave. Instead, the dominant wave aperiodically varies among several long waves. This phenomenon is known as wavenumber vacillation. Hysteresis is further examined in terms of eddy heat flux. It is shown that total eddy heat flux increases as the flow becomes more unstable, but displays a sharp drop when transition to a longer wave occurs. However, in a longwave state, the heat flux always decreases with decreasing instability even pass the threshold when wave transition first occurs.

Chou, Shih-Hung↗
NASA NTRS1980

Axially symmetric motion of a rotating stratified fluid

The spherical modeling of geophysical fluid flow is examined. In particular the extension of some previous work done in spherical geometry to the specific case of interest in the Spacelab atmospheric circulation experiment is discussed. This involves changing the boundary conditions under which the basic equations are to be solved. For simplicity the linear, axially symmetric steady state solution is sought.

Martin, B. J.↗
DOE OSTI2022

Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind

Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
NASA NTRS1988

Siphon flows in isolated magnetic flux tubes

The paper considers steady siphon flows in isolated thin magnetic flux tubes surrounded by field-free gas, with plasma beta greater than or equal to 1, appropriate for conditions in the solar photosphere. The cross-sectional area of the flux tube varies along the tube in response to pressure changes induced by the siphon flow. Consideration is also given to steady isothermal siphon flows in arched magnetic flux tubes in a stratified atmosphere. Applications of the results to intense magnetic flux tubes in the solar photosphere and to the photospheric Evershed flow in a sunspot penumbra are addressed.

Thomas, John H.↗
NASA NTRS1992

Evolution of a magnetic flux tube in two-dimensional penetrative convection

Highly supercritical compressible convection is simulated in a two-dimensional domain in which the upper half is unstable to convection while the lower half is stably stratified. This configuration is an idealization of the layers near the base of the solar convection zone. Once the turbulent flow is well developed, a toroidal magnetic field B sub tor is introduced to the stable layer. The field's evolution is governed by an advection-diffusion-type equation, and the Lorentz force does not significantly affect the flow. After many turnover times the field is stratified such that the absolute value of B sub tor/rho is approximately constant in the convective layer, where rho is density, while in the stable layer this ratio decreases linearly with depth. Consequently most of the magnetic flux is stored in the overshoot layer. The inclusion of rotation leads to travelling waves which transport magnetic flux latitudinally in a manner reminiscent of the migrations seen during the solar cycle.

Jennings, R. L.↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
NASA NTRS1992

Modeling deep convection in the Greenland Sea

The development of deep convective events in the high-latitude ocean is studied using a three-dimensional, coupled ice-ocean model. Oceanic mixing is described according to the level 2.5 turbulence closure scheme in which convection occurs in a continuous way, i.e., convective adjustment is not invoked. The model is forced by strong winds and surface cooling. Strong upwelling at the multilyear ice edge and consequent entrainment of warm Atlantic waters into the mixed layer is produced by winds parallel to the ice edge. Concomitant cooling drives deep convection and produces chimneylike structures. Inclusion of a barotropic mean flow over topography to the model provides important preconditioning and selects the location of deep convection. The most efficient preconditioning occurs at locations where the flow ascends a slope. In a stratified environment similar to the Greenland Sea with a 12 m/s wind the model simulations show that localized deep convection takes place after about 10 days to depths of 1000 m.

Hakkinen, S.↗
NASA NTRS1995

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

Mesoscale model simulations provide insight into the complex jet streak adjustments on 11-12 July 1981 that preceded the first of two significant gravity wave events to have been generated over the Rocky Mountains in Montana. Simulations employing a variety of terrain treatments indicate that prior to wave formation, geostrophic adjustment processes modified the structure of the mid-upper tropospheric jet streak by creating secondary jetlets to the southeast of the polar jet streak in proximity to the gravity wave generation region. This simulated restructuring of the mid-upper tropospheric jet streak is the result of a four stage process. During stage 1, the wind adjusts to the mass field as the jet streak exit region propagates into the inflection point between the upstream trough and downstream ridge in the height field. Stage 2 is initiated as the mass field is forced to adjust to the new ageostrophic wind field created during stage 1. Stage 3 is defined by a second geostrophic adjustment process occurring in a similar manner but to the south and east of the adjustment which occurs during stage 1. A low-level mesoscale jetlet is formed during stage 4 in response to the low-level pressure falls that are established during stage 3. The perturbation of this jetlet, caused by orographically-induced adiabatic and diabatic physical processes, is the likely mechanism responsible for the generation of the first and second episode of observed gravity waves. The dynamics responsible for this wave episode are discussed as differential surface sensible heating inducing an orographically-forced mountain-plains solenoid, resulting in the formation of additional mesoscale jetlets and internal gravity waves. Also discussed is how convective latent heating modifies the numerically simulated terrain-induced internal gravity waves, especially their amplitude and phase velocities, which provide better agreement with those wave characteristics observed in nature. Finally, the three-dimensional linear response of a zonally uniform barotropic flow in a vertically unbounded, continuously stratified, Boussinesq atmosphere which is perturbed from geostrophic equilibrium is investigated.

Lin, Yuh-Lang↗
NASA NTRS2001

MHD Instability and Turbulence in the Tachocline

In this quarter we have begun simulations on the Cray T3E at PSC and we are debugging our code on the TSC. The PSC simulations are examining stratified shear turbulence with a flow-aligned magnetic field and passive tracer particles. We have conducted analysis of neutral simulations to establish a firm basis of comparison. Second-order structure functions have been computed, fit, and compared to theoretical expressions relating the dissipation fields and the structure-function-fit parameters. Agreement with high-Reynolds number observations is excellent, giving us confidence that the lower-Re simulations are relevant to higher-Re flows. We have also evaluated the neutral layer anisotropy.

Werne, Joseph↗
NASA NTRS1990

Siphon flows in isolated magnetic flux tubes. III - The equilibrium path of the flux-tube arch

It is shown how to calculate the equilibrium path of a thin magnetic flux tube in a stratified, nonmagnetic atmosphere when the flux tube contains a steady siphon flow. The equilbrium path of a static thin flux tube in an infinite stratified atmosphere generally takes the form of a symmetric arch of finite width, with the flux tube becoming vertical at either end of the arch. A siphon flow within the flux tube increases the curvature of the arched equilibrium path in order that the net magnetic tension force can balance the inertial force of the flow, which tries to straighten the flux tube. Thus, a siphon flow reduces the width of the arched equilibrium path, with faster flows producing narrower arches. The effect of the siphon flow on the equilibrium path is generally greater for flux tubes of weaker magnetic field strength. Examples of the equilibrium are shown for both isothemal and adiabatic siphon flows in thin flux tubes in an isothermal external atmosphere.

Thomas, John H.↗
NASA NTRS1989

Analysis of rotary engine combustion processes based on unsteady, three-dimensional computations

A new computer code was developed for predicting the turbulent, and chemically reacting flows with sprays occurring inside of a stratified charge rotary engine. The solution procedure is based on an Eulerian Lagrangian approach where the unsteady, 3-D Navier-Stokes equations for a perfect gas mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite volume, Steger-Warming flux vector splitting scheme, and the liquid phase equations are solved in Lagrangian coordinates. Both the details of the numerical algorithm and the finite difference predictions of the combustor flow field during the opening of exhaust and/or intake, and also during fuel vaporization and combustion, are presented.

Raju, M. S.↗