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,009 records · Page 56

Gradient-Based Aerodynamic Shape Optimization Using ADI Method for Large-Scale Problems

A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.

Pandya, Mohagna J.↗

Laminar Heating Validation of the OVERFLOW Code

OVERFLOW, a structured finite difference code, was applied to the solution of hypersonic laminar flow over several configurations assuming perfect gas chemistry. By testing OVERFLOW's capabilities over several configurations encompassing a variety of flow physics a validated laminar heating was produced. Configurations tested were a flat plate at 0 degrees incidence, a sphere, a compression ramp, and the X-38 re-entry vehicle. This variety of test cases shows the ability of the code to predict boundary layer flow, stagnation heating, laminar separation with re-attachment heating, and complex flow over a three-dimensional body. In addition, grid resolutions studies were done to give recommendations for the correct number of off-body points to be applied to generic problems and for wall-spacing values to capture heat transfer and skin friction. Numerical results show good comparison to the test data for all the configurations.

Lillard, Randolph P.↗

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.↗

Verification of a New NOAA/NSIDC Passive Microwave Sea-Ice Concentration Climate Record

A new satellite-based passive microwave sea-ice concentration product developed for the National Oceanic and Atmospheric Administration (NOAA)Climate Data Record (CDR) programme is evaluated via comparison with other passive microwave-derived estimates. The new product leverages two well-established concentration algorithms, known as the NASA Team and Bootstrap, both developed at and produced by the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC). The sea ice estimates compare well with similar GSFC products while also fulfilling all NOAA CDR initial operation capability (IOC) requirements, including (1) self describing file format, (2) ISO 19115-2 compliant collection-level metadata,(3) Climate and Forecast (CF) compliant file-level metadata, (4) grid-cell level metadata (data quality fields), (5) fully automated and reproducible processing and (6) open online access to full documentation with version control, including source code and an algorithm theoretical basic document. The primary limitations of the GSFC products are lack of metadata and use of untracked manual corrections to the output fields. Smaller differences occur from minor variations in processing methods by the National Snow and Ice Data Center (for the CDR fields) and NASA (for the GSFC fields). The CDR concentrations do have some differences from the constituent GSFC concentrations, but trends and variability are not substantially different.

Passive Microwave↗

Postflight aerothermodynamic analysis of Pegasus(tm) using computational fluid dynamic techniques

The objective was to validate the computational capability of the NASA Ames Navier-Stokes code, F3D, for flows at high Mach numbers using comparison flight test data from the Pegasus (tm) air launched, winged space booster. Comparisons were made with temperature and heat fluxes estimated from measurements on the wing surfaces and wing-fuselage fairings. Tests were conducted for solution convergence, sensitivity to grid density, and effects of distributing grid points to provide high density near temperature and heat flux sensors. The measured temperatures were from sensors embedded in the ablating thermal protection system. Surface heat fluxes were from plugs fabricated of highly insulative, nonablating material, and mounted level with the surface of the surrounding ablative material. As a preflight design tool, the F3D code produces accurate predictions of heat transfer and other aerodynamic properties, and it can provide detailed data for assessment of boundary layer separation, shock waves, and vortex formation. As a postflight analysis tool, the code provides a way to clarify and interpret the measured results.

Kuhn, Gary D.↗

High-latitude filtering in a global grid-point model using model normal modes

A normal modes expansion technique is applied to perform high latitude filtering in the GLAS fourth order global shallow water model with orography. The maximum permissible time step in the solution code is controlled by the frequency of the fastest propagating mode, which can be a gravity wave. Numerical methods are defined for filtering the data to identify the number of gravity modes to be included in the computations in order to obtain the appropriate zonal wavenumbers. The performances of the model with and without the filter, and with a time tendency and a prognostic field filter are tested with simulations of the Northern Hemisphere winter. The normal modes expansion technique is shown to leave the Rossby modes intact and permit 3-5 day predictions, a range not possible with the other high-latitude filters.

Takacs, L. L.↗

Expanding the Time Dimension of Hyperspectral Infrared Sounding Observations: Designing the NUCAPS-Forecast System

Low Earth orbit (LEO) hyper-spectral infrared (IR) sounders have significant yet untapped potential for characterizing thermodynamic environments of convective initiation and ongoing convection. While LEO soundings are of enormous value to weather forecasters, they cannot provide the temporal resolution needed to resolve the rapidly evolving thermodynamics of the convective environment. A novel methodology was developed to extend snapshots of LEO soundings forward in time up to six hours into the future and has been made available within the NOAA National Weather Service Advanced Weather Interactive Processing System (AWIPS) for testing in a quasi-operational setting. This methodology is based on parcel forward-trajectory calculations from the satellite observing time to recreate future soundings of temperature and moisture at regularly gridded intervals in space and time. The soundings are based on NOAA-Unique Combined Atmospheric Processing System (NUCAPS) retrievals from the Suomi NPP and NOAA-20 satellite platforms to allow for rapid near-real time product processing. Science code was transitioned to operational production in 2019 and the product, NUCAPS-Forecast, was tested at the NOAA Hazardous Weather Testbed Spring Experiment in 2019 and 2021. Forecaster feedback has led to improving the product and processing system to address product limitations and improve its applicability for pre-convective forecasting. This presentation highlights the design of the NUCAPS-Forecast system. Solutions were tested and implemented to remove excess moisture and conjoin forward-trajected parcels into vertical profiles for calculating grids of stability indices. Additionally, use cases that explore the tendencies of derived convective available potential energy (CAPE) and convective inhibition (CIN) were investigated. When evaluated against gridded, hourly accumulated rainfall obtained from the Multi-Radar Multi-Sensor (MRMS) observational data set areas with elevated CAPE (reduced CIN) are shown to be associated with areas of precipitation and increase confidence the LEO soundings advected forward in time physically represent atmospheric processes. These results suggest that adiabatic parcel advection of LEO satellite sounding snapshots forward in time exhibits skill in discriminating locations and times that are more likely to lead to convective initiation. Leveraging trajectory modeling to increase the spatial and temporal resolution of LEO IR sounding observations today demonstrates the value of these observations from the future constellation of JPSS satellites and GeoXO mission.

weather↗

Numerical solution of Space Shuttle Orbiter flow field

The supersonic, viscous laminar flow around the Space Shuttle Orbiter forebody has been computed with a parabolized Navier-Stokes code using a generalized coordinate transformation. The initial solution for the nose part of the Orbiter geometry was obtained with a three-dimensional time-dependent Navier-Stokes solver. It was necessary to employ a wind axis oriented coordinate system to obtain the initial solution with the time-dependent code. The generalized PNS technique used in this study allows the solution to be marched from the given initial data surface to any desired surface downstream. A grid point clustering scheme was employed to accurately describe the body shape by clustering points at the wing tip and at the wing body juncture. The computed heat transfer coefficients, pressure coefficients, and shock shapes are compared with the available experimental data for 0 degrees and 30 degrees angle of attack.

Venkatapathy, E.↗

A point implicit unstructured grid solver for the Euler and Navier-Stokes equations

An upwind finite element technique that uses cell centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured triangular grids. A structured grid of quadrilaterals is laid out near the body surface. For inviscid flows the method is stable at Courant numbers of over 100,000. A first order basic scheme and a higher order flux corrected transport (FCT) scheme have been implemented. This technique has been applied to the problem of predicting type III and IV shock wave interactions on a cylinder, with a view of simulating the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively generated grids. The adaptive mesh generator and the upwind Navier-Stokes solver are combined in a set of programs called LARCNESS, an acronym for Langley Adaptive Remeshing Code and Navier-Stokes Solver.

Thareja, Rajiv R.↗

Development of a defect stream function, law of the wall/wake method for compressible turbulent boundary layers

The method presented is designed to improve the accuracy and computational efficiency of existing numerical methods for the solution of flows with compressible turbulent boundary layers. A compressible defect stream function formulation of the governing equations assuming an arbitrary turbulence model is derived. This formulation is advantageous because it has a constrained zero-order approximation with respect to the wall shear stress and the tangential momentum equation has a first integral. Previous problems with this type of formulation near the wall are eliminated by using empirically based analytic expressions to define the flow near the wall. The van Driest law of the wall for velocity and the modified Crocco temperature-velocity relationship are used. The associated compressible law of the wake is determined and it extends the valid range of the analytical expressions beyond the logarithmic region of the boundary layer. The need for an inner-region eddy viscosity model is completely avoided. The near-wall analytic expressions are patched to numerically computed outer region solutions at a point determined during the computation. A new boundary condition on the normal derivative of the tangential velocity at the surface is presented; this condition replaces the no-slip condition and enables numerical integration to the surface with a relatively coarse grid using only an outer region turbulence model. The method was evaluated for incompressible and compressible equilibrium flows and was implemented into an existing Navier-Stokes code using the assumption of local equilibrium flow with respect to the patching. The method has proven to be accurate and efficient.

Wahls, Richard A.↗

A study of transonic aerodynamic analysis methods for use with a hypersonic aircraft synthesis code

A means of performing routine transonic lift, drag, and moment analyses on hypersonic all-body and wing-body configurations were studied. The analysis method is to be used in conjunction with the Hypersonic Vehicle Optimization Code (HAVOC). A review of existing techniques is presented, after which three methods, chosen to represent a spectrum of capabilities, are tested and the results are compared with experimental data. The three methods consist of a wave drag code, a full potential code, and a Navier-Stokes code. The wave drag code, representing the empirical approach, has very fast CPU times, but very limited and sporadic results. The full potential code provides results which compare favorably to the wind tunnel data, but with a dramatic increase in computational time. Even more extreme is the Navier-Stokes code, which provides the most favorable and complete results, but with a very large turnaround time. The full potential code, TRANAIR, is used for additional analyses, because of the superior results it can provide over empirical and semi-empirical methods, and because of its automated grid generation. TRANAIR analyses include an all body hypersonic cruise configuration and an oblique flying wing supersonic transport.

Sandlin, Doral R.↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 5: Unsteady counterrotation ducted propfan analysis. Computer program user's manual

The primary objective of this study was the development of a time-marching three-dimensional Euler/Navier-Stokes aerodynamic analysis to predict steady and unsteady compressible transonic flows about ducted and unducted propfan propulsion systems employing multiple blade rows. The computer codes resulting from this study are referred to as ADPAC-AOACR (Advanced Ducted Propfan Analysis Codes-Angle of Attack Coupled Row). This report is intended to serve as a computer program user's manual for the ADPAC-AOACR codes developed under Task 5 of NASA Contract NAS3-25270, Unsteady Counterrotating Ducted Propfan Analysis. The ADPAC-AOACR program is based on a flexible multiple blocked grid discretization scheme permitting coupled 2-D/3-D mesh block solutions with application to a wide variety of geometries. For convenience, several standard mesh block structures are described for turbomachinery applications. Aerodynamic calculations are based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. Steady flow predictions are accelerated by a multigrid procedure. Numerical calculations are compared with experimental data for several test cases to demonstrate the utility of this approach for predicting the aerodynamics of modern turbomachinery configurations employing multiple blade rows.

Hall, Edward J.↗

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 5: Unsteady counterrotation ducted propfan analysis

The primary objective of this study was the development of a time-marching three-dimensional Euler/Navier-Stokes aerodynamic analysis to predict steady and unsteady compressible transonic flows about ducted and unducted propfan propulsion systems employing multiple blade rows. The computer codes resulting from this study are referred to as ADPAC-AOAR\CR (Advanced Ducted Propfan Analysis Codes-Angle of Attack Coupled Row). This document is the final report describing the theoretical basis and analytical results from the ADPAC-AOACR codes developed under task 5 of NASA Contract NAS3-25270, Unsteady Counterrotating Ducted Propfan Analysis. The ADPAC-AOACR Program is based on a flexible multiple blocked grid discretization scheme permitting coupled 2-D/3-D mesh block solutions with application to a wide variety of geometries. For convenience, several standard mesh block structures are described for turbomachinery applications. Aerodynamic calculations are based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. Steady flow predictions are accelerated by a multigrid procedure. Numerical calculations are compared with experimental data for several test cases to demonstrate the utility of this approach for predicting the aerodynamics of modern turbomachinery configurations employing multiple blade rows.

Hall, Edward J.↗

Full 3D Analysis of the GE90 Turbofan Primary Flowpath

The multistage simulations of the GE90 turbofan primary flowpath components have been performed. The multistage CFD code, APNASA, has been used to analyze the fan, fan OGV and booster, the 10-stage high-pressure compressor and the entire turbine system of the GE90 turbofan engine. The code has two levels of parallel, and for the 18 blade row full turbine simulation has 87.3 percent parallel efficiency with 121 processors on an SGI ORIGIN. Grid generation is accomplished with the multistage Average Passage Grid Generator, APG. Results for each component are shown which compare favorably with test data.

Turner, Mark G.↗

Helicon Plasma Source Optimization Studies for VASIMR

A helicon plasma source at Oak Ridge National Laboratory is being used to investigate operating scenarios relevant to the VASIMR (VAriable Specific Impulse Magnetoplasma Rocket). These include operation at high magnetic field (> = 0.4 T), high frequency (<= 30 MHz), high power (< = 3 kW), and with light ions (He+, H+). To date, He plasmas have been produced with n(sub e0) = 1.7 x 10(exp 19)/cu m (measured with an axially movable 4mm microwave interferometer), with Pin = I kW at f = 13.56 MHz and absolute value of B(sub 0) = 0.16 T. In the near future, diagnostics including a mass flow meter and a gridded energy analyzer array will be added to investigate fueling efficiency and the source power balance. The latest results, together with modeling results using the EMIR rf code, will be presented.

Goulding, R. H.↗

Multi-Megawatt Gas Turbine Power Systems for Lunar Colonies

A concept for development of second generation 10 MWe prototype lunar power plant utilizing a gas cooled fission reactor supplying heated helium working fluid to two parallel 5 MWe closed cycle gas turbines is presented. Such a power system is expected to supply the energy needs for an initial lunar colony with a crew of up to 50 persons engaged in mining and manufacturing activities. System performance and mass details were generated by an author developed code (BRMAPS). The proposed pilot power plant can be a model for future plants of the same capacity that could be tied to an evolutionary lunar power grid.

Juhasz, Albert J.↗

Isotopic Effects in Nuclear Fragmentation and GCR Transport Problems

Improving the accuracy of the galactic cosmic ray (GCR) environment and transport models is an important goal in preparing for studies of the projected risks and the efficiency of potential mitigations methods for space exploration. In this paper we consider the effects of the isotopic composition of the primary cosmic rays and the isotopic dependence of nuclear fragmentation cross sections on GCR transport models. Measurements are used to describe the isotopic composition of the GCR including their modulation throughout the solar cycle. The quantum multiple-scattering approach to nuclear fragmentation (QMSFRG) is used as the data base generator in order to accurately describe the odd-even effect in fragment production. Using the Badhwar and O'Neill GCR model, the QMSFRG model and the HZETRN transport code, the effects of the isotopic dependence of the primary GCR composition and on fragment production for transport problems is described for a complete GCR isotopic-grid. The principle finding of this study is that large errors ( 100%) will occur in the mass-flux spectra when comparing the complete isotopic-grid (141 ions) to a reduced isotopic-grid (59 ions), however less significant errors 30%) occur in the elemental-flux spectra. Because the full isotopic-grid is readily handled on small computer work-stations, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗

Aeroacoustic Simulations of a Nose Landing Gear Using FUN3D on Pointwise Unstructured Grids

Numerical simulations have been performed for a partially-dressed, cavity-closed (PDCC) nose landing gear configuration that was tested in the University of Florida's open-jet acoustic facility known as the UFAFF. The unstructured-grid flow solver FUN3D is used to compute the unsteady flow field for this configuration. Mixed-element grids generated using the Pointwise(TradeMark) grid generation software are used for these simulations. Particular care is taken to ensure quality cells and proper resolution in critical areas of interest in an effort to minimize errors introduced by numerical artifacts. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these simulations. Solutions are also presented for a wall function model coupled to the standard turbulence model. Time-averaged and instantaneous solutions obtained on these Pointwise grids are compared with the measured data and previous numerical solutions. The resulting CFD solutions are used as input to a Ffowcs Williams-Hawkings noise propagation code to compute the farfield noise levels in the flyover and sideline directions. The computed noise levels compare well with previous CFD solutions and experimental data.

Vatsa, Veer N.↗