SEARCH · Search NASA
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.
For a lifting body of revolution, the exact solutions to the gas dynamical equation are given for the flow field in a stratified atmosphere. The boundary condition is satisfied by using the (linearized) slender-body theory. The solution allows one to make a prediction of the magnitude and the location of the sonic boom in a stratified atmosphere as well as to handle cutoff effects.
The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).
It is shown that many characteristics of the Great Red Spot (GRS) and numerous other features that have been observed on Jupiter can be explained by solitary waves on a horizontally sheared zonal flow in a rotating, stratified atmosphere. Streamline patterns for waves corresponding to combined depression-elevation solitary waves (D-E solitrons) show a strong resemblence to the flow around the GRS. The morphology and flow pattern of the South Tropical Disturbance indicate that it was a D solitron. Numerous spot-like features situated in regions between cloud bands where horizontal shear forces might be expected have the morphology of E solitrons. Restrictions placed on the atmospheric parameters by the model are consistent with available models and observations.
The three-dimensional time-dependent turbulent flow in a neutrally stratified Ekman layer over a smooth flat surface was numerically simulated by directly solving the Navier-Stokes equations. Issues addressed using the direct numerical simulation (DNS) fields include the presence or absence of large-scale coherent structures ('longitudinal' or 'roll' vortices) in neutrally stratified Ekman-layer turbulence, the effects of the horizontal component of the angular velocity vector (i.e., latitude), and implications for models of the PBL. Experimental and DNS profiles are compared.
The flow of an incompressible, stratified fluid over an obstacle will produce an oscillation in which buoyancy is the restoring force, called a gravity wave. For disturbances of this scale, the atmosphere may be treated as incompressible; and even the linear approximation will explain many of the phenomena observed in the lee of mountains. However, nonlinearities arise in two ways: (1) through the large (scaled) size of the mountain, and (2) from dynamically singular levels in the fluid field. These produce a complicated array of phenomena that present hazards to aircraft and to lee surface areas. If there is no dynamic barrier, these waves can penetrate vertically into the middle atmosphere (30-100 km attitude), where recent observations show them to be of a length scale that must involve the Coriolis force in any modeling. At these altitudes, the amplitude of the waves is very large, and the waves are studied with a view to their potential impact on the projected National Aerospace Plane. This paper presents the results of analyses and state-of-the-art numerical simulations, validated where possible by observational data.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.