Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large-scale simulation”

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 451 records · Page 25

Local isotropy in distorted turbulent boundary layers at high Reynolds number

This is a report on the continuation of our experimental investigations of the hypothesis of local isotropy in shear flows. This hypothesis, which states that at sufficiently high Reynolds numbers the small-scale structures of turbulent motions are independent of large-scale structures and mean deformations, has been used in theoretical studies of turbulence and computational methods such as large-eddy simulation. Since Kolmogorov proposed his theory, there have been many experiments, conducted in wakes, jets, mixing layers, a tidal channel, and atmospheric and laboratory boundary layers, in which attempts have been made to verify - or refute - the local-isotropy hypothesis. However, a review of the literature over the last five decades indicated that, despite all these experiments in shear flows, there was no consensus in the scientific community regarding this hypothesis, and, therefore, it seemed worthwhile to undertake a fresh experimental investigation into this question.

Saddoughi, Seyed G.↗

Computation of large-scale statistics in decaying isotropic turbulence

We have performed large-eddy simulations of decaying isotropic turbulence to test the prediction of self-similar decay of the energy spectrum and to compute the decay exponents of the kinetic energy. In general, good agreement between the simulation results and the assumption of self-similarity were obtained. However, the statistics of the simulations were insufficient to compute the value of gamma which corrects the decay exponent when the spectrum follows a k(exp 4) wave number behavior near k = 0. To obtain good statistics, it was found necessary to average over a large ensemble of turbulent flows.

Chasnov, Jeffrey R.↗

Remote-sensing based technique to account for sub-grid scale variability of land surface properties

A method has been presented for the representation of sub-grid scale variability of surface properties within a land surface processes model. The method uses remotely-sensed data to directly or indirectly estimate probability density functions (PDF's) or key surface variables. Application of this technique in a coupled land surface-atmosphere model requires only grid-scale values of the variables of interest, obtained from low-resolution satellite imagery or surface/remote sensing data assimilation. The PDF's of each controlling surface property are superimposed on the respective grid-scale values to simulate sub-grid scale heterogeneity. Sensitivity studies will be carried out to ascertain the relative importance of the heterogeneity of several variables, and the degree to which non-linear property-process interactions impact large-scale fluxes.

Crosson, William L.↗

Small-scale behavior in distorted turbulent boundary layers at low Reynolds number

During the last three years we have conducted high- and low-Reynolds-number experiments, including hot-wire measurements of the velocity fluctuations, in the test-section-ceiling boundary layer of the 80- by 120-foot Full-Scale Aerodynamics Facility at NASA Ames Research Center, to test the local-isotropy predictions of Kolmogorov's universal equilibrium theory. This hypothesis, which states that at sufficiently high Reynolds numbers the small-scale structures of turbulent motions are independent of large-scale structures and mean deformations, has been used in theoretical studies of turbulence and computational methods such as large-eddy simulation; however, its range of validity in shear flows has been a subject of controversy. The present experiments were planned to enhance our understanding of the local-isotropy hypothesis. Our experiments were divided into two sets. First, measurements were taken at different Reynolds numbers in a plane boundary layer, which is a 'simple' shear flow. Second, experiments were designed to address this question: will our criteria for the existence of local isotropy hold for 'complex' nonequilibrium flows in which extra rates of mean strain are added to the basic mean shear?

Saddoughi, Seyed G.↗

Voids and constraints on nonlinear clustering of galaxies

Void statistics of the galaxy distribution in the Center for Astrophysics Redshift Survey provide strong constraints on galaxy clustering in the nonlinear regime, i.e., on scales R equal to or less than 10/h Mpc. Computation of high-order moments of the galaxy distribution requires a sample that (1) densely traces the large-scale structure and (2) covers sufficient volume to obtain good statistics. The CfA redshift survey densely samples structure on scales equal to or less than 10/h Mpc and has sufficient depth and angular coverage to approach a fair sample on these scales. In the nonlinear regime, the void probability function (VPF) for CfA samples exhibits apparent agreement with hierarchical scaling (such scaling implies that the N-point correlation functions for N greater than 2 depend only on pairwise products of the two-point function xi(r)) However, simulations of cosmological models show that this scaling in redshift space does not necessarily imply such scaling in real space, even in the nonlinear regime; peculiar velocities cause distortions which can yield erroneous agreement with hierarchical scaling. The underdensity probability measures the frequency of 'voids' with density rho less than 0.2 -/rho. This statistic reveals a paucity of very bright galaxies (L greater than L asterisk) in the 'voids.' Underdensities are equal to or greater than 2 sigma more frequent in bright galaxy samples than in samples that include fainter galaxies. Comparison of void statistics of CfA samples with simulations of a range of cosmological models favors models with Gaussian primordial fluctuations and Cold Dark Matter (CDM)-like initial power spectra. Biased models tend to produce voids that are too empty. We also compare these data with three specific models of the Cold Dark Matter cosmogony: an unbiased, open universe CDM model (omega = 0.4, h = 0.5) provides a good match to the VPF of the CfA samples. Biasing of the galaxy distribution in the 'standard' CDM model (omega = 1, b = 1.5; see below for definitions) and nonzero cosmological constant CDM model (omega = 0.4, h = 0.6 lambda(sub 0) = 0.6, b = 1.3) produce voids that are too empty. All three simulations match the observed VPF and underdensity probability for samples of very bright (M less than M asterisk = -19.2) galaxies, but produce voids that are too empty when compared with samples that include fainter galaxies.

Vogeley, Michael S.↗

Simulated life cycles of persistent anticyclonic anomalies over the North Pacific: Role of synoptic-scale eddies

This study examines the role of synoptic-scale eddies during the development of persistent anticyclonic height anomalies over the central North Pacific in a general circulation model under perpetual January conditions. The General Circulation Model (GCM) replicates the basic characteristics of the evolution of the anomaly patterns found in observations. The life cycle is characterized by the rapid establishment of the major anomaly center and considerably longer maintenance and decay phases, which include the development of downstream anomaly centers. The simulation also shows a realistic evolution of synoptic-scale activity beginning with enhanced activity off the east coast of Asia prior to onset, followed by a northward shift of the Pacific storm track, which lasts throughout the maintenance phase. The initial enhancement of synoptic-scale eddy activity is associated with a large-scale cyclonic anomaly that developes over Siberia several days prior to the onset of the main anticyclonic anomaly over the central North Pacific. The observations, however, show considerable interdecadel variability in the details of the composite onset behavior; it is unclear whether this variability is real or whether it reflects differences in the data assimilation systems. The role of the time mean flow and synoptic-scale eddies in the development of the persistent Pacific anomalies is studied within the context of a kinetic energy budget in which the flow is decomposed into the time-mean, low-frequency (timescales longer than 10 days), and synoptic (timescales less than 6 days) components. The budget, which is carried out for the simulation at 500 mb, shows that the initial growth of the persistent anticyclonic anomalies is associated with barotropic conversions of energy, with approximately equal contributions coming from the mean flow and the synoptic-scale eddies. After onset the barotropic conversion from the mean flow dominates, whereas the decay phase is associated with baroclinic processes within the low-frequency flow.

Higgins, R. W.↗

Topological analysis of the CfA redshift survey

We study the topology of large-scale structure in the Center for Astrophysics Redshift Survey, which now includes approximately 12,000 galaxies with limiting magnitude m(sub B) is less than or equal to 15.5. The dense sampling and large volume of this survey allow us to compute the topology on smoothing scales from 6 to 20/h Mpc; we thus examine the topology of structure in both 'nonlinear' and 'linear' regimes. On smoothing scales less than or equal to 10/h Mpc this sample has 3 times the number of resolution elements of samples examined in previous studies. Isodensity surface of the smoothed galaxy density field demonstrate that coherent high-density structures and large voids dominate the galaxy distribution. We compute the genus-threshold density relation for isodensity surfaces of the CfA survey. To quantify phase correlation in these data, we compare the CfA genus with the genus of realizations of Gaussian random fields with the power spectrum measured for the CfA survey. On scales less than or equal to 10/h Mpc the observed genus amplitude is smaller than random phase (96% confidence level). This decrement reflects the degree of phase coherence in the observed galaxy distribution. In other words the genus amplitude on these scales is not good measure of the power spectrum slope. On scales greater than 10/h Mpc, where the galaxy distribution is rougly in the 'linear' regime, the genus ampitude is consistent with the random phase amplitude. The shape of the genus curve reflects the strong coherence in the observed structure; the observed genus curve appears broader than random phase (94% confidence level for smoothing scales less than or equal to 10/h Mpc) because the topolgoy is spongelike over a very large range of density threshold. This departre from random phase consistent with a distribution like a filamentary net of 'walls with holes.' On smoothing scales approaching approximately 20/h Mpc the shape of the CfA genus curve is consistent with random phase. There is very weak evidence for a shift of the genus toward a 'bubble-like' topology. To test cosmological models, we compute the genus for mock CfA surveys drawn from large (L greater than or approximately 400/h Mpc) N-body simulations of three variants of the cold dark matter (CDM) cosmogony. The genus amplitude of the 'standard' CDM model (omega h = 0.5, b = 1.5) differs from the observations (96% confidence level) on smoothing scales is less than or approximately 10/h Mpc. An open CDM model (omega h = 0.2) and a CDM model with nonzero cosmological constant (omega h = 0.24, lambda (sub 0) = 0.6) are consistent with the observed genus amplitude over the full range of smoothing scales. All of these models fail (97% confidence level) to match the broadness of the observed genus curve on smoothing scales is less than or equal to 10/h Mpc.

Vogeley, Michael S.↗

Gravitational collapse of small-scale structure as the origin of the Lyman-alpha forest

If gravitational clustering is a hierarchical process, the present large-scale structure of the galaxy distribution implies that structures on smaller scales must have formed at high redshift. We simulate the formation of small-scale structure (average cell mass: (delta) bar m(sub b) = 10(sup 4.2) solar mass) and the evolution of photoionized gas, in the specific case of the CDM + Lambda model. The photoionized gas has a natural minimal scale of collapse, the Jeans scale (m(sub b,J) is approximately equal to 10(exp 9) solar mass). We find that low column density (N(sub HI) is less than or equal to 10(exp 14)/sq cm) lines originate in regions resembling Zel'dovich pancakes, where gas with overdensities in the range 3-30 is enclosed by two shocks but is typically reexpanding at approximately the Hubble velocity. However, higher column density lines stem from more overdense regions where the shocked gas is cooling. We show that this model can probably account for the observed number of lines, their distribution in column density and b-parameters, as well as the cloud physical sizes as observed in gravitationally lensed quasars. We find a redshift evolution that is too steep; however, this may be due to insufficient dynamical range in the simulation or because the specific model is incorrect. The model predicts that high signal-to-noise observations should find systematic deviations from Voigt profiles, mainly in the form of broad wings in the line profiles, and that a fluctuating Gunn-Peterson effect will be detected, which can be modeled as a superposition of weak lines with a wide range of b-parameters.

Cen, Renyue↗

Termination shock response to large-scale solar wind fluctuations

The analysis of data recorded by the Voyager 2 spacecraft indicates the presence of large-scale fluctuations in the solar wind ram pressure on the time scale of tens of days. The amplitude of the fluctuations is highly variable but often lies within a factor of 5 to 10 change from an average or mean value of the ram pressure. Since the spacecraft has presumably not encountered the termination shock yet, these fluctuations should eventually interact with the shock and thereby play a role in determining the shock location. Numerical solutions of the time-dependent gasdynamic equations are used to simulate the response of the termination shock to fluctuations in the solar wind ram pressure comparable to those observed. The primary result of this study is that the maximum shock excursion due to the fluctuations is of the order of 1 AU, which is much smaller than that predicted by other studies. Additional simulations show that the limited movement is due to the fact that the time scale for the termination shock response is substantially larger than the time scale of the fluctuations. It is also shown that the heliopause acts as a barrier for the fluctuations and confines them to the heliosphere.

Steinolfson, R. S.↗

Termination shock response to large-scale solar wind fluctuations

The analysis of data recorded by the Voyager 2 spacecraft indicates the presence of large-scale fluctuations in the solar wind ram pressure on the time scale of tens of days. The amplitude of the fluctuations is highly variable but often lies within a factor of 5 to 10 change from an average or mean value of the ram pressure. Since the spacecraft has presumably not encountered the termination shock yet, these fluctuations should eventually interact with the shock and thereby play a role in determining the shock location. Numerical solutions of the time-dependent gas dynamic equations are used to simulate the response of the termination shock to fluctuations in the solar wind ram pressure comparable to those observed. The primary result of this study is that the maximum shock excursion due to the fluctuations is of the order of 1 AU, which is much smaller than that predicted by other studies. Additional simulations show that the limited movement is due to the fact that the time scale for the termination shock response is substantially larger than the time scale of the fluctuations. It is also shown that the heliopause acts as a barrier for the fluctuations and confines them to the heliosphere.

Steinolfson, R. S.↗

The NASA/MSFC Global Reference Atmospheric Model-1995 version (GRAM-95)

The latest version of the Global Reference Atmospheric Model (GRAM-95) is presented and discussed. GRAM-95 uses the new Global Upper Air Climatic Atlas (GUACA) CD-ROM data set, for 0- to 27-km altitudes. As with earlier versions, GRAM-95 provides complete geographical and altitude coverage for each month of the year. Individual years 1985 to 1991 and a period-of-record (1980 to 1991) can be simulated for the GUACA height range. GRAM-95 uses a specially developed data set, based on Middle Atmosphere Program (MAP) data, for the 20- to 120-km height range, and the NASA Marshall Engineering Thermosphere (MET) model for heights above 90 km. Fairing techniques assure a smooth transition in the overlap height ranges (20 to 27 km and 90 to 120 km). In addition to the traditional GRAM variables of pressure, density, temperature and wind components, GRAM-95 now includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He, and H). A new, variable-scale perturbation model provides both large-scale and small-scale deviations from mean values for the thermodynamic variables and horizontal and vertical wind components. The perturbation model includes new features that simulate intermittency (patchiness) in turbulence and small-scale perturbation fields. The density perturbations and density gradients (density shears) computed by the new model compare favorably in their statistical characteristics with observed density perturbations and density shears from 32 space shuttle reentry profiles. GRAM-95 provides considerable improvement in wind estimates from the new GUACA data set, compared to winds calculated from the geostrophic wind relations previously used in the 0- to 25-km height range. The GRAM-95 code has been put into a more modular form, easier to incorporate as subroutines in other programs (e.g., trajectory codes). A complete user's guide for running the program, plus sample input and output, is provided.

Justus, C. G.↗

Large eddy simulations of compressible turbulent flows

An evaluation of existing models for Large Eddy Simulations (LES) of incompressible turbulent flows has been completed. LES is a computation in which the large, energy-carrying structures to momentum and energy transfer is computed exactly, and only the effect of the smallest scales of turbulence is modeled. That is, the large eddies are computed and the smaller eddies are modeled. The dynamics of the largest eddies are believed to account for most of sound generation and transport properties in a turbulent flow. LES analysis is based on an observation that pressure, velocity, temperature, and other variables are the sum of their large-scale and small-scale parts. For instance, u(i) (velocity) can be written as the sum of bar-u(i) and u(i)-prime, where bar-u(i) is the large-scale and u(i)-prime is the subgrid-scale (SGS). The governing equations for large eddies in compressible flows are obtained after filtering the continuity, momentum, and energy equations, and recasting in terms of Favre averages. The filtering operation maintains only large scales. The effects of the small-scales are present in the governing equations through the SGS stress tensor tau(ij) and SGS heat flux q(i). The mathematical formulation of the Favre-averaged equations of motion for LES is complete.

Porter-Locklear, Freda↗

What determines the direction of minimum variance of the magnetic field fluctuations in the solar wind?

The solar wind is not an isotropic medium; two symmetry axis are provided, first the radial direction (because the mean wind is radial) and second the spiral direction of the mean magnetic field, which depends on heliocentric distance. Observations show very different anisotropy directions, depending on the frequency waveband; while the large-scale velocity fluctuations are essentially radial, the smaller scale magnetic field fluctuations are mostly perpendicular to the mean field direction, which is not the expected linear (WkB) result. We attempt to explain how these properties are related, with the help of numerical simulations.

Grappin, R.↗

Direct Computation of Sound Radiation by Jet Flow Using Large-scale Equations

Jet noise is directly predicted using large-scale equations. The computational domain is extended in order to directly capture the radiated field. As in conventional large-eddy-simulations, the effect of the unresolved scales on the resolved ones is accounted for. Special attention is given to boundary treatment to avoid spurious modes that can render the computed fluctuations totally unacceptable. Results are presented for a supersonic jet at Mach number 2.1.

Mankbadi, R. R.↗

Dynamic Models for LES of Turbulent Front Propagation With a Spectral Method

Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Im, H. G.↗

Towards Numerical Simulations of Trailing-Edge Aeroacoustics

The aeroacoustics of flow-hydrofoil interactions exhibits distinct characteristics depending on the physical length scales involved. In the small-foil (relative to acoustic wavelength) limit characteristic of the noise generated by large-scale vortex shedding at low flow Mach number, the noise calculation is facilitated by the use of the Lighthill analogy in conjunction with a free-space Green's function, in the sense of Curle's formulation. A methodology for computing the vortex-shedding noise using the Curle formulation, including both surface-induced dipole sources and volume quadrupole sources, has been developed.

Wang, Meng↗

A Comparison between the TOPEX/POSEIDON Data and a Global Ocean General Circulation Model during 1992-1993

The TOPEX/POSEIDON altimetric sea level observation during 1992-1993 was used to validate the simulation made by a global ocean general circulation model (OGCM) forced by the daily wind stress and heat flux derived from the National Meteorological Center operational analysis. The OGCM is a version of the modular ocean model with a horizontal resolution of 2 deg longitude and 1 deg latitude and 22 levels in the vertical. The model simulation is compared to the observation at spatial scales of the order of 500 km and larger. Only the temporal variations are examined. The variability is composed primarily of the annual cycle and intraseasonal fluctuations (periods shorter than 100 days). The basic features of the annual cycle are simulated well by the model. Major discrepancies are found in the eastern tropical Pacific, as well as the eastern North Pacific and most of the interior of the North Atlantic. The culprit is suspected to be the inadequate heat forcing and mixing parameterizations of the model. Significant intraseasonal variability is found in the central North Pacific and the Southern Ocean. The simulation is highly correlated with the observation at periods from 20 to 100 days. The spatial scales are larger than 1000 km in many places. These variabilities are apparently the barotropic response of the ocean to wind forcing. The results of the study provide a basis for future assimilation of the data into the OGCM for improved description of the large-scale ocean variabilities.

Chao, Yi↗

The Evolution of Plane Wakes Subjected to Irrotational Strains

Three direct numerical simulations of time-evolving turbulent plane wakes developing in the presence of various irrotational plane strains have been generated. A pseudospectral numerical method with up to 26 million modes is used to solve the equations in a reference frame moving with the irrotational strain. The initial condition for each simulation is taken from a previous turbulent self-similar plane wake direct numerical simulation at a velocity deficit Reynolds number, Re, of about 2000. Three different plane strains (of the same magnitude) are imposed. In the first two simulations the strain is in a plane normal to the streamwise wake direction (the two cases having strain of opposite sign); in the third the wake is compressed in the streamwise direction and stretched in the inhomogeneous cross-stream direction. The two flows that are stretched in the cross-stream direction experience an exponential increase of Re; flow visualization indicates many small-scale vortices with little or no organized large-scale structure. In the flow that is compressed in the cross-stream direction Re decays exponentially and the layer appears to be relaminarizing. The evolution of several turbulence statistics in each of these flows is examined.

Rogers, M. R.↗