SEARCH · Engineering Papers

Results for “Grid Analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1995

TADS: A CFD-based turbomachinery and analysis design system with GUI. Volume 1: Method and results

The primary objective of this study was the development of a computational fluid dynamics (CFD) based turbomachinery airfoil analysis and design system, controlled by a graphical user interface (GUI). The computer codes resulting from this effort are referred to as the Turbomachinery Analysis and Design System (TADS). This document describes the theoretical basis and analytical results from the TADS system. TADS couples a throughflow solver (ADPAC) with a quasi-3D blade-to-blade solver (RVCQ3D) in an interactive package. Throughflow analysis capability was developed in ADPAC through the addition of blade force and blockage terms to the governing equations. A GUI was developed to simplify user input and automate the many tasks required to perform turbomachinery analysis and design. The coupling of various programs was done in a way that alternative solvers or grid generators could be easily incorporated into the TADS framework. Results of aerodynamic calculations using the TADS system are presented for a highly loaded fan, a compressor stator, a low-speed turbine blade, and a transonic turbine vane.

Topp, D. A.↗

NASA NTRS2003

Comparison of Flow-Dependent and Static Error Correlation Models in the DAO Ozone Data Assimilation System

In a data assimilation system the forecast error covariance matrix governs the way in which the data information is spread throughout the model grid. Implementation of a correct method of assigning covariances is expected to have an impact on the analysis results. The simplest models assume that correlations are constant in time and isotropic or nearly isotropic. In such models the analysis depends on the dynamics only through assumed error standard deviations. In applications to atmospheric tracer data assimilation this may lead to inaccuracies, especially in regions with strong wind shears or high gradient of potential vorticity, as well as in areas where no data are available. In order to overcome this problem we have developed a flow-dependent covariance model that is based on short term evolution of error correlations. The presentation compares performance of a static and a flow-dependent model applied to a global three- dimensional ozone data assimilation system developed at NASA s Data Assimilation Office. We will present some results of validation against WMO balloon-borne sondes and the Polar Ozone and Aerosol Measurement (POAM) III instrument. Experiments show that allowing forecast error correlations to evolve with the flow results in positive impact on assimilated ozone within the regions where data were not assimilated, particularly at high latitudes in both hemispheres and in the troposphere. We will also discuss statistical characteristics of both models; in particular we will argue that including evolution of error correlations leads to stronger internal consistency of a data assimilation ,

Wargan, K.↗

NASA NTRS2005

Seasonal-scale Observational Data Analysis and Atmospheric Phenomenology for the Cold Land Processes Experiment

The Cold Land Processes Experiment (CLPX) experiment emphasized the development of a strong synergism between process-oriented understanding, land surface models and microwave remote sensing. Our work sought to investigate which topographically- generated atmospheric phenomena are most relevant to the CLPX MSA's for the purpose of evaluating their climatic importance to net local moisture fluxes and snow transport through the use of high-resolution data assimilation/atmospheric numerical modeling techniques. Our task was to create three long-term, scientific quality atmospheric datasets for quantitative analysis (for all CLPX researchers) and provide a summary of the meteorologically-relevant phenomena of the three MSAs (see Figure) over northern Colorado. Our efforts required the ingest of a variety of CLPX datasets and the execution an atmospheric and land surface data assimilation system based on the Navier-Stokes equations (the Local Analysis and Prediction System, LAPS, and an atmospheric numerical weather prediction model, as required) at topographically- relevant grid spacing (approx. 500 m). The resulting dataset will be analyzed by the CLPX community as a part of their larger research goals to determine the relative influence of various atmospheric phenomena on processes relevant to CLPX scientific goals.

Poulos, Gregory S.↗

NASA NTRS2021

Intergenerational Inequities in Exposure to Climate Extremes: Young Generations Are Severely Threatened By Climate Change

Under continued global warming, extreme events such as heat waves will continue to rise in frequency, intensity, duration, and spatial extent over the next decades. Younger generations are therefore expected to face more such events across their lifetimes compared with older generations. This raises important issues of solidarity and fairness across generations that have fueled a surge of climate protests led by young people in recent years and that underpin issues of intergenerational equity raised in recent climate litigation. However, the standard scientific paradigm is to assess climate change in discrete time windows or at discrete levels of warming, a “period” approach that inhibits quantification of how much more extreme events a particular generation will experience over its lifetime compared with another. By developing a “cohort” perspective to quantify changes in lifetime exposure to climate extremes and compare across generations (see the first figure), we estimate that children born in 2020 will experience a two- to sevenfold increase in extreme events, particularly heat waves, compared with people born in 1960, under current climate policy pledges. Our results highlight a severe threat to the safety of young generations and call for drastic emission reductions to safeguard their future. Meteorological extremes, hazards, or climate change impacts are mostly studied as they evolve over time under varying emission scenarios and socioeconomic pathways. For example, applying a heat wave indicator (see table S1) to four bias-adjusted global climate models indicates that the land area annually affected by such heat waves will increase from ~15% around 2020 to ~22% by 2100 under a scenario compatible with limiting global warming to 1.5°C, and to ~46% under a scenario in line with current emission reduction pledges (see the first figure). Recent studies extended this approach, studying aspects of climate change as a function of global mean temperature (GMT) increments, highlighting the scenario-independence of several extreme event indicators but remaining, in essence, a comparison of time windows. By contrast, we performed a birth cohort analysis by combining a collection of multimodel extreme event projections with country-scale life expectancy information, gridded population data, and future global temperature trajectories from the Intergovernmental Panel on Climate Change (IPCC) Special Report on Global Warming of 1.5°C (see supplementary materials). By integrating the exposure of an average person in a country or region to extreme events across their lifetime, we encapsulate spatiotemporal changes in climate hazards, population density, cohort size, and life expectancy (see the first figure).

global warming↗

NASA NTRSDate not supplied

Progress on Transonic Flutter and Shock Buffet Computationsin Support of the Third Aeroelastic Prediction Workshop

This paper reports on the progress of the NASA Langley team contributions to the third Aeroelastic Pre-diction Workshop’s (AePW-3) High Angle Working Group (HAWG). The primary objectives of HAWG is to predict the fluter dynamic pressure of the NASA Benchmark Supercritical Wing (BSCW) configuration at Mach 0.8 and 5◦angle of attack. The secondary objective is to determine if a shock-buffet onset is present at or near that flow condition. The computational results are obtained using FUN3D, an unstructured grid Reynolds-averaged Navier-Stokes solver developed at the NASA Langley Research Center. The preliminary analysis results show a computationally-obtained flutter dynamic pressure of approximately 120 psf. Initial results describing unforced BSCW unsteady flow environment at flutter condition are also presented.

Pawel Chwalowski↗

NASA NTRS2002

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗

NASA NTRS2018

Novel Applications of Focused Ion Beam Technique for Planetary Sample Analyses

We are using innovative FIB techniques to prepare samples of planetary materials for different types of coordinated analyses using ion microprobes, synchrotron beamlines, and specialized transmission electron microscopy (TEM) techniques. In these cases, the FIB sample preparation is the critical step in enabling these specialized analyses. We discuss several examples below utilizing the FEI Quanta3D instrument at the NASA Johnson Space Center. Trace element analyses utilizing synchrotron x-ray fluorescence. The trace element content of mineral grains in comet dust provides important clues on their formation and processing in the early solar system. We preformed coordinated analyses of a comet dust particle that had been prepared using ultramicrotomy for TEM analysis. Following the TEM analyses, we extracted a 70 nm thick section from a region of the carbon (C) film of the TEM grid, for additional analyses. A carbon ring ~2-3 μm thick was deposited on top of the C film using the FIB. The C film on the outer rim of the ring was milled away using various patterns to uniformly release the stresses on the film, preventing rupture and collapse, and was attached to the micromanipulator needle. We then isolated the ring completely and transferred the section to a silicon sample holder for analysis using the HXN (hard X-ray nanoprobe) beamline at NSLSII at Brookhaven National Lab. Coordinated Analyses of Presolar Grains. Rare sub-m presolar grains that originate in evolved stars and supernovae, occur in primitive astromaterials and are identified by their exotic isotopic compositions. Coordinated analyses of these grains using NanoSIMS, TEM, and other techniques on the same grain is enabled by innovative FIB sample preparation. In order to obtain subsequent isotopic analyses of Mg and Fe, contributions from surrounding grains were minimized. We precisely deposited a protective cap of Pt on top of the grain to preserve the grain of interest and then milled away about 5 μm diameter of the surrounding material. Following the isotopic analyses, the spindle was extracted and thinned to electron transparency for TEM microstructural analyses. In situ heating TEM experiments on lunar samples. We extracted a FIB thin section from Apollo 17 lunar rock 76015. To avoid ion-beam damage, e-beam deposition was used to deposit the first 500 nm of the C strap, followed by ion beam-assisted deposition of ~3 μm carbon. We performed an ex situ lift-out of the section and placed the section on one of the elements of a microelectromechanical systems (MEMS) - specialized heating substrate and attached the section to the substrate by depositing small C straps with the FIB. The heating chips utilize silicon nitride windows to support the samples and provide uniform heating while enabling TEM imaging. The heating chip was loaded into a Hitachi “Blaze” heating holder and analyzed using a Hitachi HF5000 at the University of Arizona.

Rahman, Z.↗

NASA NTRS1992

Integrating aerodynamic surface modeling for computational fluid dynamics with computer aided structural analysis, design, and manufacturing

This presentation will discuss the development of a NASA Geometry Exchange Specification for transferring aerodynamic surface geometry between LeRC systems and grid generation software used for computational fluid dynamics research. The proposed specification is based on a subset of the Initial Graphics Exchange Specification (IGES). The presentation will include discussion of how the NASA-IGES standard will accommodate improved computer aided design inspection methods and reverse engineering techniques currently being developed. The presentation is in viewgraph format.

Thorp, Scott A.↗

NASA NTRS2003

An Aerodynamic Analysis of a Spinning Missile with Dithering Canards

A generic spinning missile with dithering canards is used to demonstrate the utility of an overset structured grid approach for simulating the aerodynamics of rolling airframe missile systems. The approach is used to generate a modest aerodynamic database for the generic missile. The database is populated with solutions to the Euler and Navier-Stokes equations. It is used to evaluate grid resolution requirements for accurate prediction of instantaneous missile loads and the relative aerodynamic significance of angle-of-attack, canard pitching sequence, viscous effects, and roll-rate effects. A novel analytical method for inter- and extrapolation of database results is also given.

Meakin, Robert L.↗

NASA NTRS1998

Parallel 3D Multi-Stage Simulation of a Turbofan Engine

A 3D multistage simulation of each component of a modern GE Turbofan engine has been made. An axisymmetric view of this engine is presented in the document. This includes a fan, booster rig, high pressure compressor rig, high pressure turbine rig and a low pressure turbine rig. In the near future, all components will be run in a single calculation for a solution of 49 blade rows. The simulation exploits the use of parallel computations by using two levels of parallelism. Each blade row is run in parallel and each blade row grid is decomposed into several domains and run in parallel. 20 processors are used for the 4 blade row analysis. The average passage approach developed by John Adamczyk at NASA Lewis Research Center has been further developed and parallelized. This is APNASA Version A. It is a Navier-Stokes solver using a 4-stage explicit Runge-Kutta time marching scheme with variable time steps and residual smoothing for convergence acceleration. It has an implicit K-E turbulence model which uses an ADI solver to factor the matrix. Between 50 and 100 explicit time steps are solved before a blade row body force is calculated and exchanged with the other blade rows. This outer iteration has been coined a "flip." Efforts have been made to make the solver linearly scaleable with the number of blade rows. Enough flips are run (between 50 and 200) so the solution in the entire machine is not changing. The K-E equations are generally solved every other explicit time step. One of the key requirements in the development of the parallel code was to make the parallel solution exactly (bit for bit) match the serial solution. This has helped isolate many small parallel bugs and guarantee the parallelization was done correctly. The domain decomposition is done only in the axial direction since the number of points axially is much larger than the other two directions. This code uses MPI for message passing. The parallel speed up of the solver portion (no 1/0 or body force calculation) for a grid which has 227 points axially.

Turner, Mark G.↗

NASA NTRS2013

Uncertainty of Passive Imager Cloud Optical Property Retrievals to Instrument Radiometry and Model Assumptions: Examples from MODIS

The optical and microphysical structure of clouds is of fundamental importance for understanding a variety of cloud radiation and precipitation processes. With the advent of MODIS on the NASA EOS Terra and Aqua platforms, simultaneous global-daily 1 km retrievals of cloud optical thickness (COT) and effective particle radius (CER) are provided, as well as the derived water path (WP). The cloud product (MOD06/MYD06 for MODIS Terra and Aqua, respectively) provides separate retrieval datasets for various two-channel retrievals, typically a VISNIR channel paired with a 1.6, 2.1, and 3.7 m spectral channel. The MOD06 forward model is derived from on a homogeneous plane-parallel cloud. In Collection 5 processing (completed in 2007 with a modified Collection 5.1 completed in 2010), pixel-level retrieval uncertainties were calculated for the following non-3-D error sources: radiometry, surface spectral albedo, and atmospheric corrections associated with model analysis uncertainties (water vapor only). The latter error source includes error correlation across the retrieval spectral channels. Estimates of uncertainty in 1 aggregated (Level-3) means were also provided assuming unity correlation between error sources for all pixels in a grid for a single day, and zero correlation of error sources from one day to the next. I n Collection 6 (expected to begin in late summer 2013) we expanded the uncertainty analysis to include: (a) scene-dependent calibration uncertainty that depends on new band and detector-specific Level 1B uncertainties, (b) new model error sources derived from the look-up tables which includes sensitivities associated with wind direction over the ocean and uncertainties in liquid water and ice effective variance, (c) thermal emission uncertainties in the 3.7 m band associated with cloud and surface temperatures that are needed to extract reflected solar radiation from the total radiance signal, (d) uncertainty in the solar spectral irradiance at 3.7 m, and (e) addition of stratospheric ozone uncertainty in visible atmospheric corrections. A summary of the approach and example Collection 6 results will be shown.

Cloud Properties Observations↗

NASA NTRS2017

Nearfield Summary and Statistical Analysis of the Second AIAA Sonic Boom Prediction Workshop

A summary is provided for the Second AIAA Sonic Boom Workshop held 8-9 January 2017 in conjunction with AIAA SciTech 2017. The workshop used three required models of increasing complexity: an axisymmetric body, a wing body, and a complete configuration with flow-through nacelle. An optional complete configuration with propulsion boundary conditions is also provided. These models are designed with similar nearfield signatures to isolate geometry and shock/expansion interaction effects. Eleven international participant groups submitted nearfield signatures with forces, pitching moment, and iterative convergence norms. Statistics and grid convergence of these nearfield signatures are presented. These submissions are propagated to the ground, and noise levels are computed. This allows the grid convergence and the statistical distribution of a noise level to be computed. While progress is documented since the first workshop, improvement to the analysis methods for a possible subsequent workshop are provided. The complete configuration with flow-through nacelle showed the most dramatic improvement between the two workshops. The current workshop cases are more relevant to vehicles with lower loudness and have the potential for lower annoyance than the first workshop cases. The models for this workshop with quieter ground noise levels than the first workshop exposed weaknesses in analysis, particularly in convective discretization.

Park, Michael A.↗

NASA NTRS2018

Nearfield Summary and Statistical Analysis of the Second AIAA Sonic Boom PredictionWorkshop

A summary is provided for the Second AIAA Sonic Boom Workshop held 8–9 January 2017 in conjunction with AIAA SciTech 2017. The workshop used three required models of increasing complexity: an axisymmetric body, a wing body, and a complete configuration with flow-through nacelle. An optional complete configuration with propulsion boundary conditions is also provided. These models are designed with similar nearfield signatures to isolate geometry and shock/expansion interaction e ects. Eleven international participant groups submitted nearfield signatures with forces, pitching moment, and iterative convergence norms. Statistics and grid convergence of these nearfield signatures are presented. These submissions are propagated to the ground, and noise levels are computed. This allows the grid convergence and the statistical distribution of a noise level to be computed. While progress is documented since the first workshop, improvement to the analysis methods for a possible subsequent workshop are provided. The complete configuration with flow-through nacelle showed the most dramatic improvement between the two workshops. The current workshop cases are more relevant to vehicles with lower loudness and have the potential for lower annoyance than the first workshop cases. The models for this workshop with quieter ground noise levels than the first workshop exposed weaknesses in analysis, particularly in convective discretization.

Park, Michael A.↗

NASA NTRS2000

A Tool for Automatic Data Distribution for CFD Applications on Structured Grids

Development of HPF versions of NPB and ARC3D has shown that HPF provides an efficient, concise way to express parallelism and to organize data traffic. The use of HPF, as noted in the papers, requires an intimate knowledge of the applications and a detailed analysis of data affinity, data movement, and data granularity. To simplify and accelerate the task of developing HPF versions of existing CFD applications we have designed and implemented ADAPT (Automatic Data Alignment and Placement Tool). ADAPT analyzes a CFD application working on a single structured grid and generates HPF TEMPLATE, (RE)DISTRIBUTION, ALIGNMENT, and INDEPENDENT directives. The directives can be generated on the nest level, subroutine level, application level, or on the application interface level. ADAPT annotates an existing CFD FORTRAN application, performing computations on single or multiple grids. On each grid the application is considered as a sequence of operators, each applied to a set of variables defined in a particular grid domain. ADAPT automatically detects implicit operators (i.e., having data dependences) and explicit operators (without data dependences). For parallelization of an explicit operator ADAPT creates a template for the operator domain, aligns arrays used in the operator with the template, distributes the template, and declares the loops over the distributed dimensions as INDEPENDENT. For parallelization of an implicit operator, the distribution of the operator's domain should be consistent with the operator's dependences. Any dependence between sections distributed on different processors would preclude parallelization if the compiler does not have an ability to pipeline computations. If a data distribution is "orthogonal" to the dependences of an implicit operator, then the loop which implements the operator can be declared as INDEPENDENT. ADAPT starts with an analysis of array index expressions of the loop nests. For each pair of arrays referenced in an assignment statement, it generates an arc in the alignment graph and annotates it with an affinity relation. The template, alignment, and distribution directives for a particular loop nest are then derived from a transitive closure of the affinity relation. A compromise of data distributions in different nests and subroutines is achieved by merging annotated alignment graphs for adjacent nests/stibroutine calls in the nest/call graph of the application in the process called distribution lifting. ADAPT has been implemented as a C++ program running in conjunction with a parallelization tool called CAPTools. ADAPT uses the parse tree, interprocedural analysis and application database generated by CAPTools. It also uses the Directed Graph class, initially implemented in p2d2 (parallel debugger oi distributed programs), and some other classes supporting symbolic computations. ADAPT uses data distribution techniques described. ADAPT was tested with ARC3D and the FT benchmark and has demonstrated a code performance within a factor of 1.5 of handwritten versions.

Frumkin, Michael↗

NASA NTRS2006

Multiphysics Analysis of a Solid-Core Nuclear Thermal Engine Thrust Chamber

The objective of this effort is to develop an efficient and accurate thermo-fluid computational methodology to predict environments for a hypothetical solid-core, nuclear thermal engine thrust chamber. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics methodology. Formulations for heat transfer in solids and porous media were implemented and anchored. A two-pronged approach was employed in this effort: A detailed thermo-fluid analysis on a multi-channel flow element for mid-section corrosion investigation; and a global modeling of the thrust chamber to understand the effect of hydrogen dissociation and recombination on heat transfer and thrust performance. The formulations and preliminary results on both aspects are presented.

Wang, Ten-See↗

NASA NTRS1981

The reduction, verification and interpretation of MAGSAT magnetic data over Canada

Analysis of Investigator B tapes was extended to use a statistical screening of individual data points in each grid cell over the map area. The mean value of data from each cell was computed and all data values more than 4nT from the mean were rejected. The revised mean and standard deviation was computed. The average standard deviation for all map grid cells was 1.4nT for the scalar anomaly map and about 2nT for the vector component anomaly maps. At this stage, the grid cell is larger and spatial resolution lower for the vector maps than for the scalar maps, because the vector components are considerably more disturbed and more data points were rejected. An approach for altitude adjustment is in progress so that corrections can be made and data reduction procede. To evaluate the use of scalar and vector MAGSAT data for the production of regional magnetic charts, charts of the north, east, vertical, and scalar field derived from MAGSAT and from Canadian aeromagnetic data were compared. It appears possible to retrieve crustal signals from the horizontal vector components despite external disturbance fields at high latitude.

Coles, R. L.↗

NASA NTRS1986

Comparison between the South Pacific Convergence Zone and South Atlantic Convergence Zone in January 1979

Using one of the FGGE data sets, the eddy potential and eddy kinetic energy budgets during the period of January 10-27, 1979 were derived for the regions containing the South Pacific Convergence Zone (SPCZ) and the South Atlantic Convergence Zone. A complete set of energy equations developed by Brennan and Vincent (1980) was partitioned into zonal and eddy components; the eddy forms of energy were derived from a grid point's departure from the zonal average over a 5-deg latitude belt. The results of the analysis indicate that the eddy kinetic energy content (KE) is greater than the eddy available potential energy content (AE) throughout the period. Both AE and KE show large losses beginning when the SPCZ enters a decaying state: AE shows a gradual but slow decrease until after January 22, when it decreases rapidly; KE oscillates, reaching a maximum on January 22 before a rapid decrease.

Vincent, Dayton G.↗

NASA NTRS1994

Multigrid one shot methods for optimal control problems: Infinite dimensional control

The multigrid one shot method for optimal control problems, governed by elliptic systems, is introduced for the infinite dimensional control space. ln this case, the control variable is a function whose discrete representation involves_an increasing number of variables with grid refinement. The minimization algorithm uses Lagrange multipliers to calculate sensitivity gradients. A preconditioned gradient descent algorithm is accelerated by a set of coarse grids. It optimizes for different scales in the representation of the control variable on different discretization levels. An analysis which reduces the problem to the boundary is introduced. It is used to approximate the two level asymptotic convergence rate, to determine the amplitude of the minimization steps, and the choice of a high pass filter to be used when necessary. The effectiveness of the method is demonstrated on a series of test problems. The new method enables the solutions of optimal control problems at the same cost of solving the corresponding analysis problems just a few times.

Arian, Eyal↗