Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid analysis and integration”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Annual Research Briefs, 2002: Center for Turbulence Research

Turbulent combustion remains the largest component of the CTR's core program. This program and several related activities at CTR are supported by NASA's Ultra Efficient Engine Technology Program. It is also intimately connected with the Department of Energy's ASCI program at Stanford which develops the technology for numerical simulation of realistic aircraft engines using state of the art massively parallel computers. In combustion modeling the attention has been directed to the modeling of higher levels of complexity such as spray dynamics, radiation and soot formation. Major aircraft engine manufacturers have shown considerable interest in this program; in particular, a significant active collaboration exists between CTR and the Pratt & Whitney Corporation. CTR's combustion program is essentially based on the large-eddy simulation technique, LES, which is actively being pursued at CTR for this and many other applications. Important accomplishments in LES included simulations with three-dimensional filters which result in grid independent calculations (that is why we call it "true" LES), and the development of the methodology for integration of LES and Reynolds Averaged computations. Optimization techniques are being studied and used for the important problem of wall boundary conditions for LES as well as for optimal shape design for aeroacoustic and aerodynamic performance gains.

ULTRA EFFICIENT ENGINE TECHNOLOGY↗
NASA NTRS2023

Drought and Deluge: the Recurrence of Hydroclimate Extremes During the Past 600 Years in Eastern Australia’s Natural Resource Management (NRM) Clusters

Recent extremes of flood and drought across Australia have raised questions about the recurrence of such rare events and highlighted the importance of understanding multi-decadal climate variability. However, instrumental records over the past century are too short to adequately characterise climate variability on multi-decadal and longer timescales or robustly determine extreme event frequencies and their duration. Palaeoclimate reconstructions can provide much-needed information to help address this problem. Here, we use the 600-year hydroclimate record captured in the eastern Australian and New Zealand Drought Atlas (ANZDA) to analyse drought and pluvial frequency trends for East Australian Natural Resource Management (NRM) clusters. This partitioning of the drought atlas grid points into recognised biophysical areas (i.e. NRM clusters) revealed their differences and similarities in drought intensity and pluvial events over time. We find sustained multi-decadal periods of a wet–dry geographic 'seesaw' between eastern to central and southern NRMs (e.g. 1550–1600 CE and 1700–1750 CE). In contrast, other periods reveal spatially consistent wetting (e.g. 1500–1550 CE) or drying (e.g. 1750–1800 CE). Emerging hot spot analysis further shows that some areas that appear naturally buffered from severe drought during the instrumental period have a greater exposure risk when the longer 600-year record is considered. These findings are particularly relevant to management plans when dealing with the impacts of climate extremes developed at regional scales. Our results demonstrate that integrating and extending instrumental records with palaeoclimate datasets will become increasingly important for developing robust and locally specific extreme event frequency information for managing water resources.

Australian hydroclimate↗
NASA NTRS2015

Biomaterials Out of Thin Air: in Situ, On-Demand Printing of Advanced Biocomposites

Upmass is the single most significant limitation of our current space mission capability. Although biomaterials and biocomposites have mass, strength, flexibility, and self-healing properties that could significantly reduce upmass, their use is limited by the following drawbacks: Expensive, specific production. Many biomaterials can only be produced as part of significant support ecosystem; Inaccessible functional customization. The grain of wood, the porosity of bone, and so on are an integral part of the materials' desired mechanical properties, but are not deterministic when the material is naturally grown; Limited compositions. Most biomaterials (unlike metal, plastic, etc.) cannot be easily combined or modified to produce new materials. This project builds on recent advances in: Synthetic biology. Libraries of standardized genetic parts which can be used for controlled cellular material production, delivery, and binding; 3D printing. Commercial off-the-shelf components which can be used to make of a pico- to nanoliter cell deposition system; Tissue engineering. Proven cell-compatible support hydrogels and scaffolds can be modified to bind the deposited biomaterials of interest. Objectives: Feasibility and benefit analysis. Two mission contexts span the concept's scope (see below); Proof-of-concept demonstration. A simple grid of two proteins, fluorescent for easy detection, to validate the core technology concept; Proposed implementations for follow-on work. Avenues for future work on each core component (host cell, production control, material delivery, material binding, etc.); Complementary studies exploration. A survey of other emerging areas (in situ resource utilization, protein engineering, etc.) with the potential to multiply our technology's impact. Potential Impacts: This application could dramatically expand manufacturing capabilities on Earth and in space: In situ resource utilization. A far greater range of materials and products will be available from the limited palette offered by in situ resource extraction techniques; Reduced equipment and material upmass for off-Earth habitats. Ready- to-use highly specialized construction materials (radiation hardened, compressive/tensile, light or dense) from an extremely low starting mass; Structured biomaterial production. New ready-to-use macro, micro, and molecular manufacturing techniques for traditional biomaterials such as wood, bone and shell; New and novel biocomposite creation. The ability to create completely novel material composites from any base material that cells can be engineered to produce. Suggested Mission Contexts: ISS part manufacturing. A 'minimal working example' making a finished biomaterial part aboard the International Space Station; A long-term Mars habitat. 'Cradle-to-grave' use at a hypothetical Mars habitat, covering everything from tools to construction materials. Alternate Abstract: Imagine being able to print anything from tools and composite building materials to food and human tissues. Imagine being on Mars with the ability to replace any broken part, whether it's a part of your spacesuit, your habitat, or your own body. We propose a technique that would allow just that. By printing 3D arrays of cells engineered to secrete the necessary materials, the abundant in situ resources of atmosphere and regolith become organic, inorganic, or organic-inorganic composite materials. Such materials include novel, biologically derived materials not previously possible to fabricate.

habitats↗
NASA NTRS2015

Biomaterials Out of Thin Air: in Situ, On-demand Printing of Advanced Biocomposites

Upmass is the single most significant limitation of our current space mission capability. Although biomaterials and biocomposites have mass, strength, flexibility, and self-healing properties that could significantly reduce upmass, their use is limited by the following drawbacks: Expensive, specific production. Many biomaterials can only be produced as part of significant support ecosystem; Inaccessible functional customization. The grain of wood, the porosity of bone, and so on are an integral part of the materials' desired mechanical properties, but are not deterministic when the material is naturally grown; Limited compositions. Most biomaterials (unlike metal, plastic, etc.) cannot be easily combined or modified to produce new materials. This project builds on recent advances in: Synthetic biology. Libraries of standardized genetic parts which can be used for controlled cellular material production, delivery, and binding; 3D printing. Commercial off-the-shelf components which can be used to make of a pico- to nanoliter cell deposition system; Tissue engineering. Proven cell-compatible support hydrogels and scaffolds can be modified to bind the deposited biomaterials of interest. Objectives: Feasibility and benefit analysis. Two mission contexts span the concept's scope (see below); Proof-of-concept demonstration. A simple grid of two proteins, fluorescent for easy detection, to validate the core technology concept; Proposed implementations for follow-on work. Avenues for future work on each core component (host cell, production control, material delivery, material binding, etc.); Complementary studies exploration. A survey of other emerging areas (in situ resource utilization, protein engineering, etc.) with the potential to multiply our technology's impact. Potential Impacts: This application could dramatically expand manufacturing capabilities on Earth and in space: In situ resource utilization. A far greater range of materials and products will be available from the limited palette offered by in situ resource extraction techniques; Reduced equipment and material upmass for off-Earth habitats. Ready- to-use highly specialized construction materials (radiation hardened, compressive/tensile, light or dense) from an extremely low starting mass; Structured biomaterial production. New ready-to-use macro, micro, and molecular manufacturing techniques for traditional biomaterials such as wood, bone and shell; New and novel biocomposite creation. The ability to create completely novel material composites from any base material that cells can be engineered to produce. Suggested Mission Contexts: ISS part manufacturing. A 'minimal working example' making a finished biomaterial part aboard the International Space Station; A long-term Mars habitat. 'Cradle-to-grave' use at a hypothetical Mars habitat, covering everything from tools to construction materials. Alternate Abstract: Imagine being able to print anything from tools and composite building materials to food and human tissues. Imagine being on Mars with the ability to replace any broken part, whether it's a part of your spacesuit, your habitat, or your own body. We propose a technique that would allow just that. By printing 3D arrays of cells engineered to secrete the necessary materials, the abundant in situ resources of atmosphere and regolith become organic, inorganic, or organic-inorganic composite materials. Such materials include novel, biologically derived materials not previously possible to fabricate.

habitats↗
NASA NTRSDate not supplied

The System Modeling and Analysis of Resiliency in STEReO (SMARt-STEReO)

Wildfire emergency response has remained rooted in relatively low-tech solutions for coordination between ground and aerial assets. These low-tech solutions are robust for the remote environments in which wildfires are usually fought, but limit strategic cross-organizational support and the ability to deploy and effectively utilize aerial assets. As aircraft become more advanced and new technology, including drones, become available to firefighters, a new, more modern method of asset coordination is needed. NASA is working on a project called ‘Scalable Traffic Management for Emergency Response Operations’ (STEReO) to integrate unmanned aerial systems (UAS)and UAS traffic management (UTM)into wildfire response. STEReO’s goals include simplifying the coordination of aerial assets, improving the existing UAS framework, and increasing the role of additional autonomous systems to reduce human risk and to increase system resilience. This paper describes the development of the ‘System Modeling and Analysis of Resiliency in STEReO’ (SMARt-STEReO) project, which aims to model wildfire response and to quantify the additional system resilience that STEReO technology provides firefighters. This paper verifies SMARt-STEReO and defines its scope; it includes experimental and statistical analysis of the impact that the addition of UAS has on both performance metrics and also on performance resiliency response to a given fault. SMARt-STEReO is a grid-based model of fire propagation that incorporates varying crew responses. Through the use of a Python package called ‘fmdtools’, the model easily allows for the addition of faults to the system. These faults allow analysts to investigate various response parameters. Factors including terrain, fuel type and wind speed can be modified to affect the fire propagation; additionally, the number of ground crews, engines, fixed wing aircraft, helicopters, and UAS can be changed to affect the crew response. The communication lines between actors mimic those used in real life situations. This paper explains the development of SMARt-STEReO including background research, verification and validation, and preliminary experimental analysis of system resilience to both a minor and major fault in systems with and without UAS.

Resiliency↗
NASA NTRS2024

Evaluation of 10‐m Wind Speed From ISD Meteorological Stations and the MERRA‐2 Reanalysis: Impacts on Dust Emission in the Arabian Peninsula

Mineral dust is one of the most important aerosols when studying the radiative balance and climate of the planet. There are different dust emission schemes utilized by the atmospheric modeling communities, many of which disagree on basic output quantities such as mass of dust emitted and distribution of mass among size bins. In this work, we examined mineral dust emission from a leading model scheme, the Goddard Chemistry Aerosol Radiation and Transport (GOCART), as utilized in the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) Reanalysis and compared it to dust emissions calculated using wind measurements from ground based weather stations located in the Arabian Peninsula that are included in the National Oceanic and Atmospheric Administration’s (NOAA) integrated surface database (ISD). An intercomparison of 10-m wind speed is shown for the Arabian Peninsula region, differences of the observed and modeled wind field are quantified, and impacts of differences on dust emissions are calculated. This analysis shows 10-m winds in the ISD were generally lower than MERRA-2 winds, which propagated to dust emissions errors. Our estimate of one of the most significant mass impacts in dust emission is 0.178 Tg/year/grid box with a percent change of over 200% to the recalculated dust emissions from MERRA-2. These differences in wind speed propagated to a difference in dust mass emitted by the use of a static source function which aids in scaling the mass emitted by the availability of dust in each grid. Additionally, the magnitude of these differences varies seasonally.

winds↗
NASA NTRS1985

Vector diffraction analysis of reflector antennas with mesh surfaces

Reflector antennas with mesh surfaces are used extensively in many satellite and ground antenna systems. A strip-aperture modeling of commonly used mesh surfaces is presented which provides considerable versatility in characterizing the mesh cells. The mesh transmission coefficients are constructed using a Floquet-modal expansion in conjuction with two dominant aperture modes. To account for the mesh local coordinates, the Eulerian angle transformation is invoked to obtain the total induced current on the curved reflector surface. General formulas are presented to show how the solid surface induced current is modified due to the transmission through the mesh. The effects of a variety of mesh configurations on both the co-polar and cross-polar patterns of reflector antennas are studied by numerically evaluating the vector diffraction integral using the Jacobi-Bessel expansion. For some special cases, a comparison is made with the results of the commonly used wire-grid formulation. Many of the numerical data are tailored to the dimensions of a conceptually designed mesh deployable offset reflector of the land mobile satellite system (LMSS).

Rahmat-Samii, Y.↗
NASA NTRS1989

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗
NASA NTRS2009

Spacecraft Internal Acoustic Environment Modeling

Acoustic modeling can be used to identify key noise sources, determine/analyze sub-allocated requirements, keep track of the accumulation of minor noise sources, and to predict vehicle noise levels at various stages in vehicle development, first with estimates of noise sources, later with experimental data. In FY09, the physical mockup developed in FY08, with interior geometric shape similar to Orion CM (Crew Module) IML (Interior Mode Line), was used to validate SEA (Statistical Energy Analysis) acoustic model development with realistic ventilation fan sources. The sound power levels of these sources were unknown a priori, as opposed to previous studies that RSS (Reference Sound Source) with known sound power level was used. The modeling results were evaluated based on comparisons to measurements of sound pressure levels over a wide frequency range, including the frequency range where SEA gives good results. Sound intensity measurement was performed over a rectangular-shaped grid system enclosing the ventilation fan source. Sound intensities were measured at the top, front, back, right, and left surfaces of the and system. Sound intensity at the bottom surface was not measured, but sound blocking material was placed tinder the bottom surface to reflect most of the incident sound energy back to the remaining measured surfaces. Integrating measured sound intensities over measured surfaces renders estimated sound power of the source. The reverberation time T6o of the mockup interior had been modified to match reverberation levels of ISS US Lab interior for speech frequency bands, i.e., 0.5k, 1k, 2k, 4 kHz, by attaching appropriately sized Thinsulate sound absorption material to the interior wall of the mockup. Sound absorption of Thinsulate was modeled in three methods: Sabine equation with measured mockup interior reverberation time T60, layup model based on past impedance tube testing, and layup model plus air absorption correction. The evaluation/validation was carried out by acquiring octave band microphone data simultaneously at ten fixed locations throughout the mockup. SPLs (Sound Pressure Levels) predicted by our SEA model match well with measurements for our CM mockup, with a more complicated shape. Additionally in FY09, background NC noise (Noise Criterion) simulation and MRT (Modified Rhyme Test) were developed and performed in the mockup to determine the maximum noise level in CM habitable volume for fair crew voice communications. Numerous demonstrations of simulated noise environment in the mockup and associated SIL (Speech Interference Level) via MRT were performed for various communities, including members from NASA and Orion prime-/sub-contractors. Also, a new HSIR (Human-Systems Integration Requirement) for limiting pre- and post-landing SIL was proposed.

Chu, SShao-sheng R.↗
NASA NTRS1977

Atmospheric model development in support of SEASAT. Volume 3: Prediction models

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric models are developed. The four models use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical. Conservation forms of the difference equations are integrated using either a fifteen or four minute time step on a 381 km or 127 km grid (at 60 deg N) for the 63 x 63 or 187 x 187 models respectively. A nonlinear pressure smoother is used along with momentum and temperature diffusion to help control computational noise. The horizontal boundary conditions are insulated slippery walls. Centered time differencing with time averaging of the pressure gradient force term in the momentum equations is used. Robert time filtering on of the temperature and moisture solutions is used for computational stability. The moisture and heat source/sink terms are modeled. Terms representing evaporation and large scale condensation, sensible heat exchange, parameterized cumulus convection and precipitation, and solar and terrestrial radiation are included. Initialization of the models is based on a pattern conservation technique to obtain objective analysis of the state parameter structure from the surface to 50 mb.

Wellck, R. E.↗
NASA NTRS2000

TRMM Surface Backscatter Measurements - Temporal and Spatial Analysis: Inter-Comparison with NSCAT

Global rainfall is the primary redistributor of earth's energy by the process of latent heat release. This forms the main driving force for the tropical circulation, which in turn impacts the global circulation .through transient events-such as El Nino. Hence, more precise and long-term time series of the rainfall and its variability is crucial to the understanding and prediction of the global climate and climate change. The Precipitation Radar (PR) abroad the US/Japan Tropical Rainfall Measuring Mission (TRMM) is the first radar ever launched into space that measures detailed vertical profiles of rain intensity over the tropics. One of the challenges in estimating rainfall from spaceborne radars is the presence of attenuation at frequencies, such as 14 GHz of the TRMM PR and future planned systems at this and higher frequencies. A common approach in current rainfall retrieval algorithms is to employed the path integrated attenuation (PIA) as a constraint to the retrieval, and hence overcome errors in the radar calibration or in the assumed rainfall parameters. PIA can either be derived from a radiometer or from the surface reference technique, in which a clear air radar measurement is compared with the measurement in the raining area. The current TRMM 2A21 PIA data product makes use of both a temporal and spatial clear air database for comparison to rainy measurements. In this paper we present results from analysis of TRMM surface backscatter cross-section (sigmaO) measurements from Nov 97-Feb99, and a comparison with sigmaO measurements obtained by the NASA Scatterometer (NSCAT) between Sept96-June97. Measurements for a given month from both instruments are compiled on a 1 deg. (lat.) x 1 deg. (lon.) x 1 hr. grid. This enables TRMM--NSCAT comparison and the investigation of seasonal and diurnal trends in both data sets. From preliminary analysis of TRMM sigmaO's we have decided not to treat the ocean as a single homogeneous region but to select a number of ocean sub-regions and individually analyze their trends. Likewise, and in a similar approach to previous studies of Seasat over-land data, we have selected a number of over-land regions for study. From said sigmaO maps and regional trend analysis we investigate possible sources of trends and variability. In addition, we study the effects of TRMM PR sensitivity through the PR "possible rain" class. Given NSCAT's inability to flag rain contaminated measurements we are able to gauge the impact of rain contamination on NSCAT monthly sigmaO maps, using TRMM measurements. The research described in this paper was carried out by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA), U.S.A.

Im, Eastwood↗
NASA NTRS2012

Development of a Robust and Efficient Parallel Solver for Unsteady Turbomachinery Flows

The traditional design and analysis practice for advanced propulsion systems relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. In the design of advanced propulsion systems, CFD plays a major role in defining the required performance over the entire flight regime, as well as in testing the sensitivity of the design to the different modes of operation. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment. The turbomachinery simulation capability presented here is being developed in a computational tool called Loci-STREAM [1]. It integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci [2] which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective is to be able to routinely simulate problems involving complex geometries requiring large unstructured grids and complex multidisciplinary physics. An immediate application of interest is simulation of unsteady flows in rocket turbopumps, particularly in cryogenic liquid rocket engines. The key components of the overall methodology presented in this paper are the following: (a) high fidelity unsteady simulation capability based on Detached Eddy Simulation (DES) in conjunction with second-order temporal discretization, (b) compliance with Geometric Conservation Law (GCL) in order to maintain conservative property on moving meshes for second-order time-stepping scheme, (c) a novel cloud-of-points interpolation method (based on a fast parallel kd-tree search algorithm) for interfaces between turbomachinery components in relative motion which is demonstrated to be highly scalable, and (d) demonstrated accuracy and parallel scalability on large grids (approx 250 million cells) in full turbomachinery geometries.

West, Jeff↗
NASA NTRS2024

Tail-Cone Thruster Fan Rotor Flutter Analysis

The aeroelastic stability of a tail cone thruster fan was analyzed for three structural modes and various nodal diameter patterns. Computational fluid dynamics modeling with blade vibrations was used in the analyses performed at design (peak-efficiency) and two near-stall conditions at the design rotational speed. Aerodynamic damping was used to assess flutter since this integrally bladed rotor design will have minimal structural damping. Aerodynamic damping was found to be low at the design condition. For the first mode, aerodynamic damping was seen to decrease towards stall, vanishing for a zero nodal diameter pattern. Further decrease in mass flow led to negative aerodynamic damping (flutter) for the first and second nodal diameters in the near stall condition. The flowfield at this lowest mass flow rate exhibited some additional unsteadiness at higher frequency. The sensitivity of results to time step and amplitudes of vibration was examined and found to be very satisfactory. The largest uncertainty in the results comes from the observed sensitivity of the aerodynamic damping to details of the variation of grid deformations in the blade passage mesh from the blade passage to the boundaries. For the parametric variations analyzed, the calculated aerodynamic damping is lower if the grid deformations near the blade surface are higher. Further work is needed to understand the cause of this sensitivity and find possible solutions to reduce the effects on flutter predictions.

Fan↗
NASA NTRS1986

A variational assimilation method for satellite and conventional data: Development of basic model for diagnosis of cyclone systems

A three-dimensional diagnostic model for the assimilation of satellite and conventional meteorological data is developed with the variational method of undetermined multipliers. Gridded fields of data from different type, quality, location, and measurement source are weighted according to measurement accuracy and merged using least squares criteria so that the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation are satisfied. The model is used to compare multivariate variational objective analyses with and without satellite data with initial analyses and the observations through criteria that were determined by the dynamical constraints, the observations, and pattern recognition. It is also shown that the diagnoses of local tendencies of the horizontal velocity components are in good comparison with the observed patterns and tendencies calculated with unadjusted data. In addition, it is found that the day-night difference in TOVS biases are statistically different (95% confidence) at most levels. Also developed is a hybrid nonlinear sigma vertical coordinate that eliminates hydrostatic truncation error in the middle and upper troposphere and reduces truncation error in the lower troposphere. Finally, it is found that the technique used to grid the initial data causes boundary effects to intrude into the interior of the analysis a distance equal to the average separation between observations.

Achtemeier, G. L.↗
NASA NTRS2006

Validation of a Pressure-Based Combustion Simulation Tool Using a Single Element Injector Test Problem

The traditional design and analysis practice for advanced propulsion systems, particularly chemical rocket engines, relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment by non-CFD specialists. A computational tool, called Loci-STREAM is being developed for this purpose. It is a pressure-based, Reynolds-averaged Navier-Stokes (RANS) solver for generalized unstructured grids, which is designed to handle all-speed flows (incompressible to hypersonic) and is particularly suitable for solving multi-species flow in fixed-frame combustion devices. Loci-STREAM integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective of the ongoing work is to develop a robust simulation capability for combustion problems in rocket engines. As an initial step towards validating this capability, a model problem is investigated in the present study which involves a gaseous oxygen/gaseous hydrogen (GO2/GH2) shear coaxial single element injector, for which experimental data are available. The sensitivity of the computed solutions to grid density, grid distribution, different turbulence models, and different near-wall treatments is investigated. A refined grid, which is clustered in the vicinity of the solid walls as well as the flame, is used to obtain a steady state solution which may be considered as the best solution attainable with the steady-state RANS methodology. From a design point of view, quick turnaround times are desirable; with this in mind, coarser grids are also employed and the resulting solutions are evaluated with respect to the fine grid solution.

Thakur, Siddarth↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
NASA NTRS1999

Parallel Grid Manipulations in Earth Science Calculations

The National Aeronautics and Space Administration (NASA) Data Assimilation Office (DAO) at the Goddard Space Flight Center is moving its data assimilation system to massively parallel computing platforms. This parallel implementation of GEOS DAS will be used in the DAO's normal activities, which include reanalysis of data, and operational support for flight missions. Key components of GEOS DAS, including the gridpoint-based general circulation model and a data analysis system, are currently being parallelized. The parallelization of GEOS DAS is also one of the HPCC Grand Challenge Projects. The GEOS-DAS software employs several distinct grids. Some examples are: an observation grid- an unstructured grid of points at which observed or measured physical quantities from instruments or satellites are associated- a highly-structured latitude-longitude grid of points spanning the earth at given latitude-longitude coordinates at which prognostic quantities are determined, and a computational lat-lon grid in which the pole has been moved to a different location to avoid computational instabilities. Each of these grids has a different structure and number of constituent points. In spite of that, there are numerous interactions between the grids, e.g., values on one grid must be interpolated to another, or, in other cases, grids need to be redistributed on the underlying parallel platform. The DAO has designed a parallel integrated library for grid manipulations (PILGRIM) to support the needed grid interactions with maximum efficiency. It offers a flexible interface to generate new grids, define transformations between grids and apply them. Basic communication is currently MPI, however the interfaces defined here could conceivably be implemented with other message-passing libraries, e.g., Cray SHMEM, or with shared-memory constructs. The library is written in Fortran 90. First performance results indicate that even difficult problems, such as above-mentioned pole rotation- a sparse interpolation with little data locality between the physical lat-lon grid and a pole rotated computational grid- can be solved efficiently and at the GFlop/s rates needed to solve tomorrow's high resolution earth science models. In the subsequent presentation we will discuss the design and implementation of PILGRIM as well as a number of the problems it is required to solve. Some conclusions will be drawn about the potential performance of the overall earth science models on the supercomputer platforms foreseen for these problems.

Sawyer, W.↗
NASA NTRS2013

Flight Performance Feasibility Studies for the Max Launch Abort System

In 2007, the NASA Engineering and Safety Center (NESC) initiated the Max Launch Abort System Project to explore crew escape system concepts designed to be fully encapsulated within an aerodynamic fairing and smoothly integrated onto a launch vehicle. One objective of this design was to develop a more compact launch escape vehicle that eliminated the need for an escape tower, as was used in the Mercury and Apollo escape systems and what is planned for the Orion Multi-Purpose Crew Vehicle (MPCV). The benefits for the launch vehicle of eliminating a tower from the escape vehicle design include lower structural weights, reduced bending moments during atmospheric flight, and a decrease in induced aero-acoustic loads. This paper discusses the development of encapsulated, towerless launch escape vehicle concepts, especially as it pertains to the flight performance and systems analysis trade studies conducted to establish mission feasibility and assess system-level performance. Two different towerless escape vehicle designs are discussed in depth: one with allpropulsive control using liquid attitude control thrusters, and a second employing deployable aft swept grid fins to provide passive stability during coast. Simulation results are presented for a range of nominal and off-nominal escape conditions.

Tarabini, Paul V.↗