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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS1994

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.↗

NASA NTRS1994

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↗

NASA NTRS1994

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.↗

NASA NTRS1994

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.↗

NASA NTRS1995

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.↗

NASA NTRS1995

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↗

NASA NTRS1995

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.↗

NASA NTRS1995

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.↗

NASA NTRS1997

Impact of Periodic Unsteadiness on Performance and Heat Load in Axial Flow Turbomachines

Results of an analytical and experimental investigation, directed at the understanding of the impact of periodic unsteadiness on the time-averaged flows in axial flow turbomachines, are presented. Analysis of available experimental data, from a large-scale rotating rig (LSRR) (low speed rig), shows that in the time-averaged axisymmetric equations the magnitude of the terms representing the effect of periodic unsteadiness (deterministic stresses) are as large or larger than those due to random unsteadiness (turbulence). Numerical experiments, conducted to highlight physical mechanisms associated with the migration of combustor generated hot-streaks in turbine rotors, indicated that the effect can be simulated by accounting for deterministic stress like terms in the time-averaged mass and energy conservation equations. The experimental portion of this program shows that the aerodynamic loss for the second stator in a 1-1/2 stage turbine are influenced by the axial spacing between the second stator leading edge and the rotor trailing edge. However, the axial spacing has little impact on the heat transfer coefficient. These performance changes are believed to be associated with the change in deterministic stress at the inlet to the second stator. Data were also acquired to quantify the impact of indexing the first stator relative to the second stator. For the range of parameters examined, this effect was found to be of the same order as the effect of axial spacing.

Sharma, Om P.↗

NASA NTRS1996

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.↗

NASA NTRS1996

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↗

NASA NTRS1996

Dynamic Smagorinsky model on anisotropic grids

Large Eddy Simulation (LES) of complex-geometry flows often involves highly anisotropic meshes. To examine the performance of the dynamic Smagorinsky model in a controlled fashion on such grids, simulations of forced isotropic turbulence are performed using highly anisotropic discretizations. The resulting model coefficients are compared with a theoretical prediction (Scotti et al., 1993). Two extreme cases are considered: pancake-like grids, for which two directions are poorly resolved compared to the third, and pencil-like grids, where one direction is poorly resolved when compared to the other two. For pancake-like grids the dynamic model yields the results expected from the theory (increasing coefficient with increasing aspect ratio), whereas for pencil-like grids the dynamic model does not agree with the theoretical prediction (with detrimental effects only on smallest resolved scales). A possible explanation of the departure is attempted, and it is shown that the problem may be circumvented by using an isotropic test-filter at larger scales. Overall, all models considered give good large-scale results, confirming the general robustness of the dynamic and eddy-viscosity models. But in all cases, the predictions were poor for scales smaller than that of the worst resolved direction.

Scotti, A.↗

NASA NTRS1997

The Universe at Moderate Redshift

The report covers the work done in the past year and a wide range of fields including properties of clusters of galaxies; topological properties of galaxy distributions in terms of galaxy types; patterns of gravitational nonlinear clustering process; development of a ray tracing algorithm to study the gravitational lensing phenomenon by galaxies, clusters and large-scale structure, one of whose applications being the effects of weak gravitational lensing by large-scale structure on the determination of q(0); the origin of magnetic fields on the galactic and cluster scales; the topological properties of Ly(alpha) clouds the Ly(alpha) optical depth distribution; clustering properties of Ly(alpha) clouds; and a determination (lower bound) of Omega(b) based on the observed Ly(alpha) forest flux distribution. In the coming year, we plan to continue the investigation of Ly(alpha) clouds using larger dynamic range (about a factor of two) and better simulations (with more input physics included) than what we have now. We will study the properties of galaxies on 1 - 100h(sup -1) Mpc scales using our state-of-the-art large scale galaxy formation simulations of various cosmological models, which will have a resolution about a factor of 5 (in each dimension) better than our current, best simulations. We will plan to study the properties of X-ray clusters using unprecedented, very high dynamic range (20,000) simulations which will enable us to resolve the cores of clusters while keeping the simulation volume sufficiently large to ensure a statistically fair sample of the objects of interest. The details of the last year's works are now described.

Cen, Renyue↗

NASA NTRS1998

A New View of the Dynamics of Reynolds Stress Generation in Turbulent Boundary Layers

The structure of a numerically simulated turbulent boundary layer over a flat plate at Re(theta) = 670 was studied using the invariants of the velocity gradient tensor (Q and R) and a related scalar quantity, the cubic discriminant (D = 27R(exp 2)/4 + Q(exp 3)). These invariants have previously been used to study the properties of the small-scale motions responsible for the dissipation of turbulent kinetic energy. In addition, these scalar quantities allow the local flow patterns to be unambiguously classified according to the terminology proposed by Chong et al. The use of the discriminant as a marker of coherent motions reveals complex, large-scale flow structures that are shown to be associated with the generation of Reynolds shear stress -u'v'(bar). These motions are characterized by high spatial gradients of the discriminant and are believed to be an important part of the mechanism that sustains turbulence in the near-wall region.

Cantwell, Brian J.↗

NASA NTRS1995

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↗

NASA NTRS1998

Boreal Winter Predictions with the GEOS-2 GCM: The Role of Boundary Forcing and Initial Conditions

Ensembles of atmospheric General Circulation Model (GCM) seasonal forecasts and long-term simulations (1980-94) are analyzed to assess the controlling influences of boundary forcing and memory of the initial conditions. Both the forecasts and simulations are carried out with version 2 of the Goddard Earth Observing System (GEOS-2) GCM forced with observed sea surface temperatures (SSTs). While much of the focus is on the seasonal time scale (January- March) and the Pacific North American (PNA) region, we also present results for other regions, shorter time scales, and other known modes of variability in the northern hemisphere extratropics. Forecasts of indices of some of the key large-scale modes of variability show that there is considerable variability in skill between different regions of the Northern Hemisphere. The eastern North Atlantic region has the poorest long lead forecast skill showing no skill beyond about 10 days. Skillful seasonal forecasts are primarily confined to the wave-like ENSO response emanating from the tropical Pacific. In the Northern Hemisphere, this is associated with the well-known Pacific/North American (PNA) pattern. Memory of the initial conditions is the major factor leading to skillful extratropical forecasts of lead time less than one month, while SST forcing is the only factor at the seasonal time scale. SST forcing contributes to skillful forecasts at sub- seasonal time scales only over the PNA region. The GEOS-2 GCM produces average (1980-94) signal to noise ratios which are less than one everywhere in the extratropics, except for the subtropical Pacific where they approach 1.5. When confined to the ENSO years, the maximum signal to noise ratios occur in the PNA region where they exceed three. An assessment of the sampling distribution of the forecasts suggests the model's ENSO response is very likely too weak. These results show some sensitivity to the uncertainties in the estimates of the SST forcing fields. In the North Pacific region, the sensitivity to SST forcing manifests itself primarily as changes in the variability of the PNA response, underscoring the need for an ensemble approach to the seasonal prediction problem.

Chang, Yehui↗

NASA NTRS1995

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.↗

NASA NTRS1997

Vertical Transport Processes for Inert and Scavenged Species: TRACE-A Measurements

The TRACE-A mission of the NASA DC-8 aircraft made a large-scale survey of the tropical and subtropical atmosphere in September and October of 1992. Both In-situ measurements of CO (G. Sachsen NASA Langley) and aerosol size (J. Browell group, NASA Langley) provide excellent data sets with which to constrain vertical transport by planetary boundary layer mixing and deep-cloud cumulus convection. Lidar profiles of aerosol-induced scattering and ozone (also by Bremen) are somewhat require more subtle interpretation as tracers, but the vertical information on layering largely compensates for these complexities. The reason this DC-8 dataset is so useful is that very large areas of biomass burning over Africa and South America provide surface sources of appropriate sizes with which to characterize vertical and horizontal motions; the major limitation of our source description is that biomass burning patterns move considerably every few days, and daily burning inventories are a matter of concurrent, intensive research. We use the Penn State / NCAR MM5 model in an assimilation mode on the synoptic and intercontinental scale, and assess the success it shows in vertical transport descriptions. We find that the general level of emissions suggested by the climatological approach (Will. Has, U. of Montana) appears to be approximately correct, possibly a bit low, for this October, 1992, time period. Vertical transport in planetary boundary layer mixing to 5.5 kin was observed and reproduced in our simulations. Furthermore we find evidence that Blackader "transilient" or matrix-transport scheme is needed, but may require some adaptation in our tracer model: CO seems to exhibit very high values at the top of the planetary boundary layer, a process that stretches the eddy-diffusion parameterization. We will report on progress in improving the deep convective transport of carbon monoxide: the Grail scheme as we used it at 100 kin resolution did not transport enough material to the upper troposphere. We expect to be able to attribute this to either parameterization reasons (inadequacy of this parameterization at the large 100km scale) or other reasons. Nevertheless, the qualitative nature of deep transport by clouds shows up well in the simulations. As for scavengable species, the simulations predict tens of micrograms per standard cubic meter of smoke aerosol in the boundary layer. In a straightforward illustration of our simple bulk-mass scavenging parameterization, to one or two micrograms per standard cubic meter of smoke aerosol in the free troposphere just above the source regions: very high concentrations for the free troposphere. We expect to report on comparisons of these predictions to a variety of observations.

Chatfield, Robert B.↗