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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
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 NTRS1999

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗
NASA NTRS1999

Estimates of Long Term Surface Soil Moisture in the Midwestern U.S. Derived from Satellite Microwave Observations

Soil moisture is a key component of the water and energy balances of the Earth's surface, and has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, soil moisture is often somewhat difficult to measure accurately in both space and time, especially at large spatial scales. Soil moisture is highly variable, and while point measurements are typically quite accurate, subsequent areal averaging of these measurements often leads to large errors. Since remotely sensed land surface observations are already a spatially averaged or areally integrated value, they are a logical input parameter to regional or larger scale land process models. A database of long-term soil moisture was compared to satellite microwave observations over test sites in the Midwestern United States. Ground measurements of average volumetric surface soil moisture in the top ten cm were made bimonthly at 19 locations throughout the state of Illinois. Nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and a temporal frequency over the test area of about 3 nighttime orbits per week. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Because the time of satellite coverage does not always coincide with the ground measurements of soil moisture, the existing ground data were used to calibrate a water balance for the top IO cm surface layer in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area average soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture, Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.

Owe, M.↗
NASA NTRS1999

Horizontal Wavenumber Spectra of Winds, Temperature and Trace Gases During the Pacific Exploratory Missions: Gravity Waves, Quasi-Two-Dimensional Turbulence, and Vortical Modes - 2

We examine the horizontal wavenumber spectra of horizontal velocity and potential temperature collected by aircraft above the Pacific Ocean to determine whether gravity waves, quasi-two-dimensional (Q-2-D) turbulence, or vortical modes dominate atmospheric fluctuations at scale sizes of 1-100 km and altitudes of 2-12 km. We conclude from the study of Doppler-shifting effects that Q-2-D turbulence and/or vortical modes are more prevalent than gravity waves over the ocean, except in the equatorial zone. The results are consistent with recent numerical simulations of Q-2-D turbulence, which show that the characteristic inverse cascade of energy is greatly facilitated by the presence of background rotation. Furthermore, a Stokes-parameter analysis reveals the general paucity of coherent wavelike motions, although specific cases of gravity-wave propagation are observed. Finally, a case study of a long flight segment displays a k(sup -3) horizontal velocity variance spectrum at scales longer than about 100 km. A Stokes-parameter analysis indicates that these large-scale fluctuations were likely due to vortical modes rather than inertio-gravity waves.

Cho, John Y. N.↗
NASA NTRS1999

Rocket-in-a-Duct Performance Analysis

An axisymmetric, 110 N class, rocket configured with a free expansion between the rocket nozzle and a surrounding duct was tested in an altitude simulation facility. The propellants were gaseous hydrogen and gaseous oxygen and the hardware consisted of a heat sink type copper rocket firing through copper ducts of various diameters and lengths. A secondary flow of nitrogen was introduced at the blind end of the duct to mix with the primary rocket mass flow in the duct. This flow was in the range of 0 to 10% of the primary massflow and its effect on nozzle performance was measured. The random measurement errors on thrust and massflow were within +/-1%. One dimensional equilibrium calculations were used to establish the possible theoretical performance of these rocket-in-a-duct nozzles. Although the scale of these tests was small, they simulated the relevant flow expansion physics at a modest experimental cost. Test results indicated that lower performance was obtained at higher free expansion area ratios and longer ducts, while, higher performance was obtained with the addition of secondary flow. There was a discernable peak in specific impulse efficiency at 4% secondary flow. The small scale of these tests resulted in low performance efficiencies, but prior numerical modeling of larger rocket-in-a-duct engines predicted performance that was comparable to that of optimized rocket nozzles. This remains to be proven in large-scale, rocket-in-a-duct tests.

Schneider, Steven J.↗
NASA NTRS1999

Global Aerosol Radiative Forcing Derived from Sea WiFS-Inferred Aerosol Optical Properties

Aerosol optical properties inferred from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) radiance measurements are used to compute the aerosol shortwave radiative forcing using a radiative transfer model. The aerosol optical thickness at the wavelength of 865-nm is taken from the SeaWIFS archive. It is found that the nominal optical thickness over oceans ranges from 0.1 to 0.2. Using a maritime aerosol model and the radiances measured at the various SeaWiFS channels, the Angstrom exponent is determined to be 0.2174, the single-scattering albedo to be 0.995, and the asymmetry factor to be 0.786. The radiative transfer model has eight bands in the visible and ultraviolet spectral regions and three bands in the near infrared. It includes the absorption due to aerosols, water vapor, carbon dioxide, and oxygen, and the scattering due to aerosols and gases (Rayleigh scattering). The radiative forcing is computed over global oceans for four months (January, April, July, and October, 1998) to represent four seasons. It is found that the aerosol radiative forcing is large and changes significantly with seasons near the continents with large-scale forest fires and desert dust. Averaged over oceans and the four months, the aerosol radiative forcing is approximately 7 W/sq m at the top of the atmosphere. This large radiative forcing is expected to have a significant cooling effect on the Earth's climate as implied from simulations of a number of general circulation models.

Chou, Ming-Dah↗