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 1,153 records · Page 64

Three-dimensional multigrid Navier-Stokes computations for turbomachinery applications

The fully three-dimensional, time-dependent compressible Navier-Stokes equations in cylindrical coordinates are presently used, in conjunction with the multistage Runge-Kutta numerical integration scheme for solution of the governing flow equations, to simulate complex flowfields within turbomechanical components whose pertinent effects encompass those of viscosity, compressibility, blade rotation, and tip clearance. Computed results are presented for selected cascades, emphasizing the code's capabilities in the accurate prediction of such features as airfoil loadings, exit flow angles, shocks, and secondary flows. Computations for several test cases have been performed on a Cray-YMP, using nearly 90,000 grid points.

Subramanian, S. V.↗

Experimental and computational studies of hovering rotor flows

HELIX-I, an essentially standard full-potential CFD helicopter rotor code, is unique in its use of the Clebisch kinematical flow description to specify a freely convecting wake and its capacity for predicting hover performance. A study is presently performed to assess the method's sensitivity to grids and solution-starting techniques. The effects of these parameters on thrust, power, load distribution, and wake geometry are ascertained and compared with an extensive rotor data base. The use of a fairly accurate starting solution yields no obvious advantage over the use of a novel starting method which employs a succession of diminishing artificial flows.

Nsi Mba, M.↗

Computational study of generic hypersonic vehicle flow fields

The geometric data of the generic hypersonic vehicle configuration included body definitions and preliminary grids for the forebody (nose cone excluded), midsection (propulsion system excluded), and afterbody sections. This data was to be augmented by the nose section geometry (blunt conical section mated with the noncircular cross section of the forebody initial plane) along with a grid and a detailed supersonic combustion ramjet (scramjet) geometry (inlet and combustor) which should be merged with the nozzle portion of the afterbody geometry. The solutions were to be obtained by using a Navier-Stokes (NS) code such as TUFF for the nose portion, a parabolized Navier-Stokes (PNS) solver such as the UPS and STUFF codes for the forebody, a NS solver with finite rate hydrogen-air chemistry capability such as TUFF and SPARK for the scramjet and a suitable solver (NS or PNS) for the afterbody and external nozzle flows. The numerical simulation of the hypersonic propulsion system for the generic hypersonic vehicle is the major focus of this entire work. Supersonic combustion ramjet is such a propulsion system, hence the main thrust of the present task has been to establish a solution procedure for the scramjet flow. The scramjet flow is compressible, turbulent, and reacting. The fuel used is hydrogen and the combustion process proceeds at a finite rate. As a result, the solution procedure must be capable of addressing such flows.

Narayan, Johnny R.↗

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 1: Code validation and stationary accretors

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing stationary spheres of varying sizes (from 10.0 down to 0.02 Bondi radii) accrete matter from a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3 or 1.2). To accommodate the long-range gravitational forces, the extent of the computational volume is typically a factor of 100 larger than the radius of the accretor. We compare the numerical mass accretion rates with the theoretical predictions of Bondi, to assess the validity of the code. The hydrodynamics is modeled by the piecewise parabolic method. No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (6-8) grids around the stationary sphere, each finer grid being a factor of 2 smaller spatially than the next coarser grid. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid while at the same time evolving the gas on the coarser grids. The accretion rates derived numerically are in in very good agreement (to about 10% over several orders of magnitude) with the values given by Bondi for a stationary accretor within a hydrodynamic medium. However, the equations have to be changed in order to include the finite size of the accretor (in some cases very large compared to the sonic point or even to the Bondi radius).

Ruffert, Maximilian↗

Forecasting of Storm Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite element based codes. It affords a capability for simulating tidal circulation and storm surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate.

Valenti, Elizabeth↗

NASA Tech Briefs, September 2005

Topivs include: Diamond-Coated Carbon Nanotubes for Efficient Field Emission; Improved Anode Coatings for Direct Methanol Fuel Cells; Advanced Ablative Insulators and Methods of Making Them; PETIs as High-Temperature Resin-Transfer-Molding Materials; Stable Polyimides for Terrestrial and Space Uses; Low-Density, Aerogel-Filled Thermal-Insulation Tiles; High-Performance Polymers Having Low Melt Viscosities; Nonflammable, Hydrophobic Aerogel Composites for Insulation; Front-Side Microstrip Line Feeding a Raised Antenna Patch; Medium-Frequency Pseudonoise Georadar; Facilitating Navigation Through Large Archives; Program for Weibull Analysis of Fatigue Data; Comprehensive Micromechanics-Analysis Code - Version 4.0; Component-Based Visualization System; Software for Engineering Simulations of a Spacecraft; LabVIEW Interface for PCI-SpaceWire Interface Card; Path Following with Slip Compensation for a Mars Rover; International Space Station Electric Power System Performance Code-SPACE; Software for Automation of Real-Time Agents, Version 2; Software for Optimizing Plans Involving Interdependent Goals; Computing Gravitational Fields of Finite-Sized Bodies; Custom Sky-Image Mosaics from NASA's Information Power Grid; ANTLR Tree Grammar Generator and Extensions; Generic Kalman Filter Software; Alignment Stage for a Cryogenic Dilatometer; Rugged Iris Mechanism; Treatments To Produce Stabilized Aluminum Mirrors for Cryogenic Uses; Making AlNx Tunnel Barriers Using a Low-Energy Nitrogen-Ion Beam; Making Wide-IF SIS Mixers with Suspended Metal-Beam Leads; Sol-Gel Glass Holographic Light-Shaping Diffusers; Automated Counting of Particles To Quantify Cleanliness; Phase Correction for GPS Antenna with Nonunique Phase Center; Compact Infrasonic Windscreen; Broadband External-Cavity Diode Laser; High-Efficiency Solar Cells Using Photonic-Bandgap Materials; Generating Solid Models from Topographical Data; Computationally Lightweight Air-Traffic-Control Simulation; Spool Valve for Switching Air Flows Between Two Beds; Partial Model of Insulator/ Insulator Contact Charging; Asymmetric Electrostatic Radiation Shielding for Spacecraft; and Reusable Hybrid Propellant Modules for Outer-Space Transport.

Source record↗

Performance Tests of Snow-Related Variables Over the Tibetan Plateau and Himalayas Using a New Version of NASA GEOS-5 Land Surface Model that Includes the Snow Darkening Effect

The snow darkening effect, i.e. the reduction of snow albedo, is caused by absorption of solar radiation by absorbing aerosols (dust, black carbon, and organic carbon) deposited on the snow surface. This process is probably important over Himalayan and Tibetan glaciers due to the transport of highly polluted Atmospheric Brown Cloud (ABC) from the Indo-Gangetic Plain (IGP). This effect has been incorporated into the NASA Goddard Earth Observing System model, version 5 (GEOS-5) atmospheric transport model. The Catchment land surface model (LSM) used in GEOS-5 considers 3 snow layers. Code was developed to track the mass concentration of aerosols in the three layers, taking into account such processes as the flushing of the compounds as liquid water percolates through the snowpack. In GEOS-5, aerosol emissions, transports, and depositions are well simulated in the Goddard Chemistry Aerosol Radiation and Transport (GO CART) module; we recently made the connection between GOCART and the GEOS-5 system fitted with the revised LSM. Preliminary simulations were performed with this new system in "replay" mode (i.e., with atmospheric dynamics guided by reanalysis) at 2x2.5 degree horizontal resolution, covering the period 1 November 2005 - 31 December 2009; we consider the final three years of simulation here. The three simulations used the following variants of the LSM: (1) the original Catchment LSM with a fixed fresh snowfall density of 150 kg m-3 ; (2) the LSM fitted with the new snow albedo code, used here without aerosol deposition but with changes in density formulation and melting water effect on snow specific surface area, (3) the LSM fitted with the new snow albedo code as same as (2) but with fixed aerosol deposition rates (computed from GOCART values averaged over the Tibetan Plateau domain [Ion.: 60-120E; lat.: 20-50N] during March-May 2008) applied to all grid points at every time step. For (2) and (3), the same setting on the fresh snowfall density as in (1) was used.

Yasunari, Tppei J.↗

Computational and Experimental Study of Supersonic Nozzle Flow and Aft-Deck Interactions

NASA has been conducting research into reducing sonic boom and changing FAA regulations to allow for supersonic commercial transport over land in the United States. This particular study looks at a plume passing through a shock generated from an aft deck on a nacelle; the aft deck is meant to represent the trailing edge of a wing. NASA Langley Research Center USM3D CFD code results are compared to the experimental data taken at the NASA Glenn Research Center 1-foot by 1-foot Supersonic Wind Tunnel. This study included examining two turbulence models along with different volume sourcing methods for grid generation. The results show that using the k-epsilon turbulence model within USM3D produced shock signatures that closely follow the experimental data at a variety of nozzle pressure ratio settings.

Bruce, Walter E., IV↗

Implementation of a Body Force Model in OVERFLOW for Propulsor Simulations

We present an implementation of a propulsor model based on body force method into the OVERFLOW computational fluid dynamics code to model turbofan engines and propulsors of similar type. The model estimates the forces imparted on the fluid by the blade camber surfaces as a body force source terms in the momentum and energy equations over grids that represent the rotor and stator of a fan stage. We tested the implementation on three cases: 1) Source Diagnostics Test (SDT) fan with R4 rotors, which has extensive test data on aerodynamic performance and rotor wake surveys, 2) A stand-alone Aeronaut TF8000 propulsor which is an off-the-shelf propulsor used on model aircraft and 3) The D8 aircraft model with TF8000 propulsors placed in a wind tunnel. Despite missing some of the features in the velocity profiles near the endwalls (i.e. hub and casing surfaces), the OVERFLOW simulations with body force model predicted area-averaged flow speed and total pressure rise through the SDT propulsor within a few percent of the LDV (Laser Doppler Velocimetry) measurements. In the case of TF8000 propulsor on the D8 airframe, the model under-predicted mechanical power coefficient by several percent of the wind tunnel test results when the horizontal force balance condition over the airframe is targeted by tuning rotor speed. By investigating the upstream influence of the rotor swirl, it was found out that the induced swirl velocity effects upstream were relatively small and they rapidly vanished before reaching one fan diameter upstream of the fan face. The body force model provided insights on aerothermodynamics and aeromechanics of boundary layer ingesting propulsor; these insights could not be obtained by using the uniform pressure jump model.

Propulsor↗

Isolated Nozzle Analysis Using OVERFLOW

LaRC conducted a code validation study for the OVERFLOW code to ascertain its accuracy for boattail drag prediction. The OVERFLOW results compared favorably with the LaRC 16-ft. Transonic Wind Tunnel (TWT) data, and prior CFD solutions from PAB3D and CFL3D. The ultimate goal is to investigate the installation drag of the nacelle boattails with powered nozzles at transonic mach numbers. The OVERFLOW solver was chosen because of its ability to accept volume overlapping structured grid for very complex airframe configurations. Structured grid components for representing the transonic nozzle boattail can be added to the BCAG grid for a TCA airframe with 2D bifurcated inlet and flow through nacelle without alteration. The focus of this research was to determine the suitability of the OVERFLOW solver for accomplishing this ultimate goal. This presentation will first introduce the transonic nozzle boattail wind-tunnel model geometry, followed by an examination of aerodynamic features based on the current OVERFLOW solutions and the solutions obtained previously using PAB3D, comparisons of Cp on the flap surface between the OVERFLOW solutions, wind tunnel data, and solutions from other CFD codes, an assessment of boattail drag count prediction, and a work plan for FY99.

Wong, Tin-Chee↗

Developing Information Power Grid Based Algorithms and Software

This exploratory study initiated our effort to understand performance modeling on parallel systems. The basic goal of performance modeling is to understand and predict the performance of a computer program or set of programs on a computer system. Performance modeling has numerous applications, including evaluation of algorithms, optimization of code implementations, parallel library development, comparison of system architectures, parallel system design, and procurement of new systems. Our work lays the basis for the construction of parallel libraries that allow for the reconstruction of application codes on several distinct architectures so as to assure performance portability. Following our strategy, once the requirements of applications are well understood, one can then construct a library in a layered fashion. The top level of this library will consist of architecture-independent geometric, numerical, and symbolic algorithms that are needed by the sample of applications. These routines should be written in a language that is portable across the targeted architectures.

Dongarra, Jack↗

Extended Conformal Mappings for Supersonic Aircraft Calculations

A gasdynamical analysis is presented for a three-dimensional, supersonic, inviscid, steady, shockless flow past an arbitrary airframe, using computational grids. The analysis, which includes special treatments for body points and bow-shock points, relies on the equations of motion written in terms of logarithm of pressure entropy and two angles reprising the velocity vector. A FORTRAN code was employed (sample geometry is shown). The flow analysis is considered reliable for aircraft cross-sections that are elliptical in shape.

Moretti, G.↗

Some recent applications of XTRAN3S

A time marching finite difference code, XTRAN3S that solves the three dimensional transonic small perturbation equation for flow over isolated wings was developed. During initial applications of the program, problems were encountered in the prediction of unsteady forces. The use of a revised grid and force calculation scheme improved those predictions. Comparisons are made between predicted and experimental pressure data for a rectangular supercritical wing. Comparisons of steady and unsteady data at M sub infinity = 0.700 show good agreement between calculated and experimental values. A comparison of steady data at M sub infinity 0.825 shows poor agreement between calculations and experiment. Program difficulties were encountered with swept and tapered configurations.

Seidel, D. A.↗

Development of direct-inverse 3-D methods for applied transonic aerodynamic wing design and analysis

Progress in the direct-inverse wing design method in curvilinear coordinates has been made. This includes the remedying of a spanwise oscillation problem and the assessment of grid skewness, viscous interaction, and the initial airfoil section on the final design. It was found that, in response to the spanwise oscillation problem that designing at every other spanwise station produced the best results for the cases presented, a smoothly varying grid is especially needed for the accurate design at the wing tip, the boundary layer displacement thicknesses must be included in a successful wing design, the design of high and medium aspect ratio wings is possible with this code, and the final airfoil section designed is fairly independent of the initial section.

Carlson, Leland A.↗

Three-dimensional Navier-Stokes heat transfer predictions for turbine blade rows

Results are shown for a three-dimensional Navier-Stokes analysis of both the flow and the surface heat transfer for turbine applications. Heat transfer comparisons are made with the experimental shock-tunnel data of Dunn and Kim, and with the data of Blair for the rotor of the large scale rotating turbine. The analysis was done using the steady-state, three-dimensional, thin-layer Navier-Stokes code developed by Chima, which uses a multistage Runge-Kutta scheme with implicit residual smoothing. An algebraic mixing length turbulence model is used to calculate turbulent eddy viscosity. The variation in heat transfer due to variations in grid parameters is examined. The effects of rotation, tip clearance, and inlet boundary layer thickness variation on the predicted blade and endwall heat transfer are examined.

Boyle, R. J.↗

Three-dimensional Navier-Stokes heat transfer predictions for turbine blade rows

Results are shown for a three-dimensional Navier-Stokes analysis of both the flow and the surface heat transfer for turbine applications. Heat transfer comparisons are made with the experimental shock-tunnel data of Dunn and Kim, and with the data of Blair for the rotor of the large scale rotating turbine. The analysis was done using the steady-state, three-dimensional, thin-layer Navier-Stokes code developed by Chima, which uses a multistage Runge-Kutta scheme with implicit residual smoothing. An algebraic mixing length turbulence model is used to calculate turbulent eddy viscosity. The variation in heat transfer due to variations in grid parameters is examined. The effects of rotation, tip clearance, and inlet boundary layer thickness variation on the predicted blade and endwall heat transfer are examined.

Boyle, Robert J.↗

Computational Fluid Dynamic Modeling of Rocket Based Combined Cycle Engine Flowfields

Computational Fluid Dynamic techniques are used to study the flowfield of a fixed geometry Rocket Based Combined Cycle engine operating in rocket ejector mode. Heat addition resulting from the combustion of injected fuel causes the subsonic engine flow to choke and go supersonic in the slightly divergent combustor-mixer section. Reacting flow computations are undertaken to predict the characteristics of solutions where the heat addition is determined by the flowfield. Here, adaptive gridding is used to improve resolution in the shear layers. Results show that the sonic speed is reached in the unheated portions of the flow first, while the heated portions become supersonic later. Comparison with results from another code show reasonable agreement. The coupled solutions show that the character of the combustion-based thermal choking phenomenon can be controlled reasonably well such that there is opportunity to optimize the length and expansion ratio of the combustor-mixer.

Daines, Russell L.↗