SEARCH · Search NASA

Results for “subgrid”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Turbulence - The filtering approach

The paper describes a new operatorial approach to the study of turbulence, based on the general algebraic properties of the filtered representations of a turbulence field at different levels. The main results of this analysis is the averaging invariance of the filtered Navier-Stokes eaquations in terms of the generalized central moments, and an algebraic identity that relates the turbulent stresses at different levels. The resolved turbulence is defined, the algebraic consistency rules that relate these resolved quantities to the turbulent stresses at different levels are derived, and their possible uses in subgrid modeling is discussed.

Germano, M.↗

NASA NTRS1992

Turbulent convection with overshooting - Reynolds stress approach

The Reynolds stress formalism is adopted to treat turbulent convection. This methodology is reviewed, and it is suggested that it may prove very useful to treat stellar and accretion disk turbulent convection, as well as in the construction of the subgrid models needed in large eddy simulations. A set of differential equations that yield the mean and turbulent quantities such as convective fluxes, turbulent viscosity, and turbulent conductivity is presented.

Canuto, V. M.↗

NASA NTRS1993

Nonhydrostatic effects in numerical modeling of mesoscale convective systems and baroclinic waves

The present investigation is concerned with the role of convection upon mesoscale modeling results, particularly when the grid resolution becomes small enough that there is not a clear scale separation between the explicitly resolved circulations and the parameterized clouds. In those situations, the vertical accelerations in explicitly resolved circulations become strong enough that the hydrostatic assumption may no longer be valid. These concerns arise from interests in improving mesoscale modeling per se and in improving the subgrid-scale parameterizations in global models. The hydrostatic and the nonhydrostatic options of the Colorado State University Regional Atmospheric Modeling System were used to simulate dry gravity currents in two dimensions, using several different horizontal grid sizes. With horizontal grid intervals of 10 km or less, nonhydrostatic simulations produce wider and colder heads and weaker but wider forced updrafts than do the hydrostatic simulations. Comparing the hydrostatic and nonhydrostatic models show that the difference between the vertical mass fluxes is much less than the difference between the vertical velocities. When the grid is fine enough to resolve the head of the gravity current, horizontal convergence at the gust front extends upwards almost to the head of the cold air. Vertical mass flux in the forced updraft at the front varies with horizontal grid size mainly as a function of the height of the simulated head. For coarser grids, which do not resolve the head, vertical mass flux at all heights decreases with increasing horizontal grid size. A comparison on nonhydrostatic simulations with horizontal grid intervals of 1 km and 2 km illustrates how decreasing the grid size does not necessarily increase the intensity of the resolved circulation. The smaller grid enables the simulated gravity current to entrain a bubble of warm air behind the head, which results in a weaker circulation with a shorter head and weaker updraft.

Cohen, Charles↗

NASA NTRS1993

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The primary objective of this research is to extend current capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first two years of this research have been concentrated on a priori investigations of single-point Probability Density Function (PDF) methods for providing subgrid closures in reacting turbulent flows. In the efforts initiated in the third year, our primary focus has been on performing actual LES by means of PDF methods. The approach is based on assumed PDF methods and we have performed extensive analysis of turbulent reacting flows by means of LES. This includes simulations of both three-dimensional (3D) isotropic compressible flows and two-dimensional reacting planar mixing layers. In addition to these LES analyses, some work is in progress to assess the extent of validity of our assumed PDF methods. This assessment is done by making detailed companions with recent laboratory data in predicting the rate of reactant conversion in parallel reacting shear flows. This report provides a summary of our achievements for the first six months of the third year of this program.

Givi, P.↗

NASA NTRS1993

The structure of supersonic jet flow and its radiated sound

Large-eddy simulation of a supersonic jet is presented with emphasis on capturing the unsteady features of the flow pertinent to sound emission. A high-accuracy numerical scheme is used to solve the filtered, unsteady, compressible Navier-Stokes equations while modelling the subgrid-scale turbulence. For random inflow disturbance, the wave-like feature of the large-scale structure is demonstrated. The large-scale structure was then enhanced by imposing harmonic disturbances to the inflow. The limitation of using the full Navier-Stokes equation to calculate the far-field sound is discussed. Application of Lighthill's acoustic analogy is given with the objective of highlighting the difficulties that arise from the non-compactness of the source term.

Mankbadi, Reda R.↗

NASA NTRS1993

Inertial transfers in the helical decomposition

The directions of the triadic energy transfers assumed in a previous paper agree with the spectral closures, in a similarity range, if the exponent of the power-law energy spectrum is less than unity. The helical interactions showing a strong local energy transfer when the triad is nonlocal, sum up to a reverse cascade unless the spectrum falls off faster than a -7/3 power of the wave number. The energy cascades from each type of helical interaction are calculated for a -5/3 inertial range using the eddy damped quasinormal Markovian model. One type of interaction is responsible for 86 percent of the cascade. The contributions of the two classes of helical interactions to the subgrid-scale eddy viscosity are presented, together with the contributions from the forward and reverse cascading interactions. An application of the assumption on the triadic energy transfers to turbulence under rapid rotation gives a simple argument for the tendency toward nonlinear two-dimensionalization of the flow.

Waleffe, Fabian↗

NASA NTRS1993

Large-eddy simulations of compressible convection on massively parallel computers

We report preliminary implementation of the large-eddy simulation (LES) technique in 2D simulations of compressible convection carried out on the CM-2 massively parallel computer. The convective flow fields in our simulations possess structures similar to those found in a number of direct simulations, with roll-like flows coherent across the entire depth of the layer that spans several density scale heights. Our detailed assessment of the effects of various subgrid scale (SGS) terms reveals that they may affect the gross character of convection. Yet, somewhat surprisingly, we find that our LES solutions, and another in which the SGS terms are turned off, only show modest differences. The resulting 2D flows realized here are rather laminar in character, and achieving substantial turbulence may require stronger forcing and less dissipation.

Xie, Xin↗

NASA NTRS1991

Recipes for degrees of freedom of frequency stability estimators

The Allan variance for an averaging time tau can be estimated either from all available phase samples or from a subgrid of samples with spacing tau. This paper gives a set of computational recipes that yield the variance of both estimators, with less than 2 percent error, for the five power-law components of the classical continuous-time clock noise model.

Greenhall, Charles A.↗

NASA NTRS1993

A description of rainfall interception over large areas

A new scheme is developed for describing interception at spatial scales comparable to the typical resolution of climate models. The scheme is based on the Rutter model of interception and statistical description of the subgrid-scale spatial variability of canopy storage and rainfall. The interception loss simulated by the new scheme is significantly smaller than those simulated by other schemes that do not include considerations for spatial variability. The explanation of this result is partly in the enhancement of spatially averaged canopy drainage due to the large local drainage from the few buckets of large canopy storage. The relative reduction in interception loss simulated by the new scheme may explain the overestimation of interception loss by climate models that do not include the effects of spatial variability on interception processes.

Eltahir, E. A. B.↗

NASA NTRS1993

The implementation and validation of improved land-surface hydrology in an atmospheric general circulation model

New land-surface hydrologic parameterizations are implemented into the NASA Goddard Institute for Space Studies (GISS) General Circulation Model (GCM). These parameterizations are: 1) runoff and evapotranspiration functions that include the effects of subgrid-scale spatial variability and use physically based equations of hydrologic flux at the soil surface and 2) a realistic soil moisture diffusion scheme for the movement of water and root sink in the soil column. A one-dimensional climate model with a complete hydrologic cycle is used to screen the basic sensitivities of the hydrological parameterizations before implementation into the full three-dimensional GCM. Results of the final simulation with the GISS GCM and the new land-surface hydrology indicate that the runoff rate, especially in the tropics, is significantly improved. As a result, the remaining components of the heat and moisture balance show similar improvements when compared to observations. The validation of model results is carried from the large global (ocean and land-surface) scale to the zonal, continental, and finally the regional river basin scales.

Johnson, Kevin D.↗

NASA NTRS1993

Generalized Smagorinsky model for anisotropic grids

The paper revises the Smagorinsky subgrid model to properly account for grid anisotropy on the basis of energy equilibrium considerations in isotropic turbulence. For moderate resolution anisotropies, Deardorff's (1970) estimate involving an equivalent grid scale Delta(eq) = (Delta1 Delta2 Delta3) exp 1/3 is given a rigorous basis. It is argued that the results obtained should be used in conjunction with the dynamic model of Germano et al. (1991) whenever the anisotropy of the test filter differs significantly from that of the basic grid.

Scotti, Alberto↗

NASA NTRS1993

Direct simulation of compressible wall-bounded turbulence

When analyzing many turbulent flows, the effects of compressibility can be neglected. Even some relatively high-speed flows, such as boundary layers generated by a supersonic aircraft, produce turbulent statistics that are similar to those found for the incompressible case. There are other situations, however, in which the non-zero divergence of the turbulence leads to behavior that is fundamentally different from that found at constant density. Examples include flows created by internal combustion engines, hypersonic flight, and supersonic combustion. It is with instances such as these that this project is concerned. In particular, we are interested in the effects of compressibility on turbulence near a smooth solid constant temperature surface; our primary objective is an increased physical understanding that can be used to improve turbulence models of wall-bounded compressible flows. With this in mind, we have begun a direct numerical simulation (DNS) study of turbulence in a plane channel. Because all of the relevant spatial and temporal scales are to be resolved, the simulations require no subgrid scale parameterization. The DNS code developed by Buell to study compressible plane Couette flow was modified to solve the compressible Navier-Stokes equations in the plane channel. The channel was chosen over the Couette flow for two reasons: (1) to avoid using the very large streamwise domains needed to adequately capture the large Couette vortical structures; and (2) to make use of previous experience by considering the compressible version of a well established case and isolate finite Mach number effects by comparing to the incompressible channel.

Coleman, Gary N.↗

NASA NTRS1993

Numerical simulation using vorticity-vector potential formulation

An accurate and efficient computational method is needed for three-dimensional incompressible viscous flows in engineering applications. On solving the turbulent shear flows directly or using the subgrid scale model, it is indispensable to resolve the small scale fluid motions as well as the large scale motions. From this point of view, the pseudo-spectral method is used so far as the computational method. However, the finite difference or the finite element methods are widely applied for computing the flow with practical importance since these methods are easily applied to the flows with complex geometric configurations. However, there exist several problems in applying the finite difference method to direct and large eddy simulations. Accuracy is one of most important problems. This point was already addressed by the present author on the direct simulations on the instability of the plane Poiseuille flow and also on the transition to turbulence. In order to obtain high efficiency, the multi-grid Poisson solver is combined with the higher-order, accurate finite difference method. The formulation method is also one of the most important problems in applying the finite difference method to the incompressible turbulent flows. The three-dimensional Navier-Stokes equations have been solved so far in the primitive variables formulation. One of the major difficulties of this method is the rigorous satisfaction of the equation of continuity. In general, the staggered grid is used for the satisfaction of the solenoidal condition for the velocity field at the wall boundary. However, the velocity field satisfies the equation of continuity automatically in the vorticity-vector potential formulation. From this point of view, the vorticity-vector potential method was extended to the generalized coordinate system. In the present article, we adopt the vorticity-vector potential formulation, the generalized coordinate system, and the 4th-order accurate difference method as the computational method. We present the computational method and apply the present method to computations of flows in a square cavity at large Reynolds number in order to investigate its effectiveness.

Tokunaga, Hiroshi↗

NASA NTRS1993

The helical decomposition and the instability assumption

Direct numerical simulations show that the triadic transfer function T(k,p,q) peaks sharply when q (or p) is much smaller than k. The triadic transfer function T(k,p,q) gives the rate of energy input into wave number k from all interactions with modes of wave number p and q, where k, p, q form a triangle. This observation was thought to suggest that energy is cascaded downscale through non-local interactions with local transfer and that there was a strong connection between large and small scales. Both suggestions were in contradiction with the classical Kolmogorov picture of the energy cascade. The helical decomposition was found useful in distinguishing between kinematically independent interactions. That analysis has gone beyond the question of non-local interaction with local transfer. In particular, an assumption about the statistical direction of triadic energy transfer in any kinematically independent interaction was introduced (the instability assumption). That assumption is not necessary for the conclusions about non-local interactions with local transfer recalled above. In the case of turbulence under rapid rotation, the instability assumption leads to the prediction that energy is transferred in spectral space from the poles of the rotation axis toward the equator. The instability assumption is thought to be of general validity for any type of triad interactions (e.g. internal waves). The helical decomposition and the instability assumption offer detailed information about the homogeneous statistical dynamics of the Navier-Stokes equations. The objective was to explore the validity of the instability assumption and to study the contributions of the various types of helical interactions to the energy cascade and the subgrid-scale eddy-viscosity. This was done in the context of spectral closures of the Direct Interaction or Quasi-Normal type.

Waleffe, Fabian A.↗

NASA NTRS1992

LES versus DNS: A comparative study

We have performed Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of forced isotropic turbulence at moderate Reynolds numbers. The subgrid scale model used in the LES is based on an eddy viscosity which adjusts instantaneously the energy spectrum of the LES to that of the DNS. The statistics of the large scales of the DNS (filtered DNS field or fDNS) are compared to that of the LES. We present results for the transfer spectra, the skewness and flatness factors of the velocity components, the PDF's of the angle between the vorticity and the eigenvectors of the rate of strain, and that between the vorticity and the vorticity stretching tensor. The above LES statistics are found to be in good agreement with those measured in the fDNS field. We further observe that in all the numerical measurements, the trend was for the LES field to be more gaussian than the fDNS field. Future research on this point is planned.

Shtilman, L.↗

NASA NTRS1993

Parallel spatial direct numerical simulations on the Intel iPSC/860 hypercube

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube is documented. The direct numerical simulation approach is used to compute spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows. The feasibility of using the PSDNS on the hypercube to perform transition studies is examined. The results indicate that the direct numerical simulation approach can effectively be parallelized on a distributed-memory parallel machine. By increasing the number of processors nearly ideal linear speedups are achieved with nonoptimized routines; slower than linear speedups are achieved with optimized (machine dependent library) routines. This slower than linear speedup results because the Fast Fourier Transform (FFT) routine dominates the computational cost and because the routine indicates less than ideal speedups. However with the machine-dependent routines the total computational cost decreases by a factor of 4 to 5 compared with standard FORTRAN routines. The computational cost increases linearly with spanwise wall-normal and streamwise grid refinements. The hypercube with 32 processors was estimated to require approximately twice the amount of Cray supercomputer single processor time to complete a comparable simulation; however it is estimated that a subgrid-scale model which reduces the required number of grid points and becomes a large-eddy simulation (PSLES) would reduce the computational cost and memory requirements by a factor of 10 over the PSDNS. This PSLES implementation would enable transition simulations on the hypercube at a reasonable computational cost.

Joslin, Ronald D.↗

NASA NTRS1993

Evapotranspiration and runoff from large land areas: Land surface hydrology for atmospheric general circulation models

A land surface hydrology parameterization for use in atmospheric GCM's is presented. The parameterization incorporates subgrid scale variability in topography, soils, soil moisture and precipitation. The framework of the model is the statistical distribution of a topography-soils index, which controls the local water balance fluxes, and is therefore taken to represent the large land area. Spatially variable water balance fluxes are integrated with respect to the topography-soils index to yield our large topography-soils distribution, and interval responses are weighted by the probability of occurrence of the interval. Grid square averaged land surface fluxes result. The model functions independently as a macroscale water balance model. Runoff ratio and evapotranspiration efficiency parameterizations are derived and are shown to depend on the spatial variability of the above mentioned properties and processes, as well as the dynamics of land surface-atmosphere interactions.

Famiglietti, J. S.↗

NASA NTRS1993

A regional analysis of cloudy mean spherical albedo over the marine stratocumulus region and the tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. The relationships between cloud properties and cloud fraction are studied in order to supplement grid scale parameterizations. The satellite data used is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.5 deg x 2.5 deg latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution, and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakley's theory of reflection for uniform and broken layered cloud and to Kedem, et al.'s findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 m) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al's theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.↗