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

Time-dependent Navier-Stokes computations of separated flows over airfoils

The application of both a central-difference and an upwind implicit approximate factorization Navier-Stokes algorithm to highly separated flow is described. Using the thin-layer approximation, both algorithms are employed to solve the low Reynolds number laminar flow around a circular cylinder with periodic shedding. Results agree with experimental Strouhal numbers to within about 10 percent. A grid refinement study is conducted using the upwind code for both the thin-layer and complete Navier-Stokes approximations. The circular cylinder solution is shown to be highly dependent upon both grid density outward from the body and grid extent. Finally, the thin-layer central-difference code is used to predict the time-accurate unsteady flow about two stalled airfoils. The NACA 0012 airfoil turbulent solutions are periodic with Strouhal numbers of about 0.1; a strong vortex is periodically shed from the leading edge. The LS(1)-0417 airfoil turbulent solution shows a high frequency oscillatory trailing edge separation. The 20 degree NACA 0012 airfoil solution is massively separated, but shows no evidence of periodicity.

Rumsey, C. L.↗
NASA NTRS2002

Flowfield Comparisons from Three Navier-Stokes Solvers for an Axisymmetric Separate Flow Jet

To meet new noise reduction goals, many concepts to enhance mixing in the exhaust jets of turbofan engines are being studied. Accurate steady state flowfield predictions from state-of-the-art computational fluid dynamics (CFD) solvers are needed as input to the latest noise prediction codes. The main intent of this paper was to ascertain that similar Navier-Stokes solvers run at different sites would yield comparable results for an axisymmetric two-stream nozzle case. Predictions from the WIND and the NPARC codes are compared to previously reported experimental data and results from the CRAFT Navier-Stokes solver. Similar k-epsilon turbulence models were employed in each solver, and identical computational grids were used. Agreement between experimental data and predictions from each code was generally good for mean values. All three codes underpredict the maximum value of turbulent kinetic energy. The predicted locations of the maximum turbulent kinetic energy were farther downstream than seen in the data. A grid study was conducted using the WIND code, and comments about convergence criteria and grid requirements for CFD solutions to be used as input for noise prediction computations are given. Additionally, noise predictions from the MGBK code, using the CFD results from the CRAFT code, NPARC, and WIND as input are compared to data.

Koch, L. Danielle↗
NASA NTRS1988

Accurate modeling of parallel scientific computations

Scientific codes are usually parallelized by partitioning a grid among processors. To achieve top performance it is necessary to partition the grid so as to balance workload and minimize communication/synchronization costs. This problem is particularly acute when the grid is irregular, changes over the course of the computation, and is not known until load time. Critical mapping and remapping decisions rest on the ability to accurately predict performance, given a description of a grid and its partition. This paper discusses one approach to this problem, and illustrates its use on a one-dimensional fluids code. The models constructed are shown to be accurate, and are used to find optimal remapping schedules.

Nicol, David M.↗
NASA NTRS1989

The simulation of hypersonic flow using a zonal approach

The development of a Navier-Stokes code for the simulation of hypersonic flow will be described. The code uses a zonal philosophy which helps alleviate the grid generation problem for complex configurations and allows an efficient solution procedure. The code has been used to simulate equilibrium flow past a simple generic hypersonic configuration. Validation of the code on a blunt cone has also been conducted.

Flores, Jolen↗
NASA NTRS2013

The Impact of Back-Sputtered Carbon on the Accelerator Grid Wear Rates of the NEXT and NSTAR Ion Thrusters

A study was conducted to quantify the impact of back-sputtered carbon on the downstream accelerator grid erosion rates of the NEXT (NASA's Evolutionary Xenon Thruster) Long Duration Test (LDT1). A similar analysis that was conducted for the NSTAR (NASA's Solar Electric Propulsion Technology Applications Readiness Program) Life Demonstration Test (LDT2) was used as a foundation for the analysis developed herein. A new carbon surface coverage model was developed that accounted for multiple carbon adlayers before complete surface coverage is achieved. The resulting model requires knowledge of more model inputs, so they were conservatively estimated using the results of past thin film sputtering studies and particle reflection predictions. In addition, accelerator current densities across the grid were rigorously determined using an ion optics code to determine accelerator current distributions and an algorithm to determine beam current densities along a grid using downstream measurements. The improved analysis was applied to the NSTAR test results for evaluation. The improved analysis demonstrated that the impact of back-sputtered carbon on pit and groove wear rate for the NSTAR LDT2 was negligible throughout most of eroded grid radius. The improved analysis also predicted the accelerator current density for transition from net erosion to net deposition considerably more accurately than the original analysis. The improved analysis was used to estimate the impact of back-sputtered carbon on the accelerator grid pit and groove wear rate of the NEXT Long Duration Test (LDT1). Unlike the NSTAR analysis, the NEXT analysis was more challenging because the thruster was operated for extended durations at various operating conditions and was unavailable for measurements because the test is ongoing. As a result, the NEXT LDT1 estimates presented herein are considered preliminary until the results of future posttest analyses are incorporated. The worst-case impact of carbon back-sputtering was determined to be the full power operating condition, but the maximum impact of back-sputtered carbon was only a four percent reduction in wear rate. As a result, back-sputtered carbon is estimated to have an insignificant impact on the first failure mode of the NEXT LDT at all operating conditions.

Soulas, George C.↗
NASA NTRS2013

The Impact of Back-Sputtered Carbon on the Accelerator Grid Wear Rates of the NEXT and NSTAR Ion Thrusters

A study was conducted to quantify the impact of back-sputtered carbon on the downstream accelerator grid erosion rates of the NASA's Evolutionary Xenon Thruster (NEXT) Long Duration Test (LDT1). A similar analysis that was conducted for the NASA's Solar Electric Propulsion Technology Applications Readiness Program (NSTAR) Life Demonstration Test (LDT2) was used as a foundation for the analysis developed herein. A new carbon surface coverage model was developed that accounted for multiple carbon adlayers before complete surface coverage is achieved. The resulting model requires knowledge of more model inputs, so they were conservatively estimated using the results of past thin film sputtering studies and particle reflection predictions. In addition, accelerator current densities across the grid were rigorously determined using an ion optics code to determine accelerator current distributions and an algorithm to determine beam current densities along a grid using downstream measurements. The improved analysis was applied to the NSTAR test results for evaluation. The improved analysis demonstrated that the impact of back-sputtered carbon on pit and groove wear rate for the NSTAR LDT2 was negligible throughout most of eroded grid radius. The improved analysis also predicted the accelerator current density for transition from net erosion to net deposition considerably more accurately than the original analysis. The improved analysis was used to estimate the impact of back-sputtered carbon on the accelerator grid pit and groove wear rate of the NEXT Long Duration Test (LDT1). Unlike the NSTAR analysis, the NEXT analysis was more challenging because the thruster was operated for extended durations at various operating conditions and was unavailable for measurements because the test is ongoing. As a result, the NEXT LDT1 estimates presented herein are considered preliminary until the results of future post-test analyses are incorporated. The worst-case impact of carbon back-sputtering was determined to be the full power operating condition, but the maximum impact of back-sputtered carbon was only a 4 percent reduction in wear rate. As a result, back-sputtered carbon is estimated to have an insignificant impact on the first failure mode of the NEXT LDT1 at all operating conditions.

Soulas, George C.↗
NASA NTRS2001

Efficient Cache use for Stencil Operations on Structured Discretization Grids

We derive tight bounds on the cache misses for evaluation of explicit stencil operators on structured grids. Our lower bound is based on the isoperimetrical property of the discrete octahedron. Our upper bound is based on a good surface to volume ratio of a parallelepiped spanned by a reduced basis of the interference lattice of a grid. Measurements show that our algorithm typically reduces the number of cache misses by a factor of three, relative to a compiler optimized code. We show that stencil calculations on grids whose interference lattice have a short vector feature abnormally high numbers of cache misses. We call such grids unfavorable and suggest to avoid these in computations by appropriate padding. By direct measurements on a MIPS R10000 processor we show a good correlation between abnormally high numbers of cache misses and unfavorable three-dimensional grids.

Frumkin, Michael↗
NASA NTRS2003

Numerical Simulation of Ion Thruster Optics

A three-dimensional simulation code (ffx) designed to analyze ion thruster optics is described. It is an extension of an earlier code and includes special features like the ability to model a wide range of grid geometries, cusp details, and mis-aligned aperture pairs to name a few. However, the principle reason for advancing the code was in the study of ion optics erosion. Ground based testing of ion thruster optics, essential to the understanding of the processes of grid erosion, can be time consuming and costly. Simulation codes that can accurately predict grid lifetimes and the physical mechanisms of grid erosion can be of great utility in the development of future ion thruster optics designed for more ambitious applications. Results of simulations are presented that describe wear profiles for several standard and nonstandard aperture geometries, such as those grid sets with square- or slotted-hole layout patterns. The goal of this paper will be to introduce the methods employed in the ffx code and to briefly demonstrate their use.

Rawlin, Vincent K.↗
NASA NTRS1991

Multitasking the INS3D-LU code on the Cray Y-MP

This paper presents the results of multitasking the INS3D-LU code on eight processors. The code is a full Navier-Stokes solver for incompressible fluid in three dimensional generalized coordinates using a lower-upper symmetric-Gauss-Seidel implicit scheme. This code has been fully vectorized on oblique planes of sweep and parallelized using autotasking with some directives and minor modifications. The timing results for five grid sizes are presented and analyzed. The code has achieved a processing rate of over one Gflops.

Fatoohi, Rod↗
NASA NTRS1988

Body conforming grids for general unsteady airfoil motion

A method for computing a new body conforming grid at each time step of an unsteady flow calculation is presented. The method is based on a spring system analogy. It simply updates the grid from the previous time step, hence it requires that an initial grid be created using another method. The method is designed to be independent of the flow equations used, so it can be easily transferred between aerodynamic codes with the same type of initial grid (i.e., O-type). The initial development was done using an O-type grid. It was then extended to C-type grids which required only minor changes to the routine. Results are presented for both O-type and C-type grids. The quality of the updated grids were examined by observing their movement and the flow results they produced. The effect due to the shape of the prescribed grid lines for the O-type grids was examined for a NACA 0012 airfoil undergoing harmonic pitch. The time history of the pressure coefficients of this case was compared with the interpolation scheme and experimental data.

Levine, M.↗
NASA NTRS1981

Black box multigrid

The black box multigrid (BOXMG) code, which only needs specification of the matrix problem for application in the multigrid method was investigated. It is contended that a major problem with the multigrid method is that each new grid configuration requires a major programming effort to develop a code that specifically handles that grid configuration. The SOR and ICCG methods only specify the matrix problem, no matter what the grid configuration. It is concluded that the BOXMG does everything else necessary to set up the auxiliary coarser problems to achieve a multigrid solution.

Dendy, J. E., Jr.↗
NASA NTRS1991

PARCEQ2D heat transfer grid sensitivity analysis

The material presented in this paper is an extension of two-dimensional Aeroassist Flight Experiment (AFE) results shown previously. This study has focused on the heating rate calculations to the AFE obtained from an equilibrium real gas code, with attention placed on the sensitivity of grid dependence and wall temperature. Heat transfer results calculated by the PARCEQ2D code compare well with those computed by other researchers. Temperature convergence in the case of kinetic transport has been accomplished by increasing the wall temperature gradually from 300 K to the wall temperature of 1700 K.

Saladino, Anthony J.↗
NASA NTRS1998

User Manual for Beta Version of TURBO-GRD: A Software System for Interactive Two-Dimensional Boundary/ Field Grid Generation, Modification, and Refinement

TURBO-GRD is a software system for interactive two-dimensional boundary/field grid generation. modification, and refinement. Its features allow users to explicitly control grid quality locally and globally. The grid control can be achieved interactively by using control points that the user picks and moves on the workstation monitor or by direct stretching and refining. The techniques used in the code are the control point form of algebraic grid generation, a damped cubic spline for edge meshing and parametric mapping between physical and computational domains. It also performs elliptic grid smoothing and free-form boundary control for boundary geometry manipulation. Internal block boundaries are constructed and shaped by using Bezier curve. Because TURBO-GRD is a highly interactive code, users can read in an initial solution, display its solution contour in the background of the grid and control net, and exercise grid modification using the solution contour as a guide. This process can be called an interactive solution-adaptive grid generation.

Choo, Yung K.↗
NASA NTRS1997

Creating a Simple Single Computational Approach to Modeling Rarefied and Continuum Flow About Aerospace Vehicles

We proposed to create a single computational code incorporating methods that can model both rarefied and continuum flow to enable the efficient simulation of flow about space craft and high altitude hypersonic aerospace vehicles. The code was to use a single grid structure that permits a smooth transition between the continuum and rarefied portions of the flow. Developing an appropriate computational boundary between the two regions represented a major challenge. The primary approach chosen involves coupling a four-speed Lattice Boltzmann model for the continuum flow with the DSMC method in the rarefied regime. We also explored the possibility of using a standard finite difference Navier Stokes solver for the continuum flow. With the resulting code we will ultimately investigate three-dimensional plume impingement effects, a subject of critical importance to NASA and related to the work of Drs. Forrest Lumpkin, Steve Fitzgerald and Jay Le Beau at Johnson Space Center. Below is a brief background on the project and a summary of the results as of the end of the grant.

Goldstein, David B.↗
NASA NTRS2001

Coarse Grid Modeling of Turbine Film Cooling Flows Using Volumetric Source Terms

The recent trend in numerical modeling of turbine film cooling flows has been toward higher fidelity grids and more complex geometries. This trend has been enabled by the rapid increase in computing power available to researchers. However, the turbine design community requires fast turnaround time in its design computations, rendering these comprehensive simulations ineffective in the design cycle. The present study describes a methodology for implementing a volumetric source term distribution in a coarse grid calculation that can model the small-scale and three-dimensional effects present in turbine film cooling flows. This model could be implemented in turbine design codes or in multistage turbomachinery codes such as APNASA, where the computational grid size may be larger than the film hole size. Detailed computations of a single row of 35 deg round holes on a flat plate have been obtained for blowing ratios of 0.5, 0.8, and 1.0, and density ratios of 1.0 and 2.0 using a multiblock grid system to resolve the flows on both sides of the plate as well as inside the hole itself. These detailed flow fields were spatially averaged to generate a field of volumetric source terms for each conservative flow variable. Solutions were also obtained using three coarse grids having streamwise and spanwise grid spacings of 3d, 1d, and d/3. These coarse grid solutions used the integrated hole exit mass, momentum, energy, and turbulence quantities from the detailed solutions as volumetric source terms. It is shown that a uniform source term addition over a distance from the wall on the order of the hole diameter is able to predict adiabatic film effectiveness better than a near-wall source term model, while strictly enforcing correct values of integrated boundary layer quantities.

Heidmann, James D.↗
NASA NTRS2002

Incompressible Navier-Stokes Solvers in Primative Variables and their Applications to Steady and Unsteady Flow Simulations

This paper reviews recent progress made in incompressible Navier-Stokes simulation procedures and their application to problems of engineering interest. Discussions are focused on the methods designed for complex geometry applications in three dimensions, and thus are limited to primitive variable formulation. A summary of efforts in flow solver development is given followed by numerical studies of a few example problems of current interest. Both steady and unsteady solution algorithms and their salient features are discussed. Solvers discussed here are based on a structured-grid approach using either a finite -difference or a finite-volume frame work. As a grand-challenge application of these solvers, an unsteady turbopump flow simulation procedure has been developed which utilizes high performance computing platforms. In the paper, the progress toward the complete simulation capability of the turbo-pump for a liquid rocket engine is reported. The Space Shuttle Main Engine (SSME) turbo-pump is used as a test case for evaluation of two parallel computing algorithms that have been implemented in the INS3D code. The relative motion of the grid systems for the rotorstator interaction was obtained using overact grid techniques. Unsteady computations for the SSME turbo-pump, which contains 114 zones with 34.5 million grid points, are carried out on SCSI Origin 3000 systems at NASA Ames Research Center. The same procedure has been extended to the development of NASA-DeBakey Ventricular Assist Device (VAD) that is based on an axial blood pump. Computational, and clinical analysis of this device are presented.

Kiris, Cetin C.↗
NASA NTRS1997

Advances in Computational Capabilities for Hypersonic Flows

The paper reviews the growth and advances in computational capabilities for hypersonic applications over the period from the mid-1980's to the present day. The current status of the code development issues such as surface and field grid generation, algorithms, physical and chemical modeling, and validation is provided. A brief description of some of the major codes being used at NASA Langley Research Center for hypersonic continuum and rarefied flows is provided, along with their capabilities and deficiencies. A number of application examples are presented, and future areas of research to enhance accuracy, reliability, efficiency, and robustness of computational codes are discussed.

Kumar, Ajay↗
NASA NTRS1992

IGGy: An interactive environment for surface grid generation

A graphically interactive derivative of the EAGLE boundary code is presented. This code allows the user to interactively build and execute commands and immediately see the results. Strong ties with a batch oriented script language are maintained. A generalized treatment of grid definition parameters allows a more generic definition of the grid generation process and allows the generation of command scripts which can be applied to topologically similar configurations. The use of the graphical user interface is outlined and example applications are presented.

Prewitt, Nathan C.↗