Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic modeling”

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

Stressed Oxidation Life Prediction for C/SiC Composites

The residual strength and life of C/SiC is dominated by carbon interface and fiber oxidation if seal coat and matrix cracks are open to allow oxygen ingress. Crack opening is determined by the combination of thermal, mechanical and thermal expansion mismatch induced stresses. When cracks are open, life can be predicted by simple oxidation based models with reaction controlled kinetics at low temperature, and by gas phase diffusion controlled kinetics at high temperatures. Key life governing variables in these models include temperature, stress, initial strength, oxygen partial pressure, and total pressure. These models are described in this paper.

Levine, Stanley R.↗

Decoupled Method for Reconstruction of Surface Conditions From Internal Temperatures On Ablative Materials With Uncertain Recession Model

Obtaining measurements of flight environments on ablative heat shields is both critical for spacecraft development and extremely challenging due to the harsh heating environment and surface recession. Thermocouples installed several millimeters below the surface are commonly used to measure the heat shield temperature response, but an ill-posed inverse heat conduction problem must be solved to reconstruct the surface heating environment from these measurements. Ablation can contribute substantially to the measurement response making solutions to the inverse problem strongly dependent on the recession model, which is often poorly characterized. To enable efficient surface reconstruction for recession model sensitivity analysis, a method for decoupling the surface recession evaluation from the inverse heat conduction problem is presented. The decoupled method is shown to provide reconstructions of equivalent accuracy to the traditional coupled method but with substantially reduced computational effort. These methods are applied to reconstruct the environments on the Mars Science Laboratory heat shield using diffusion limit and kinetically limited recession models.

Oliver, A. Brandon↗

Collisionless relaxation in spiral galaxy models

The increase in random kinetic energy of stars by rapidly fluctuating gravitational fields (collisionless or violent relaxation) in disk galaxy models is investigated for three interaction potentials of the stars corresponding to (1) point stars, (2) rod stars of length 2 kpc, and (3) uniform density spherical stars of radius 2 kpc. To stabilize the galaxy against the large scale bar forming instability, a fixed field corresponding to a central core or halo component of stars was added with the stars containing at most 20 percent of the total mass of the galaxy. Considerable heating occurred for both the point stars and the rod stars, whereas the use of spherical stars resulted in a very low heating rate. The use of spherical stars with the resulting low heating rate will be desirable for the study of large scale galactic stability or density wave propagation, since collective heating effects will no longer mask the phenomena under study.

Hohl, F.↗

Exact solutions and low-frequency instability of the adiabatic auroral arc model

The adiabatic auroral arc model couples a kinetic theory parallel current driven by mirror forces to horizontal ionospheric currents; the resulting equations are nonlinear. Some exact stationary solutions to these equations, some of them based on the Liouville equation, are developed, with both latitudinal and longitudinal spatial variations. These Liouville equation exact solutions are related to stability boundaries of low-frequency instabilities such as Kelvin-Helmholtz, as shown by a study of a simplified model.

Cornwall, John M.↗

Calcium kinetics with microgram stable isotope doses and saliva sampling

Studies of calcium kinetics require administration of tracer doses of calcium and subsequent repeated sampling of biological fluids. This study was designed to develop techniques that would allow estimation of calcium kinetics by using small (micrograms) doses of isotopes instead of the more common large (mg) doses to minimize tracer perturbation of the system and reduce cost, and to explore the use of saliva sampling as an alternative to blood sampling. Subjects received an oral dose (133 micrograms) of 43Ca and an i.v. dose (7.7 micrograms) of 46Ca. Isotopic enrichment in blood, urine, saliva and feces was well above thermal ionization mass spectrometry measurement precision up to 170 h after dosing. Fractional calcium absorptions determined from isotopic ratios in blood, urine and saliva were similar. Compartmental modeling revealed that kinetic parameters determined from serum or saliva data were similar, decreasing the necessity for blood samples. It is concluded from these results that calcium kinetics can be assessed with micrograms doses of stable isotopes, thereby reducing tracer costs and with saliva samples, thereby reducing the amount of blood needed.

Clinical Trial↗

Advances in Studies of Electrode Kinetics and Mass Transport in AMTEC Cells

Previous work reported from JPL has included characterization of electrode kinetics and alkali atom transport from electrodes including Mo, W, WRh x (Mn), in sodium AMTEC cells and vapor exposure cells; and Mo in potassium vapor exposure cells. These studies were generally performed in cells with small area electrodes (about 1 to 5 cm 2 ), and device geometry had little effect on transport. Alkali metal diffusion coefficients through these electrodes have been characterized, and approximate surface diffusion coefficients were derived in cases of activated transport. A basic model of electrode kinetics and transport at the alkali metal vapor/pourous metal electrode/alkali beta"-alumina solid electrolyte (BASE) three phase boundary has been proposed which accounts for electrochemical reaction rates with a collision frequency near the three phase boundary and tunneling from the porous electrode partially covered with adsorbed alkali metal atoms. The small electrode effect in AMTEC cells has been discussed in several papers, but quantitative investigations have described only the overall effect and the important contribution of electrolyte resistance. The quantitative characterization of transport losses in cells with large area electrodes has been limited to simulations of large area electrode effects, or characterization of transport losses from large area electrodes with significant longitudinal temperature gradients. This paper describes new investigations of electrochemical kinetics and transport, with four 14.4 cm 2 WPt 3.5 electrodes, including the influence of electrode size on the mass transport loss in the AMTEC cell. These electrodes exhibit very slow sintering, as well as excellent sodium transport properties, making them attractive candidates for AMTEC power conversion use. However, the facile sodium transport in WPt 3.5 electrodes makes characterization of the transport process difficult.

electrode↗

Advances in Studies of Electrode Kinetics and Mass Transport in AMTEC Cells (abstract)

Previous work reported from JPL has included characterization of electrode kinetics and alkali atom transport from electrodes including Mo, W, WRh(sub x), WPt(sub x)(Mn), in sodium AMTEC cells and vapor exposure cells, and Mo in potassium vapor exposure cells. These studies were generally performed in cells with small area electrodes (about 1 to 5 cm(sup 2)), and device geometry had little effect on transport. Alkali diffusion coefficients through these electrodes have been characterized, and approximate surface diffusion coefficients derived in cases of activated transport. A basic model of electrode kinetic at the alkali metal vapor/porous metal electrode/alkali beta'-alumina solid electrolyte three phase boundary has been proposed which accounts for electrochemical reaction rates with a collision frequency near the three phase boundary and tunneling from the porous electrode partially covered with adsorbed alkali metal atoms. The small electrode effect in AMTEC cells has been discussed in several papers, but quantitative investigations have described only the overall effect and the important contribution of electrolyte resistance. The quantitative characterization of transport losses in cells with large area electrodes has been limited to simulations of large area electrode effects, or characterization of transport losses from large area electrodes with significant longitudinal temperature gradients. This paper describes new investigations of electrochemical kinetics and transport, particularily with WPt(sub 3.5) electrodes, including the influence of electrode size on the mass transport loss in the AMTEC cell. These electrodes possess excellent sodium transport properties making verification of device limitations on transport much more readily attained.

electrode kinetics alkali atom transport transport↗

Active chlorine and nitric oxide formation from chemical rocket plume afterburning

Chlorine and oxides of nitrogen (NO(x)) released into the atmosphere contribute to acid rain (ground level or low-altitude sources) and ozone depletion from the stratosphere (high-altitude sources). Rocket engines have the potential for forming or activating these pollutants in the rocket plume. For instance, H2/O2 rockets can produce thermal NO(x) in their plumes. Emphasis, in the past, has been placed on determining the impact of chlorine release on the stratosphere. To date, very little, if any, information is available to understand what contribution NO(x) emissions from ground-based engine testing and actual rocket launches have on the atmosphere. The goal of this work is to estimate the afterburning emissions from chemical rocket plumes and determine their local stratospheric impact. Our study focuses on the space shuttle rocket motors, which include both the solid rocket boosters (SRB's) and the liquid propellant main engines (SSME's). Rocket plume afterburning is modeled employing a one-dimensional model incorporating two chemical kinetic systems: chemical and thermal equilibria with overlayed nitric oxide chemical kinetics (semi equilibrium) and full finite-rate chemical kinetics. Additionally, the local atmospheric impact immediately following a launch is modeled as the emissions diffuse and chemically react in the stratosphere.

Leone, D. M.↗

Multilayer Insulation Ascent Venting Model

The thermal and venting transient experienced by tank-applied multilayer insulation (MLI) in the Earth-to-orbit environment is very dynamic and not well characterized. This new predictive code is a first principles-based engineering model which tracks the time history of the mass and temperature (internal energy) of the gas in each MLI layer. A continuum-based model is used for early portions of the trajectory while a kinetic theory-based model is used for the later portions of the trajectory, and the models are blended based on a reference mean free path. This new capability should improve understanding of the Earth-to-orbit transient and enable better insulation system designs for in-space cryogenic propellant systems.

Tramel, R. W.↗

Modeling Self-Pressurization and Spray Bar Pressure Control of a Cryogenic Storage Tank in Normal Gravity

This paper presents Computational Fluid Dynamics (CFD) models for simulating self-pressurization and spray bar cooling processes in a large-scale liquid hydrogen storage tank under normal gravity. The self-pressurization model utilizes the kinetics-based Schrage equation with the Volume-Of-Fluid (VOF) model to account for the evaporative and condensing interfacial mass fluxes. It uses a laminar approach to model natural convection and compares the pressures and temperatures predicted by the CFD model with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) self-pressurization experiment. Additionally, the predicted interfacial mass flux is presented. A CFD model for simulating pressure control using a spray bar is also presented. An Eulerian-Lagrangian approach models the interactions between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by tracking particles in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. The T-sat model is used to calculate the droplet-ullage heat and mass transfer. In this model, droplets warm up to the saturation temperature corresponding to the ullage vapor pressure thereafter evaporating while remaining at the saturation temperature. The tank pressure, vapor, and liquid temperature evolutions are validated against data from the MHTB spray bar mixing experiment.

Self-Pressurization↗

Cycling phosphorus on the Archean Earth: Part I. Continental weathering and riverine transport of phosphorus

Phosphorus (P) is the key nutrient thought to limit primary productivity on geological timescales. Phosphate levels in Archean marine sediments are low, but quantification of the P cycle and how it changed through a billion years of recorded Archean history remain a challenge, hindering our understanding of the role played by P in biosphere/geosphere co-evolution on the early Earth. Here, we design kinetic and thermodynamic models to quantitatively assess one key component of the early P cycle – continental weathering – by considering the emergence and elevation of continents, as well as the evolution of climate, the atmosphere, and the absence of macroscopic vegetation during the Archean Eon. Our results suggest that the weathering rate of apatite, the major P-hosting mineral in the rocks, was at least five times higher in the Archean Eon than today, attributable to high levels of pCO2,g. Despite this, the weathering flux of P to the oceans was negligible in the early Archean Eon, increasing to a level comparable to or greater than the modern by the end of eon, a consequence of accelerating continental emergence. Furthermore, our thermodynamic calculations indicate high solubilities of primary and secondary P-hosting minerals in the acidic weathering fluids on land, linking to high Archean pCO2,g. Thus, weathering of P was both kinetically and thermodynamically favorable on the Archean Earth, and river water could transport high levels of dissolved P to the oceans, as also supported by the observed P-depletion in our new compilation of Archean paleosols. Lastly, we evaluated the relative rates of physical erosion and chemical weathering of silicates during the Archean Eon. The results suggest that continental weathering on the early and middle Archean Earth might have been transport-limited due to low erosion rates associated with limited subaerial emergence and low plateau elevations; by the late Archean, however, continental weathering would have transited to kinetically-limited state because of continental emergence, increased plateau elevation, and weakening weathering rates. Overall, our weathering calculations together with paleosol evidence indicate an increasing flux of bioavailable P to the oceans through time, associated with late Archean continental emergence, reaching levels comparable to or higher than modern values by the end of the eon. Increased P fluxes could have fueled increasing rates of primary production, including oxygenic photosynthesis, through time, facilitating the irreversible oxidation of the Earth’s atmosphere early in the Proterozoic Eon.

Phosphorus↗

Errors in Sounding of the Atmosphere Using Broadband Emission Radiometry (SABER) Kinetic Temperature Caused by Non-Local Thermodynamic Equilibrium Model Parameters

The vast set of near global and continuous atmospheric measurements made by the SABER instrument since 2002, including daytime and nighttime kinetic temperature (T(sub k)) from 20 to 105 km, is available to the scientific community. The temperature is retrieved from SABER measurements of the atmospheric 15 micron CO2 limb emission. This emission separates from local thermodynamic equilibrium (LTE) conditions in the rarefied mesosphere and thermosphere, making it necessary to consider the CO2 vibrational state non-LTE populations in the retrieval algorithm above 70 km. Those populations depend on kinetic parameters describing the rate at which energy exchange between atmospheric molecules take place, but some of these collisional rates are not well known. We consider current uncertainties in the rates of quenching of CO2 (v2 ) by N2 , O2 and O, and the CO2 (v2 ) vibrational-vibrational exchange to estimate their impact on SABER T(sub k) for different atmospheric conditions. The T(sub k) is more sensitive to the uncertainty in the latter two and their effects depend on altitude. The T(sub k) combined systematic error due to non-LTE kinetic parameters does not exceed +/- 1.5 K below 95 km and +/- 4-5 K at 100 km for most latitudes and seasons (except for polar summer) if the Tk profile does not have pronounced vertical structure. The error is +/- 3 K at 80 km, +/- 6 K at 84 km and +/- 18 K at 100 km under the less favourable polar summer conditions. For strong temperature inversion layers, the errors reach +/- 3 K at 82 km and +/- 8 K at 90 km. This particularly affects tide amplitude estimates, with errors of up to +/- 3 K.

Garcia-Comas, Maya↗

Test Outline for Flutter Analysis of Rectangular Panels in Rarefied Flow Conditions

Jet plume impingement forces acting on large flexible space structures may precipitate dynamically unstable behavior during space flights. Typical operating conditions in space involve rarefied gas flow regimes which are intrinsically distinct from continuum gas flow and are normally modeled using the kinetic theory of gas flow. Docking and undocking operations of the Space Shuttle with the Russian Mir space laboratory represent a scenario in which the stability boundaries of solar panels may be of interest. Extensive literature review of research work on the dynamic stability of rectangular panels in rarefied gas flow conditions indicated the lack of published reports dealing with this phenomenon. A recently completed preliminary study for NASA JSC dealing with the mathematical analysis of the stability of two-degree-of-freedom elastically supported rigid panels under the effect of rarefied gas flow was reviewed. A test plan outline is prepared for the purpose of conducting a series of experiments on four rectangular rigid test articles in a vacuum chamber under the effect of continuous and pulsating Nitrogen jet plumes. The purpose of the test plan is to gather enough data related to a number of key parameters to allow the validation of the two-degree-of-freedom mathematical model. The hardware required careful design to select a very lightweight material while satisfying rigidity and frequency requirements within the constraints of the test environment. The data to be obtained from the vacuum chamber tests can be compared with the predicted behavior of the theoretical two-degree-of-freedom model. Using the data obtained in this study, further research can identify the limitations of the mathematical model. In addition modifications to the mathematical model can be made, if warranted, to accurately predict the behavior of rigid panels under rarefied gas flow regimes.

Akl, Fred A.↗

Modeling the Uptake of Hydrogen Chloride onto Interior Spacecraft Materials

Hydrogen chloride (HCl) is a major combustion product from the pyrolysis of polyvinyl chloride (PVC) insulated electrical wires, a common spacecraft fire safety concern. Models at two different scales were developed to predict HCl uptake on anodized, chromate conversion coated (Iridite), and bare aluminum surfaces, as well as on Nomex fabric: a macroscopic one-step global surface reaction model where all the active sites are on the exterior surface, and a pore model where the interior active sites deeper into the oxide layer can also be accessed by HCl. Experiments were performed to calibrate kinetic and diffusion constants in the models. A cast acrylic test cell was used to measure the differences between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. For the materials with a thin (< 200 µm) or no oxide layer, the macroscopic surface reaction model adequately predicts the experimental measurements. For the anodized aluminum with a thicker oxide layer, the pore model provided a better match to experimental results. The results will be discussed with respect to the spacecraft fire safety project (Saffire).

modeling↗

Modeling the internal combustion engine

A flexible and computationally economical model of the internal combustion engine was developed for use on large digital computer systems. It is based on a system of ordinary differential equations for cylinder-averaged properties. The computer program is capable of multicycle calculations, with some parameters varying from cycle to cycle, and has restart capabilities. It can accommodate a broad spectrum of reactants, permits changes in physical properties, and offers a wide selection of alternative modeling functions without any reprogramming. It readily adapts to the amount of information available in a particular case because the model is in fact a hierarchy of five models. The models range from a simple model requiring only thermodynamic properties to a complex model demanding full combustion kinetics, transport properties, and poppet valve flow characteristics. Among its many features the model includes heat transfer, valve timing, supercharging, motoring, finite burning rates, cycle-to-cycle variations in air-fuel ratio, humid air, residual and recirculated exhaust gas, and full combustion kinetics.

Zeleznik, F. J.↗

RTM user's guide

RTM is a FORTRAN '77 computer code which simulates the infiltration of textile reinforcements and the kinetics of thermosetting polymer resin systems. The computer code is based on the process simulation model developed by the author. The compaction of dry, woven textile composites is simulated to describe the increase in fiber volume fraction with increasing compaction pressure. Infiltration is assumed to follow D'Arcy's law for Newtonian viscous fluids. The chemical changes which occur in the resin during processing are simulated with a thermo-kinetics model. The computer code is discussed on the basis of the required input data, output files and some comments on how to interpret the results. An example problem is solved and a complete listing is included.

Claus, Steven J.↗

Developing a gas rocket performance prediction technique

A simple, semi-empirical performance correlation/prediction technique applicable to gaseous and liquid propellant rocket engines is presented. Excellent correlations were attained for over 100 test firings by adjusting the computation of the gaseous mixing of an unreactive, coaxial jet using a correlation factor, F, which resulted in prediction of the experimental combustion efficiency for each firing. Static pressure, mean velocity and turbulence intensity in the developing region of non-reactive coaxial jets, typical of those of coaxial injector elements were determined. Detailed profiles were obtained at twelve axial locations (extending from the nozzle exit for a distance of five diameters) downstream from a single element of the Bell Aerospace H2/O2 19-element coaxial injector. These data are compared with analytical predictions made using both eddy viscosity and turbulence kinetic energy mixing models and available computer codes. Comparisons were disappointing, demonstrating the necessity of developing improved turbulence models and computational techniques before detailed predictions of practical coaxial free jet flows are attempted.

Morgenthaler, J. H.↗