Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Solar transition region response to variations in the heating rate

The response of a numerical model for the upper chromosphere, transition region, and corona to variations in the energy input has been examined. The numerical model solves the set of one-dimensional two-fluid hydrodynamic equations in a simple vertical magnetic flux tube. The atmosphere responds to both the increase and decrease in energy deposition by smoothly readjusting the temperature gradient and the amount of material in the region of peak radiating efficiency to radiate away energy being deposited. At no time during this readjustment is a departure from a thin laminar transition region structure seen. In addition, a time-dependent description of the nonequilibrium ionization of all of the ionization stages of oxygen has been included. This calculation is coupled with the self-consistent calculations of the dynamical variables. It is found that the nonequilibrium ionization balance calculations for both heating and cooling small loops in the quiet sun predict relative ionic abundances which differ substantially from those which would be predicted by an equilibrium calculation

Mariska, J. T.↗
NASA NTRS2005

Computing Thermal Effects of Cavitation in Cryogenic Liquids

A computer program implements a numerical model of thermal effects of cavitation in cryogenic fluids. The model and program were developed for use in designing and predicting the performances of turbopumps for cryogenic fluids. Prior numerical models used for this purpose do not account for either the variability of properties of cryogenic fluids or the thermal effects (especially, evaporative cooling) involved in cavitation. It is important to account for both because in a cryogenic fluid, the thermal effects of cavitation are substantial, and the cavitation characteristics are altered by coupling between the variable fluid properties and the phase changes involved in cavitation. The present model accounts for both thermal effects and variability of properties by incorporating a generalized representation of the properties of cryogenic fluids into a generalized compressible-fluid formulation for a cavitating pump. The model has been extensively validated for liquid nitrogen and liquid hydrogen. Using the available data on the properties of these fluids, the model has been shown to predict accurate temperature-depression values.

Hosangadi, Ashvin↗
NASA NTRS1985

Coating life prediction

The investigation combines both experimental studies and numerical modeling to predict coating life in an oxidizing environment. The experimental work provides both input to and verification of two numerical models. The coatings being examined are an aluminide coating on Udimet 700 (U-700), a low-pressure plasma spray (LPPS) Ni-18Co-17Cr-24Al-0.2Y overlay coating also on U- 700, and bulk deposits of the LPPS NiCoCrAlY coating.

Nesbitt, James A.↗
NASA NTRS2016

The Impacts of Dry Dynamic Cores on Asymmetric Hurricane Intensification

The fundamental pathways for tropical cyclone (TC) intensification are explored by considering axisymmetric and asymmetric impulsive thermal perturbations to balanced, TC-like vortices using the dynamic cores of three different nonlinear numerical models. Attempts at reproducing the results of previous work, which used the community WRF Model, revealed a discrepancy with the impacts of purely asymmetric thermal forcing. The current study finds that thermal asymmetries can have an important, largely positive role on the vortex intensification, whereas other studies find that asymmetric impacts are negligible. Analysis of the spectral energetics of each numerical model indicates that the vortex response to asymmetric thermal perturbations is significantly damped in WRF relative to the other models. Spectral kinetic energy budgets show that this anomalous damping is primarily due to the increased removal of kinetic energy from the vertical divergence of the vertical pressure flux, which is related to the flux of inertia-gravity wave energy. The increased kinetic energy in the other two models is shown to originate around the scales of the heating and propagate upscale with time from nonlinear effects. For very large thermal amplitudes (50 K), the anomalous removal of kinetic energy due to inertia-gravity wave activity is much smaller, resulting in good agreement between models. The results of this paper indicate that the numerical treatment of small-scale processes that project strongly onto inertia-gravity wave energy can lead to significant differences in asymmetric TC intensification. Sensitivity tests with different time integration schemes suggest that diffusion entering into the implicit solution procedure is partly responsible for the anomalous damping of energy.

Guimond, Stephen R.↗
NASA NTRS1999

Development of a Linearized Unsteady Euler Analysis with Application to Wake/Blade-Row Interactions

A three-dimensional, linearized, Euler analysis is being developed to provide a comprehensive and efficient unsteady aerodynamic analysis for predicting the aeroacoustic and aeroelastic responses of axial-flow turbomachinery blading. The mathematical models needed to describe nonlinear and linearized, inviscid, unsteady flows through a blade row operating within a cylindrical annular duct are presented in this report. A numerical model for linearized inviscid unsteady flows, which couples a near-field, implicit, wave-split, finite volume analysis to far-field eigen analyses, is also described. The linearized aerodynamic and numerical models have been implemented into the three-dimensional unsteady flow code, LINFLUX. This code is applied herein to predict unsteady subsonic flows driven by wake or vortical excitations. The intent is to validate the LINFLUX analysis via numerical results for simple benchmark unsteady flows and to demonstrate this analysis via application to a realistic wake/blade-row interaction. Detailed numerical results for a three-dimensional version of the 10th Standard Cascade and a fan exit guide vane indicate that LINFLUX is becoming a reliable and useful unsteady aerodynamic prediction capability that can be applied, in the future, to assess the three-dimensional flow physics important to blade-row, aeroacoustic and aeroelastic responses.

Verdon, Joseph M.↗
NASA NTRS2008

FAA Development of Reliable Modeling Methodologies for Fan Blade Out Containment Analysis: Ballistic Impact Testing - Part 2

This report summarizes the ballistic impact testing that was conducted to provide validation data for the development of numerical models of blade out events in fabric containment systems. The ballistic impact response of two different fiber materials - Kevlar 49 (E.I. DuPont Nemours and Company) and Zylon AS (Toyobo Co., Ltd.) was studied by firing metal projectiles into dry woven fabric specimens using a gas gun. The shape, mass, orientation and velocity of the projectile were varied and recorded. In most cases the tests were designed such that the projectile would perforate the specimen, allowing measurement of the energy absorbed by the fabric. The results for both Zylon and Kevlar presented here represent a useful set of data for the purposes of establishing and validating numerical models for predicting the response of fabrics under conditions simulating those of a jet engine blade release situations. In addition some useful empirical observations were made regarding the effects of projectile orientation and the relative performance of the different materials.

Revilock, Duane M.↗
NASA NTRS2003

Particle Engulfment and Pushing By Solidifying Interfaces

The study of particle behavior at solid/liquid interfaces (SLI s) is at the center of the Particle Engulfment and Pushing (PEP) research program. Interactions of particles with SLI s have been of interest since the 1960 s, starting with geological observations, i.e., frost heaving. Ever since, this field of research has become significant to such diverse areas as metal matrix composite materials, fabrication of superconductors, and inclusion control in steels. The PEP research effort is geared towards understanding the fundamental physics of the interaction between particles and a planar SLI. Experimental work including 1-g and mu-g experiments accompany the development of analytical and numerical models. The experimental work comprised of substantial groundwork with aluminum (Al) and zinc (Zn) matrices containing spherical zirconia particles, mu-g experiments with metallic Al matrices and the use of transparent organic metal-analogue materials. The modeling efforts have grown from the initial steady-state analytical model to dynamic models, accounting for the initial acceleration of a particle at rest by an advancing SLI. To gain a more comprehensive understanding, numerical models were developed to account for the influence of the thermal and solutal field. Current efforts are geared towards coupling the diffusive 2-D front tracking model with a fluid flow model to account for differences in the physics of interaction between 1-g and -g environments. A significant amount of this theoretical investigation has been and is being performed by co-investigators at NASA MSFC.

Source record↗
NASA NTRS2005

Fluid Distribution for In-space Cryogenic Propulsion

The ultimate goal of this task is to enable the use of a single supply of cryogenic propellants for three distinct spacecraft propulsion missions: main propulsion, orbital maneuvering, and attitude control. A fluid distribution system is sought which allows large propellant flows during the first two missions while still allowing control of small propellant flows during attitude control. Existing research has identified the probable benefits of a combined thermal management/power/fluid distribution system based on the Solar Integrated Thermal Management and Power (SITMAP) cycle. Both a numerical model and an experimental model are constructed in order to predict the performance of such an integrated thermal management/propulsion system. This research task provides a numerical model and an experimental apparatus which will simulate an integrated thermal/power/fluid management system based on the SITMAP cycle, and assess its feasibility for various space missions. Various modifications are done to the cycle, such as the addition of a regeneration process that allows heat to be transferred into the working fluid prior to the solar collector, thereby reducing the collector size and weight. Fabri choking analysis was also accounted for. Finally the cycle is to be optimized for various space missions based on a mass based figure of merit, namely the System Mass Ratio (SMR). -. 1 he theoretical and experimental results from these models are be used to develop a design code (JETSIT code) which is able to provide design parameters for such a system, over a range of cooling loads, power generation, and attitude control thrust levels. The performance gains and mass savings will be compared to those of existing spacecraft systems.

Lear, William↗
NASA NTRS1996

A Three-Dimensional Linearized Unsteady Euler Analysis for Turbomachinery Blade Rows

A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic response characteristics of axial-flow turbomachinery blading. The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. In addition, a numerical model for linearized inviscid unsteady flow, which is based upon an existing nonlinear, implicit, wave-split, finite volume analysis, is described. These aerodynamic and numerical models have been implemented into an unsteady flow code, called LINFLUX. A preliminary version of the LINFLUX code is applied herein to selected, benchmark three-dimensional, subsonic, unsteady flows, to illustrate its current capabilities and to uncover existing problems and deficiencies. The numerical results indicate that good progress has been made toward developing a reliable and useful three-dimensional prediction capability. However, some problems, associated with the implementation of an unsteady displacement field and numerical errors near solid boundaries, still exist. Also, accurate far-field conditions must be incorporated into the FINFLUX analysis, so that this analysis can be applied to unsteady flows driven be external aerodynamic excitations.

Montgomery, Matthew D.↗
NASA NTRS2017

Uncertainty Aware Structural Topology Optimization Via a Stochastic Reduced Order Model Approach

This work presents a stochastic reduced order modeling strategy for the quantification and propagation of uncertainties in topology optimization. Uncertainty aware optimization problems can be computationally complex due to the substantial number of model evaluations that are necessary to accurately quantify and propagate uncertainties. This computational complexity is greatly magnified if a high-fidelity, physics-based numerical model is used for the topology optimization calculations. Stochastic reduced order model (SROM) methods are applied here to effectively 1) alleviate the prohibitive computational cost associated with an uncertainty aware topology optimization problem; and 2) quantify and propagate the inherent uncertainties due to design imperfections. A generic SROM framework that transforms the uncertainty aware, stochastic topology optimization problem into a deterministic optimization problem that relies only on independent calls to a deterministic numerical model is presented. This approach facilitates the use of existing optimization and modeling tools to accurately solve the uncertainty aware topology optimization problems in a fraction of the computational demand required by Monte Carlo methods. Finally, an example in structural topology optimization is presented to demonstrate the effectiveness of the proposed uncertainty aware structural topology optimization approach.

Aguilo, Miguel A.↗
NASA NTRS1999

The Effect of Velocity on the Extinction Behavior of a Diffusion Flame during Transient Depressurization

Current fire suppression plans for the International Space Station include the use of venting (depressurization) as a method for extinguishing a fire. Until recently this process had only been examined as part of a material flammability experiment performed on Skylab in the early 1970's. Due to the low initial pressure (0.35 Atm) and high oxygen concentration (65%), the Skylab experimental results are not applicable for understanding the effects of venting on a fire in a space station environment (21%O2, 1 Atm). Recent research examined the extinction behavior of a diffusion flame over a polymethyl methacrylate (PMMA) cylinder during a transient depressurization in low-gravity. The numerical model was used to examine extinction limits as a function of depressurization rate, forced flow velocity, and initial solid phase temperature. The experimental and numerically predicted extinction data indicated that as the solid phase temperature increased the pressure required to extinguish the flame decreased. The numerical model was also used to examine conditions not obtainable in the low-gravity experiments. From these simulations, a series of extinction boundaries were generated that showed a region of increased flammability existed at a forced flow of 10 cm/s. Analysis of these extinction boundaries indicated that they were quasi-steady in nature, and that the final extinction conditions were independent of the transient process. The velocity range in the previous study was limited and thus the results did not examine the effects of velocities less than 1 cm/s or greater than 20 cm/s. This study utilized low-gravity experiments performed on NASA's Reduced-gravity Research Aircraft Laboratory and numerical simulations to examine conditions applicable to the Space Station environment. This paper extends the analysis of the previous study to a comprehensive examination of the effect of increased velocity on extinction behavior and extinction limits during a transient depressurization in low-gravity. This is achieved by examining extinction data from buoyant (normal-gravity) and low-buoyant (low-gravity) depressurization. experiments, as well as from numerical predictions of flame behavior during depressurization in a non-buoyant (zero-gravity) environment.

Goldmeer, Jeffrey S.↗
NASA NTRS1987

Transport models for numerical forecast

The explosive growth of computing power, coupled with scientific and technological emphasis on a national scale, has led to significant major advances in operational numerical weather prediction (NWP) during the last two decades. There are about half a dozen major centers around the world running global NWP models operationally. Many more countries have operational hemispheric or limited-area models which provide weather forecasts. The global models typically have several hundred kilometer resolution, while the limited-area models usually have horizontal spacing of 50 to 100 km. Given the pace of burgeoning growth in this area, it seems warranted to occasionally take an overview of aspects of the field common to all modelers. Here, a brief look is taken at the nature of subgrid scale turbulence transport parameterization, and some of the difficulties pertaining thereto, with particular emphasis on operational NWP models.

Burk, Stephen D.↗
NASA NTRSDate not supplied

Determining the Cause of Reduced Concurrent Flame Spread over Thin Solid Fuels in Low Pressure and Low Gravity

The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.

Fire Safety↗
NASA NTRSDate not supplied

Determining the Cause of Reduced Concurrent Flame Spread over Thin Solid Fuels in Low Pressure and Low Gravity

The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.

microgravity↗
NASA NTRS1988

Analysis and modeling of summertime convective cloud and precipitation structure over the southeastern United States

Described is work performed under NASA Grant NAG8-654 for the period 15 March to 15 September 1988. This work entails primarily data analysis and numerical modeling efforts related to the 1986 Satellite Precipitation and Cloud Experiment (SPACE). In the following, the SPACE acronym is used along with the acronym COHMEX, which represents the encompassing Cooperative Huntsville Meteorological Experiment. Progress made during the second half of the first year of the study included: (1) installation and testing of the RAMS numerical Modeling system on the Alabama CRAY X-MP/24; (2) a start on the analysis of the mesoscale convection system (MCS) of 13 July 1986 COHMEX case; and (3) a cursory examination of a small MCS that formed over the COHMEX region on 15 July 1986. Details of each of these individual tasks are given.

Knupp, Kevin R.↗
NASA NTRS2009

Explicit Finite Element Modeling of Multilayer Composite Fabric for Gas Turbine Engine Containment Systems, Phase II: Ballistic Impact Testing - Part 2

This report summarizes the ballistic impact testing that was conducted to provide validation data for the development of numerical models of blade-out events in fabric containment systems. The ballistic impact response of two different fiber materials - Kevlar(TradeName) 49 and Zylon(TradeName) AS (as spun) was studied by firing metal projectiles into dry woven fabric specimens using a gas gun. The shape, mass, orientation, and velocity of the projectile were varied and recorded. In most cases, the tests were designed so the projectile would perforate the specimen, allowing measurement of the energy absorbed by the fabric. The results for both Zylon and Kevlar presented here represent a useful set of data for the purposes of establishing and validating numerical models to predict the response of fabrics under conditions that simulate those of a jet engine blade-release situation. In addition, some useful empirical observations were made regarding the effects of projectile orientation and the relative performance of the different fabric materials.

Zylon↗
NASA NTRS2006

The Use of RGPS Kinematic Data to Estimate Nonlinear Sea Ice Motion

In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motion In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motio

Pruis, M.↗
NASA NTRSDate not supplied

Numerical Prediction of Heat Transfer Coefficients during Xenon Tank Fill in Microgravity

NASA’s Gateway will serve as an orbital outpost to enable sustained human lunar exploration and facilitate access to destinations beyond the Moon. To maintain sustainable deep space missions, on-orbit propellant transfer and refueling must be realized. One of Gateway’s critical components, the Power and Propulsion Element (PPE), will use solar electric propulsion to perform attitude control maneuvers and orbit transfers. Xenon will be used as the propellant for the PPE and stored in a composite overwrapped pressure vessel (COPV). During the propellant refueling process, xenon will compress causing the fluid and tank walls to warm. There is concern that the bond between the tank liner and composite may degrade at elevated temperatures, along with the possibility of the tank exceeding design pressure. While on-orbit testing is being proposed, numerical models are being used to predict the thermal response and reduce risk during a refueling operation in microgravity where natural convection is diminished. Accurate numerical models may be used to inform decisions for follow-on testing, or design and operation of refueling architectures. Challenges of modeling xenon tank fill in microgravity include non-ideal gas behavior, operating near the fluid’s critical point, and lack of experimental flight data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict averaged heat transfer coefficients during tank filling. At the anticipated injection flow rates, the CFD results showed that forced convection dominates in microgravity. Initial results from a multi-node Thermal Desktop model that assumed a purely conducting fluid (Nu=1) over-predicted the CFD fluid temperature and pressure, while a similar model allowing for natural convection under-predicted. The heat transfer coefficients predicted by CFD were implemented into the multi-node model, and the COPV thermal response from the two computational approaches were compared with good agreement.

Heat Transfer↗