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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

DAC↗
NASA NTRS2002

High-Fidelity Micromechanics Model Developed for the Response of Multiphase Materials

A new high-fidelity micromechanics model has been developed under funding from the NASA Glenn Research Center for predicting the response of multiphase materials with arbitrary periodic microstructures. The model's analytical framework is based on the homogenization technique, but the method of solution for the local displacement and stress fields borrows concepts previously employed in constructing the higher order theory for functionally graded materials. The resulting closed-form macroscopic and microscopic constitutive equations, valid for both uniaxial and multiaxial loading of periodic materials with elastic and inelastic constitutive phases, can be incorporated into a structural analysis computer code. Consequently, this model now provides an alternative, accurate method.

Aboudi, Jacob↗
NASA NTRS2006

Accuracy Quantification of the Loci-CHEM Code for Chamber Wall Heat Fluxes in a G02/GH2 Single Element Injector Model Problem

A robust rocket engine combustor design and development process must include tools which can accurately predict the multi-dimensional thermal environments imposed on solid surfaces by the hot combustion products. Currently, empirical methods used in the design process are typically one dimensional and do not adequately account for the heat flux rise rate in the near-injector region of the chamber. Computational Fluid Dynamics holds promise to meet the design tool requirement, but requires accuracy quantification, or validation, before it can be confidently applied in the design process. This effort presents the beginning of such a validation process for the Loci- CHEM CPD code. The model problem examined here is a gaseous oxygen (GO2)/gaseous hydrogen (GH2) shear coaxial single element injector operating at a chamber pressure of 5.42 MPa. The GO2/GH2 propellant combination in this geometry represents one the simplest rocket model problems and is thus foundational to subsequent validation efforts for more complex injectors. Multiple steady state solutions have been produced with Loci-CHEM employing different hybrid grids and two-equation turbulence models. Iterative convergence for each solution is demonstrated via mass conservation, flow variable monitoring at discrete flow field locations as a function of solution iteration and overall residual performance. A baseline hybrid grid was used and then locally refined to demonstrate grid convergence. Solutions were also obtained with three variations of the k-omega turbulence model.

West, Jeff↗
NASA NTRS1986

Progress of Stirling cycle analysis and loss mechanism characterization

An assessment of Stirling engine thermodynamic modeling and design codes shows a general deficiency; this deficiency is due to poor understanding of the fluid flow and heat transfer phenomena that occur in the oscillating flow and pressure level environment within the engines. Stirling engine thermodynamic loss mechanisms are listed. Several experimental and computational research efforts now underway to characterize various loss mechanisms are reviewed. The need for additional experimental rigs and rig upgrades is discussed. Recent developments and current efforts in Stirling engine thermodynamic modeling are also reviewed.

Tew, R. C., Jr.↗
NASA NTRS2006

Accuracy Quantification of the Loci-CHEM Code for Chamber Wall Heat Transfer in a GO2/GH2 Single Element Injector Model Problem

A robust rocket engine combustor design and development process must include tools which can accurately predict the multi-dimensional thermal environments imposed on solid surfaces by the hot combustion products. Currently, empirical methods used in the design process are typically one dimensional and do not adequately account for the heat flux rise rate in the near-injector region of the chamber. Computational Fluid Dynamics holds promise to meet the design tool requirement, but requires accuracy quantification, or validation, before it can be confidently applied in the design process. This effort presents the beginning of such a validation process for the Loci-CHEM CFD code. The model problem examined here is a gaseous oxygen (GO2)/gaseous hydrogen (GH2) shear coaxial single element injector operating at a chamber pressure of 5.42 MPa. The GO2/GH2 propellant combination in this geometry represents one the simplest rocket model problems and is thus foundational to subsequent validation efforts for more complex injectors. Multiple steady state solutions have been produced with Loci-CHEM employing different hybrid grids and two-equation turbulence models. Iterative convergence for each solution is demonstrated via mass conservation, flow variable monitoring at discrete flow field locations as a function of solution iteration and overall residual performance. A baseline hybrid was used and then locally refined to demonstrate grid convergence. Solutions were obtained with three variations of the k-omega turbulence model.

West, Jeff↗
NASA NTRS2009

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗
NASA NTRS2002

Numerical Simulation of Flow Field Within Parallel Plate Plastometer

Parallel Plate Plastometer (PPP) is a device commonly used for measuring the viscosity of high polymers at low rates of shear in the range 10 4 to 10 9 poises. This device is being validated for use in measuring the viscosity of liquid glasses at high temperatures having similar ranges for the viscosity values. PPP instrument consists of two similar parallel plates, both in the range of 1 inch in diameter with the upper plate being movable while the lower one is kept stationary. Load is applied to the upper plate by means of a beam connected to shaft attached to the upper plate. The viscosity of the fluid is deduced from measuring the variation of the plate separation, h, as a function of time when a specified fixed load is applied on the beam. Operating plate speeds measured with the PPP is usually in the range of 10 -3 cm/s or lower. The flow field within the PPP can be simulated using the equations of motion of fluid flow for this configuration. With flow speeds in the range quoted above the flow field between the two plates is certainly incompressible and laminar. Such flows can be easily simulated using numerical modeling with computational fluid dynamics (CFD) codes. We present below the mathematical model used to simulate this flow field and also the solutions obtained for the flow using a commercially available finite element CFD code.

Basil N. Antar↗
NASA NTRS2011

A Short Review of Ablative-Material Response Models and Simulation Tools

A review of the governing equations and boundary conditions used to model the response of ablative materials submitted to a high-enthalpy flow is proposed. The heritage of model-development efforts undertaken in the 1960s is extremely clear: the bases of the models used in the community are mathematically equivalent. Most of the material-response codes implement a single model in which the equation parameters may be modified to model different materials or conditions. The level of fidelity of the models implemented in design tools only slightly varies. Research and development codes are generally more advanced but often not as robust. The capabilities of each of these codes are summarized in a color-coded table along with research and development efforts currently in progress.

Lachaud, Jean↗
NASA NTRS2002

Propulsion System Modeling and Simulation

The Aerospace Systems Design Laboratory at the School of Aerospace Engineering in Georgia Institute of Technology has developed a core competency that enables propulsion technology managers to make technology investment decisions substantiated by propulsion and airframe technology system studies. This method assists the designer/manager in selecting appropriate technology concepts while accounting for the presence of risk and uncertainty as well as interactions between disciplines. This capability is incorporated into a single design simulation system that is described in this paper. This propulsion system design environment is created with a commercially available software called iSIGHT, which is a generic computational framework, and with analysis programs for engine cycle, engine flowpath, mission, and economic analyses. iSIGHT is used to integrate these analysis tools within a single computer platform and facilitate information transfer amongst the various codes. The resulting modeling and simulation (M&S) environment in conjunction with the response surface method provides the designer/decision-maker an analytical means to examine the entire design space from either a subsystem and/or system perspective. The results of this paper will enable managers to analytically play what-if games to gain insight in to the benefits (and/or degradation) of changing engine cycle design parameters. Furthermore, the propulsion design space will be explored probabilistically to show the feasibility and viability of the propulsion system integrated with a vehicle.

Tai, Jimmy C. M.↗
NASA NTRS1995

Investigation of infra-red and nonequilibrium air radiation

This report describes progress on the first year of a research program on the infrared radiation of air plasmas conducted in the High Temperature Gasdynamics Laboratory at Stanford University. This program is intended to investigate the masking of infrared signatures by the air plasma formed behind the bow shock of high velocity missiles. To this date, the radiative emission of air plasmas in the infrared has been the object of few experimental investigations, and although several infrared systems are already modeled in radiation codes such as NEQAIR, measurements are required to validate numerical predictions and indicate whether all transitions of importance are accounted for. The present program is motivated by the fact that 9 excited states (A, B, C, D, B', F, H, and H') of NO radiate in the infrared, especially between 1 and 1.5 microns where at least 9 transitions involving can be observed. Because these IR transitions are relatively well separated from each other, excited NO states concentrations can be easily measured, thus providing essential information on excited-state chemistry for use in optical diagnostics or in electronic excitation model validation. Developing accurate collisional-radiative models for these excited NO states is of importance as the UV-VUV transitions of NO (beta, gamma, epsilon, beta prime, gamma prime) produce a major, if not dominant, fraction of the radiation emitted by air plasmas. During the first year of the program, research has focused on the spectral range 1.0 to 1.5 microns, as detailed in Section 2 of this report. The measurements, conducted in a 50 kW radio-frequency inductively coupled plasma torch operating on air at atmospheric pressure, extend previous shock tube investigations by Wray to a wider spectral range (1.0 to 1.5 microns vs 0.9 to 1.2 microns) and higher temperatures (7600 K in the plasma torch versus 6700 K in the shock-tube). These higher temperatures in the present experiment have made it possible to observe high-lying excited NO states that were previously undetectable. These measurements are currently being extended up to 5 microns, with particular attention paid to the rovibronic bands of ground state NO, molecular continua, CO transitions, and other systems of importance. Publications and presentations resulting from or related to this work are cited in Section 3, and Section 4 lists the personnel who contributed to this report.

Kruger, Charles H.↗
NASA NTRS1992

Modeling of solid/porous wall boundary conditions for the validation of computational fluid dynamics codes

A computational study has been undertaken to investigate method of modeling solid and porous wall boundary conditions in computational fluid dynamics (CFD) codes. The procedure utilizes experimental measurements at the walls to develop a flow field solution based on the method of singularities. This flow field solution is then imposed as a boundary condition in a CFD simulation of the internal flow field. The effectiveness of this method in describing the boundary conditions at the wind tunnel walls using only sparse experimental measurements has been investigated. Position and refinement of experimental measurement locations required to describe porous wall boundary conditions has also been considered.

Beutner, Thomas J.↗
NASA NTRS1992

Analysis of film cooling in rocket nozzles

Computational Fluid Dynamics (CFD) programs are customarily used to compute details of a flow field, such as velocity fields or species concentrations. Generally they are not used to determine the resulting conditions at a solid boundary such as wall shear stress or heat flux. However, determination of this information should be within the capability of a CFD code, as the code supposedly contains appropriate models for these wall conditions. Before such predictions from CFD analyses can be accepted, the credibility of the CFD codes upon which they are based must be established. This report details the progress made in constructing a CFD model to predict the heat transfer to the wall in a film cooled rocket nozzle. Specifically, the objective of this work is to use the NASA code FDNS to predict the heat transfer which will occur during the upcoming hot-firing of the Pratt & Whitney 40K subscale nozzle (1Q93). Toward this end, an M = 3 wall jet is considered, and the resulting heat transfer to the wall is computed. The values are compared against experimental data available in Reference 1. Also, FDNS's ability to compute heat flux in a reacting flow will be determined by comparing the code's predictions against calorimeter data from the hot firing of a 40K combustor. The process of modeling the flow of combusting gases through the Pratt & Whitney 40K subscale combustor and nozzle is outlined. What follows in this report is a brief description of the FDNS code, with special emphasis on how it handles solid wall boundary conditions. The test cases and some FDNS solution are presented next, along with comparison to experimental data. The process of modeling the flow through a chamber and a nozzle using the FDNS code will also be outlined.

Woodbury, Keith A.↗
NASA NTRS1979

Initial comparison of single cylinder Stirling engine computer model predictions with test results

A NASA developed digital computer code for a Stirling engine, modelling the performance of a single cylinder rhombic drive ground performance unit (GPU), is presented and its predictions are compared to test results. The GPU engine incorporates eight regenerator/cooler units and the engine working space is modelled by thirteen control volumes. The model calculates indicated power and efficiency for a given engine speed, mean pressure, heater and expansion space metal temperatures and cooler water inlet temperature and flow rate. Comparison of predicted and observed powers implies that the reference pressure drop calculations underestimate actual pressure drop, possibly due to oil contamination in the regenerator/cooler units, methane contamination in the working gas or the underestimation of mechanical loss. For a working gas of hydrogen, the predicted values of brake power are from 0 to 6% higher than experimental values, and brake efficiency is 6 to 16% higher, while for helium the predicted brake power and efficiency are 2 to 15% higher than the experimental.

Tew, R. C., Jr.↗
NASA NTRS1991

Theoretical and experimental studies relevant to interpretation of auroral emissions

The accomplishments achieved over the past year are detailed with emphasis on the interpretation or auroral emissions and studies of potential spacecraft-induced contamination effects. Accordingly, the research was divided into two tasks. The first task is designed to add to the understanding of space vehicle induced external contamination. An experimental facility for simulation of the external environment for a spacecraft in low earth orbit was developed. The facility was used to make laboratory measurements of important phenomena required for improving the understanding of the space vehicle induced external environment and its effect on measurement of auroral emissions from space-based platforms. A workshop was sponsored to provide a forum for presentation of the latest research by nationally recognized experts on space vehicle contamination and to discuss the impact of this research on future missions involving space-based platforms. The second task is to add an ab initio auroral calculation to the extant ionospheric/thermospheric global modeling capabilities. Once the addition of the code was complete, the combined model was to be used to compare the relative intensities and behavior of various emission sources (dayglow, aurora, etc.). Such studies are essential to an understanding of the types of vacuum ultraviolet (VUV) auroral images which are expected to be available within two years with the successful deployment of the Ultraviolet Imager (UVI) on the ISTP POLAR spacecraft. In anticipation of this, the second task includes support for meetings of the science working group for the UVI to discuss operational and data analysis needs. Taken together, the proposed tasks outline a course of study designed to make significant contributions to the field of space-based auroral imaging.

Keffer, Charles E.↗
NASA NTRS2002

A Procedure for Structural Weight Estimation of Single Stage to Orbit Launch Vehicles (Interim User's Manual)

This is an interim user's manual for current procedures used in the Vehicle Analysis Branch at NASA Langley Research Center, Hampton, Virginia, for launch vehicle structural subsystem weight estimation based on finite element modeling and structural analysis. The process is intended to complement traditional methods of conceptual and early preliminary structural design such as the application of empirical weight estimation or application of classical engineering design equations and criteria on one dimensional "line" models. Functions of two commercially available software codes are coupled together. Vehicle modeling and analysis are done using SDRC/I-DEAS, and structural sizing is performed with the Collier Research Corp. HyperSizer program.

Martinovic, Zoran N.↗
NASA NTRS2005

Heat Transfer on a Flat Plate with Uniform and Step Temperature Distributions

Heat transfer associated with turbulent flow on a step-heated or cooled section of a flat plate at zero angle of attack with an insulated starting section was computationally modeled using the GASP Navier-Stokes code. The algebraic eddy viscosity model of Baldwin-Lomax and the turbulent two-equation models, the K- model and the Shear Stress Turbulent model (SST), were employed. The variations from uniformity of the imposed experimental temperature profile were incorporated in the computations. The computations yielded satisfactory agreement with the experimental results for all three models. The Baldwin- Lomax model showed the closest agreement in heat transfer, whereas the SST model was higher and the K-omega model was yet higher than the experiments. In addition to the step temperature distribution case, computations were also carried out for a uniformly heated or cooled plate. The SST model showed the closest agreement with the Von Karman analogy, whereas the K-omega model was higher and the Baldwin-Lomax was lower.

Bahrami, Parviz A.↗
NASA NTRS2010

Concerning the Perturbations Generated by Flux Transfer Events

The interaction of the solar wind with the Earth's magnetosphere is often highly unsteady. Bursts of magnetic reconnection at multiple locations on the dayside equatorial magnetopause generate flux transfer events or FTEs: twisted ropes of interconnected magnetosheath and magnetospheric magnetic field lines. Once formed, the events move antisunward, displacing and perturbing the ambient media. This talk explores the perturbations predicted by both global numerical simulations and analytical models. Results from global hybrid code simulation confirm the predictions of analytical models indicating that the events generate standing forward slow mode waves as their speeds relative to the magnetosheath flow approach the Alfven velocity. Geometric considerations lead to the conclusion that events generated by component reconnection on the dayside magnetopause move poleward and exhibit strong signatures during intervals of southward IMF orientation, but move towards the flanks and exhibit weak signatures during intervals of northward IMF orientation. Changing event orientations and magnetosheath/magnetospheric magnetic field orientation can enhance the amplitudes of events reaching the flanks. Although the orientations of events on the flanks inferred from multispacecraft timing techniques are consistent with the predictions of the component reconnection model, occurrence patterns versus latitude and IMF orientation require an explanation in terms of both the component and antiparallel reconnection models.

Sibeck, David G.↗
NASA NTRS2013

Vertical Settling and Radial Segregation of Large Dust Grains in the Circumstellar Disk of the Butterfly Star

Circumstellar disks are considered as the environment for the formation of planets. The growth of dust grains in these disks is the first step in the core accretion-gas capture planet formation scenario. Indicators as well as evidence of disk evolution can be traced in spatially resolved images and the spectral energy distribution (SED) of these objects. Aims. We aim at developing a model for the dust phase of the edge-on orientated circumstellar disk of the Butterfly Star which allows one to fit observed multi-wavelength images as well as the SED simultaneously. Our model is based on spatially resolved high angular resolution observations at 1.3 mm, 894 μm, 2.07 μm, 1.87 μm, 1.60 μm and 1.13 μm as well as an extensively covered SED ranging from 12 μm to 2.7 mm, including a detailed spectrum obtained with the Spitzer Space Telescope in the range from 12 μm to 38 μm. A parameter study based on a grid search method, involving the detailed analysis of every parameter, was performed to constrain the disk parameters and find the best-fit model for the independent observations. The individual observations were modeled simultaneously, using our continuum radiative transfer code MC3D. We derived a model that is capable of reproducing all of the observations of the disk at the same time. We find quantitative evidence for grain growth up to ~100 μm sized particles, vertical settling of larger dust grains toward the disk midplane, and radial segregation of the latter toward the central star. Our results are consistent with current theoretical models for the evolution of circumstellar disks and the early stages of planet formation.

Radiative transfer↗