Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid size”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Renormalization group analysis of reduced magnetohydrodynamics with application to subgrid modeling

The technique for obtaining a subgrid model for Navier-Stokes turbulence, based on renormalization group analysis (RNG), is extended to the reduced magnetohydrodynamic (RMND) equations. It is shown that a RNG treatment of the Alfven turbulence supported by the RMHD equations leads to effective values of the viscosity and resistivity at large scales, k yields 0, dependent on the amplitude of turbulence. The effective viscosity and resistivity become independent of the molecular quantities when the RNG analysis is augmented by the Kolmogorov argument for energy cascade. A self-contained system of equations is derived for the range of scales, k = 0-K, where K = pi/Delta is the maximum wave number for a grid size Delta. Differential operators, whose coefficients depend upon the amplitudes of the large-scale quantities, represent in this system the resistive and viscous dissipation.

Longcope, D. W.↗
NASA NTRS1992

Applications of multigrid software in the atmospheric sciences

Elliptic partial differential equations from different areas in the atmospheric sciences are efficiently and easily solved utilizing the multigrid software package named MUDPACK. It is demonstrated that the multigrid method is more efficient than other commonly employed techniques, such as Gaussian elimination and fixed-grid relaxation. The efficiency relative to other techniques, both in terms of storage requirement and computational time, increases quickly with grid size.

Adams, J.↗
NASA NTRS1993

Numerical simulations of galactic wakes

We have performed a series of numerical simulations which examine the simple but realistic case of a galaxy with a smooth potential corresponding to a King profile, and with no interstellar medium. This corresponds to the case of an early-type cluster galaxy that has previously been stripped of its gas. The simulations provide a numerical basis for future work involving galaxies with interstellar media. We use a time-explicit finite difference code (Lufkin & Hawley 1992, ApJ submitted) to obtain solutions to the nonlinear equations for adiabatic flow in axisymmetry. The computational grid covers a region extending 500 kpc from a stationary galaxy radially and along the symmetry axis. We use a graded mesh, allowing for full resolution near the galaxy. Because the assumed potential has a finite depth, no artificial inner boundary condition is necessary; reflection symmetry at the axis is assumed. With a total grid size of 64 radial x 128 vertical, each run requires approximately 5 minutes of cpu time on the NCSA Cray-2.

Lufkin, Eric A.↗
NASA NTRS1993

Parallel solution of high-order numerical schemes for solving incompressible flows

A new parallel numerical scheme for solving incompressible steady-state flows is presented. The algorithm uses a finite-difference approach to solving the Navier-Stokes equations. The algorithms are scalable and expandable. They may be used with only two processors or with as many processors as are available. The code is general and expandable. Any size grid may be used. Four processors of the NASA LeRC Hypercluster were used to solve for steady-state flow in a driven square cavity. The Hypercluster was configured in a distributed-memory, hypercube-like architecture. By using a 50-by-50 finite-difference solution grid, an efficiency of 74 percent (a speedup of 2.96) was obtained.

Milner, Edward J.↗
NASA NTRS1993

EFEDA - European field experiment in a desertification-threatened area

During June 1991 more than 30 scientific teams worked in Castilla-La Mancha, Spain, studying the energy and water transfer processes between soil, vegetation, and the atmosphere in semiarid conditions within the coordinated European research project EFEDA (European Field Experiment in Desertification-threatened Areas). Measurements were made from the microscale (e.g., measurements on single plants) up to a scale compatible with the grid size of global models. For this purpose three sites were selected 70 km apart and heavily instrumented at a scale in the order of 30 sq km. Aircraft missions, satellite data, and movable equipment were deployed to provide a bridge to the larger scale. This paper gives a description of the experimental design along with some of the preliminary results of this successful experiment.

Bolle, H.-J.↗
NASA NTRS1991

Grid-refinement study of hypersonic laminar flow over a 2-D ramp

Computations were made for those test cases of Problem 3 which were designated as laminar flows, viz., test cases 3.1, 3.2, 3.4, and 3.5. These test cases corresponded to flows over a flat plate and a compression ramp at high Mach number and at high Reynolds number. The computations over the compression ramps indicate a substantial streamwise extent of separation. Based on previous experience with separated laminar flows at high Mach numbers which indicated a substantial effect with spatial grid refinement, a series of computations with different grid sizes were performed. Also, for the flat plate, comparisons of the results for two different algorithms were made.

Thomas, James L.↗
NASA NTRS1993

Numerical investigations in three-dimensional internal flows

In the present reporting period, the 3D version of the OVERFLOW code was used to solve the flow within the internal portion of the supersonic inlet. The internal portion of this inlet is bounded by an inflow plane containing the leading edge of the sidewalls, the sidewalls, the ramp and cowl surfaces and an outflow plane just downstream of the minimum geometric area of the inlet. Boundary layer bleed was used in the two-dimensional calculations discussed in the previous progress report and that same bleed was applied in the present study. For reference, this bleed corresponds to locations designated as R2 and R3 in the Mach 5 inlet model test. Using the GRIDGEN code, a three dimensional grid was generated that accounted for the viscous effects expected to occur on the sidewall, as well as those known to occur on the ramp and cowl surfaces. The internal flow grid size was 141 streamwise by 101 cross stream by 71 in the lateral direction between sidewalls. Since the flow entering the inlet was not symmetrical, the inlet was solved from sidewall to sidewall (without using a symmetry plane). In addition to the short sidewalls proposed in the Langley geometry database, a set of shorter sidewalls was also investigated in the present study and was shown to have beneficial effects with respect to the flow distortion exiting the supersonic inlet. In addition to these calculations, additional 3D solutions using the OVERFLOW code were obtained for the flow downstream of the throat of the supersonic inlet, including a terminal shock wave system produced by a backpressured subsonic diffuser.

Rose, William C.↗
NASA NTRS1993

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 323 simulation (Run B). Other than the grid size, the primary difference between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. The numerical results indicate that gravitational collapse can proceed in a variety of ways. In the Run A, collapse led to an elongated tube-like structure, while in the Run B, collapse led to a flatter, disklike structure.

Shebalin, John V.↗
NASA NTRS1994

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗
NASA NTRS1993

Parameterization of sub-grid scale convection

The following topics are discussed: an overview of the cumulus parameterization problem; interactions between explicit and implicit processes in mesoscale models; effects of model grid size on the cumulus parameterization problem; parameterizing convective effects on momentum fields in mesoscale models; differences between slantwise and vertical cumulus parameterization; experiments with different closure hypotheses; and coupling cumulus parameterizations to boundary layer, stable cloud, and radiation schemes.

Frank, William↗
NASA NTRS1993

Simulation of wave propagation in three-dimensional random media

Quantitative error analysis for simulation of wave propagation in three dimensional random media assuming narrow angular scattering are presented for the plane wave and spherical wave geometry. This includes the errors resulting from finite grid size, finite simulation dimensions, and the separation of the two-dimensional screens along the propagation direction. Simple error scalings are determined for power-law spectra of the random refractive index of the media. The effects of a finite inner scale are also considered. The spatial spectra of the intensity errors are calculated and compared to the spatial spectra of intensity. The numerical requirements for a simulation of given accuracy are determined for realizations of the field. The numerical requirements for accurate estimation of higher moments of the field are less stringent.

Coles, William A.↗
NASA NTRS1994

The delayed-detonation model of a type Ia supernovae. 1: The deflagration phase

The nature of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is investigated by using two-dimensional numerical hydrodynamics simulations. A new algorithm is used to treat the deflagration front. Assuming that it propagates locally at the laminar flame speed, the deflagration is insufficient to unbind the star. Expansion shuts of the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. The burning front does become wrinkled, but the wavelength of the instability is much larger than the computational grid size and is resolved; this is consistent with previous analysis. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. This paper deals with the deflagration phase, from the onset of burning, through expansion and quenching of the flame, to the first contraction.

Arnett, David↗
NASA NTRS1994

Application of a distributed network in computational fluid dynamic simulations

A general-purpose 3-D, incompressible Navier-Stokes algorithm is implemented on a network of concurrently operating workstations using parallel virtual machine (PVM) and compared with its performance on a CRAY Y-MP and on an Intel iPSC/860. The problem is relatively computationally intensive, and has a communication structure based primarily on nearest-neighbor communication, making it ideally suited to message passing. Such problems are frequently encountered in computational fluid dynamics (CDF), and their solution is increasingly in demand. The communication structure is explicitly coded in the implementation to fully exploit the regularity in message passing in order to produce a near-optimal solution. Results are presented for various grid sizes using up to eight processors.

Deshpande, Manish↗
NASA NTRS1994

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 32(exp 3) simulation (Run A) and a 64(exp 3) simulation (Run B). Other than the grid size, the primary differences between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. ('Compressible motion' is that part of the velocity which has zero curl, but non-zero divergence). The numerical results indicate that gravitational collapse can proceed in a variety of ways. In Run A, collapse led to an elongated tube-like structure, while in Run B, collapse led to a flatter, disk-like structure.

Shebalin, J. V.↗
NASA NTRS1995

The role and status of Euler solvers in impulsive rotor noise computations

Several recent applications (in the last five years) of Euler solvers in the computation of impulsive noise from rotor blades emphasize their emerging role in complementing other methods and experimental work. In the area of high-speed impulsive noise the use of Euler solvers as research tools has become fairly mature with very favorable comparisons with experimental data, especially in hover. The grid sizes and resulting computational times are reasonable when compared to those required for accurate surface aerodynamics alone. Furthermore, Euler solvers have provided a rich database with the resolution and accuracy needed for input to Kirchhoff and acoustic analogy methods for predicting the far-field noise. On the other hand, the application of Euler solvers to calculate blade-vortex interaction noise is still far from mature. The computational resources required for accurate calculations away from the blade are much larger than for high-speed impulsive noise. Current calculations help improve the basic understanding of the phenomena involved, but to date no comparisons with experiment have been made. Fortunately, the use of coupled Euler solver/Kirchhoff methods seems to offer promise for a robust and efficient technique for predicting both high-speed impulsive noise and blade-vortex interaction noise. Finally, a simple model problem of an isolated vortex interacting with an arbitrarily prescribed pitching airfoil demonstrates the feasibility of using Euler solvers to examine noise reduction techniques. The use of simple aerodynamic quasi-static theory and the computed lift time history as feedback to determine the required pitching motion appears sufficient to significantly dampen the unsteady loading and subsequent acoustics by an order of magnitude within a few blade passages.

Baeder, James D.↗
NASA NTRS1995

Ground-based grasslands data to support remote sensing and ecosystem modeling of terrestrial primary production

Estimating terrestrial net primary production (NPP) using remote-sensing tools and ecosystem models requires adequate ground-based measurements for calibration, parameterization, and validation. These data needs were strongly endorsed at a recent meeting of ecosystem modelers organized by the International Geosphere-Biosphere Program's (IGBP's) Data and Information System (DIS) and its Global Analysis, Interpretation, and Modelling (GAIM) Task Force. To meet these needs, a multinational, multiagency project is being coordinated by the IGBP DIS to compile existing NPP data from field sites and to regionalize NPP point estimates to various-sized grid cells. Progress at Oak Ridge National Laboratory (ORNL) on compiling NPP data for grasslands as part of the IGBP DIS data initiative is described. Site data and associated documentation from diverse field studies are being acquired for selected grasslands and are being reviewed for completeness, consistency, and adequacy of documentation, including a description of sampling methods. Data are being compiled in a database with spatial, temporal, and thematic characteristics relevant to remote sensing and global modeling. NPP data are available from the ORNL Distributed Active Archive Center (DAAC) for biogeochemical dynamics. The ORNL DAAC is part of the Earth Observing System Data and Information System, of the US National Aeronautics and Space Administration.

Olson, R. J.↗
NASA NTRS1996

A Structured Grid Based Solution-Adaptive Technique for Complex Separated Flows

The objective of this work was to enhance the predictive capability of widely used computational fluid dynamic (CFD) codes through the use of solution adaptive gridding. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. In order to study the accuracy and efficiency improvements due to the grid adaptation, it is necessary to quantify grid size and distribution requirements as well as computational times of non-adapted solutions. Flow fields about launch vehicles of practical interest often involve supersonic freestream conditions at angle of attack exhibiting large scale separate vortical flow, vortex-vortex and vortex-surface interactions, separated shear layers and multiple shocks of different intensity. In this work, a weight function and an associated mesh redistribution procedure is presented which detects and resolves these features without user intervention. Particular emphasis has been placed upon accurate resolution of expansion regions and boundary layers. Flow past a wedge at Mach=2.0 is used to illustrate the enhanced detection capabilities of this newly developed weight function.

Thornburg, Hugh↗
NASA NTRS1996

Aerodynamic Optimization of the High Speed Civil Transport

The aerodynamic optimization program used for wing and fuselage optimization of the High Speed Civil Transport was modified for utilization of multiple processors on parallel processor computers. The modified version uses multiple processors to simultaneously conduct three-dimensional flow solutions of different wing and fuselage geometries for calculations of the gradient functions and for directional searches to minimize an objective function. Demonstrations have shown the parallel program to be useful for coarse grid optimization, however memory problems for processors on the IBM SP2 were encountered when finer grid sizes were used.

Brawley, Stephen C.↗