Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model code”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Computational Work to Support FAP/SRW Variable-Speed Power-Turbine Development

The purpose of this report is to document the work done to enable a NASA CFD code to model the transition on a blade. The purpose of the present work is to down-select a transition model that would allow the flow simulation of a Variable-Speed Power-Turbine (VSPT) to be accurately performed. The modeling is to be ultimately performed to also account for the blade row interactions and effect on transition and therefore accurate accounting for losses. The present work is limited to steady flows. The low Reynolds number k-omega model of Wilcox and a modified version of same will be used for modeling of transition on experimentally measured blade pressure and heat transfer. It will be shown that the k-omega model and its modified variant fail to simulate the transition with any degree of accuracy. A case is therefore made for more accurate transition models. Three-equation models based on the work of Mayle on Laminar Kinetic Energy were explored and the Walters and Leylek model which was thought to be in a more mature state of development is introduced and implemented in the Glenn-HT code. Two-dimensional flat plate results and three-dimensional results for flow over turbine blades and the resulting heat transfer and its transitional behavior are reported. It is shown that the transition simulation is much improved over the baseline k-omega model.

Ameri, Ali A.↗
NASA NTRS1986

Hydroburst test of a carbon-carbon involute exit cone

A hydroburst test of the aft portion of the PAM-D exit cone and the test procedure are described in detail. The hydrostatic pressure required to buckle the cone was 9.75 psi. Meanwhile, the PAM-D exit cone was modeled using the finite element method and a theoretical bucking pressure (8.76 psi) was predicted using the SPAR finite element code. The modeling technique employed is discussed. By comparing the theoretical to predicted critical pressures, this report verifies the modeling technique and calculates a material knockdown factor for the carbon-carbon exit cone.

Sullivan, Roy M.↗
NASA NTRS1993

Advanced turbulence models for turbomachinery

Development and assessment of the single-time-scale k-epsilon turbulence model with different near-wall treatments and the multi-scale turbulence model for rotating flows are presented. These turbulence models are coded as self contained module decks that can be interfaced with a number of CFD main flow solvers. For each model, a stand-alone module deck with its own formulation, discretization scheme, solver and boundary condition implementations is presented. These satellite decks will take as input (from a main flow solver) the velocity field, grid, boundary condition specifications and will deliver turbulent quantities as output. These modules were tested as separate entities and, although many logical and programming problems were overcome, only wider use and further testing can render the modules sufficiently 'fool proof'.

Hadid, Ali H.↗
NASA NTRS2010

The OMPS Limb Profiler Instrument: Two-Dimensional Retrieval Algorithm

The upcoming Ozone Mapper and Profiler Suite (OMPS), which will be launched on the NPOESS Preparatory Project (NPP) platform in early 2011, will continue monitoring the global distribution of the Earth's middle atmosphere ozone and aerosol. OMPS is composed of three instruments, namely the Total Column Mapper (heritage: TOMS, OMI), the Nadir Profiler (heritage: SBUV) and the Limb Profiler (heritage: SOLSE/LORE, OSIRIS, SCIAMACHY, SAGE III). The ultimate goal of the mission is to better understand and quantify the rate of stratospheric ozone recovery. The focus of the paper will be on the Limb Profiler (LP) instrument. The LP instrument will measure the Earth's limb radiance (which is due to the scattering of solar photons by air molecules, aerosol and Earth surface) in the ultra-violet (UV), visible and near infrared, from 285 to 1000 nm. The LP simultaneously images the whole vertical extent of the Earth's limb through three vertical slits, each covering a vertical tangent height range of 100 km and each horizontally spaced by 250 km in the cross-track direction. Measurements are made every 19 seconds along the orbit track, which corresponds to a distance of about 150km. Several data analysis tools are presently being constructed and tested to retrieve ozone and aerosol vertical distribution from limb radiance measurements. The primary NASA algorithm is based on earlier algorithms developed for the SOLSE/LORE and SAGE III limb scatter missions. All the existing retrieval algorithms rely on a spherical symmetry assumption for the atmosphere structure. While this assumption is reasonable in most of the stratosphere, it is no longer valid in regions of prime scientific interest, such as polar vortex and UTLS regions. The paper will describe a two-dimensional retrieval algorithm whereby the ozone distribution is simultaneously retrieved vertically and horizontally for a whole orbit. The retrieval code relies on (1) a forward 2D Radiative Transfer code (to model limb radiances within a non-uniform atmosphere and evaluate 2D analytical partial derivatives) and (2) an optimal estimator inversion routine. The algorithm uses the typically sparse nature of the kernel matrices as well as fast matrix inversion techniques to allow for fast inversion of limb data with efficient memory management (as was done for MIPAS data processing). While the method has so far only been developed in the context of Single Scatter, the paper will show how the CPU intensive Multiple Scatter modeling can be implemented using parallel CPU processing. Initial results will be presented in terms of retrieved ozone profiles and code performance.

Rault, Didier F.↗
NASA NTRS2000

Solar Power System Evaluated for the Human Exploration of Mars

The electric power system is a crucial element of any mission for the human exploration of the Martian surface. The bulk of the power generated will be delivered to crew life support systems, extravehicular activity suits, robotic vehicles, and predeployed in situ resource utilization (ISRU) equipment. In one mission scenario, before the crew departs for Mars, the ISRU plant operates for 435 days producing liquefied methane and oxygen for ascent-stage propellants and water for crew life support. About 200 days after ISRU production is completed, the crew arrives for a 500-day surface stay. In this scenario, the power system must operate for a total of 1130 days (equivalent to 1100 Martian "sols"), providing 400 MW-hr of energy to the ISRU plant and up to 18 kW of daytime user power. A photovoltaic power-generation system with regenerative fuel cell (RFC) energy storage has been under study at the NASA Glenn Research Center at Lewis Field. The conceptual power system is dominated by the 4000- m2 class photovoltaic array that is deployed orthogonally as four tent structures, each approximately 5 m on a side and 100-m long. The structures are composed of composite members deployed by an articulating mast, an inflatable boom, or rover vehicles, and are subsequently anchored to the ground. Array panels consist of thin polymer membranes with thin-film solar cells. The array is divided into eight independent electrical sections with solar cell strings operating at 600 V. Energy storage is provided by regenerative fuel cells based on hydrogen-oxygen proton exchange membrane technology. Hydrogen and oxygen reactants are stored in gaseous form at 3000 psi, and the water produced is stored at 14.7 psi. The fuel cell operating temperature is maintained by a 40-m2 deployable pumped-fluid loop radiator that uses water as the working fluid. The power management and distribution (PMAD) architecture features eight independent, regulated 600-Vdc channels. Power management and distribution power cables use various gauges of copper conductors with ethylene tetrafluoroethylene insulation. To assess power system design options and sizing, we developed a dedicated Fortran code to predict detailed power system performance and estimate system mass. This code also modeled the requisite Mars surface environments: solar insolation, Sun angles, dust storms, dust deposition, and thermal and ultraviolet radiation. Using this code, trade studies were performed to assess performance and mass sensitivities to power system design parameters (photovoltaic array geometry and orientation) and mission parameters (landing date and landing site latitude, terrain slope, and dust storm activity). Mission analysis cases were also run. Power results are shown in this graph for an analysis case with a September 1, 2012, landing date; 18.95 North latitude landing site; two seasonal dusts storms; and tent arrays. To meet user load requirements and the ISRU energy requirement, an 8-metric ton (MT) power system and 4000-m2 photovoltaic array area were required for the assumed advanced CuInS2 thin-film solar cell technology. In this figure, the top curve is the average daytime photovoltaic array power, the middle curve is average daytime user load power, and the bottom curve is nighttime power. At mission day 1, daytime user power exceeds 120 kW before falling off to 80 kW at the end of the mission. Throughout the mission, nighttime user power is set to the nighttime power requirement. In this analysis, "nighttime" is defined as the 13- to 15-hr period when array power output is below the daytime power requirement. During dust storms, power system capability falls off dramatically so that by mission day 900, a daily energy balance cannot be maintained. Under these conditions, the ISRU plant is placed in standby mode, and the regenerative fuel cell energy storage is gradually discharged to meet user loads.

Kerslake, Thomas W.↗
NASA NTRS2014

A System of Conservative Regridding for Ice-Atmosphere Coupling in a General Circulation Model (GCM)

The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere-ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.

Ice-atmosphere↗
NASA NTRS2018

Finite Element Simulations of Two Vertical Drop Tests of F-28 Fuselage Sections

In March 2017, a vertical drop test of a forward fuselage section of a Fokker F-28 MK4000 aircraft was conducted as part of a joint NASA/FAA project to investigate the performance of transport aircraft under realistic crash conditions. In June 2017, a vertical drop test was conducted of a wing-box fuselage section of the same aircraft. Both sections were configured with two rows of aircraft seats, in a triple-double configuration. A total of ten Anthropomorphic Test Devices (ATDs) were secured in seats using standard lap belt restraints. The forward fuselage section was also configured with luggage in the cargo hold. Both sections were outfitted with two hat racks, each with added ballast mass. The drop tests were performed at the Landing and Impact Research facility located at NASA Langley Research Center in Hampton, Virginia. The measured impact velocity for the forward fuselage section was 346.8-in/s onto soil. The wing-box section was dropped with a downward facing pitch angle onto a sloping soil surface in order to create an induced forward acceleration in the airframe. The vertical impact velocity of the wing-box section was 349.2-in/s. A second objective of this project was to assess the capabilities of finite element simulations to predict the test responses. Finite element models of both fuselage sections were developed for execution in LS-DYNA(Registered Trademark), a commercial explicit nonlinear transient dynamic code. The models contained accurate representations of the airframe structure, the hat racks and hat rack masses, the floor and seat tracks, the luggage in the cargo hold for the forward section, and the detailed under-floor structure in the wing-box section. Initially, concentrated masses were used to represent the inertial properties of the seats, restraints, and ATD occupants. However, later simulations were performed that included finite element representations of the seats, restraints, and ATD occupants. These models were developed to more accurately replicate the seat loading of the floor and to enable prediction of occupant impact responses. Models were executed to generate analytical predictions of airframe responses, which were compared with test data to validate the model. Comparisons of predicted and experimental structural deformation and failures were made. Finally, predicted and experimental soil deformation and crater depths were also compared for both drop test configurations.

Jackson, Karen E.↗
NASA NTRS2005

NESTEM Probabilistic Analysis Used to Study Mistuned Bladed Disks and Blisks With Aerodynamic and Structural Coupling

This document summarizes the initial results from a research effort at the NASA Glenn Research Center on blisk and bladed disk mistuning, including both structural and aerodynamic coupling. The structural coupling model is based on the Fundamental Mistuning Model (FMM, developed by Feiner and Griffin). This effort extends the FMM technique to accept aerodynamic coupling coefficients from computational fluid dynamic codes. The model was applied to a representative modern front compressor stage. Flutter stability and forced response were determined with structural coupling only, with aerodynamic coupling only, and with both structural and aerodynamic coupling. Tuned, randomly mistuned, and near alternately mistuned rotors were considered.

Kielb, Robert E.↗
NASA NTRS1987

Three-dimensional computations of transverse hydrogen jet combustion in a supersonic airstream

A computational fluid dynamics (CFD) code is being developed to compute the mixing and combustion of hydrogen fuel in the turbulent flow fields of supersonic combustion ramjets (scramjet). The code solves the three-dimensional Reynolds time-averaged complete Navier-Stokes equations including transport equations for a four species, two reaction, global finite rate chemistry model. The code was applied to the case of transverse injection of hydrogen from a sonic circular orifice into a supersonic airstream. The equations were numerically integrated using MacCormack's explicit method, and the algebraic eddy viscosity model of Baldwin-Lomax was used to model the turbulence. In the species transport and energy equations, diffusion coefficients based on Fick's Law and an assumption of unit Lewis number were applied. Computed features of the three-dimensional flow field are depicted by static pressure, static temperature, mass fraction of species, and velocity vectors. For engineering interest, mixing and combustion parameters were examined to assess the effect of injector diameter, injected fuel pressure, fuel-air ratio, and spacing of fuel injectors. The objective of the present paper is to demonstrate the capability of the present three-dimensional spatially elliptic, CFD code for turbulent, reacting flow. Application of the code to specific supersonic combustion configurations is planned.

Uenishi, K.↗
NASA NTRS1991

Low altitude plume impingement handbook

Plume Impingement modeling is required whenever an object immersed in a rocket exhaust plume must survive or remain undamaged within specified limits, due to thermal and pressure environments induced by the plume. At high altitudes inviscid plume models, Monte Carlo techniques along with the Plume Impingement Program can be used to predict reasonably accurate environments since there are usually no strong flowfield/body interactions or atmospheric effects. However, at low altitudes there is plume-atmospheric mixing and potential large flowfield perturbations due to plume-structure interaction. If the impinged surface is large relative to the flowfield and the flowfield is supersonic, the shock near the surface can stand off the surface several exit radii. This results in an effective total pressure that is higher than that which exists in the free plume at the surface. Additionally, in two phase plumes, there can be strong particle-gas interaction in the flowfield immediately ahead of the surface. To date there have been three levels of sophistication that have been used for low altitude plume induced environment predictions. Level 1 calculations rely on empirical characterizations of the flowfield and relatively simple impingement modeling. An example of this technique is described by Piesik. A Level 2 approach consists of characterizing the viscous plume using the SPF/2 code or RAMP2/LAMP and using the Plume Impingement Program to predict the environments. A Level 3 analysis would consist of using a Navier-Stokes code such as the FDNS code to model the flowfield and structure during a single calculation. To date, Level 1 and Level 2 type analyses have been primarily used to perform environment calculations. The recent advances in CFD modeling and computer resources allow Level 2 type analysis to be used for final design studies. Following some background on low altitude impingement, Level 1, 2, and 3 type analysis will be described.

Smith, Sheldon D.↗
NASA NTRS1993

Al-26 production profile and model comparisons in Canyon Diablo

The large preatmospheric size of the Canyon Diablo meteorite, a radius of about 15 m, makes it especially suitable for systematic studies of cosmogenic nuclide production rates of iron objects in a 2 pi geometry. To reconstruct the exposure history of the meteoroid, Heymann et al. investigated several fragments recovered from known geographic locations around the crater for their shock features and cosmogenic nobel gases. They applied the Signer-Nier noble gas production rate model to establish the preatmospheric depth of the specimens in the meteoroid. Cosmic ray exposure ages suggested a multi-episodic irradiation, with 170 or 540 Ma being inferred for most of the samples studied while two anomalous specimens indicated a possible third exposure age at 940 Ma. Be-10 and Cl-36 have been measured in a number of these same samples by accelerator mass spectrometry (AMS), with use being made of the preatmospheric depths determined in Heymann et al. to construct production profiles. The present study extends the cosmogenic radionuclide data to Al-26 and compares the results with both the production rate model of Reedy and Arnold and production rates determined from the cross sections used by the Reedy-Arnold model (for the major nuclear reactions making Al-26) in combination with differential fluxes calculated using the Los Alamos High Energy Transport (LAHET) Code System. Model calculations for Be-10 and Cl-36 have also been obtained.

Michlovich, E.↗
NASA NTRS2013

Modeling Analysis for NASA GRC Vacuum Facility 5 Upgrade

A model of the VF5 test facility at NASA Glenn Research Center was developed using the direct simulation Monte Carlo Hypersonic Aerothermodynamics Particle (HAP) code. The model results were compared to several cold flow and thruster hot fire cases. The main uncertainty in the model is the determination of the effective sticking coefficient -- which sets the pumping effectiveness of the cryopanels and oil diffusion pumps including baffle transmission. An effective sticking coefficient of 0.25 was found to provide generally good agreement with the experimental chamber pressure data. The model, which assumes a cold diffuse inflow, also fared satisfactorily in predicting the pressure distribution during thruster operation. The model was used to assess other chamber configurations to improve the local effective pumping speed near the thruster. A new configuration of the existing cryopumps is found to show more than 2x improvement over the current baseline configuration.

cryopumping↗
NASA NTRS1994

Analysis of subgrid models using direct and large-eddy simulations of isotropic turbulence

Direct and large eddy simulations of forced and decaying isotropic turbulence have been performed using a pseudospectral and a finite-difference code. Subgrid models that include a one-equation subgrid kinetic energy model with and without a stochastic backscatter forcing term and a new scale similarity model have been analyzed in both Fourier space and physical space. The Fourier space analysis showed that the energy transfer across the cutoff wavenumber k(sub c) is dominated by local interaction. The correlation between the exact and the modeled (by a spectral eddy viscosity) nonlinear terms and the subgrid energy transfer in physical space was found to be quite low. In physical space, a similar correlation analysis was carried out using top hat filtering. Results show that the subgrid stress and the energy flux predicted by the subgrid models correlates very well with the exact data. The scale similarity model showed very high correlation for reasonable grid resolution. However, with decrease in grid resolution, the scale similarity model became more uncorrelated, when compared to the kinetic energy subgrid model. The subgrid models were then used for large-eddy simulations for a range of Reynolds number. It was determined that the dissipation was modeled poorly and that the correlation with the exact results was quite low for all the models. In general, for coarse grid resolution, the scale similarity model consistently showed very low correlation while the kinetic energy model showed a relatively higher correlation. These results suggest that to use the scale similarity model relatively fine grid resolution may be required, whereas, the kinetic energy model could be used even in coarse grid.

Menon, S.↗
NASA NTRS1983

A one-dimensional unsteady model of dual mode scramjet operation

A one-dimensional unsteady model and code have been developed which can handle internal flows involving area change, heat release, transients and coupled inviscid core-boundary layer flows. The model has been shown to work well in diverging channels with a supersonic entrance Mach number and unseparated boundary layers. Transient flows have also been studied. The model needs no prior knowledge of shock (or shocks) location and can be extended to two or three dimensions. The results show that the one-dimensional unsteady model developed here can be used for analysis of internal flows as occur in scramjet combustors.

Bussing, T. R. A.↗
NASA NTRS1988

Automatic mathematical modeling for real time simulation system

A methodology for automatic mathematical modeling and generating simulation models is described. The models will be verified by running in a test environment using standard profiles with the results compared against known results. The major objective is to create a user friendly environment for engineers to design, maintain, and verify their model and also automatically convert the mathematical model into conventional code for conventional computation. A demonstration program was designed for modeling the Space Shuttle Main Engine Simulation. It is written in LISP and MACSYMA and runs on a Symbolic 3670 Lisp Machine. The program provides a very friendly and well organized environment for engineers to build a knowledge base for base equations and general information. It contains an initial set of component process elements for the Space Shuttle Main Engine Simulation and a questionnaire that allows the engineer to answer a set of questions to specify a particular model. The system is then able to automatically generate the model and FORTRAN code. The future goal which is under construction is to download the FORTRAN code to VAX/VMS system for conventional computation. The SSME mathematical model will be verified in a test environment and the solution compared with the real data profile. The use of artificial intelligence techniques has shown that the process of the simulation modeling can be simplified.

Wang, Caroline↗
NASA NTRS2010

Comparisons of Flutter Analyses for an Experimental Fan

Two propulsion aeroelasticity codes were used to model the aeroelastic characteristics of an experimental forward-swept fan that encountered flutter during wind tunnel testing. Both of these three-dimensional codes model the unsteady flowfield due to blade vibrations using the Navier-Stokes equations. In the first approach, the unsteady flow equations are solved using an implicit time-marching approach. In the second approach, the unsteady flow equations are converted to a harmonic balance form and solved using a pseudo-time marching method. This paper describes the flutter calculations and compares the results to experimental measurements.

Bakhle, Milind A.↗
NASA NTRS2022

Application of Rotor Disk Model in the Overflow CFD Code

To reduce the complexity and to reduce the cost of turnaround time of simulating unsteady flows of rotating blades with moving overset grids, amid-fidelity option is provided in the OVERFLOW computational fluid dynamics analysis to model rotors and propellers as infinitely thin disks. The effect of the flow through the disk is represented by the momentum source terms in the governing Navier-Stokes equations. The model is loosely coupled with the rotorcraft comprehensive analysis CAMRADII. A brief description of the coupling algorithm is provided. Simulations of a single rotor and multiple rotors are shown, including comparisons with measurements and with the full rotor where all blades are modeled, and unsteady simulation performed with moving grids for rotating geometry,

ARMD↗