Engineering Papers⌕ Search

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.

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

At least 199 records · Page 11

Atmospheric water vapor transport: Estimation of continental precipitation recycling and parameterization of a simple climate model

The advective transport of atmospheric water vapor and its role in global hydrology and the water balance of continental regions are discussed and explored. The data set consists of ten years of global wind and humidity observations interpolated onto a regular grid by objective analysis. Atmospheric water vapor fluxes across the boundaries of selected continental regions are displayed graphically. The water vapor flux data are used to investigate the sources of continental precipitation. The total amount of water that precipitates on large continental regions is supplied by two mechanisms: (1) advection from surrounding areas external to the region; and (2) evaporation and transpiration from the land surface recycling of precipitation over the continental area. The degree to which regional precipitation is supplied by recycled moisture is a potentially significant climate feedback mechanism and land surface-atmosphere interaction, which may contribute to the persistence and intensification of droughts. A simplified model of the atmospheric moisture over continents and simultaneous estimates of regional precipitation are employed to estimate, for several large continental regions, the fraction of precipitation that is locally derived. In a separate, but related, study estimates of ocean to land water vapor transport are used to parameterize an existing simple climate model, containing both land and ocean surfaces, that is intended to mimic the dynamics of continental climates.

Brubaker, Kaye L.↗

Design of Optimum Ducts Using an Efficient 3-D Viscous Computational Flow Analysis

Design of fluid dynamically efficient ducts is addressed through the combination of an optimization analysis with a three-dimensional viscous fluid dynamic analysis code. For efficiency, a parabolic fluid dynamic analysis was used. Since each function evaluation in an optimization analysis is a full three-dimensional viscous flow analysis requiring 200,000 grid points, it is important to use both an efficient fluid dynamic analysis and an efficient optimization technique. Three optimization techniques are evaluated on a series of test functions. The Quasi-Newton (BFGS, eta = .9) technique was selected as the preferred technique. A series of basic duct design problems are performed. On a two-parameter optimization problem, the BFGS technique is demonstrated to require half as many function evaluations as a steepest descent technique.

Madabhushi, Ravi K.↗

Electric propulsion system technology

The work performed on the Ion Propulsion System Technology Task in FY90 is described. The objectives of this work fall under two broad categories. The first of these deals with issues associated with the application of xenon ion thrusters for primary propulsion of planetary spacecraft, and the second with the investigation of technologies which will facilitate the development of larger, higher power ion thrusters to support more advanced mission applications. Most of the effort was devoted to investigation of the critical issues associated with the use of ion thrusters for planetary spacecraft. These issues may be succinctly referred to as life time, system integration, and throttling. Chief among these is the engine life time. If the engines do not have sufficient life to perform the missions of interest, then the other issues become unimportant. Ion engine life time was investigated through two experimental programs: an investigation into the reduction of ion engine internal sputter erosion through the addition of small quantities of nitrogen, and a long duration cathode life test. In addition, a literature review and analysis of accelerator grid erosion were performed. The nitrogen addition tests indicated that the addition of between 0.5 and 1.0 percent of nitrogen by mass to the xenon propellant results in a reduction in the sputter erosion of discharge chamber components by a factor of between 20 and 50, with negligible reduction in thruster performance. The long duration test of a 6.35-mm dia. xenon hollow cathode is still in progress, and has accumulated more than 4,000 hours of operation at an emission current of 25 A at the time of this writing. One of the major system integration issues concerns possible interactions of the ion thruster produced charge exchange plasma with the spacecraft. A computer model originally developed to describe the behavior of mercury ion thruster charge exchange plasmas was resurrected and modified for xenon propellant. This model enables one to calculate the flow direction and local density of the charge exchange plasma, and indicates the degree to which this plasma can flow upstream of the thruster exhaust plane. A continuing effort to investigate the most desirable throttling technique for noble gas ion thrusters concentrated this year on experimentally determining the fixed flow rate throttling range of a 30-cm dia. thruster with a two-grid accelerator system. These experiments demonstrated a throttling capability which covers a 2.8 to 1 variation in input power. This throttling range is 55 percent greater than expected, and is due to better accelerator system performance at low net-to-total voltage ratios than indicated in the literature. To facilitate the development of large, higher power ion thrusters several brief studies were performed. These include the development of a technique which simulates ion thruster operation without beam extraction, the development of an optical technique to measure ion thruster grid distortion due to thermal expansion, tests of a capacitance measurement technique to quantify the accelerator system grid separation, and the development of a segmented thruster geometry which enables near term development of ion thrusters at power levels greater than 100 kW. Finally, a paper detailing the benefits of electric propulsion for the Space Exploration Initiative was written.

Brophy, John R.↗

A User's Manual for ROTTILT Solver: Tiltrotor Fountain Flow Field Prediction

A CFD solver has been developed to provide the time averaged details of the fountain flow typical for tiltrotor aircraft in hover. This Navier-Stokes solver, designated as ROTTILT, assumes the 3-D fountain flowfield to be steady and incompressible. The theoretical background is described in this manual. In order to enable the rotor trim solution in the presence of tiltrotor aircraft components such as wing, nacelle, and fuselage, the solver is coupled with a set of trim routines which are highly efficient in CPU and suitable for CFD analysis. The Cartesian grid technique utilized provides the user with a unique capability for insertion or elimination of any components of the bodies considered for a given tiltrotor aircraft configuration. The flowfield associated with either a semi or full-span configuration can be computed through user options in the ROTTILT input file. Full details associated with the numerical solution implemented in ROTTILT and assumptions are presented. A description of input surface mesh topology is provided in the appendices along with a listing of all preprocessor programs. Input variable definitions and default values are provided for the V22 aircraft. Limited predicted results using the coupled ROTTILT/WOPWOP program for the V22 in hover are made and compared with measurement. To visualize the V22 aircraft and predictions, a preprocessor graphics program GNU-PLOT3D was used. This program is described and example graphic results presented.

Tadghighi, Hormoz↗

Microshutter Arrays for the NGST NIRSpec

We are developing a programmable multiobject field selector for the NGST NIRSpec. This device is a microshutter array, a close-packed 2000 x 2000 array of 100 micrometer x 100 micrometer shutters fabricated with Micro-ElectroMechanical Systems (MEMS) technologies. The shutters are opened using a magnet scanned over the array. The shutters are held open by electrostatic potentials on adjacent vertical electrodes. Modeling of the magnetic actuation of the microshutters allowed optimization of the magnetic field configuration, the distance between magnet and shutters, and the force and torque produced by the magnet. The results of laboratory tests are consistent with our modeling. We have demonstrated actuation, latching and addressing. We are able to rotate shutters out of the plane up to 90 degree and more, to latch them electrostatically to the walls and release them selectively. During the last six month all critical elements of the array have been demonstrated. Fabrication processes were developed to deposit vertical electrodes on support grid side walls and insulated light-shields that block the gaps between shutters and the support grid. Physical optics analysis was performed to assess the expected optical performance of the microshutters. A preliminary Fourier optics analysis on the microshutter array performance has been carried out. More detailed analysis of wave propagation through 3D structures is a next step. An optical test station for verifying optical parameters of the shutter array at room temperature has been developed. It allows to measure parameters, such as transmission, scattering, contrast ratio, and diffraction of the array structure. Initial data have been acquired and are being analyzed and compared to theoretical predictions.

Moseley, Harvey↗

RTF Analysis-Propulsion/Flowliner

The objective of the current task is to provide a computational framework for design and analysis of the entire fuel supply system of a liquid rocket engine, including high fidelity unsteady turbopump flow analysis. Time-accurate results obtained from shuttle fuel flowliner analysis using 66 Million grid points with 262 overlapped zones will be shown. We will present analysis results along with performance data of the simulation runs on Supercomputers such as Columbia.

Kiris, Centin↗

The Global Precipitation Climatology Project (GPCP): Results, Status and Future

The Global Precipitation Climatology Project (GPCP) is one of a number of long-term, satellite-based, global analyses routinely produced under the auspices of the World Climate Research Program (WCRP) and its Global Energy and Watercycle EXperiment (GEWEX) program. The research quality analyses are produced a few months after real-time through the efforts of scientists at various national agencies and universities in the U.S., Europe and Japan. The primary product is a monthly analysis of surface precipitation that is globally complete and spans the period 1979-present. There are also pentad analyses for the same period and a daily analysis for the 1997-present period. Although generated with somewhat different data sets and analysis schemes, the pentad and daily data sets are forced to agree with the primary monthly analysis on a grid box by grid box basis. The primary input data sets are from low-orbit passive microwave observations, geostationary infrared observations and surface raingauge information. Examples of research with the data sets are discussed, focusing on tropical (25N-25s) rainfall variations and possible long-term changes in the 28-year (1979-2006) monthly dataset. Techniques are used to discriminate among the variations due to ENSO, volcanic events and possible long-term changes for rainfall over both land and ocean. The impact of the two major volcanic eruptions over the past 25 years is estimated to be about a 5% maximum reduction in tropical rainfall during each event. Although the global change of precipitation in the data set is near zero, a small upward linear change over tropical ocean (0.06 mm/day/l0yr) and a slight downward linear change over tropical land (-0.03 mm/day/l0yr) are examined to understand the impact of the inhomogeneity in the data record and the length of the data set. These positive changes correspond to about a 5% increase (ocean) and 3% increase (ocean plus land) during this time period. Relations between variations in surface temperature and precipitation are analyzed on seasonal to inter-decadal time scales. A new, version 3 of GPCP is being planned to incorporate new satellite information (e.g., TRMM) and provide higher spatial and temporal resolution for at least part of the data record. The goals and plans for that GPCP re-processing will be outlined.

Adler, Robert F.↗

A Computational Investigation of Gear Windage

A CFD method has been developed for application to gear windage aerodynamics. The goals of this research are to develop and validate numerical and modeling approaches for these systems, to develop physical understanding of the aerodynamics of gear windage loss, including the physics of loss mitigation strategies, and to propose and evaluate new approaches for minimizing loss. Absolute and relative frame CFD simulation, overset gridding, multiphase flow analysis, and sub-layer resolved turbulence modeling were brought to bear in achieving these goals. Several spur gear geometries were studied for which experimental data are available. Various shrouding configurations and free-spinning (no shroud) cases were studied. Comparisons are made with experimental data from the open literature, and data recently obtained in the NASA Glenn Research Center Gear Windage Test Facility. The results show good agreement with experiment. Interrogation of the validative and exploratory CFD results have led, for the first time, to a detailed understanding of the physical mechanisms of gear windage loss, and have led to newly proposed mitigation strategies whose effectiveness is computationally explored.

Hill, Matthew J.↗

Air Traffic Management TestBed: Weather Visualization using Map Tiles

Running realistic simulations and developing algorithms for decision support tools usually require weather data access. This document describes the weather data service capability implemented in the Air Traffic Management (ATM) TestBed. The capability provides ways for querying weather data from the weather products provided by the National Aeronautics and Space Administration’s (NASA) Sherlock ATM Data Warehouse. Depending on use cases, weather data can be queried using a Java application programming interface, web data access, and map tile image access. Map tile images are generated using triangular meshes obtained using the Matching Squares contour algorithm. An analysis determines a grid size that will produce both high image quality and fast image rendering. In addition, the visualization of weather data as map tiles on a TestBed visualization tool called Traffic Viewer is presented.

Chok Fung Lai↗

Asteroseismic Inversions of Mixed Acoustic-Gravity Modes to Probe the Stellar Core Structure

The discovery of mixed acoustic-gravity modes of oscillations of moderate mass stars opens a unique opportunity to infer the structure of the inner energy-generating cores and thus test the stellar evolution theory. The mixed modes have properties of internal gravity waves (g-modes) in the convectively stable helium core and properties of acoustic modes outside the core. We select several sets of the oscillation mode frequencies in the mass range from about 1.3 to 1.6 solar masses from the Kepler Legacy database, and apply the optimally localized averaging inversion technique previously developed for low-degree helioseismology. The inversion technique takes into account the uncertainties in the determination of the mass and radius of the stars, as well as the surface effects. The methodology provides sensitivity kernels for various structure properties, including the sound speed, density, and Ledoux parameter of convective stability, and, thus, the direct relationship between the stellar properties and the deviation of observed frequencies from the reference models. The inversion results reveal significant deviations in the core structure from the reference models calculated using the MESA evolutionary code for the stellar parameters obtained by the asteroseismic model grid fitting. Our analysis shows that the best resolution of the inner helium core and surrounding shell is achieved in inversions for the Ledoux parameter.

SMD↗

Boundary-Layer Stability Analysis of the Mean Flows Obtained Using Unstructured Grids

Boundary-layer stability analyses of mean flows extracted from unstructured-grid Navier- Stokes solutions have been performed. A procedure has been developed to extract mean flow profiles from the FUN3D unstructured-grid solutions. Extensive code-to-code validations have been performed by comparing the extracted mean ows as well as the corresponding stability characteristics to the predictions based on structured-grid solutions. Comparisons are made on a range of problems from a simple at plate to a full aircraft configuration-a modified Gulfstream-III with a natural laminar flow glove. The future aim of the project is to extend the adjoint-based design capability in FUN3D to include natural laminar flow and laminar flow control by integrating it with boundary-layer stability analysis codes, such as LASTRAC.

Liao, Wei↗

A computational design system for rapid CFD analysis

A computation design system (CDS) is described in which these tools are integrated in a modular fashion. This CDS ties together four key areas of computational analysis: description of geometry; grid generation; computational codes; and postprocessing. Integration of improved computational fluid dynamics (CFD) analysis tools through integration with the CDS has made a significant positive impact in the use of CFD for engineering design problems. Complex geometries are now analyzed on a frequent basis and with greater ease.

Ascoli, E. P.↗

Aeroacoustic sensitivity analysis and optimal aeroacoustic design of turbomachinery blades

During the first year of the project, we have developed a theoretical analysis - and wrote a computer code based on this analysis - to compute the sensitivity of unsteady aerodynamic loads acting on airfoils in cascades due to small changes in airfoil geometry. The steady and unsteady flow though a cascade of airfoils is computed using the full potential equation. Once the nominal solutions have been computed, one computes the sensitivity. The analysis takes advantage of the fact that LU decomposition is used to compute the nominal steady and unsteady flow fields. If the LU factors are saved, then the computer time required to compute the sensitivity of both the steady and unsteady flows to changes in airfoil geometry is quite small. The results to date are quite encouraging, and may be summarized as follows: (1) The sensitivity procedure has been validated by comparing the results obtained by 'finite difference' techniques, that is, computing the flow using the nominal flow solver for two slightly different airfoils and differencing the results. The 'analytic' solution computed using the method developed under this grant and the finite difference results are found to be in almost perfect agreement. (2) The present sensitivity analysis is computationally much more efficient than finite difference techniques. We found that using a 129 by 33 node computational grid, the present sensitivity analysis can compute the steady flow sensitivity about ten times more efficiently that the finite difference approach. For the unsteady flow problem, the present sensitivity analysis is about two and one-half times as fast as the finite difference approach. We expect that the relative efficiencies will be even larger for the finer grids which will be used to compute high frequency aeroacoustic solutions. Computational results show that the sensitivity analysis is valid for small to moderate sized design perturbations. (3) We found that the sensitivity analysis provided important insight into how airfoils should be modified to improve their aeroelastic stability. Using this insight, we redesigned an aeroelastically unstable cascade to produce a stable cascade.

Hall, Kenneth C.↗

Current Grid Generation Strategies and Future Requirements in Hypersonic Vehicle Design, Analysis and Testing

Recent advances in computational power enable computational fluid dynamic modeling of increasingly complex configurations. A review of grid generation methodologies implemented in support of the computational work performed for the X-38 and X-33 are presented. In strategizing topological constructs and blocking structures factors considered are the geometric configuration, optimal grid size, numerical algorithms, accuracy requirements, physics of the problem at hand, computational expense, and the available computer hardware. Also addressed are grid refinement strategies, the effects of wall spacing, and convergence. The significance of grid is demonstrated through a comparison of computational and experimental results of the aeroheating environment experienced by the X-38 vehicle. Special topics on grid generation strategies are also addressed to model control surface deflections, and material mapping.

Papadopoulos, Periklis↗

Grid Resolution Study over Operability Space for a Mach 1.7 Low Boom External Compression Inlet

This paper presents a statistical methodology whereby the probability limits associated with CFD grid resolution of inlet flow analysis can be determined which provide quantitative information on the distribution of that error over the specified operability range. The objectives of this investigation is to quantify the effects of both random (accuracy) and systemic (biasing) errors associated with grid resolution in the analysis of the Lockheed Martin Company (LMCO) N+2 Low Boom external compression supersonic inlet. The study covers the entire operability space as defined previously by the High Speed Civil Transport (HSCT) High Speed Research (HSR) program goals. The probability limits in terms of a 95.0% confidence interval on the analysis data were evaluated for four ARP1420 inlet metrics, namely (1) total pressure recovery (PFAIP), (2) radial hub distortion (DPH/P), (3) ) radial tip distortion (DPT/P), and (4) ) circumferential distortion (DPC/P). In general, the resulting +/-0.95 delta Y interval was unacceptably large in comparison to the stated goals of the HSCT program. Therefore, the conclusion was reached that the "standard grid" size was insufficient for this type of analysis. However, in examining the statistical data, it was determined that the CFD analysis results at the outer fringes of the operability space were the determining factor in the measure of statistical uncertainty. Adequate grids are grids that are free of biasing (systemic) errors and exhibit low random (precision) errors in comparison to their operability goals. In order to be 100% certain that the operability goals have indeed been achieved for each of the inlet metrics, the Y+/-0.95 delta Y limit must fall inside the stated operability goals. For example, if the operability goal for DPC/P circumferential distortion is ≤0.06, then the forecast Y for DPC/P plus the 95% confidence interval on DPC/P, i.e. +/-0.95 delta Y, must all be less than or equal to 0.06.

Grid↗

Parallel Anisotropic Unstructured Grid Adaptation

Computational Fluid Dynamics (CFD) has become critical to the design and analysis of aerospace vehicles. Parallel grid adaptation that resolves multiple scales with anisotropy is identified as one of the challenges in the CFD Vision 2030 Study to increase the capacity and capability of CFD simulation. The Study also cautions that computer architectures are undergoing a radical change and dramatic increases in algorithm concurrency will be required to exploit full performance. This paper reviews four different methods to parallel anisotropic grid generation. They cover both ends of the spectrum: (i) using existing state-of-the-art software optimized for a single core and modifying it for parallel platforms and (ii) designing and implementing scalable software with incomplete, but rapidly maturating functionality. A brief overview for each grid adaptation system is presented in the context of a telescopic approach for multilevel concurrency. These methods employ different approaches to enable parallel execution, which provides a unique opportunity to illustrate the relative behavior of each approach. Qualitative and quantitative metric evaluations are used to draw lessons for future developments in this critical area for parallel CFD simulation.

Tsolakis, Christos↗

Validating a large geophysical data set: Experiences with satellite-derived cloud parameters

We are validating the global cloud parameters derived from the satellite-borne HIRS2 and MSU atmospheric sounding instrument measurements, and are using the analysis of these data as one prototype for studying large geophysical data sets in general. The HIRS2/MSU data set contains a total of 40 physical parameters, filling 25 MB/day; raw HIRS2/MSU data are available for a period exceeding 10 years. Validation involves developing a quantitative sense for the physical meaning of the derived parameters over the range of environmental conditions sampled. This is accomplished by comparing the spatial and temporal distributions of the derived quantities with similar measurements made using other techniques, and with model results. The data handling needed for this work is possible only with the help of a suite of interactive graphical and numerical analysis tools. Level 3 (gridded) data is the common form in which large data sets of this type are distributed for scientific analysis. We find that Level 3 data is inadequate for the data comparisons required for validation. Level 2 data (individual measurements in geophysical units) is needed. A sampling problem arises when individual measurements, which are not uniformly distributed in space or time, are used for the comparisons. Standard 'interpolation' methods involve fitting the measurements for each data set to surfaces, which are then compared. We are experimenting with formal criteria for selecting geographical regions, based upon the spatial frequency and variability of measurements, that allow us to quantify the uncertainty due to sampling. As part of this project, we are also dealing with ways to keep track of constraints placed on the output by assumptions made in the computer code. The need to work with Level 2 data introduces a number of other data handling issues, such as accessing data files across machine types, meeting large data storage requirements, accessing other validated data sets, processing speed and throughput for interactive graphical work, and problems relating to graphical interfaces.

Kahn, Ralph↗