Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid codes”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

SSME Turbopump Turbine Computations

A two-dimensional viscous code was developed to be used in the prediction of the flow in the SSME high-pressure turbopump blade passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme to solve the two-dimensional, thin-layer Navier-Stokes equations. The Baldwin-Lomax eddy-viscosity model is used for these turbulent flow calculations. A viable method was developed to use the relative exit conditions from an upstream blade row as the inlet conditions to the next blade row. The blade loading diagrams are compared with the meridional values obtained from an in-house quasithree-dimensional inviscid code. Periodic boundary conditions are imposed on a body-fitted C-grid computed by using the GRAPE GRids about Airfoils using Poisson's Equation (GRAPE) code. Total pressure, total temperature, and flow angle are specified at the inlet. The upstream-running Riemann invariant is extrapolated from the interior. Static pressure is specified at the exit such that mass flow is conserved from blade row to blade row, and the conservative variables are extrapolated from the interior. For viscous flows the noslip condition is imposed at the wall. The normal momentum equation gives the pressure at the wall. The density at the wall is obtained from the wall total temperature.

Jorgenson, P. G. E.↗
NASA NTRS1999

A CFD Assessment of Several High-Lift Reference H Configuration Using Structured Grids

The objective of this study is to calibrate a Navier-Stokes code for a high-lift Reference H configuration using structured grids. The outline of this presentation will first include a brief description of the grids used and the flow solver. Next the results will be presented in terms of convergence and resources used on the C-90. Predicted force and moment and surface pressure results are compared to experiment and off- and on-surface flow viz. is discussed.

Lessard, Wendy B.↗
DOE OSTI2020

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Characterizing Artificial Viscosity Parameters with Approximate Symmetries

We are often faced with trying to capture the physics of compressible shocks, which are governed by the Euler equations. However, Euler shocks are formally discontinuous at the shock front (translating to a step-function behavior of rel evant flow variables). This poses a practical problem for codes with finite-sized grid elements. As a result, one must make a concession in simulating the be havior of shocks within a discrete framework. In particular, we must blur, or ‘regularize’ Euler shocks so that they may be captured on a finite grid.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

Numerical solutions of the Navier-Stokes equations for transonic afterbody flows

The time dependent Navier-Stokes equations in mass averaged variables are solved for transonic flow over axisymmetric boattail plume simulator configurations. Numerical solution of these equations is accomplished with the unsplit explict finite difference algorithm of MacCormack. A grid subcycling procedure and computer code vectorization are used to improve computational efficiency. The two layer algebraic turbulence models of Cebeci-Smith and Baldwin-Lomax are employed for investigating turbulence closure. Two relaxation models based on these baseline models are also considered. Results in the form of surface pressure distribution for three different circular arc boattails at two free stream Mach numbers are compared with experimental data. The pressures in the recirculating flow region for all separated cases are poorly predicted with the baseline turbulence models. Significant improvements in the predictions are usually obtained by using the relaxation models.

Swanson, R. C., Jr.↗
NASA NTRS1982

Characterization of silicon-gate CMOS/SOS integrated circuits processed with ion implantation

The procedure used to generate MEBES masks and produce test wafers from the 10X Mann 1600 Pattern Generator Tape using existing CAD utility programs and the MEBES machine in the RCA Solid State Technology Center are described. The test vehicle used is the MSFC-designed SC102 Solar House Timing Circuit. When transforming the Mann 1600 tapes into MEBES tapes, extreme care is required in order to obtain accurate minimum linewidths when working with two different coding systems because the minimum grid sizes may be different for the two systems. The minimum grid sizes are 0.025 mil for MSFC Mann 1600 and 0.02 mil for MEBES. Some snapping to the next grid is therefore inevitable, and the results of this snapping effect are significant when submicron lines are present. However, no problem was noticed in the SC102 circuit because its minimum linewidth is 0.3 mil (7.6 microns). MEBES masks were fabricated and wafers were processed using the silicon-gate CMOS/SOS and aluminum-gate COS/MOS processing.

Woo, D. S.↗
NASA NTRS1984

Calculations of turbulent mass transport in a bluff-body diffusion-flame combustor

Experimental measurements of turbulent mass transit in a bluff-body diffusion-flame combustor (using CO2 instead of fuel) are analyzed to evaluate the accuracy of physical models used in flow computations. The data of Lightman and Magill (1981), Lightman et al. (1983), and Roquemore et al. (1983) are used to calculate apparent turbulent Schmidt numbers (TSN) for a series of flow conditions by a modified TEACH viscous-flow code. The modeling principles, calculation grid, and boundary conditions are discussed, and the results are presented in graphs comparing calculated and measured values. It is shown that models such as the two-equation (K - epsilon) model which use a single value of the TSN are inappropriate for conditions typical of gas-turbine combustors: local TSN variations due to turbulence behavior must be taken into account in improved models.

Sturgess, G. J.↗
NASA NTRS1989

Hypersonic CFD applications at NASA Langley using CFL3D and CFL3DE

The CFL3D/CFL3DE CFD codes and the industrial use status of the codes are outlined. Comparison of grid density, pressure, heat transfer, and aerodynamic coefficience are presented. Future plans related to the National Aerospace Plane Program are briefly outlined.

Richardson, Pamela F.↗
NASA NTRS1991

Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation

Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in a high gradient region or the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational costs. A detailed description is given of the enrichment and coarsening procedures and comparisons with alternative results and experimental data are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.

Rausch, Russ D.↗
NASA NTRS1991

Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation

Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with alternative results and experimental data to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.

Rausch, Russ D.↗
NASA NTRS1991

The design/analysis of flows through turbomachinery - A viscous/inviscid approach

A new design/analysis system for the flows through turbomachinery is currently being developed for studying turbomachinery problems with an axisymmetric viscous/inviscid 'average-passage' throughflow code. The advantage of this approach, compared to streamline curvature codes, is that the solutions obtained simulate some of the unsteadiness, compressibility and viscous effects of a multistage turbomachine. The design/analysis system consists of three elemental parts, the axisymmetric block grid generator, the blade surface element code, and the axisymmetric flow code. Each element of the system will be discussed and the flow solutions for three axisymmetric geometries will be shown compared to experimental data where available. The computations are shown to be in very good agreement with test data for SR7 spinner body and transonic boattail geometry obtained in the wind tunnels at NASA Lewis Research Center. The VIADAC Rotor 67 Fan results were compared to PARC2D calculated results and shown to be in very good agreement.

Miller, D. P.↗
NASA NTRS1993

Spatial adaptation procedures on tetrahedral meshes for unsteady aerodynamic flow calculations

Spatial adaptation procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaptation procedures were developed and implemented within a three-dimensional, unstructured-grid, upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. A detailed description of the enrichment and coarsening procedures are presented and comparisons with experimental data for an ONERA M6 wing and an exact solution for a shock-tube problem are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady results, obtained using spatial adaptation procedures, are shown to be of high spatial accuracy, primarily in that discontinuities such as shock waves are captured very sharply.

Rausch, Russ D.↗
NASA NTRS1993

Spatial adaptation procedures on tetrahedral meshes for unsteady aerodynamic flow calculations

Spatial adaptation procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaptation procedures were developed and implemented within a three-dimensional, unstructured-grid, upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with experimental data for an ONERA M6 wing and an exact solution for a shock-tube problem to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady results, obtained using spatial adaptation procedures, are shown to be of high spatial accuracy, primarily in that discontinuities such as shock waves are captured very sharply.

Rausch, Russ D.↗
NASA NTRS1995

Ion propulsion

An ion engine is a plasma thruster which produces thrust by extracting ions from the plasma and accelerating them to high velocity with an electrostatic field. The ions are then neutralized and leave the engine as high velocity neutral particles. The advantages of ion engines are high specific impulse and efficiency and their ability to operate over a wide range of input powers. In comparison with other electric thrusters, the ion engine has higher efficiency and specific impulse than thermal electric devices such as the arcjet, microwave, radiofrequency and laser heated thrusters and can operate at much lower current levels than the MPD thruster. However, the thrust level for an ion engine may be lower than a thermal electric thruster of the same operating power, consistent with its higher specific impulse, and therefore ion engines are best suited for missions which can tolerate longer duration propulsive phases. The critical issue for the ion engine is lifetime, since the prospective missions may require operation for several thousands of hours. The critical components of the ion engine, with respect to engine lifetime, are the screen and accelerating grid structures. Typically, these are large metal screens that must support a large voltage difference and maintain a small gap between them. Metallic whisker growth, distortion and vibration can lead to arcing, and over a long period of time ion sputtering will erode the grid structures and change their geometry. In order to study the effects of long time operation of the grid structure, we are developing computer codes based on the Particle-In-Cell (PIC) technique and Laser Induced Fluorescence (LIF) diagnostic techniques to study the physical processes which control the performance and lifetime of the grid structures.

Meserole, J. S.↗
NASA NTRS1989

Ion propulsion

An ion engine is a plasma thruster which produces thrust by extracting ions from the plasma and accelerating them to high velocity with an electrostatic field. The ions are then neutralized and leave the engine as high velocity neutral particles. The advantages of ion engines are high specific impulse and efficiency and their ability to operate over a wide range of input powers. In comparison with other electric thrusters, the ion engine has higher efficiency and specific impulse than thermal electric devices such as the arcjet, microwave, radio frequency, and laser heated thrusters and can operate at much lower current levels than the MPD thruster. However, the thrust level for an ion engine may be lower than a thermal electric thruster of the same operating power, consistent with its higher specific impulse, and therefore ion engines are best suited for missions which can tolerate longer duration propulsive phases. The critical issue for the ion engine is lifetime, since the prospective missions may require operation for several thousands of hours. The critical components of the ion engine, with respect to engine lifetime, are the screen and accelerating grid structures. Typically, these are large metal screens that must support a large voltage difference and maintain a small gap between them. Metallic whisker growth, distortion, and vibration can lead to arcing, and over a long period of time ion sputtering will erode the grid structures and change their geometry. In order to study the effects of long time operation of the grid structure, we are developing computer codes based on the particle-in-cell (PIC) technique and laser induced fluorescence (LIF) diagnostic techniques to study the physical processes which control the performance and lifetime of the grid structures.

Meserole, J. S.↗
NASA NTRS1995

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗
NASA NTRS2005

Navier-Stokes Simulation of Airconditioning Facility of a Large Modem Computer Room

NASA recently assembled one of the world's fastest operational supercomputers to meet the agency's new high performance computing needs. This large-scale system, named Columbia, consists of 20 interconnected SGI Altix 512-processor systems, for a total of 10,240 Intel Itanium-2 processors. High-fidelity CFD simulations were performed for the NASA Advanced Supercomputing (NAS) computer room at Ames Research Center. The purpose of the simulations was to assess the adequacy of the existing air handling and conditioning system and make recommendations for changes in the design of the system if needed. The simulations were performed with NASA's OVERFLOW-2 CFD code which utilizes overset structured grids. A new set of boundary conditions were developed and added to the flow solver for modeling the roomls air-conditioning and proper cooling of the equipment. Boundary condition parameters for the flow solver are based on cooler CFM (flow rate) ratings and some reasonable assumptions of flow and heat transfer data for the floor and central processing units (CPU) . The geometry modeling from blue prints and grid generation were handled by the NASA Ames software package Chimera Grid Tools (CGT). This geometric model was developed as a CGT-scripted template, which can be easily modified to accommodate any changes in shape and size of the room, locations and dimensions of the CPU racks, disk racks, coolers, power distribution units, and mass-storage system. The compute nodes are grouped in pairs of racks with an aisle in the middle. High-speed connection cables connect the racks with overhead cable trays. The cool air from the cooling units is pumped into the computer room from a sub-floor through perforated floor tiles. The CPU cooling fans draw cool air from the floor tiles, which run along the outside length of each rack, and eject warm air into the center isle between the racks. This warm air is eventually drawn into the cooling units located near the walls of the room. One major concern is that the hot air ejected to the middle isle might recirculate back into the cool rack side and cause thermal short-cycling. The simulations analyzed and addressed the following important elements of the computer room: 1) High-temperature build-up in certain regions of the room; 2) Areas of low air circulation in the room; 3) Potential short-cycling of the computer rack cooling system; 4) Effectiveness of the perforated cooling floor tiles; 5) Effect of changes in various aspects of the cooling units. Detailed flow visualization is performed to show temperature distribution, air-flow streamlines and velocities in the computer room.

Source record↗
NASA NTRS2007

A Perspective on Computational Aerothermodynamics at NASA

The evolving role of computational aerothermodynamics (CA) within NASA over the past 20 years is reviewed. The presentation highlights contributions to understanding the Space Shuttle pitching moment anomaly observed in the first shuttle flight, prediction of a static instability for Mars Pathfinder, and the use of CA for damage assessment in post-Columbia mission support. In the view forward, several current challenges in computational fluid dynamics and aerothermodynamics for hypersonic vehicle applications are discussed. Example simulations are presented to illustrate capabilities and limitations. Opportunities to advance the state-of-art in algorithms, grid generation and adaptation, and code validation are identified.

Gnoffo, Peter A.↗