Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structured Grids”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Lunar crater chains of non-impact origin

Apollo 15, 16, and 17 photographs were scanned for lunar crater chains consisting of three or more aligned craters with similar states of degradation, and the origin of these chains is considered. The mode of origin for mare crater chains appears to be associated with either collapsed lava tubes or cinder cones, while the origin of highland crater chains generally seems to involve highland volcanoes. Crater chains from mare regions range from 20 to 40 km long and appear to have no preferred structural control of their orientation. Highland crater chains range from 1 to 113 km long and may be structurally controlled by the lunar grid system.

Eppler, D.↗
NASA NTRS1981

Transonic flow computations in cascades using finite element method

Analysis of transonic flow through a cascade of airfoils is investigated using the finite element method. Development of a computational grid suitable for complex flow structures and different types of boundary conditions is presented. An efficient pseudo-time integration scheme is developed for the solution of equations. Modeling of the shock and the convergence characteristics of the developed scheme are discussed. Numerical results include a 45 deg staggered cascade of NACA 0012 airfoils with inlet flow Mach number of 0.8 and angles of attack 1, 0, and 1 deg.

Akay, H. U.↗
NASA NTRS1990

Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

Improved algorithms for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration shceme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. The paper presents a description of the Euler solvers along with results and comparisons which assess the capability.

Batina, John T.↗
NASA NTRS1990

Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes

Improved algorithms for the solution of the three-dimensional time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multi-stage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons which assess the capability.

Batina, John T.↗
NASA NTRS1990

Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

Improved algorithm for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗
NASA NTRS1990

Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes

Improved algorithms for the solution of the 3-D time dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time integration scheme using a multistage Runge-Kutta procedure or an implicit time integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗
NASA NTRS1991

Accuracy of an unstructured-grid upwind-Euler algorithm for the ONERA M6 wing

Improved algorithms for the solution of the three-dimensional, time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured-grid flow solvers. The spatial discretization involves a flux-split approach that is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multistage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure, which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the Office National d'Etudes et de Recherches Aerospatiales M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons that assess the capability.

Batina, John T.↗
NASA NTRS1991

Euler/experiment correlations of sonic boom pressure signatures

The ability of inviscid computational fluid dynamics (CFD) codes to compute sonic boom pressure signatures is examined using three different codes that solve the Euler equations of fluid flow on structured hexahedral and unstructured tetrahedral grids. The results of these Euler codes were evaluated by comparing the computed pressure signatures with near-field experimental data. The computational pressure signatures were determined at distances of one body length or less below the configuration in the plane of symmetry and extrapolated to experimental distances. The extrapolated CFD pressure signatures gave acceptable correlations with experimental data, provided that fine grids were used near the surface and downstream of the configuration.

Cliff, Susan E.↗
NASA NTRS1995

A 3DGRAPE/AL: The Ames/Langley technology upgrade

This paper describes a new three-dimensional structured multiple-block volume grid generator called 3DGRAPE/AL. It is a significantly improved version of the previously-released and widely-distributed program 3DGRAPE, with many of the improvements taken from the grid-generator program 3DMAGGS. It generates volume grids by iteratively solving the Poisson Equations in three-dimensions. The right-hand-side terms are designed so that user-specified grid cell heights and user-specified grid cell skewness near boundary surfaces result automatically, with little user intervention. Versatility was a high priority in this code's development, and as a result it can generate grids in almost any three-dimensional physical domain. Improvements include added kinds of forcing functions, improved control of cell skewness, improved initial conditions, convergence acceleration, the ability to take as input the output from GRIDGEN, and a simple but powerful graphical user interface (GUI).

Sorenson, Reese L.↗
NASA NTRS2000

Interference Lattice-based Loop Nest Tilings for Stencil Computations

A common method for improving performance of stencil operations on structured multi-dimensional discretization grids is loop tiling. Tile shapes and sizes are usually determined heuristically, based on the size of the primary data cache. We provide a lower bound on the numbers of cache misses that must be incurred by any tiling, and a close achievable bound using a particular tiling based on the grid interference lattice. The latter tiling is used to derive highly efficient loop orderings. The total number of cache misses of a code is the sum of (necessary) cold misses and misses caused by elements being dropped from the cache between successive loads (replacement misses). Maximizing temporal locality is equivalent to minimizing replacement misses. Temporal locality of loop nests implementing stencil operations is optimized by tilings that avoid data conflicts. We divide the loop nest iteration space into conflict-free tiles, derived from the cache miss equation. The tiling involves the definition of the grid interference lattice an equivalence class of grid points whose images in main memory map to the same location in the cache-and the construction of a special basis for the lattice. Conflicts only occur on the boundaries of the tiles, unless the tiles are too thin. We show that the surface area of the tiles is bounded for grids of any dimensionality, and for caches of any associativity, provided the eccentricity of the fundamental parallelepiped (the tile spanned by the basis) of the lattice is bounded. Eccentricity is determined by two factors, aspect ratio and skewness. The aspect ratio of the parallelepiped can be bounded by appropriate array padding. The skewness can be bounded by the choice of a proper basis. Combining these two strategies ensures that pathologically thin tiles are avoided. They do not, however, minimize replacement misses per se. The reason is that tile visitation order influences the number of data conflicts on the tile boundaries. If two adjacent tiles are visited successively, there will be no replacement misses on the shared boundary. The iteration space may be covered with pencils larger than the size of the cache while avoiding data conflicts if the pencils are traversed by a scanning-face method. Replacement misses are incurred only on the boundaries of the pencils, and the number of misses is minimized by maximizing the volume of the scanning face, not the volume of the tile. We present an algorithm for constructing the most efficient scanning face for a given grid and stencil operator. In two dimensions it is based on a continued fraction algorithm. In three dimensions it follows Voronoi's successive minima algorithm. We show experimental results of using the scanning face, and compare with canonical loop orderings.

VanderWijngaart, Rob F.↗
NASA NTRS2010

NASA's Evolutionary Xenon Thruster (NEXT) Project Qualification Propellant Throughput Milestone: Performance, Erosion, and Thruster Service Life Prediction After 450 kg

The NASA s Evolutionary Xenon Thruster (NEXT) program is tasked with significantly improving and extending the capabilities of current state-of-the-art NSTAR thruster. The service life capability of the NEXT ion thruster is being assessed by thruster wear test and life-modeling of critical thruster components, such as the ion optics and cathodes. The NEXT Long-Duration Test (LDT) was initiated to validate and qualify the NEXT thruster propellant throughput capability. The NEXT thruster completed the primary goal of the LDT; namely to demonstrate the project qualification throughput of 450 kg by the end of calendar year 2009. The NEXT LDT has demonstrated 28,500 hr of operation and processed 466 kg of xenon throughput--more than double the throughput demonstrated by the NSTAR flight-spare. Thruster performance changes have been consistent with a priori predictions. Thruster erosion has been minimal and consistent with the thruster service life assessment, which predicts the first failure mode at greater than 750 kg throughput. The life-limiting failure mode for NEXT is predicted to be loss of structural integrity of the accelerator grid due to erosion by charge-exchange ions.

Herman, Daniel A.↗
NASA NTRS2010

Status of the NASA's Evolutionary Xenon Thruster (NEXT) Long-Duration Test After 30,352 Hours of Operation

The NASA s Evolutionary Xenon Thruster (NEXT) program is tasked with significantly improving and extending the capabilities of current state-of-the-art NSTAR thruster. The service life capability of the NEXT ion thruster is being assessed by thruster wear test and life-modeling of critical thruster components, such as the ion optics and cathodes. The NEXT Long-Duration Test (LDT) was initiated to validate and qualify the NEXT thruster propellant throughput capability. The NEXT thruster completed the primary goal of the LDT; namely to demonstrate the project qualification throughput of 450 kg by the end of calendar year 2009. The NEXT LDT has demonstrated 30,352 hr of operation and processed 490 kg of xenon throughput--surpassing the NSTAR Extended Life Test hours demonstrated and more than double the throughput demonstrated by the NSTAR flight-spare. Thruster performance changes have been consistent with a priori predictions. Thruster erosion has been minimal and consistent with the thruster service life assessment, which predicts the first failure mode at greater than 750 kg throughput. The life-limiting failure mode for NEXT is predicted to be loss of structural integrity of the accelerator grid due to erosion by charge-exchange ions.

Herman, Daniel A.↗
NASA NTRS2003

Prognostics Methodology for Complex Systems

An automatic method to schedule maintenance and repair of complex systems is produced based on a computational structure called the Informed Maintenance Grid (IMG). This method provides solutions to the two fundamental problems in autonomic logistics: (1) unambiguous detection of deterioration or impending loss of function and (2) determination of the time remaining to perform maintenance or other corrective action based upon information from the system. The IMG provides a health determination over the medium-to-longterm operation of the system, from one or more days to years of study. The IMG is especially applicable to spacecraft and both piloted and autonomous aircraft, or industrial control processes.

Gulati, Sandeep↗
NASA NTRS2011

Modeling of Supersonic Film Cooling on the J-2X Nozzle Extension

Supersonic film cooling (SSFC) of nozzles has been used in several liquid rocket engine designs, and is being applied to the nozzle extension (NE) of the J-2X upper stage engine currently under development. Turbine exhaust gas (TEG) is injected tangentially from a manifold along the NE, and provides a thermal barrier from the core nozzle flow for the NE. As the TEG stream mixes with the nozzle flow, the effectiveness of the thermal barrier is reduced. This paper documents computational fluid dynamics (CFD) analysis work performed by NASA Marshall Space Flight Center (MSFC) to model the flow of the TEG through the manifold, into the nozzle, and the subsequent mixing of the TEG stream with the core flow. The geometry and grid of the TEG manifold, structural support ribs, and the NE wall will be shown, and the CFD boundary conditions described. The Loci-CHEM CFD code used in this work will also be briefly described. A unique approach to modeling the combined TEG manifold/thrust chamber assembly (TCA) was employed, as it was not practical to model the entire 360 circumferential range in one simulation. Prior CFD validation work modeling Calspan SSFC experiments in the early 1990s, documented in a previous AIAA paper, will also be briefly discussed. The fluid dynamics of the TEG flow through the manifold, into and between the structural support ribs, and into the nozzlette that feeds the TCA will be described. Significant swirl and non-uniformities are present, which along with the wakes from the ribs, act to degrade the film cooling effectiveness compared to idealized injection of TEG gas. The effect of these flow characteristics on the adiabatic wall temperature profile on the NE will be discussed.

Ruf, Joseph H.↗
NASA NTRS2013

Investigations of Fluid-Structure-Coupling and Turbulence Model Effects on the DLR Results of the Fifth AIAA CFD Drag Prediction Workshop

The accurate calculation of aerodynamic forces and moments is of significant importance during the design phase of an aircraft. Reynolds-averaged Navier-Stokes (RANS) based Computational Fluid Dynamics (CFD) has been strongly developed over the last two decades regarding robustness, efficiency, and capabilities for aerodynamically complex configurations. Incremental aerodynamic coefficients of different designs can be calculated with an acceptable reliability at the cruise design point of transonic aircraft for non-separated flows. But regarding absolute values as well as increments at off-design significant challenges still exist to compute aerodynamic data and the underlying flow physics with the accuracy required. In addition to drag, pitching moments are difficult to predict because small deviations of the pressure distributions, e.g. due to neglecting wing bending and twisting caused by the aerodynamic loads can result in large discrepancies compared to experimental data. Flow separations that start to develop at off-design conditions, e.g. in corner-flows, at trailing edges, or shock induced, can have a strong impact on the predictions of aerodynamic coefficients too. Based on these challenges faced by the CFD community a working group of the AIAA Applied Aerodynamics Technical Committee initiated in 2001 the CFD Drag Prediction Workshop (DPW) series resulting in five international workshops. The results of the participants and the committee are summarized in more than 120 papers. The latest, fifth workshop took place in June 2012 in conjunction with the 30th AIAA Applied Aerodynamics Conference. The results in this paper will evaluate the influence of static aeroelastic wing deformations onto pressure distributions and overall aerodynamic coefficients based on the NASA finite element structural model and the common grids.

Keye, Stefan↗
NASA NTRS2009

Commercial off-the-Shelf (COTS) Refractory Material Evaluation

Fondu Fyre (FF) is currently the only refractory material qualified for use in the flame trench at KSC's Shuttle Launch Pads 39A and 39B. However, the material is not used as it was qualified and has undergone increasingly frequent and severe degradation due to the launch blasts. This degradation is costly as well as dangerous for launch infrastructure, crew and vehicle. FF is applied at the pad via the gunnite process, where wetted refractory material is sprayed onto a steel grid mounted on a support structure. The water content in this process can be manually adjusted by operators, causing distinct visual and physical discrepancies among repair areas. Since the application process is unlikely to change for new refractory materials, it is important to understand the effects of water content on commercial off-the-shelf (COTS) refractory materials. The purpose of this study was to evaluate the performance of the FF with respect to various water contents as well as heat treatments, to simulate aging and exposure to the blast. Initial results indicated that different water contents and heat treatments result in distinct differences in crushing strength, apparent porosity and bulk density. However, water content became an insignificant factor in both crush strength and porosity when FF was cured to at least 1500 deg. Additionally, inspection of the material's surface microstructure by scanning electron microscopy indicated distinguishable characteristics for different heat treatment levels. Results from this study will help guide future studies on the development and identification of new refractory materials.

Calle, Luz Marina↗
NASA NTRS2014

Computational Simulations of Convergent Nozzles for the AIAA 1st Propulsion Aerodynamics Workshop

Computational Fluid Dynamics (CFD) simulations were completed for a series of convergent nozzles in participation of the American Institute of Aeronautics and Astronautics (AIAA) 1st Propulsion Aerodynamics Workshop. The simulations were performed using the Wind-US flow solver. Discharge and thrust coefficients were computed for four axisymmetric nozzles with nozzle pressure ratios (NPR) ranging from 1.4 to 7.0. The computed discharge coefficients showed excellent agreement with available experimental data; the computed thrust coefficients captured trends observed in the experimental data, but over-predicted the thrust coefficient by 0.25 to 1.0 percent. Sonic lines were computed for cases with NPR >= 2.0 and agreed well with experimental data for NPR >= 2.5. Simulations were also performed for a 25 deg. conic nozzle bifurcated by a flat plate at NPR = 4.0. The jet plume shock structure was compared with and without the splitter plate to the experimental data. The Wind-US simulations predicted the shock structure well, though lack of grid resolution in the plume reduced the sharpness of the shock waves. Unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations and Detached Eddy Simulations (DES) were performed at NPR = 1.6 for the 25 deg conic nozzle with splitter plate. The simulations predicted vortex shedding from the trailing edge of the splitter plate. However, the vortices of URANS and DES solutions appeared to dissipate earlier than observed experimentally. It is believed that a lack of grid resolution in the region of the vortex shedding may have caused the vortices to break down too soon

DES↗