Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structured Grids”

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 847 records · Page 47

Two and three-dimensional shock-shock interactions on the blunt leading edges of the hypersonic inlets

The effect of shock impingement on the blunt leading edges of the top and sidewall compression type inlet of a scramjet engine is studied numerically. The impinging shock is caused by the vehicle forebody. The interaction of this forebody shock with the inlet leading edge shock results in a very complex flowfield containing local regions of high pressure and intense heating. This complex flowfield in calculated by solving the Navier-Stokes equations using a finite volume flux splitting technique due to van Leer. To resolve the finer details of the flow structure as well as to predict the surface heat transfer accurately, adaptive grid technique is used in the analysis. Results of the present numerical study are compared with available experimental results.

Singh, D. J.↗

Partitioning of unstructured problems for parallel processing

Many large-scale computational problems are based on unstructured computational domains. Primary examples are unstructured grid calculations based on finite volume methods in computational fluid dynamics, or structural analysis problems based on finite element approximations. The question of how to distribute such unstructured computational domains over a large number of processors in a MIMD machine with distributed memory is addressed. A graph theoretical framework for these problems is established. Based on this framework three decomposition algorithms are introduced. In particular a new decomposition algorithm is discussed, which is based on the computation of an eigenvector of the Laplacian matrix associated with the graph. Numerical comparisons on large-scale two- and three-dimensional problems demonstrate the superiority of the new spectral bisection algorithm.

Simon, H. D.↗

Methods of testing parameterizations: Vertical ocean mixing

The ocean's velocity field is characterized by an exceptional variety of scales. While the small-scale oceanic turbulence responsible for the vertical mixing in the ocean is of scales a few centimeters and smaller, the oceanic general circulation is characterized by horizontal scales of thousands of kilometers. In oceanic general circulation models that are typically run today, the vertical structure of the ocean is represented by a few tens of discrete grid points. Such models cannot explicitly model the small-scale mixing processes, and must, therefore, find ways to parameterize them in terms of the larger-scale fields. Finding a parameterization that is both reliable and plausible to use in ocean models is not a simple task. Vertical mixing in the ocean is the combined result of many complex processes, and, in fact, mixing is one of the less known and less understood aspects of the oceanic circulation. In present models of the oceanic circulation, the many complex processes responsible for vertical mixing are often parameterized in an oversimplified manner. Yet, finding an adequate parameterization of vertical ocean mixing is crucial to the successful application of ocean models to climate studies. The results of general circulation models for quantities that are of particular interest to climate studies, such as the meridional heat flux carried by the ocean, are quite sensitive to the strength of the vertical mixing. We try to examine the difficulties in choosing an appropriate vertical mixing parameterization, and the methods that are available for validating different parameterizations by comparing model results to oceanographic data. First, some of the physical processes responsible for vertically mixing the ocean are briefly mentioned, and some possible approaches to the parameterization of these processes in oceanographic general circulation models are described in the following section. We then discuss the role of the vertical mixing in the physics of the large-scale ocean circulation, and examine methods of validating mixing parameterizations using large-scale ocean models.

Tziperman, Eli↗

Performance and Applications of ENSAERO-MPI on Scalable Components

The latest improvements and results generated by ENSAERO-MPI are presented in this paper. ENSAERO-MPI is a parallelized, high-fidelity, multi-block code with fluids, structures and controls capabilities developed at NASA Ames Research Center under the support of HPCC. It is capable of multidisciplinary simulations by simultaneously integrating the Navier-Stokes equations, the finite element structural equations as well as control dynamics equations using aeroelastically adaptive, patched grids. Improvements have been made to the code's robustness, moving grid capabilities and performance.

Farhangnia, Mehrdad↗

Discovery of a New Cassini-Size Basin on Mars from MOLA Topographic Data

MOLA profile data collected during the Science Phasing Operations have revealed a previously unknown, 450 km wide, 2 km deep basin on Mars centered at 30N, 312W near the Phison Rupes. This structure, as large as and somewhat deeper than the very obvious Cassini impact basin located 1000 km to the SW, is not apparent in the existing good quality Viking imagery. Gridded MOLA data show the basin as a closed depression, but elevation contours show only weak correlation with what little structure exists in the area and with mapped geologic units. From analysis of slope breaks readily visible in two MOLA profiles we suggest this is at least a three-ring basin. Portions of rings from concentric fits to slope breaks align with some of the linear ridge-like structures of the Phison Rupes, and outline a region of lower crater density and smoother inter-crater plains. The discovery of such a pronounced topographic depression which lacks obvious visible structure suggests more such previously unknown features may exist on Mars, and that MOLA data may be useful in finding them.

Frey, H.↗

Numerical Simulation of Inlet Bleed with Circular Holes on Plate Under Shock-Wave/Boundary-Layer Interactions

A numerical study was performed to investigate the shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. The bleed-hole patterns considered for the study include in-line multiple holes and staggered multiple-row holes that are configured to simulate the patterns used in inlet bleed systems of high performance aircraft. The focus of the study was to examine how the bleed through multiple holes affect bleed rate and the pressure and Mach number distributions. Since the bleed performance was found sensitive to the change in bleed conditions, a computational procedure was developed to give a good turnaround computational time for parametric studies involving changes in bleed hole geometry and the structure of shock-wave/boundary-layer flowfield. The procedure includes the grid-generation methodology and the flow simulation with solutions from the Navier-Stokes equations. The computational techniques permit analysis of complex bleed systems and make possible the investigation of a broader range of design variables associated with inlet bleed operation.

Chyu, Wei J.↗

Fabricating Sub-collimating Grids for an X-ray Solar Imaging Spectrometer using LIGA Techniques

This paper describes fabrication of sub-collimating X-ray grids that are used in an instrument for the High Energy Solar Spectroscopic Imager ( HESSI ), a proposed NASA mission that will promote understanding of solar particle acceleration and explosive energy release in the magnetized plasmas at the Sun. The HESSI comprises 12 rotating pairs of high-aspect-ratio, high-Z grids; each pair is separated by 1.7 meters and backed by a single Ge detector.

x-ray↗

Fine Pitch Grids for an X-Ray Solar Imaging Spectrometer Fabricated by Optical Lithography and XeF2 Etching

Fine-pitch, sub-collimating X-ray grids have been developed for an instrument in the High Energy Solar Spectroscopic Imager , a proposed NASA mission. In addition to high-energy x- rays, the instrument requires collimation of photons with energies of less than 4 keV such that free-standing grids are required that have no material between the grid slats. Fabrication of these grids is described.

x-ray↗

Conical-Domain Model for Estimating GPS Ionospheric Delays

The conical-domain model is a computational model, now undergoing development, for estimating ionospheric delays of Global Positioning System (GPS) signals. Relative to the standard ionospheric delay model described below, the conical-domain model offers improved accuracy. In the absence of selective availability, the ionosphere is the largest source of error for single-frequency users of GPS. Because ionospheric signal delays contribute to errors in GPS position and time measurements, satellite-based augmentation systems (SBASs) have been designed to estimate these delays and broadcast corrections. Several national and international SBASs are currently in various stages of development to enhance the integrity and accuracy of GPS measurements for airline navigation. In the Wide Area Augmentation System (WAAS) of the United States, slant ionospheric delay errors and confidence bounds are derived from estimates of vertical ionospheric delay modeled on a grid at regularly spaced intervals of latitude and longitude. The estimate of vertical delay at each ionospheric grid point (IGP) is calculated from a planar fit of neighboring slant delay measurements, projected to vertical using a standard, thin-shell model of the ionosphere. Interpolation on the WAAS grid enables estimation of the vertical delay at the ionospheric pierce point (IPP) corresponding to any arbitrary measurement of a user. (The IPP of a given user s measurement is the point where the GPS signal ray path intersects a reference ionospheric height.) The product of the interpolated value and the user s thin-shell obliquity factor provides an estimate of the user s ionospheric slant delay. Two types of error that restrict the accuracy of the thin-shell model are absent in the conical domain model: (1) error due to the implicit assumption that the electron density is independent of the azimuthal angle at the IPP and (2) error arising from the slant-to-vertical conversion. At low latitudes or at mid-latitudes under disturbed conditions, the accuracy of SBAS systems based upon the thin-shell model suffers due to the presence of complex ionospheric structure, high delay values, and large electron density gradients. Interpolation on the vertical delay grid serves as an additional source of delay error. The conical-domain model permits direct computation of the user s slant delay estimate without the intervening use of a vertical delay grid. The key is to restrict each fit of GPS measurements to a spatial domain encompassing signals from only one satellite. The conical domain model is so named because each fit involves a group of GPS receivers that all receive signals from the same GPS satellite (see figure); the receiver and satellite positions define a cone, the satellite position being the vertex. A user within a given cone evaluates the delay to the satellite directly, using (1) the IPP coordinates of the line of sight to the satellite and (2) broadcast fit parameters associated with the cone. The conical-domain model partly resembles the thin-shell model in that both models reduce an inherently four-dimensional problem to two dimensions. However, unlike the thin-shell model, the conical domain model does not involve any potentially erroneous simplifying assumptions about the structure of the ionosphere. In the conical domain model, the initially four-dimensional problem becomes truly two-dimensional in the sense that once a satellite location has been specified, any signal path emanating from a satellite can be identified by only two coordinates; for example, the IPP coordinates. As a consequence, a user s slant-delay estimate converges to the correct value in the limit that the receivers converge to the user s location (or, equivalently, in the limit that the measurement IPPs converge to the user s IPP).

Sparks, Lawrence↗

Infrared-Bolometer Arrays with Reflective Backshorts

Integrated circuits that incorporate square arrays of superconducting-transition- edge bolometers with optically reflective backshorts are being developed for use in image sensors in the spectral range from far infrared to millimeter wavelengths. To maximize the optical efficiency (and, thus, sensitivity) of such a sensor at a specific wavelength, resonant optical structures are created by placing the backshorts at a quarter wavelength behind the bolometer plane. The bolometer and backshort arrays are fabricated separately, then integrated to form a single unit denoted a backshort-under-grid (BUG) bolometer array. In a subsequent fabrication step, the BUG bolometer array is connected, by use of single-sided indium bump bonding, to a readout device that comprises mostly a superconducting quantum interference device (SQUID) multiplexer circuit. The resulting sensor unit comprising the BUG bolometer array and the readout device is operated at a temperature below 1 K. The concept of increasing optical efficiency by use of backshorts at a quarter wavelength behind the bolometers is not new. Instead, the novelty of the present development lies mainly in several features of the design of the BUG bolometer array and the fabrication sequence used to implement the design. Prior to joining with the backshort array, the bolometer array comprises, more specifically, a square grid of free-standing molybdenum/gold superconducting-transition-edge bolometer elements on a 1.4- m-thick top layer of silicon that is part of a silicon support frame made from a silicon-on-insulator wafer. The backshort array is fabricated separately as a frame structure that includes support beams and contains a correspond - ing grid of optically reflective patches on a single-crystal silicon substrate. The process used to fabricate the bolometer array includes standard patterning and etching steps that result in the formation of deep notches in the silicon support frame. These notches are designed to interlock with the support beams on the backshort-array structure to provide structural support and precise relative positioning. The backshort-array structure is inserted in the silicon support frame behind the bolometer array, and the notches in the frame serve to receive the support beams of the backshort-array structure and thus determine the distance between the backshort and bolometer planes. The depth of the notches and, thus, the distance between the backshort and bolometer planes, can be tailored to a value between 25 to 300 m adjusting only a few process steps. The backshort array is designed so as not to interfere with the placement of indium bumps for subsequent indium bump-bonding to the multiplexing readout circuitry

Miller, Timothy M.↗

Sensitivity Analysis of Multidisciplinary Rotorcraft Simulations

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. An unstructured sensitivity-enabled Navier-Stokes solver, FUN3D, and a nonlinear flexible multibody dynamics solver, DYMORE, are coupled to predict the aerodynamic loads and structural responses of helicopter rotor blades. A discretely-consistent adjoint-based sensitivity analysis available in FUN3D provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute DYMORE structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Numerical results verify accuracy of the FUN3D/DYMORE system by conducting simulations for a benchmark rotorcraft test model and comparing solutions with established analyses and experimental data. Complex-variable implementation of sensitivity analysis of DYMORE and the coupled FUN3D/DYMORE system is verified by comparing with real-valued analysis and sensitivities. Correctness of adjoint formulations for FUN3D/DYMORE interfaces is verified by comparing adjoint-based and complex-variable sensitivities. Finally, sensitivities of the lift and drag functions obtained by complex-variable FUN3D/DYMORE simulations are compared with sensitivities computed by the multidisciplinary sensitivity analysis, which couples adjoint-based flow and grid sensitivities of FUN3D and FUN3D/DYMORE interfaces with complex-variable sensitivities of DYMORE structural responses.

Wang, Li↗

Fluid Structure Interaction in a Turbine Blade

An unsteady, three dimensional Navier-Stokes solution in rotating frame formulation for turbomachinery applications is presented. Casting the governing equations in a rotating frame enabled the freezing of grid motion and resulted in substantial savings in computer time. The turbine blade was computationally simulated and probabilistically evaluated in view of several uncertainties in the aerodynamic, structural, material and thermal variables that govern the turbine blade. The interconnection between the computational fluid dynamics code and finite element structural analysis code was necessary to couple the thermal profiles with the structural design. The stresses and their variations were evaluated at critical points on the Turbine blade. Cumulative distribution functions and sensitivity factors were computed for stress responses due to aerodynamic, geometric, mechanical and thermal random variables.

Gorla, Rama S. R.↗

Numerical simulation of channel flow transition, resolution requirements and structure of the hairpin vortex

Three-dimensional, nonlinear numerical simulations are presented for the K-type and H-type transitions for channel flow. There are two objectives. The first is to establish firmly the resolution requirements for the various stages in the transition process. Comparisons between calculations on various grids suggest a set of guidelines for maintaining a physically meaningful calculation. The second objective is to map out the structure of the hairpin vortices which arise in K-type and H-type transitions in channel flow, to the latest stage currently feasible. Flow field details are presented for both a subcritical Reynolds number of 1500 and a supercritical Reynolds number of 8000. The diagnostics include illustrations of the vertical shear, streamwise and spanwise vorticity, helicity, vortex stretching, and vortex diffusion fields.

Krist, Steven E.↗

User's Manual and Final Report for Hot-SMAC GUI Development

A new software package called Higher Order Theory-Structural/Micro Analysis Code (HOT-SMAC) has been developed as an effective alternative to the finite element approach for Functionally Graded Material (FGM) modeling. HOT-SMAC is a self-contained package including pre- and post-processing through an intuitive graphical user interface, along with the well-established Higher Order Theory for Functionally Graded Materials (HOTFGM) thermomechanical analysis engine. This document represents a Getting Started/User's Manual for HOT-SMAC and a final report for its development. First, the features of the software are presented in a simple step-by-step example where a HOT-SMAC model representing a functionally graded material is created, mechanical and thermal boundary conditions are applied, the model is analyzed and results are reviewed. In a second step-by-step example, a HOT-SMAC model of an actively cooled metallic channel with ceramic thermal barrier coating is built and analyzed. HOT-SMAC results from this model are compared to recently published results (NASA/TM-2001-210702) for two grid densities. Finally, a prototype integration of HOTSMAC with the commercially available HyperSizer(R) structural analysis and sizing software is presented. In this integration, local strain results from HyperSizer's structural analysis are fed to a detailed HOT-SMAC model of the flange-to-facesheet bond region of a stiffened panel. HOT-SMAC is then used to determine the peak shear and peel (normal) stresses between the facesheet and bonded flange of the panel and determine the "free edge" effects.

Yarrington, Phil↗

The grid array antenna

Attention is called to the grid array as a possible useful antenna design for UHF clear air radars. This type of antenna integrates radiating elements and the feed network into a single structure so that a fairly large array can be driven from a single feed point. Figures are given which demonstrate the basic principle of the grid array, which are adapted from Conti et al. (1981). Conductors are arranged above a ground plane in a repeating, staggered array of connected rectangles. Each rectangular element is approximately one by one-half wavelength in size. The currents on the conductors at resonance is indicated. The grid array illustrated can be expanded in both vertical and horizontal directions about the feed point by adding additional rectangular conductors. The design eliminates the feed network and would provide a thin, panel-like antenna that could be easily built and transported.

Ecklund, W. L.↗

NASTRAN data generation of helicopter fuselages using interactive graphics

The development and implementation of a preprocessor system for the finite element analysis of helicopter fuselages is described. The system utilizes interactive graphics for the generation, display, and editing of NASTRAN data for fuselage models. It is operated from an IBM 2250 cathode ray tube (CRT) console driven by an IBM 370/145 computer. Real time interaction plus automatic data generation reduces the nominal 6 to 10 week time for manual generation and checking of data to a few days. The interactive graphics system consists of a series of satellite programs operated from a central NASTRAN Systems Monitor. Fuselage structural models including the outer shell and internal structure may be rapidly generated. All numbering systems are automatically assigned. Hard copy plots of the model labeled with GRID or elements ID's are also available. General purpose programs for displaying and editing NASTRAN data are included in the system. Utilization of the NASTRAN interactive graphics system has made possible the multiple finite element analysis of complex helicopter fuselage structures within design schedules.

Sainsbury-Carter, J. B.↗

A preliminary study of numerical simulation of thermosolutal convection of interest to crystal growth

Calculations were performed with computer models using three types of finite difference methods of thermosolutal convection: horizontal heating of a container filled with a stably stratified solution, finger convection in a container, and finger convection in a horizontally infinite channel. The importance of including thermosolutal convection in models of crystal growth is emphasized, and the difficulties in doing so are demonstrated. It is pointed out that these difficulties, due primarily to the fine structure of the convection, may be partly overcome by the use of fine grids and implicit time stepping methods.

Miller, T. L.↗

A 3-D chimera grid embedding technique

A three-dimensional (3-D) chimera grid-embedding technique is described. The technique simplifies the construction of computational grids about complex geometries. The method subdivides the physical domain into regions which can accommodate easily generated grids. Communication among the grids is accomplished by interpolation of the dependent variables at grid boundaries. The procedures for constructing the composite mesh and the associated data structures are described. The method is demonstrated by solution of the Euler equations for the transonic flow about a wing/body, wing/body/tail, and a configuration of three ellipsoidal bodies.

Benek, J. A.↗