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 901 records · Page 50

Minimizing distortion in truss structures via Tabu search

The shape control of large flexible space structures is of great interest to structural designers. A related problem is to seek ways to minimize the need for active controls by careful design and construction of the space structure. A tetrahedral truss structure that is used to support a precision segmented reflector or antenna surface is considered. The structure has a hexagonal platform and is characterized by the number of rings of members in the truss. For simplicity it is assumed that a flat truss geometry exists. Hence, all structural members and ball joints are required to have the same nominal length and diameter, respectively. Inaccuracies in the length of member or diameters of joints may produce unacceptable levels of surface distortion and internal forces. In the case of a truss structure supporting an antenna, surface distortions may cause unacceptable gain loss or pointing errors. The focus is solely on surface distortion, however, internal forces may be treated in a similar manner. To test the Tabu search code for DSQRMS the appropriate influence matrices are used for a flat, two-ring tetrahedral reflector truss generated by Green and Haftka (1989). In this example there are 102 members (NMEMB) and 31 ball joints (NJOINT) of the same nominal length, respectively. Hence, all the members may be interchanged and all the joints may be interchanged. In addition, 19 positions on the surface of the truss (NNODES) were used to measure error influences. After a variety of experiments a set of good parameters was choosen for Tabu search. The sample size at each iteration is 10*NMEMB and the short term memory size is 40. In addition four pruning rules were used to accelerate the search..

Kincaid, Rex K.↗

Hypersonic Navier-Stokes Comparisons to Orbiter Flight Data

During the STS-119 flight of Space Shuttle Discovery, two sets of surface temperature measurements were made. Under the HYTHIRM program3 quantitative thermal images of the windward side of the Orbiter with a were taken. In addition, the Boundary Layer Transition Flight Experiment 4 made thermocouple measurements at discrete locations on the Orbiter wind side. Most of these measurements were made downstream of a surface protuberance designed to trip the boundary layer to turbulent flow. In this paper, we use the US3D computational fluid dynamics code to simulate the Orbiter flow field at conditions corresponding to the STS-119 re-entry. We employ a standard two-temperature, five-species finite-rate model for high-temperature air, and the surface catalysis model of Stewart.1 This work is similar to the analysis of Wood et al . 2 except that we use a different approach for modeling turbulent flow. We use the one-equation Spalart-Allmaras turbulence model8 with compressibility corrections 9 and an approach for tripping the boundary layer at discrete locations. In general, the comparison between the simulations and flight data is remarkably good

Candler, Graham V.↗

Visualization of particle flux in the human body on the surface of Mars

For a given galactic cosmic ray (GCR) environment, information on the particle flux of protons, alpha particles, and heavy ions, that varies with respect to the topographical altitude on the Martian surface, are needed for planning exploration missions to Mars. The Mars Global Surveyor (MGS) mission with its Mars Orbiter Laser Altimeter (MOLA) instrument has been providing precise topographical surface map of the Mars. With this topographical data, the particle flux at the Martian surface level through the CO2 atmospheric shielding for solar minimum and solar maximum conditions are calculated. These particle flux calculations are then transported first through an anticipated shielding of a conceptual shelter with several water equivalent shield values (up to 50 g/cm2 of water in steps of 5 g/cm2) considered to represent a surface habitat, and then into the human body. Model calculations are accomplished utilizing the HZETRN, QMSFRG, and SUM-MARS codes. Particle flux calculations for 12 different locations in the human body were considered from skin depth to the internal organs including the blood-forming organs (BFO). Visualization of particle flux in the human body at different altitudes on the Martian surface behind a known shielding is anticipated to provide guidance for assessing radiation environment risk on the Martian surface for future human missions.

NASA Center JSC↗

A three-dimensional spacecraft-charging computer code

A computer code is described which simulates the interaction of the space environment with a satellite at geosynchronous altitude. Employing finite elements, a three-dimensional satellite model has been constructed with more than 1000 surface cells and 15 different surface materials. Free space around the satellite is modeled by nesting grids within grids. Applications of this NASA Spacecraft Charging Analyzer Program (NASCAP) code to the study of a satellite photosheath and the differential charging of the SCATHA (satellite charging at high altitudes) satellite in eclipse and in sunlight are discussed. In order to understand detector response when the satellite is charged, the code is used to trace the trajectories of particles reaching the SCATHA detectors. Particle trajectories from positive and negative emitters on SCATHA also are traced to determine the location of returning particles, to estimate the escaping flux, and to simulate active control of satellite potentials.

Rubin, A. G.↗

Numerical simulations of fuel droplet flows using a Lagrangian triangular mesh

The incompressible, Lagrangian, triangular grid code, SPLISH, was converted for the study of flows in and around fuel droplets. This involved developing, testing and incorporating algorithms for surface tension and viscosity. The major features of the Lagrangian method and the algorithms are described. Benchmarks of the algorithms are given. Several calculations are presented for kerosene droplets in air. Finally, extensions which make the code compressible and three dimensional are discussed.

Fritts, M. J.↗

Solar Thermal Concept Evaluation

Concentrated solar thermal energy can be utilized in a variety of high temperature applications for both terrestrial and space environments. In each application, knowledge of the collector and absorber's heat exchange interaction is required. To understand this coupled mechanism, various concentrator types and geometries, as well as, their relationship to the physical absorber mechanics were investigated. To conduct experimental tests various parts of a 5,000 watt, thermal concentrator, facility were made and evaluated. This was in anticipation at a larger NASA facility proposed for construction. Although much of the work centered on solar thermal propulsion for an upper stage (less than one pound thrust range), the information generated and the facility's capabilities are applicable to material processing, power generation and similar uses. The numerical calculations used to design the laboratory mirror and the procedure for evaluating other solar collectors are presented here. The mirror design is based on a hexagonal faceted system, which uses a spherical approximation to the parabolic surface. The work began with a few two dimensional estimates and continued with a full, three dimensional, numerical algorithm written in FORTRAN code. This was compared to a full geometry, ray trace program, BEAM 4, which optimizes the curvatures, based on purely optical considerations. Founded on numerical results, the characteristics of a faceted concentrator were construed. The numerical methodologies themselves were evaluated and categorized. As a result, the three-dimensional FORTRAN code was the method chosen to construct the mirrors, due to its overall accuracy and superior results to the ray trace program. This information is being used to fabricate and subsequently, laser map the actual mirror surfaces. Evaluation of concentrator mirrors, thermal applications and scaling the results of the 10 foot diameter mirror to a much larger concentrator, were studied. Evaluations, recommendations and pit falls regarding the structure, materials and facility design are presented.

Hawk, Clark W.↗

Common radiation analysis model for 75,000 pound thrust NERVA engine (1137400E)

The mathematical model and sources of radiation used for the radiation analysis and shielding activities in support of the design of the 1137400E version of the 75,000 lbs thrust NERVA engine are presented. The nuclear subsystem (NSS) and non-nuclear components are discussed. The geometrical model for the NSS is two dimensional as required for the DOT discrete ordinates computer code or for an azimuthally symetrical three dimensional Point Kernel or Monte Carlo code. The geometrical model for the non-nuclear components is three dimensional in the FASTER geometry format. This geometry routine is inherent in the ANSC versions of the QAD and GGG Point Kernal programs and the COHORT Monte Carlo program. Data are included pertaining to a pressure vessel surface radiation source data tape which has been used as the basis for starting ANSC analyses with the DASH code to bridge into the COHORT Monte Carlo code using the WANL supplied DOT angular flux leakage data. In addition to the model descriptions and sources of radiation, the methods of analyses are briefly described.

Warman, E. A.↗

The implementation of the graphics of program EAGLE: A numerical grid generation code on NASA Langley SNS computer system

Program EAGLE (Eglin Arbitrary Geometry Implicit Euler) Numerical Grid Generation System is a composite (multi-block) algebraic or elliptic grid generation system designed to discretize the domain in and/or around any arbitrarily shaped three dimensional regions. This system combines a boundary conforming surface generation scheme and includes plotting routines designed to take full advantage of the DISSPLA Graphics Package (Version 9.0). Program EAGLE is written to compile and execute efficiently on any Cray machine with or without solid state disk (SSD) devices. Also, the code uses namelist inputs which are supported by all Cray machines using the FORTRAN compiler CFT77. The namelist inputs makes it easier for the user to understand the inputs and operation of Program EAGLE. EAGLE's numerical grid generator is constructed in the following form: main program, EGG (executive routine); subroutine SURFAC (surface generation routine); subroutine GRID (grid generation routine); and subroutine GRDPLOT (grid plotting routines). The EAGLE code was modified to use on the NASA-LaRC SNS computer (Cray 2S) system. During the modification a conversion program was developed for the output data of EAGLE's subroutine GRID to permit the data to be graphically displayed by IRIS workstations, using Plot3D. The code of program EAGLE was modified to make operational subroutine GRDPLOT (using DI-3000 Graphics Software Packages) on the NASA-LaRC SNS Computer System. How to implement graphically, the output data of subroutine GRID was determined on any NASA-LaRC graphics terminal that has access to the SNS Computer System DI-300 Graphics Software Packages. A Quick Reference User Guide was developed for the use of program EAGLE on the NASA-LaRC SNS Computer System. One or more application program(s) was illustrated using program EAGLE on the NASA LaRC SNS Computer System, with emphasis on graphics illustrations.

Houston, Johnny L.↗

Developing New Tools for Modeling Rocket Plume-Surface Interactions

With NASA’s goal to land the next human on the lunar surface in the next few years, it has become vitally important that we have a better understanding of how future landing spacecraft will interact with the unique properties of regolith¬¬––the layer of loose, unconsolidated dust and rock on the lunar surface¬¬––which can cause hazards like visual obstructions, particulate clouds, and cratering of the landing zone. Researchers from the Fluid Dynamics Branch at NASA’s Marshall Space Flight Center are performing plume-surface interaction (PSI) simulations between lander engine plumes and unprepared regolith surfaces, and have developed new tools to provide predictive PSI environments for various NASA projects and missions, including the Human Lander System (HLS), Commercial Lunar Payload Services (CLPS), and future Mars landers. These tools allow the researchers to determine how to best meet the simulation and time requirements for each project by varying model fidelity. The highest fidelity tool is the Gas Granular Flow Solver (Loci/GGFS) that models gas-particle multi-phase interactions to predict regolith cratering and ejection of particles into the immediate surroundings of the lander. At its highest fidelity, it can model microscopic regolith particle interactions with a particle size/shape distribution that statistically replicates actual regolith, however, to be most effective with today’s computing resources, it is currently run using only one to three equivalent particle sizes/shapes. The team also incorporated engineering models into their software suite to create production-ready hybrid tools with reduced fidelity. At the lowest fidelity, the computational fluid dynamics (CFD) code Loci/CHEM+DIGGEM can predict crater depth over time by relating local CFD-predicted surface shear stresses to a model of erosion mass flux.

plume surface interaction↗

Effects of emitting surfaces and trapped ions on the sheath physics and current flow in multidimensional plasma systems

Recent one-dimensional simulations of planar sheaths with strong electron emission have shown that trapping of charge-exchange ions causes transitions from space-charge limited (SCL) to inverse sheaths. However, multidimensional emitting sheath phenomena with collisions remained unexplored, due in part to high computational cost. As such, we developed a novel continuum kinetic code to study the sheath physics, current flow and potential distributions in two-dimensional unmagnetized configurations with emitting surfaces. For small negatively biased thermionic cathodes in a plasma, the cathode sheath can exist in an equilibrium SCL state. The SCL sheath carries an immense density of trapped ions, neutralized by thermoelectrons, within the potential well of the virtual cathode. For further increases of emitted flux, the trapped ion cloud expands in space. The trapped ion space charge causes an increase of thermionic current far beyond the saturation limit predicted by conventional collisionless SCL sheath models without ion trapping. For sufficiently strong emission, the trapped ion cloud consumes the entire 2D plasma domain, forming a mode with globally confined ions and an inverse sheath at the cathode. In situations where the emitted flux is fixed and the bias is swept (e.g. emissive probe), the trapped ions cause a large thermionic current to escape for all biases below the plasma potential. Strong suppression of the thermionic emission, required for the probe to float, only occurs when the probe is above the plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Contributions to HiLiftPW-3 Using Structured, Overset Grid Methods

The High-Lift Common Research Model (HL-CRM) and the JAXA Standard Model (JSM) were analyzed computationally using both the OVERFLOW and LAVA codes for the third AIAA High-Lift Prediction Workshop. Geometry descriptions and the test cases simulated are described. With the HL-CRM, the effects of surface smoothness during grid projection and the effect of partially sealing a flap gap were studied. Grid refinement studies were performed at two angles of attack using both codes. For the JSM, simulations were performed with and without the nacelle/pylon. Without the nacelle/pylon, evidence of multiple solutions was observed when a quadratic constitutive relation is used in the turbulence modeling; however, using time-accurate simulation seemed to alleviate this issue. With the nacelle/pylon, no evidence of multiple solutions was observed. Laminar-turbulent transition modeling was applied to both JSM configuration, and had an overall favorable impact on the lift predictions.

HiLiftPW-3↗

Influence of Ply Waviness on Fatigue Life of Tapered Composite Flexbeam Laminates

Nonlinear tapered flexbeam laminates, with significant ply waviness, were cut from a full-size composite rotor hub flexbeam. The specimens were tested under combined axial tension and cyclic bending loads. All of the specimens had wavy plies through the center and near the surfaces (termed marcelled areas), although for some of the specimens the surface marcels were very obvious, and for others they were much smaller. The specimens failed by first developing cracks through the marcels at the surfaces, and then delaminations grew from those cracks, in both directions. Delamination failure occurred in these specimens at significantly shorter fatigue lives than similar specimens without waviness, tested in ref. 2. A 2D finite element model was developed which closely approximated the flexbeam geometry, boundary conditions, and loading. In addition, the FE model duplicated the waviness observed in one of the test specimens. The model was analyzed using a geometrically nonlinear FE code. Modifications were made to the original model to reduce the amplitude of the marcels near the surfaces. The analysis was repeated for each modification. Comparisons of the interlaminar normal stresses, sigma(sub n), in the various models showed that under combined axial-tension and cyclic-bending loading, for marcels of the same aspect ratio, sigma(sub n) stresses increased as the distance along the taper, from thick to thin end, increased. For marcels of the same aspect ratio and at the same X-location along the taper, sigma(sub n) stresses decreased as the distance from the surface into the flexbeam interior increased. A technique was presented for determining the smallest acceptable marcel aspect ratio at various locations in the flexbeam.

Murri, Gretchen B.↗

Prediction of Unshsrouded Rotor Blade Tip Heat Transfer

The rate of heat transfer on the tip of a turbine rotor blade and on the blade surface in the vicinity of the tip, was successfully predicted. The computations were performed with a multiblock computer code which solves the Reynolds Averaged Navier-Stokes equations using an efficient multigrid method. The case considered for the present calculations was the Space Shuttle Main Engine (SSME) high pressure fuel side turbine. The predictions of the blade tip heat transfer agreed reasonably well with the experimental measurements using the present level of grid refinement. On the tip surface, regions with high rate of heat transfer was found to exist close to the pressure side and suction side edges. Enhancement of the heat transfer was also observed on the blade surface near the tip. Further comparison of the predictions was performed with results obtained from correlations based on fully developed channel flow.

Ameri, A. A.↗

Estimation of Surface Temperature and Heat Flux by Inverse Heat Transfer Methods Using Internal Temperatures Measured While Radiantly Heating a Carbon/Carbon Specimen up to 1920 F

The ability to solve the heat conduction equation is needed when designing materials to be used on vehicles exposed to extremely high temperatures; e.g. vehicles used for atmospheric entry or hypersonic flight. When using test and flight data, computational methods such as finite difference schemes may be used to solve for both the direct heat conduction problem, i.e., solving between internal temperature measurements, and the inverse heat conduction problem, i.e., using the direct solution to march forward in space to the surface of the material to estimate both surface temperature and heat flux. The completed research first discusses the methods used in developing a computational code to solve both the direct and inverse heat transfer problems using one dimensional, centered, implicit finite volume schemes and one dimensional, centered, explicit space marching techniques. The developed code assumed the boundary conditions to be specified time varying temperatures and also considered temperature dependent thermal properties. The completed research then discusses the results of analyzing temperature data measured while radiantly heating a carbon/carbon specimen up to 1920 F. The temperature was measured using thermocouple (TC) plugs (small carbon/carbon material specimens) with four embedded TC plugs inserted into the larger carbon/carbon specimen. The purpose of analyzing the test data was to estimate the surface heat flux and temperature values from the internal temperature measurements using direct and inverse heat transfer methods, thus aiding in the thermal and structural design and analysis of high temperature vehicles.

Pizzo, Michelle↗

Enhancements to the GRIDGEN structured grid generation system for internal and external flow applications

GRIDGEN is a government domain software package for interactive generation of multiple block grids around general configurations. Though it has been freely available since 1989, it has not been widely embraced by the internal flow community due to a misconception that it was designed for external flow use only. In reality GRIDGEN has always worked for internal flow applications, and GRIDGEN ongoing enhancements are increasing the quality of and efficiency with which grids for external and internal flow problems may be constructed. The software consists of four codes used to perform the four steps of the grid generation process. GRIDBLOCK is first used to decompose the flow domain into a collection of component blocks and then to establish interblock connections and flow solver boundary conditions. GRIDGEN2D is then used to generate surface grids on the outer shell of each component block. GRIDGEN3D generates grid points on the interior of each block, and finally GRIDVUE3D is used to inspect the resulting multiple block grid. Three of these codes (GRIDBLOCK, GRIDGEN2D, and GRIDVUE3D) are highly interactive and graphical in nature, and currently run on Silicon Graphics, Inc., and IBM RS/6000 workstations. The lone batch code (GRIDGEN3D) may be run on any of several Unix based platforms. Surface grid generation in GRIDGEN2D is being improved with the addition of higher order surface definitions (NURBS and parametric surfaces input in IGES format and bicubic surfaces input in PATRAN Neutral File format) and double precision mathematics. In addition, two types of automation have been added to GRIDGEN2D that reduce the learning curve slope for new users and eliminate work for experienced users. Volume grid generation using GRIDGEN3D has been improved via the addition of an advanced hybrid control function formulation that provides both orthogonality and clustering control at the block faces and clustering control on the block interior.

Steinbrenner, John P.↗

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implementation of Combined Feather and Surface-Normal Ice Growth Models in LEWICE/X

Experimental observations have shown that discrete rime ice growths called feathers, which grow in approximately the direction of water droplet impingement, play an important role in the growth of ice on accreting surfaces for some thermodynamic conditions. An improved physical model of ice accretion has been implemented in the LEWICE 2D panel-based ice accretion code maintained by the NASA Lewis Research Center. The LEWICE/X model of ice accretion explicitly simulates regions of feather growth within the framework of the LEWICE model. Water droplets impinging on an accreting surface are withheld from the normal LEWICE mass/energy balance and handled in a separate routine; ice growth resulting from these droplets is performed with enhanced convective heat transfer approximately along droplet impingement directions. An independent underlying ice shape is grown along surface normals using the unmodified LEWICE method. The resulting dual-surface ice shape models roughness-induced feather growth observed in icing wind tunnel tests. Experiments indicate that the exact direction of feather growth is dependent on external conditions. Data is presented to support a linear variation of growth direction with temperature and cloud water content. Test runs of LEWICE/X indicate that the sizes of surface regions containing feathers are influenced by initial roughness element height. This suggests that a previous argument that feather region size is determined by boundary layer transition may be incorrect. Simulation results for two typical test cases give improved shape agreement over unmodified LEWICE.

Matthew T Velazquez↗

High-Lift Propeller Noise Prediction for a Distributed Electric Propulsion Flight Demonstrator

Over the past several years, the use of electric propulsion technologies within aircraft design has received increased attention. The characteristics of electric propulsion systems open up new areas of the aircraft design space, such as the use of distributed electric propulsion (DEP). In this approach, electric motors are placed in many different locations to achieve increased efficiency through integration of the propulsion system with the airframe. Under a project called Scalable Convergent Electric Propulsion Technology Operations Research (SCEPTOR), NASA is designing a flight demonstrator aircraft that employs many "high-lift propellers" distributed upstream of the wing leading edge and two cruise propellers (one at each wingtip). As the high-lift propellers are operational at low flight speeds (take-off/approach flight conditions), the impact of the DEP configuration on the aircraft noise signature is also an important design consideration. This paper describes efforts toward the development of a mulit-fidelity aerodynamic and acoustic methodology for DEP high-lift propeller aeroacoustic modeling. Specifically, the PAS, OVERFLOW 2, and FUN3D codes are used to predict the aerodynamic performance of a baseline high-lift propeller blade set. Blade surface pressure results from the aerodynamic predictions are then used with PSU-WOPWOP and the F1A module of the NASA second generation Aircraft NOise Prediction Program to predict the isolated high-lift propeller noise source. Comparisons of predictions indicate that general trends related to angle of attack effects at the blade passage frequency are captured well with the various codes. Results for higher harmonics of the blade passage frequency appear consistent for the CFD based methods. Conversely, evidence of the need for a study of the effects of increased azimuthal grid resolution on the PAS based results is indicated and will be pursued in future work. Overall, the results indicate that the computational approach is acceptable for fundamental assessment of low-noise high-lift propeller designs. The extent to which the various approaches may be used in a complementary manner will be further established as measured data becomes available for validation. Ultimately, it is anticipated that this combined approach may be used to provide realistic incident source fields for acoustic shielding/scattering studies on various aircraft configurations.

Nark, Douglas M.↗