Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface codes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,513 records · Page 84

Improvements to the Unstructured Mesh Generator MESH3D

The AIRPLANE process starts with an aircraft geometry stored in a CAD system. The surface is modeled with a mesh of triangles and then the flow solver produces pressures at surface points which may be integrated to find forces and moments. The biggest advantage is that the grid generation bottleneck of the CFD process is eliminated when an unstructured tetrahedral mesh is used. MESH3D is the key to turning around the first analysis of a CAD geometry in days instead of weeks. The flow solver part of AIRPLANE has proven to be robust and accurate over a decade of use at NASA. It has been extensively validated with experimental data and compares well with other Euler flow solvers. AIRPLANE has been applied to all the HSR geometries treated at Ames over the course of the HSR program in order to verify the accuracy of other flow solvers. The unstructured approach makes handling complete and complex geometries very simple because only the surface of the aircraft needs to be discretized, i.e. covered with triangles. The volume mesh is created automatically by MESH3D. AIRPLANE runs well on multiple platforms. Vectorization on the Cray Y-MP is reasonable for a code that uses indirect addressing. Massively parallel computers such as the IBM SP2, SGI Origin 2000, and the Cray T3E have been used with an MPI version of the flow solver and the code scales very well on these systems. AIRPLANE can run on a desktop computer as well. AIRPLANE has a future. The unstructured technologies developed as part of the HSR program are now targeting high Reynolds number viscous flow simulation. The pacing item in this effort is Navier-Stokes mesh generation.

Thomas, Scott D.↗

Thermal Infrared Hot Spot and Dependence on Canopy Geometry

We perform theoretical calculations of the canopy thermal infrared (TIR) hot spot using a first principles 3-D model described earlier. Various theoretical canopies of varying leaf size and for differing canopy height are used to illustrate the magnitude of the TIR effect. Our results are similar to predicted behavior in the reflective hot spot as a function of canopy geometry and comparable to TIR measurements from the literature and our own simple ground experiments. We apply the MODTRAN atmospheric code to estimate the at-sensor variation in brightness temperature with view direction in the solar principal plane. For simple homogeneous canopies, we predict canopy thermal infrared hot spot variations of 2 degrees C at the surface with respect to nadir viewing. Dependence on leaf size is weak as long as the ratio of leaf size to canopy height is maintained. However, the angular width of the hot spot increases as the ratio of leaf diameter to canopy height increases. Atmospheric effects minimize but do not eliminate the TIR hot spot at satellite altitudes.

Smith, James A.↗

A New Modular Approach for Tightly Coupled Fluid/Structure Analysis

Static aeroelastic computations are made using a C++ executive suitable for closely coupled fluid/structure interaction studies. The fluid flow is modeled using the Euler/Navier Stokes equations and the structure is modeled using finite elements. FORTRAN based fluids and structures codes are integrated under C++ environment. The flow and structural solvers are treated as separate object files. The data flow between fluids and structures is accomplished using I/O. Results are demonstrated for transonic flow over partially flexible surface that is important for aerospace vehicles. Use of this development to accurately predict flow induced structural failure will be demonstrated.

Guruswamy, Guru↗

CFD-Based Design Optimization Tool Developed for Subsonic Inlet

The traditional approach to the design of engine inlets for commercial transport aircraft is a tedious process that ends with a less-than-optimum design. With the advent of high-speed computers and the availability of more accurate and reliable computational fluid dynamics (CFD) solvers, numerical optimization processes can effectively be used to design an aerodynamic inlet lip that enhances engine performance. The designers' experience at Boeing Corporation showed that for a peak Mach number on the inlet surface beyond some upper limit, the performance of the engine degrades excessively. Thus, our objective was to optimize efficiency (minimize the peak Mach number) at maximum cruise without compromising performance at other operating conditions. Using a CFD code NPARC, the NASA Lewis Research Center, in collaboration with Boeing, developed an integrated procedure at Lewis to find the optimum shape of a subsonic inlet lip and a numerical optimization code, ADS. We used a GRAPE-based three-dimensional grid generator to help automate the optimization procedure. The inlet lip shape at the crown and the keel was described as a superellipse, and the superellipse exponents and radii ratios were considered as design variables. Three operating conditions: cruise, takeoff, and rolling takeoff, were considered in this study. Three-dimensional Euler computations were carried out to obtain the flow field. At the initial design, the peak Mach numbers for maximum cruise, takeoff, and rolling takeoff conditions were 0.88, 1.772, and 1.61, respectively. The acceptable upper limits on the takeoff and rolling takeoff Mach numbers were 1.55 and 1.45. Since the initial design provided by Boeing was found to be optimum with respect to the maximum cruise condition, the sum of the peak Mach numbers at takeoff and rolling takeoff were minimized in the current study while the maximum cruise Mach number was constrained to be close to that at the existing design. With this objective, the optimum design satisfied the upper limits at takeoff and rolling takeoff while retaining the desirable cruise performance. Further studies are being conducted to include static and cross-wind operating conditions in the design optimization procedure. This work was carried out in collaboration with Dr. E.S. Reddy of NYMA, Inc.

Source record↗

A Novel Porosity Model for Use in Hydrocode Simulations

Introduction: Numerical modeling of impact cratering has reached a high degree of sophistication; however, the treatment of porous materials still poses a large problem in hydrocode calculations. Porosity plays only a minor role in the formation of large craters on most planetary objects, but impacts on comets are believed to be highly affected by the presence of porosity, which may be as much as 80%. The upcoming Deep Impact Mission (launched January 2005) will provide more detailed data about the composition of a comet (Tempel 1) by shooting a approx.370 kg projectile onto the surface of its nucleus. The numerical simulations of such impact events requires an appropriate model for how pore space in the comet is crushed out during the violent initial stage of the impact event. Most hydro-codes compute the pressure explicitly using an "equation of state" (EOS) for each material, which relates changes in density and internal energy to changes in pressure. The added complication introduced by porosity is that changes in a material s density are due to both the closing of pore space (compaction) and compression of the matrix. The amount of resistance to volume change and the amount of irreversible work done during these two processes is very different; it is far easier to compact a porous material sample than to compress a non-porous sample of the same material. As an alternative to existing porosity models, like the Pdot(alpha) model [1], we present a novel approach for dealing with the compaction of porosity in hydrocode calculations.

Wuennemann, K.↗

Comparison of Martian Surface Radiation Predictions to the Measurements of Mars Science Laboratory Radiation Assessment Detector (MSL/RAD)

For the analysis of radiation risks to astronauts and planning exploratory space missions, detailed knowledge of particle spectra is an important factor. Detailed measurements of the energetic particle radiation environment on the surface of Mars have been made by the Mars Science Laboratory Radiation Assessment Detector (MSL-RAD) on the Curiosity rover since August 2012, and particle fluxes for a wide range of ion species (up to several hundred MeV/u) and high energy neutrons (8 - 1000 MeV) have been available for the first 200 sols. Although the data obtained on the surface of Mars for 200 sols are limited in the narrow energy spectra, the simulation results using the Badhwar-O'Neill galactic cosmic ray (GCR) environment model and the high-charge and energy transport (HZETRN) code are compared to the data. For the nuclear interactions of primary GCR through Mars atmosphere and Curiosity rover, the quantum multiple scattering theory of nuclear fragmentation (QMSFRG) is used, which includes direct knockout, evaporation and nuclear coalescence. Daily atmospheric pressure measurements at Gale Crater by the MSL Rover Environmental Monitoring Station are implemented into transport calculations for describing the daily column depth of atmosphere. Particles impinging on top of the Martian atmosphere reach the RAD after traversing varying depths of atmosphere that depend on the slant angles, and the model accounts for shielding of the RAD by the rest of the instrument. Calculations of stopping particle spectra are in good agreement with the RAD measurements for the first 200 sols by accounting changing heliospheric conditions and atmospheric pressure. Detailed comparisons between model predictions and spectral data of various particle types provide the validation of radiation transport models, and thus increase the accuracy of the predictions of future radiation environments on Mars. These contributions lend support to the understanding of radiation health risks to astronauts for the planning of various mission scenarios.

Kim, Myung-Hee Y.↗

Mineralogy and Astrobiology Detection Using Laser Remote Sensing Instrument

A multispectral instrument based on Raman, laser-induced fluorescence (LIF), laser-induced breakdown spectroscopy (LIBS), and a lidar system provides high-fidelity scientific investigations, scientific input, and science operation constraints in the context of planetary field campaigns with the Jupiter Europa Robotic Lander and Mars Sample Return mission opportunities. This instrument conducts scientific investigations analogous to investigations anticipated for missions to Mars and Jupiter's icy moons. This combined multispectral instrument is capable of performing Raman and fluorescence spectroscopy out to a >100 m target distance from the rover system and provides single-wavelength atmospheric profiling over long ranges (>20 km). In this article, we will reveal integrated remote Raman, LIF, and lidar technologies for use in robotic and lander-based planetary remote sensing applications. Discussions are focused on recently developed Raman, LIF, and lidar systems in addition to emphasizing surface water ice, surface and subsurface minerals, organics, biogenic, biomarker identification, atmospheric aerosols and clouds distributions, i.e., near-field atmospheric thin layers detection for next robotic-lander based instruments to measure all the above-mentioned parameters. OCIS codes: (120.0280) Remote sensing and sensors; (130.0250) Optoelectronics; (280.3640) Lidar; (300.2530) Fluorescence, laser-induced; (300.6450) Spectroscopy, Raman; (300.6365) Spectroscopy, laser induced breakdown

Abedin, M. Nurul↗

Progress Towards Modeling the Ablation Response of NuSil-Coated PICA

The Mars Science Laboratory (MSL) Entry, Descent and Landing Instrumentation (MEDLI) collected in-flight data largely used by the ablation community to verify and validate physics-based models for the response of the Phenolic Impregnated Carbon Ablator (PICA) material [1-4]. MEDLI data were recently used to guide the development of NASA’s high-fidelity material response models for PICA, implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software [5-6]. A follow-up instrumentation suite, MEDLI2, is planned for the upcoming Mars 2020 mission [7] after the large scientific impact of MEDLI. Recent analyses performed as part of MEDLI2 development draw the attention to significant effects of a protective coating to the aerothermal response of PICA. NuSil, a silicone-based overcoat sprayed onto the MSL heatshield as contamination control, is currently neglected in PICA ablation models. To mitigate the spread of phenolic dust from PICA, NuSil was applied to the entire MSL heatshield, including the MEDLI plugs. NuSil is a space grade designation of the siloxane copolymer, primarily used to protect against atomic oxygen erosion in the Low Earth Orbit environment. Ground testing of PICA-NuSil (PICA-N) models all exhibited surface temperature jumps of the order of 200 K due to oxide scale formation and subsequent NuSil burn-off. It is therefore critical to include a model for the aerothermal response of the coating in ongoing code development and validation efforts.

Meurisse, Jeremie B.↗

Optimization of Elastodynamic Finite Integration Technique on Intel Xeon Phi Knights Landing Processors

This work describes the development and optimization of an implementation of an isotropic elastodynamic finite integration technique (EFIT) code for parallelized computation on Intel Knights Landing (KNL) hardware. EFIT is a numerical approach resulting in standard staggered-grid finite difference equations for the elastodynamic equations of motion to simulate bulk waves is solids. The computationally efficient simulation of elastodynamic wave propagation and interactions in aerospace materials is of high-interest in the fields of nondestructive evaluation (NDE) and structural health monitoring (SHM). Ultrasonic inspection uses an ultrasonic signal, generated at the surface of the material/structure via use of a piezoelectric transducer, to propagate sound waves into the material where it interacts with any existing defects, as well as with structural boundaries and any material inhomogeneity. Reflections from defects and boundaries are then measured by a transducer. Realistic ultrasound simulation tools can significantly aid the development and optimization of inspection techniques and can assist in the interpretation of experimental data. The optimization of an elastodynamics simulation code for the KNL Many Integrated Core processor was performed. The optimization focused on data locality and vectorization. Results show that tiling of the data to exploit the cache behavior and allow for significant utilization of the KNL hardware. The MPI implementation allows for a scalable implementation enabling large problems to be simulated. The model results were validated against theoretical dispersion curves to within 2% of the group velocity, and within 0.5% of the phase velocity of the A0 mode. Aggressive use of tiling, threading, and vectorization techniques allowed for dramatically improved time to solution.

William C. Schneck↗

A flow visualization study of the leading edge separation bubble on a NACA 0012 airfoil with simulated glaze ice

As a part of an ongoing research in aircraft icing, the leading edge separation bubble on a NACA 0012 model with a 5 min simulated glaze ice was investigated. The flow visualization methods using oil, tuft, splitter plate, smoke, and liquid crystals were employed to get reattachment line data for the leading edge separation bubble on both surfaces of the airfoil. On the upper surface, the separation bubble was found to grow larger with increasing angles of attack and reduce in size with decreasing angle of attack. On the lower surface, the bubble was found to grow larger with increasing negative angles of attack and reduce in size with increasing angles of attack. The separated flow fails to reattach beyond 6° for the upper surface and -5° for the loser surface. Their results of this study compared well with those of other experiments and computational results.

Aircraft icing↗

TranAir: A full-potential, solution-adaptive, rectangular grid code for predicting subsonic, transonic, and supersonic flows about arbitrary configurations. User's manual

The TranAir computer program calculates transonic flow about arbitrary configurations at subsonic, transonic, and supersonic freestream Mach numbers. TranAir solves the nonlinear full potential equations subject to a variety of boundary conditions modeling wakes, inlets, exhausts, porous walls, and impermeable surfaces. Regions with different total temperature and pressure can be represented. The user's manual describes how to run the TranAir program and its graphical support programs.

Johnson, F. T.↗

Computational procedure for evaluation of space shuttle TPS requirements

The rudiments of a computational procedure are presented for predicting heat transfer rates to the windward pitchplane of a shuttle vehicle under conditions of nonequilibrium flow. The procedure involves computations in three regions of the flow field, namely, the stagnation region, inviscid flow region, and boundary layer region. The thermochemical state in the inviscid flow was calculated with the aid of correlations developed from shadowgraphs of shuttle configurations. The flow in the boundary layer and stagnation regions was computed with a nonequilibrium viscous layer code. For boundary layers, the outer edge thermochemical state was determined by matching the mass flow in the boundary layer with the mass flow represented by each inviscid streamline. Typical solutions are presented to show the effects of surface pressure, catalycity, and homogeneous chemical reactions on flow field variables and surface heat transfer rates.

Tong, H.↗

Analysis of high-velocity and real-gas effects on the shock-interference pattern for delta-wing orbiters

Numerical solutions have been generated to study the two-dimensional flow field which results when supersonic flow encounters double-wedge configurations whose angles are such that a Type VI pattern occurs. Two numerical codes were developed: one which used the perfect-gas relations (including the ability to vary gamma to 'simulate' real-gas effects) and a second which incorporated a Mollier table to define equilibrium air properties. Theoretical surface-pressure and heat-transfer distributions in the interaction perturbed region along the simulated 'wing-leading edge' were obtained for velocities from 1167 m/sec to 7610 m/sec.

Bertin, J. J.↗

The decrease in effective photocurrents due to saddle points in electrostatic potentials near differentially charged spacecraft

The reported investigation had the objective to illustrate the presence of important multidimensional effects in spacecraft charging. Two-dimensional codes have been under development by Parker (1976). A description is presented of a calculation which was performed using the three-dimensional NASA Charging Analyzer Program (NASCAP). NASCAP was run to calculate the electrostatic potentials on the surface of, and in the space surrounding, a sunlit Teflon-coated sphere. Currents to the sunlit surfaces were determined on the basis of an approximate photosheath model for strong differential charging.

Mandell, M. J.↗

An initial investigation into methods of computing transonic aerodynamics sensitivity coefficients

Based upon the program outlined in the modified proposal, the initial effort has concentrated on developing the quasi-analytical approach for determining aerodynamic influence coefficients. In addition, in order to keep these initial efforts computationally efficient and straightforward, it was decided to first consider only two-dimensional flow and to model the problem using the transonic small perturbation equation with surface boundary conditions applied along the slit. The small perturbation equation was selected instead of the full potential equation because it permits the capturing of important nonlinear physical phenomena while being by comparison easy to handle from a coding standpoint. Consequently, the initial objectives were: (1) investigate the feasibility of using the quasi-analytical approach for determining aerodynamic sensitivity coefficients in the transonic regime; (2) determine the accuracy of the quasi-analytical approach by comparing results to those computed using finite difference techniques; and (3) investigate efficient methods for solving the large set of equations associated with the quasi-analytical approach.

Carlson, Leland A.↗

Three-dimensional elliptic grid generation technique with application to turbomachinery cascades

Described is a numerical method for generating 3-D grids for turbomachinery computational fluid dynamic codes. The basic method is general and involves the solution of a quasi-linear elliptic partial differential equation via pointwise relaxation with a local relaxation factor. It allows specification of the grid point distribution on the boundary surfaces, the grid spacing off the boundary surfaces, and the grid orthogonality at the boundary surfaces. A geometry preprocessor constructs the grid point distributions on the boundary surfaces for general turbomachinery cascades. Representative results are shown for a C-grid and an H-grid for a turbine rotor. Two appendices serve as user's manuals for the basic solver and the geometry preprocessor.

Chen, S. C.↗

Aeroelastic tailoring for oblique wing lateral trim

Composite material aeroelastic tailoring is presently explored as a means for the correction of the roll trim imbalance of oblique-wing aircraft configurations. The concept is demonstrated through the analysis of a realistic oblique wing by a static aeroelastic computational procedure encompassing the full potential transonic aerodynamic code FLO22 and a Ritz structural plate program that models the stiffness due to symmetrical-but-unbalanced composite wing skins. Results indicate that asymetric composite tailoring reduces the aileron deflection needed for roll equilibrium, and reduces control surface hinge moment and drag. Wing skin stresses are, however, very high.

Bohlmann, Jonathan D.↗

Heat transfer, velocity-temperature correlation, and turbulent shear stress from Navier-Stokes computations of shock wave/turbulent boundary layer interaction flows

The properties of 2-D shock wave/turbulent boundary layer interaction flows were calculated by using a compressible turbulent Navier-Stokes numerical computational code. Interaction flows caused by oblique shock wave impingement on the turbulent boundary layer flow were considered. The oblique shock waves were induced with shock generators at angles of attack less than 10 degs in supersonic flows. The surface temperatures were kept at near-adiabatic (ratio of wall static temperature to free stream total temperature) and cold wall (ratio of wall static temperature to free stream total temperature) conditions. The computational results were studied for the surface heat transfer, velocity temperature correlation, and turbulent shear stress in the interaction flow fields. Comparisons of the computational results with existing measurements indicated that (1) the surface heat transfer rates and surface pressures could be correlated with Holden's relationship, (2) the mean flow streamwise velocity components and static temperatures could be correlated with Crocco's relationship if flow separation did not occur, and (3) the Baldwin-Lomax turbulence model should be modified for turbulent shear stress computations in the interaction flows.

Wang, C. R.↗