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.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1973

Derivation of a three dimensional numerical water quality model for estuary and continental shelf application

A derivation is given for a three dimensional mass transport equation which is appropriate for numerical modeling of estuary and continental shelf water quality variations for both the time dependent and steady state cases. A finite difference approximation to the derived equation is presented and a solution scheme for the resulting equations outlined. Preliminary results are obtained using the model for the extremely simple problems which have analytical solutions. The numerical model, as presented, will provide a scheme to study water quality problems in coastal waters for both steady state and time dependent cases.

Spaulding, M.↗

NASA NTRS2000

The fluid dynamic and shear environment in the NASA/JSC rotating-wall perfused-vessel bioreactor

The rotating-wall perfused-vessel (RWPV) bioreactor, used for both microgravity and Earth-based cell science experiments, is characterized in terms of the fluid dynamic and fluid shear stress environment. A numerical model of the flow field is developed and verified with laser Doppler velocimeter measurements. The effects of changes in operating conditions, including rotation rates and fluid perfusion rates, are investigated with the numerical model. The operating conditions typically used for ground-based experiments (equal rotation of the inner and outer cylinders) leads to flow patterns with relatively poor mass distribution characteristics. Approximately 50% of the inlet-perfused fluid bypasses the bulk of the fluid volume and flows to the perfusion exit. For operating conditions typical in microgravity, small differential rotation rates between the inner and outer cylinders lead to greatly improved flow distribution patterns and very low fluid shear stress levels over a large percentage of the fluid volume. Differences in flow patterns for the different operating conditions are explored. Large differences in the hydrodynamic environments for operating conditions typical of true microgravity and ground-based "microgravity simulations" are demonstrated.

Bioreactors↗

NASA NTRS1996

Computer Modeling of Non-Isothermal Crystallization

A realistic computer model for simulating isothermal and non-isothermal phase transformations proceeding by homogeneous and heterogeneous nucleation and interface-limited growth is presented. A new treatment for particle size effects on the crystallization kinetics is developed and is incorporated into the numerical model. Time-dependent nucleation rates, size-dependent growth rates, and surface crystallization are also included. Model predictions are compared with experimental measurements of DSC/DTA peak parameters for the crystallization of lithium disilicate glass as a function of particle size, Pt doping levels, and water content. The quantitative agreement that is demonstrated indicates that the numerical model can be used to extract key kinetic data from easily obtained calorimetric data. The model can also be used to probe nucleation and growth behavior in regimes that are otherwise inaccessible. Based on a fit to data, an earlier prediction that the time-dependent nucleation rate in a DSC/DTA scan can rise above the steady-state value at a temperature higher than the peak in the steady-state rate is demonstrated.

Kelton, K. F.↗

NASA NTRS1986

Modelling the pole tide and its effect on the earth's rotation

The pole tide is the response of the ocean to incremental centrifugal forces associated with the Chandler wobble. The tide has a potentially important effect on the period and damping of the wobble, but it is at present not well constrained by observations. Here, both analytical and numerical models for the pole tide are constructed. The analytical models consider the tide first in a global ocean and then in an enclosed basin on a beta-plane. The results are found to approach equilibrium linearly with decreasing frequency and inversely with increasing basin depth. The numerical models solve Laplace's tidal equations over the world's oceans using realistic continental boundaries and bottom topography. The results indicate that the effects of non-equilibrium portion of the deep ocean tide on the Chandler wobble period and damping are negligible.

Carton, J. A.↗

NASA NTRS2010

Rotorcraft Transmission Noise Path Model, Including Distributed Fluid Film Bearing Impedance Modeling

A computational approach for simulating the effects of rolling element and journal bearings on the vibration and sound transmission through gearboxes has been demonstrated. The approach, using ARL/Penn State s CHAMP methodology, uses Component Mode Synthesis of housing and shafting modes computed using Finite Element (FE) models to allow for rapid adjustment of bearing impedances in gearbox models. The approach has been demonstrated on NASA GRC s test gearbox with three different bearing configurations: in the first condition, traditional rolling element (ball and roller) bearings were installed, and in the second and third conditions, the traditional bearings were replaced with journal and wave bearings (wave bearings are journal bearings with a multi-lobed wave pattern on the bearing surface). A methodology for computing the stiffnesses and damping in journal and wave bearings has been presented, and demonstrated for the journal and wave bearings used in the NASA GRC test gearbox. The FE model of the gearbox, along with the rolling element bearing coupling impedances, was analyzed to compute dynamic transfer functions between forces applied to the meshing gears and accelerations on the gearbox housing, including several locations near the bearings. A Boundary Element (BE) acoustic model was used to compute the sound radiated by the gearbox. Measurements of the Gear Mesh Frequency (GMF) tones were made by NASA GRC at several operational speeds for the rolling element and journal bearing gearbox configurations. Both the measurements and the CHAMP numerical model indicate that the journal bearings reduce vibration and noise for the second harmonic of the gear meshing tones, but show no clear benefit to using journal bearings to reduce the amplitudes of the fundamental gear meshing tones. Also, the numerical model shows that the gearbox vibrations and radiated sound are similar for journal and wave bearing configurations.

Hambric, Stephen A.↗

NASA NTRS2017

Integrated Model of the Eye/Optic Nerve Head Biomechanical Environment

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a concern for long-duration space flight. Previously, it has been suggested that ocular changes observed in VIIP syndrome are related to the cephalad fluid shift that results in altered fluid pressures [1]. We are investigating the impact of changes in intracranial pressure (ICP) using a combination of numerical models, which simulate the effects of various environment conditions, including finite element (FE) models of the posterior eye. The specific interest is to understand how altered pressures due to gravitational changes affect the biomechanical environment of tissues of the posterior eye and optic nerve sheath. METHODS: Additional description of the numerical modeling is provided in the IWS abstract by Nelson et al. In brief, to simulate the effects of a cephalad fluid shift on the cardiovascular and ocular systems, we utilized a lumped-parameter compartment model of these systems. The outputs of this lumped-parameter model then inform boundary conditions (pressures) for a finite element model of the optic nerve head (Figure 1). As an example, we show here a simulation of postural change from supine to 15 degree head-down tilt (HDT), with primary outcomes being the predicted change in strains at the optic nerve head (ONH) region, specifically in the lamina cribrosa (LC), retrolaminar optic nerve, and prelaminar neural tissue (PLNT). The strain field can be decomposed into three orthogonal components, denoted as the first, second and third principal strains. We compare the peak tensile (first principal) and compressive (third principal) strains, since elevated strain alters cell phenotype and induces tissue remodeling. RESULTS AND CONCLUSIONS: Our lumped-parameter model predicted an IOP increase of c. 7 mmHg after 21 minutes of 15 degree HDT, which agreed with previous reports of IOP in HDT [1]. The corresponding FEM simulations predicted a relative increase in the magnitudes of the peak tensile and compressive strains in the lamina cribrosa of 42 and 43, respectively (Fig. 2). The corresponding changes in the optic nerve strains were 17 and 39, while in the PLNT they were 47 and 43. These magnitudes of relative elevations in peak strains may induce a phenotypic response in resident mechano-responsive resident cells [2]. This approach may be expanded to investigate other environmental changes (e.g. parabolic flight). Through our VIIP SCHOLAR project, we will validate and improve these integrated models by measuring patient-specific changes in optic nerve sheath geometry in patients with idiopathic intracranial hypertension before and after lumbar puncture and CSF removal.

biomechanics↗

NASA NTRS1984

Solute transport and the prediction of breakaway oxidation in gamma + beta Ni-Cr-Al alloys

The Al transport and the condition leading to breakaway oxidation during the cyclic oxidation of gamma + beta NiCrAl alloys have been studied. The Al concentration/distance profiles were measured after various cyclic oxidation exposures at 1200 C. It was observed that cyclic oxidation results in a decreasing Al concentration at the oxide/metal interface, maintaining a constant flux of Al to the Al2O3 scale. It was also observed that breakaway oxidation occurs when the Al concentration at the oxide/metal interface approaches zero. A numerical model was developed to simulate the diffusional transport of Al and to predict breakaway oxidation in gamma + beta NiCrAl alloys undergoing cyclic oxidation. In a comparison of two alloys with similar oxide spalling characteristics, the numerical model was shown to predict correctly the onset of breakaway oxidation in the higher Al-content alloy.

Nesbitt, J. A.↗

NASA NTRS1990

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) eV and 1 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass less than about 10(exp -2) eV, axions produced deep inside the neutron star simply stream out; in a previous paper this case has been addressed. Remarkably, for an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become 'trapped' and radiated from an axion sphere. In this paper the trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a leakage approximation scheme for axion-energy transport. The axion opacity is computed due to inverse nucleon-nucleon, axion bremsstrahlung, and numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande 2 (K2) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The earlier estimate is confirmed and refined of the axion mass above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration--the most sensitive barometer of axion cooling--it is concluded that for an axion mass of greater than about 0.3 eV, axion emission would not have had a significant effect on the neutrino bursts detected by K2 and IMB. The present work, together with the previous work, strongly suggests that an axion with mass in the interval 10(exp -3) eV to 0.3 eV is excluded by SN 1987A.

Burrows, Adam↗

NASA NTRS1994

On-Orbit Propellant Motion Resulting from an Impulsive Acceleration

In-space docking and separation maneuvers of spacecraft that have large fluid mass fractions may cause undesirable spacecraft motion in response to the impulsive-acceleration-induced fluid motion. An example of this potential low-gravity fluid management problem arose during the development of the shuttle/Centaur vehicle. Experimentally verified numerical modeling techniques were developed to establish the propellant dynamics, and subsequent vehicle motion, associated with the separation of the Centaur vehicle from the shuttle cargo bay. Although the shuttle/Centaur development activity has been suspended, the numerical modeling techniques are available to predict on-orbit liquid motion resulting from impulsive accelerations for other missions and spacecraft.

Hochstein, John I.↗

NASA NTRS1995

The NASA Ames 16-Inch Shock Tunnel Nozzle Simulations and Experimental Comparison

The 16-Inch Shock Tunnel at NASA Ames Research Center is a unique test facility used for hypersonic propulsion testing. To provide information necessary to understand the hypersonic testing of the combustor model, computational simulations of the facility nozzle were performed and results are compared with available experimental data, namely static pressure along the nozzle walls and pitot pressure at the exit of the nozzle section. Both quasi-one-dimensional and axisymmetric approaches were used to study the numerous modeling issues involved. The facility nozzle flow was examined for three hypersonic test conditions, and the computational results are presented in detail. The effects of variations in reservoir conditions, boundary layer growth, and parameters of numerical modeling are explored.

TokarcikPolsky, S.↗

NASA NTRS1986

Methods for computing comet core temperatures

The temperature profile within the comet nucleus provides the key to an understanding of the history of the volatiles within a comet. Certain difficulties arise in connection with current cometary temperature models. It is shown that the constraint of zero net heat flow can be used to derive general analytical expressions which will allow for the determination of comet core temperature for a spherically symmetric comet, taking into account information about the surface temperature and the thermal conductivity. The obtained results are compared with the expression for comet core temperatures considered by Klinger (1981). Attention is given to analytical results, an example case, and numerical models. The formalization developed makes it possible to determine the core temperature on the basis of the numerical models of the surface temperature.

Mckay, C. P.↗

NASA NTRS1990

Axions and SN 1987A: Axion trapping

If an axion of mass between about 10(exp -3) and 10 eV exists, axion emission would have significantly affected the cooling of the nascent neutron star associated with SN 1987A. For an axion of mass greater than about 10(exp -2) eV axions would, like neutrinos, have a mean-free path that is smaller than the size of a neutron star, and thus would become trapped and radiated from an axion sphere. The trapping regime is treated by using numerical models of the initial cooling of a hot neutron star that incorporate a diffusion approximation for axion-energy transport. The axion opacity due to inverse nucleon-nucleon, axion bremsstrahlung is computed; and then the numerical models are used to calculate the integrated axion luminosity, the temperature of the axion sphere, and the effect of axion emission on the neutrino bursts detected by the Kamiokande II (KII) and Irvine-Michigan-Brookhaven (IMB) water-Cherenkov detectors. The larger the axion mass, the stronger the trapping and the smaller the axion luminosity. The estimate of the axion mass is confirmed above which trapping is so strong that axion emission does not significantly affect the neutrino burst. Based upon the neutrino-burst duration - the most sensitive barometer of axion cooling - it is concluded that for an axion mass greater than about 3 eV axion emission would not have had a significant effect on the neutrino bursts detected by KII and IMB. It is strongly suggested that an axion with mass in the interval 10(exp -3) to 3 eV is excluded by the observation of neutrinos from SN 1987A.

Burrows, Adam↗

NASA NTRSDate not supplied

Progress on the Reduced Gravity Cryogenic Transfer (RGCT) Project

Reduced gravity transfer of cryogenic propellants is crucial to being able to perform planned extended-duration space exploration missions. However, the transfer of cryogenic propellants under reduced gravity conditions has never been demonstrated in a mass-efficient fashion. Efficient cryogenic fluid transfer methods will reduce the transfer time or amount of propellant consumed for chilldown of transfer lines and tanks and ensure successful engine restart or fill of a customer receiver tank. The purpose of the Reduced Gravity Cryogenic Transfer (RGCT) project is to enhance line chilldown, tank chilldown, and tank fill/transfer in a reduced gravity environment through (1) ground and reduced gravity cryogenic testing and (2) numerical model development and validation, which includes empirical, lumped capacitance, and computational fluid dynamics modeling. Technologies developed, data gathered, and models developed and validated under RGCT have played a critical role in enabling longer duration in-space missions. This presentation will provide a description of the ongoing cryogenic propellant transfer ground and reduced gravity testing, numerical modeling, and technology development to-date.

Boiling↗

NASA NTRS1999

1-D Numerical Analysis of ABCC Engine Performance

ABCC engine combines air breathing and rocket engine into a single engine to increase the specific impulse over an entire flight trajectory. Except for the heat source, the basic operation of the ABCC is similar to the basic operation of the RBCC engine. The ABCC is intended to have a higher specific impulse than the RBCC for single stage Earth to orbit vehicle. Computational fluid dynamics (CFD) is a useful tool for the analysis of complex transport processes in various components in ABCC propulsion system. The objective of the present research was to develop a transient 1-D numerical model using conservation of mass, linear momentum, and energy equations that could be used to predict flow behavior throughout a generic ABCC engine following a flight path. At specific points during the development of the 1-D numerical model a myriad of tests were performed to prove the program produced consistent, realistic numbers that follow compressible flow theory for various inlet conditions.

Holden, Richard↗

NASA NTRS1991

Numerical simulations of a nonequilibrium argon plasma in a shock-tube experiment

A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.

Cambier, Jean-Luc↗

NASA NTRSDate not supplied

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗

NASA NTRS1997

A Study of Fundamental Shock Noise Mechanisms

This paper investigates two mechanisms fundamental to sound generation in shocked flows: shock motion and shock deformation. Shock motion is modeled numerically by examining the interaction of a sound wave with a shock. This numerical approach is validated by comparison with results obtained by linear theory for a small-disturbance case. Analysis of the perturbation energy with Myers' energy corollary demonstrates that acoustic energy is generated by the interaction of acoustic disturbances with shocks. This analysis suggests that shock motion generates acoustic and entropy disturbance energy. Shock deformation is modeled numerically by examining the interaction of a vortex ring with a shock. These numerical simulations demonstrate the generation of both an acoustic wave and contact surfaces. The acoustic wave spreads cylindrically. The sound intensity is highly directional and the sound pressure increases with increasing shock strength. The numerically determined relationship between the sound pressure and the Mach number is found to be consistent with experimental observations of shock noise. This consistency implies that a dominant physical process in the generation of shock noise is modeled in this study.

Meadows, Kristine R.↗

NASA NTRS2015

Impact of Asian Aerosols on Precipitation Over California: An Observational and Model Based Approach

Dust and pollution emissions from Asia are often transported across the Pacific Ocean to over the western United States. Therefore, it is essential to fully understand the impact of these aerosols on clouds and precipitation forming over the eastern Pacific and western United States, especially during atmospheric river events that account for up to half of California's annual precipitation and can lead to widespread flooding. In order for numerical modeling simulations to accurately represent the present and future regional climate of the western United States, we must account for the aerosol-cloud-precipitation interactions associated with Asian dust and pollution aerosols. Therefore, we have constructed a detailed study utilizing multi-sensor satellite observations, NOAA-led field campaign measurements, and targeted numerical modeling studies where Asian aerosols interacted with cloud and precipitation processes over the western United States. In particular, we utilize aerosol optical depth retrievals from the NASA Moderate Resolution Imaging Spectroradiometer (MODIS), NOAA Geostationary Operational Environmental Satellite (GOES-11), and Japan Meteorological Agency (JMA) Multi-functional Transport Satellite (MTSAT) to effectively detect and monitor the trans-Pacific transport of Asian dust and pollution. The aerosol optical depth (AOD) retrievals are used in assimilating the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) in order to provide the model with an accurate representation of the aerosol spatial distribution across the Pacific. We conduct WRF-Chem model simulations of several cold-season atmospheric river events that interacted with Asian aerosols and brought significant precipitation over California during February-March 2011 when the NOAA CalWater field campaign was ongoing. The CalWater field campaign consisted of aircraft and surface measurements of aerosol and precipitation processes that help extensively validate our WRF-Chem model simulations. After validating the capability of the WRF-Chem in realistically simulating the aerosol-cloud precipitation interactions, we conduct sensitivity studies where the AOD is doubled to diagnose whether an increasing concentration of Asian aerosols over the western United States will lead to further impacts on the cloud and precipitation processes over California. We also perform sensitivity studies where the aerosols will be partitioned into dust-only and pollution-only in order to separate the impacts of the differing Asian aerosol species. The results of our WRF-Chem model simulations aim to show that the trans-Pacific transport of Asian aerosols influence the precipitation associated with atmospheric river events that can ultimately impact the regional climate of the western United States. 1 University

Naeger, Aaron R.↗