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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Analytic modeling of a spray diffusion flame

A detailed model for a spray diffusion flame is described. The model is based on the boundary layer form of the equations of motion, with droplet transport accounted for using a discretized droplet size distribution function. Interphase transport of mass and energy are accounted for, with a flame-sheet model used to describe the combustion process on a droplet scale. Near dynamic equilibrium is assumed for the description of droplet transport; droplets can diffuse relative to the gas phase. Gas-phase mixing is accounted for using a two-equation turbulence model; buoyancy effects are included, with a temperature fluctuation equation used to account for buoyancy effects on turbulence structure. Thermal radiation from gas-phase CO2 and H2O is included. Gas-phase chemical kinetics are modeled using a 20-reaction, 10-species version of the advanced quasi-global chemical kinetics formulation. Results are compared with data for a vaporizing Freon spray and a pentane spray flame. It is shown that the computational approach provides a reasonably valid picture of the overall development of a spray diffusion flame, and, furthermore, provides a useful tool for the parametric examination of the spray combustion process.

Harsha, P. T.↗

NASA NTRS1987

Hydrogen oxidation mechanism with applications to (1) the chaperon efficiency of carbon dioxide and (2) vitiated air testing

Ignition delay times for the hydrogen/oxygen/carbon dioxide/argon system were obtained behind reflected shock waves. A detailed kinetic mechanism modeled the experimental hydrogen/oxygen data, Skinner and Ringrose's high-pressure data, and Slack and Grillo's hydrogen/air data. A carbon dioxide chaperon efficiency of 7.0 +/- 0.2 was determined. The reaction pathway H2O yields H2O2 yields OH yields H was required to model the high-pressure data. It is suggested that some of the lowest temperature data points (1.0 and 0.5 atm) for Slack and Grillo's hydrogen/air experiments are in error. It was found that the technique of simplifying a detailed kinetic mechanism for a limited range of experimental data may render the model useless for other test conditions.

Brabbs, Theodore A.↗

NASA NTRS2005

AutoChem

AutoChem is a suite of Fortran 90 computer programs for the modeling of kinetic reaction systems. AutoChem performs automatic code generation, symbolic differentiation, analysis, and documentation. It produces a documented stand-alone system for the modeling and assimilation of atmospheric chemistry. Given databases of chemical reactions and a list of constituents defined by the user, AutoChem automatically does the following: 1. Selects the subset of reactions that involve a user-defined list of constituents and automatically prepares a document listing the reactions; 2. Constructs the ordinary differential equations (ODEs) that describe the reactions as functions of time and prepares a document containing the ODEs; 3. Symbolically differentiates the time derivatives to obtain the Jacobian and prepares a document containing the Jacobian; 4. Symbolically differentiates the Jacobian to obtain the Hessian and prepares a document containing the Hessian; and 5. Writes all the required Fortran 90 code and datafiles for a stand-alone chemical modeling and assimilation system (implementation of steps 1 through 5). Typically, the time taken for steps 1 through 5 is about 3 seconds. The modeling system includes diagnostic components that automatically analyze each ODE at run time, the relative importance of each term, time scales, and other attributes of the model.

Lary, David John↗

NASA NTRS2010

Defining Nitrogen Kinetics for Air Break in Prebreath

Actual tissue nitrogen (N2) kinetics are complex; the uptake and elimination is often approximated with a single half-time compartment in statistical descriptions of denitrogenation [prebreathe(PB)] protocols. Air breaks during PB complicate N2 kinetics. A comparison of symmetrical versus asymmetrical N2 kinetics was performed using the time to onset of hypobaric decompression sickness (DCS) as a surrogate for actual venous N2 tension. METHODS: Published results of 12 tests involving 179 hypobaric exposures in altitude chambers after PB, with and without airbreaks, provide the complex protocols from which to model N2 kinetics. DCS survival time for combined control and airbreaks were described with an accelerated log logistic model where N2 uptake and elimination before, during, and after the airbreak was computed with a simple exponential function or a function that changed half-time depending on ambient N2 partial pressure. P1N2-P2 = (Delta)P defined decompression dose for each altitude exposure, where P2 was the test altitude and P1N2 was computed N2 pressure at the beginning of the altitude exposure. RESULTS: The log likelihood (LL) without decompression dose (null model) was -155.6, and improved (best-fit) to -97.2 when dose was defined with a 240 min half-time for both N2 elimination and uptake during the PB. The description of DCS survival time was less precise with asymmetrical N2 kinetics, for example, LL was -98.9 with 240 min half-time elimination and 120 min half-time uptake. CONCLUSION: The statistical regression described survival time mechanistically linked to symmetrical N2 kinetics during PBs that also included airbreaks. The results are data-specific, and additional data may change the conclusion. The regression is useful to compute additional PB time to compensate for an airbreak in PB within the narrow range of tested conditions.

Conkin, Johnny↗

NASA NTRS2009

Defining Nitrogen Kinetics for Air Break in Prebreathe

Actual tissue nitrogen (N2) kinetics are complex; the uptake and elimination is often approximated with a single half-time compartment in statistical descriptions of denitrogenation [prebreathe (PB)] protocols. Air breaks during PB complicate N2 kinetics. A comparison of symmetrical versus asymmetrical N2 kinetics was performed using the time to onset of hypobaric decompression sickness (DCS) as a surrogate for actual venous N2 tension. Published results of 12 tests involving 179 hypobaric exposures in altitude chambers after PB, with and without air breaks, provide the complex protocols from which to model N2 kinetics. DCS survival time for combined control and air breaks were described with an accelerated log logistic model where N2 uptake and elimination before, during, and after the air break was computed with a simple exponential function or a function that changed half-time depending on ambient N2 partial pressure. P1N2-P2 = delta P defined DCS dose for each altitude exposure, where P2 was the test altitude and P1N2 was computed N2 pressure at the beginning of the altitude exposure. The log likelihood (LL) without DCS dose (null model) was -155.6, and improved (best-fit) to -97.2 when dose was defined with a 240 min half-time for both N2 elimination and uptake during the PB. The description of DCS survival time was less precise with asymmetrical N2 kinetics, for example, LL was -98.9 with 240 min half-time elimination and 120 min half-time uptake. The statistical regression described survival time mechanistically linked to symmetrical N2 kinetics during PBs that also included air breaks. The results are data-specific, and additional data may change the conclusion. The regression is useful to compute additional PB time to compensate for an air break in PB within the narrow range of tested conditions.

Conkin, Johnny↗

NASA NTRS1998

A Gas-Kinetic Scheme for Reactive Flows

In this paper, the gas-kinetic BGK scheme for the compressible flow equations is extended to chemical reactive flow. The mass fraction of the unburnt gas is implemented into the gas kinetic equation by assigning a new internal degree of freedom to the particle distribution function. The new variable can be also used to describe fluid trajectory for the nonreactive flows. Due to the gas-kinetic BGK model, the current scheme basically solves the Navier-Stokes chemical reactive flow equations. Numerical tests validate the accuracy and robustness of the current kinetic method.

Lian,Youg-Sheng↗

NASA NTRS1996

Kinetic Space Weather: Toward a Global Hybrid Model of the Polar Ionosphere-Lower Magnetosphere Plasma Transport

During the indicated period of performance, we had a number of publications concerned with kinetic polar ionosphere-lower magnetosphere plasma transport. For the IUGG 1991-4 Quadrennial Report, we reviewed aspects of U.S. accomplishments concerned with polar plasma transport, among other issues. In another review, we examined the computer simulations of multiple-scale processes in space plasmas, including polar plasma outflow and transport. We also examined specifically multiscale processes in ionospheric outflows. We developed a Generalized Semi-Kinetic(GSK) model for the topside-lower magnetosphere which explored the synergistic action of wave heating and electric potentials in the formation of auroral Ion conics, in particular the "pressure cooker" mechanism. We extended the GSK model all the way down to 120 km and applied this code to illustrate the response of the ionosphere- magnetosphere to soft-electron precipitation and convection-driven frictional ion heating, respectively. Later, the convection-driven heating work was extended to a paper for the Journal of Geophysical Research. In addition to the above full published papers, we also presented the first developments of the coupled fluid-semikinetic model for polar plasma transport during this period. The results from a steady-state treatment were presented, with the second presentation being concerned with the effects of photo-electrons on the polar wind, and the first garnering an outstanding student paper award from the American Geophysical Union. We presented the first results from a time-dependent version of this coupled fluid-semikinetic model.

Horwitz, James L.↗

NASA NTRS1993

Stratospheric aircraft exhaust plume and wake chemistry

Progress to date in an ongoing study to analyze and model emissions leaving a proposed High Speed Civil Transport (HSCT) from when the exhaust gases leave the engine until they are deposited at atmospheric scales in the stratosphere is documented. A kinetic condensation model was implemented to predict heterogeneous condensation in the plume regime behind an HSCT flying in the lower stratosphere. Simulations were performed to illustrate the parametric dependence of contrail droplet growth on the exhaust condensation nuclei number density and size distribution. Model results indicate that the condensation of water vapor is strongly dependent on the number density of activated CN. Incorporation of estimates for dilution factors into a Lagrangian box model of the far-wake regime with scale-dependent diffusion indicates negligible decrease in ozone and enhancement of water concentrations of 6-13 times background, which decrease rapidly over 1-3 days. Radiative calculations indicate a net differential cooling rate of the plume about 3K/day at the beginning of the wake regime, with a total subsidence ranging between 0.4 and 1 km. Results from the Lagrangian plume model were used to estimate the effect of repeated superposition of aircraft plumes on the concentrations of water and NO(y) along a flight corridor. Results of laboratory studies of heterogeneous chemistry are also described. Kinetics of HCl, N2O5 and ClONO2 uptake on liquid sulfuric acid were measured as a function of composition and temperature. Refined measurements of the thermodynamics of nitric acid hydrates indicate that metastable dihydrate may play a role in the nucleation of more stable trihydrates PSC's.

Miake-Lye, R. C.↗

NASA NTRS1992

Evolution of a terrestrial magma ocean: Thermodynamics, kinetics, rheology, convection, differentiation

The evolution of an initially totally molten magma ocean is constrained on the basis of analysis of various physical problems in the magma ocean. First of all an equilibrium thermodynamics of the magma ocean is developed in the melting temperature range. The equilibrium thermodynamical parameters are found as functions only of temperature and pressure and are used in the subsequent models of kinetics and convection. Kinematic processes determine the crystal size and also determine a non-equilibrium thermodynamics of the system. Rheology controls all dynamical regimes of the magma ocean. The thermal convection models for different rheological laws are developed for both the laminar convection and for turbulent convection in the case of equilibrium thermodynamics of the multiphase system. The evolution is estimated on the basis of all the above analysis.

Solomatov, V. S.↗

NASA NTRS2001

Alfven Wave Generated Electron Time Dispersion

The results from a model of kinetic Alfven waves which includes varying magnetic field and density show that time-dispersed bursts of auroral electrons can be accelerated by Alfven, wave pulses propagating from the magnetosphere to the ionosphere. The modeled electron signatures have similar energy range and temporal structure to those observed on sounding rockets and satellites suggesting that electron time dispersion is generated by Alfven waves.

Kletzing, C. A.↗

NASA NTRS1983

A study of the kinetic energy generation with general circulation models

The history data of winter simulation by the GLAS climate model and the NCAR community climate model are used to examine the generation of atmospheric kinetic energy. The contrast between the geographic distributions of the generation of kinetic energy and divergence of kinetic energy flux shows that kinetic energy is generated in the upstream side of jets, transported to the downstream side and destroyed there. The contributions from the time-mean and transient modes to the counterbalance between generation of kinetic energy and divergence of kinetic energy flux are also investigated. It is observed that the kinetic energy generated by the time-mean mode is essentially redistributed by the time-mean flow, while that generated by the transient flow is mainly responsible for the maintenance of the kinetic energy of the entire atmospheric flow.

Chen, T.-C.↗

NASA NTRS2023

Multiphysics Modeling of Electrode Dissolution under Hydrodynamic Conditions in Low Conductivity Water

Efficient plug-flow electrochemical reactors are characterized by their ability to effectively convert species with significant ionic strength and/or with the assistance of a well-behaved electrolyte. Because of the classic electroconductivity requirement, electrochemical reactors are not particularly attractive options for applications with low-conductivity conditions, especially for low-pressure feeds. However, these reactors can still find application in dedicated processes such as ultra-clean water conditioning and biological-coupled separation systems. Consequently, a multi-physics model was developed using COMSOL Multiphysics to characterize different reactor-design conditions based on the Tertiary Current Distribution which accounts for the effects of expected variations in electrolyte composition and ionic strength on the electrochemical process, as well as solution resistance and electrode kinetics. The modelling framework presented here employs single-phase laminar fluid flow, the Nernst-Planck equation, the water-based electroneutrality condition, and concentration-dependent overpotentials. This investigation particularly explores the electrolysis process of “anodic dissolution,” in which elemental species dissociate from the solid matrix of the electrode and enter the liquid phase of the electrolyte as soluble ions. The reactor configuration of interest has a rectangular parallel-electrode design and is operated galvanostatically. A series of parametric studies are carried out to inspect the response of the system when hydrodynamic, kinetic, and geometric conditions are changed. The parametric sweep resulted in variable effluent concentrations and system voltage, which reveal underlaying optimization patterns for electrode dissolution in low conductivity water.

Multiphysics↗

NASA NTRS2012

Optimization of KINETICS Chemical Computation Code

NASA JPL has been creating a code in FORTRAN called KINETICS to model the chemistry of planetary atmospheres. Recently there has been an effort to introduce Message Passing Interface (MPI) into the code so as to cut down the run time of the program. There has been some implementation of MPI into KINETICS; however, the code could still be more efficient than it currently is. One way to increase efficiency is to send only certain variables to all the processes when an MPI subroutine is called and to gather only certain variables when the subroutine is finished. Therefore, all the variables that are used in three of the main subroutines needed to be investigated. Because of the sheer amount of code that there is to comb through this task was given as a ten-week project. I have been able to create flowcharts outlining the subroutines, common blocks, and functions used within the three main subroutines. From these flowcharts I created tables outlining the variables used in each block and important information about each. All this information will be used to determine how to run MPI in KINETICS in the most efficient way possible.

planetary atmospheres↗

NASA NTRS1978

Prediction of hydrodynamics and chemistry of confined turbulent methane-air frames in a two concentric tube combustor

A formulation of the governing partial differential equations for fluid flow and reacting chemical species in a two-concentric-tube combustor is presented. A numerical procedure for the solution of the governing differential equations is described and models for chemical-equilibrium and chemical-kinetics calculations are presented. The chemical-equilibrium model is used to characterize the hydrocarbon reactions. The chemical-kinetics model is used to predict the concentrations of the oxides of nitrogen. The combustor considered consists of two coaxial ducts. Concentric streams of gaseous fuel and air enter the inlet duct at one end; the flow then reverses and flows out through the outer duct. Two sample cases with specified inlet and boundary conditions are considered and the results are discussed.

Markatos, N. C.↗

NASA NTRS1989

A multiple-time-scale turbulence model based on variable partitioning of the turbulent kinetic energy spectrum

A multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method is presented. In the model, the effect of the ratio of the production rate to the dissipation rate on eddy viscosity is modeled by use of the multiple-time-scales and a variable partitioning of the turbulent kinetic energy spectrum. The concept of a variable partitioning of the turbulent kinetic energy spectrum and the rest of the model details are based on the previously reported algebraic stress turbulence model. Example problems considered include: a fully developed channel flow, a plane jet exhausting into a moving stream, a wall jet flow, and a weakly coupled wake-boundary layer interaction flow. The computational results compared favorably with those obtained by using the algebraic stress turbulence model as well as experimental data. The present turbulence model, as well as the algebraic stress turbulence model, yielded significantly improved computational results for the complex turbulent boundary layer flows, such as the wall jet flow and the wake boundary layer interaction flow, compared with available computational results obtained by using the standard kappa-epsilon turbulence model.

Kim, S.-W.↗

NASA NTRS2002

Simultaneous Raman-Rayleigh-LIF Measurements and Numerical Modeling Results of a Lifted H2/N2 Turbulent Jet Flame in a Vitiated Coflow

An experimental and numerical investigation is presented of a H2/N2 turbulent jet flame burner that has a novel vitiated coflow. The vitiated coflow emulates the recirculation region of most combustors, such as gas turbines or furnaces. Additionally, since the vitiated gases are coflowing, the burner allows for exploration of recirculation chemistry without the corresponding fluid mechanics of recirculation. Thus the vitiated coflow burner design facilitates the development of chemical kinetic combustion models without the added complexity of recirculation fluid mechanics. Scalar measurements are reported for a turbulent jet flame of H2/N2 in a coflow of combustion products from a lean ((empty set) = 0.25) H2/Air flame. The combination of laser-induced fluorescence, Rayleigh scattering, and Raman scattering is used to obtain simultaneous measurements of the temperature, major species, as well as OH and NO. Laminar flame calculation with equal diffusivity do agree when the premixing and preheating that occurs prior to flame stabilization is accounted for in the boundary conditions. Also presented is an exploratory pdf model that predicts the flame's axial profiles fairly well, but does not accurately predict the lift-off height.

Cabra, R.↗

NASA NTRS2024

Assessing Coulomb Collision Contributions to Photoelectron Deceleration above Earth’s Polar Caps

NASA’s Endurance rocket mission resolved a subtle energy shift in photoelectrons escaping from Earth’s polar caps, from which the first successful measurement of Earth’s ambipolar electrostatic field was made. Preflight models predicted that coulomb collisions with thermal electrons could add an extra energy shift of up to 1 eV. However, no such additional shift was reported in the initial analysis. Here, we identify a previously unreported effect related to the spacecraft wake, which introduced a systematic +0.25eV energy shift that countered and masked the shift from coulomb collisions. After applying an in-flight calibration for this wake effect, we find a total photopeak energy shift consistent with the combined influence of the ambipolar field and the subtler -0.37eV coulomb shift, predicted by the fully kinetic STET model. Further analysis of these recalibrated measurements yields an ambipolar field estimate in close agreement with the initial analysis.

Glyn A Collinson↗

NASA NTRS2003

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗