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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

An Engergetic Electron Flux Dropout Due to Magnetopause Shadowing on 1 June 2013

We examine the mechanisms responsible for the dropout of energetic electron flux during 31 May to 1 June 2013 using Van Allen Probe (Radiation Belt Storm Probes (RBSP)) electron flux data and simulations with the Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. During the storm main phase, L-shells at RBSP locations are greater than ~8, which are connected to open drift shells. Consequently, diminished electron fluxes were observed over a wide range of energies. The combination of drift shell splitting, magnetopause shadowing, and drift loss all results in butterfly electron pitch angle distributions (PADs) at the nightside. During storm sudden commencement, RBSP observations display electron butterfly PADs over a wide range of energies. However, it is difficult to determine whether there are butterfly PADs during the storm main phase since the maximum observable equatorial pitch angle from RBSP is not larger than ~40° during this period. To investigate the causes of the dropout, the CIMI model is used as a global 4-D kinetic inner magnetosphere model. The CIMI model reproduces the dropout with very similar timing and flux levels and PADs along the RBSP trajectory for 593 keV. Furthermore, the CIMI simulation shows butterfly PADs for 593 keV during the storm main phase. Based on comparison of observations and simulations, we suggest that the dropout during this event mainly results from magnetopause shadowing.

Kang, Suk-Bin↗
NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A Liever↗
NASA NTRSDate not supplied

Demonstration of Capability to Simulate Particle Irregular Shape and Poly-Disperse Mixtures Within Lunar Lander Plume-Surface Interaction

Plume-Surface Interaction (PSI) between lander engine plumes and regolith soil creates hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch is developing simulation tools to offer a predictive PSI capability to NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). The Gas-Granular Flow Solver (GGFS) is the main application tool for coupled gas-particle two-phase flow simulations to predict the range of PSI effects from onset of surface erosion to deep crater formation. GGFS features an Eulerian-Eulerian modeling approach, treating both gas and granular material as interacting continuum phases. Modeling the lunar regolith granular material fluidic characteristics poses special challenges due to complex particle shapes and mixture composition. The lunar regolith is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking jagged particle shapes. Eulerian granular material flow modeling requires closure formulations for the granular material constitutive models (stress, friction, collisional and kinetic energy dissipation, drag, etc.). While closure models for spherical particles are available from particle kinetic theory, closure models for realistic non-spherical particles must be extracted from unit physics Discrete Element Model (DEM) particle interaction simulations and provided in the form of tabular datasets. The effects of particle irregular shape (non-spherical shape factors, angular particle surface roughness, and interlocking features) are simulated by approximating the particle features in the form of grouped elemental spheres to form composite particles in the DEM simulations. The effects of the wide range of regolith mixture particle sizes and the strong effects of the presence of the small particle sizes results in high cohesion and low porosity of the regolith mixture. The range of particle sizes is simulated by binning the particle sizes into an appropriate finite number of particle-size species and solving the problem as a species mixture. Combining these two modeling approaches enables simulations to capture both, the contributions of the irregular particle shape and the particle size distribution. The integration and maturation of the DEM-based constitutive model database generation process and poly-disperse mixture binning approach into the GGFS simulation framework are proceeding under funding by the NASA Game Changing Development program. The status of current capabilities will be presented in comparisons of crater characteristics resulting for spherical and irregular shape particles, and for mono-, bi-, and tri-disperse mixture simulations of Apollo LM plume-surface interaction. The computational results confirm the significance of including the particle shape and mixture effects. Going forward plans for the full implementation of the general poly-disperse regolith modeling capability and maturation towards NASA project application readiness under the GCD program will be presented.

Peter A. Liever↗
NASA NTRS1993

A thermal NO(x) prediction model - Scalar computation module for CFD codes with fluid and kinetic effects

A thermal NO(x) prediction model is developed to interface with a CFD, k-epsilon based code. A converged solution from the CFD code is the input to the postprocessing model for prediction of thermal NO(x). The model uses a decoupled analysis to estimate the equilibrium level of (NO(x))e which is the constant rate limit. This value is used to estimate the flame (NO(x)) and in turn predict the rate of formation at each node using a two-step Zeldovich mechanism. The rate is fixed on the NO(x) production rate plot by estimating the time to reach equilibrium by a differential analysis based on the reaction: O + N2 = NO + N. The rate is integrated in the nonequilibrium time space based on the residence time at each node in the computational domain. The sum of all nodal predictions yields the total NO(x) level.

Mcbeath, Giorgio↗
NASA NTRS1993

Computation of the flow field in an annular gas turbine combustor

The KIVA-II code was modified to calculate the 3D flow field in a typical annular gas turbine combustor. The airblast fuel nozzle, cooling baffle, cooling slots, primary and dilution jets, and effusion cooling (bleed) pads were accounted for in this calculation. The turbulence and combustion were modeled using the k-epsilon model and laminar Arrhenius kinetics, respectively. The fuel was modeled as an evaporating liquid spray. The results illustrate the complicated flow fields present in such combustors. From the results obtained to date it appears that the modified KIVA-II code can be used to study the effects of different annular combustor designs and operating conditions.

Cline, Michael C.↗
NASA NTRS2017

Modeling for Battery Prognostics

For any battery-powered vehicles (be it unmanned aerial vehicles, small passenger aircraft, or assets in exoplanetary operations) to operate at maximum efficiency and reliability, it is critical to monitor battery health as well performance and to predict end of discharge (EOD) and end of useful life (EOL). To fulfil these needs, it is important to capture the battery's inherent characteristics as well as operational knowledge in the form of models that can be used by monitoring, diagnostic, and prognostic algorithms. Several battery modeling methodologies have been developed in last few years as the understanding of underlying electrochemical mechanics has been advancing. The models can generally be classified as empirical models, electrochemical engineering models, multi-physics models, and molecular/atomist. Empirical models are based on fitting certain functions to past experimental data, without making use of any physicochemical principles. Electrical circuit equivalent models are an example of such empirical models. Electrochemical engineering models are typically continuum models that include electrochemical kinetics and transport phenomena. Each model has its advantages and disadvantages. The former type of model has the advantage of being computationally efficient, but has limited accuracy and robustness, due to the approximations used in developed model, and as a result of such approximations, cannot represent aging well. The latter type of model has the advantage of being very accurate, but is often computationally inefficient, having to solve complex sets of partial differential equations, and thus not suited well for online prognostic applications. In addition both multi-physics and atomist models are computationally expensive hence are even less suited to online application An electrochemistry-based model of Li-ion batteries has been developed, that captures crucial electrochemical processes, captures effects of aging, is computationally efficient, and is of suitable accuracy for reliable EOD prediction in a variety of operational profiles. The model can be considered an electrochemical engineering model, but unlike most such models found in the literature, certain approximations are done that allow to retain computational efficiency for online implementation of the model. Although the focus here is on Li-ion batteries, the model is quite general and can be applied to different chemistries through a change of model parameter values. Progress on model development, providing model validation results and EOD prediction results is being presented.

Prognostics↗
NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗
NASA NTRS2011

Modeling the Inner Magnetosphere: Radiation Belts, Ring Current, and Composition

The space environment is a complex system defined by regions of differing length scales, characteristic energies, and physical processes. It is often difficult, or impossible, to treat all aspects of the space environment relative to a particular problem with a single model. In our studies, we utilize several models working in tandem to examine this highly interconnected system. The methodology and results will be presented for three focused topics: 1) Rapid radiation belt electron enhancements, 2) Ring current study of Energetic Neutral Atoms (ENAs), Dst, and plasma composition, and 3) Examination of the outflow of ionospheric ions. In the first study, we use a coupled MHD magnetosphere - kinetic radiation belt model to explain recent Akebono/RDM observations of greater than 2.5 MeV radiation belt electron enhancements occurring on timescales of less than a few hours. In the second study, we present initial results of a ring current study using a newly coupled kinetic ring current model with an MHD magnetosphere model. Results of a dst study for four geomagnetic events are shown. Moreover, direct comparison with TWINS ENA images are used to infer the role that composition plays in the ring current. In the final study, we directly model the transport of plasma from the ionosphere to the magnetosphere. We especially focus on the role of photoelectrons and and wave-particle interactions. The modeling methodology for each of these studies will be detailed along with the results.

Glocer, Alex↗
NASA NTRS2009

Analysis of Electromagnetic Wave Propagation in a Magnetized Re-Entry Plasma Sheath Via the Kinetic Equation

Based on a theoretical model of the propagation of electromagnetic waves through a hypersonically induced plasma, it has been demonstrated that the classical radiofrequency communications blackout that is experienced during atmospheric reentry can be mitigated through the appropriate control of an external magnetic field of nominal magnitude. The model is based on the kinetic equation treatment of Vlasov and involves an analytical solution for the electric and magnetic fields within the plasma allowing for a description of the attendant transmission, reflection and absorption coefficients. The ability to transmit through the magnetized plasma is due to the magnetic windows that are created within the plasma via the well-known whistler modes of propagation. The case of 2 GHz transmission through a re-entry plasma is considered. The coefficients are found to be highly sensitive to the prevailing electron density and will thus require a dynamic control mechanism to vary the magnetic field as the plasma evolves through the re-entry phase.

Manning, Robert M.↗
NASA NTRS2018

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

This work presents a general framework for model reduction of non-equilibrium energy transfer and dissociation processes. The multi-group maximum-entropy method is coupled with the quasi-classical trajectory method to directly construct a reduced order model for chemical non-equilibrium. Kinetic data is calculated by applying detailed balance at a microscopic level, overcoming the limitation of calculating recombination kinetic data. This approach enables the construction of a reduced order model for kinetics which bypasses the need to compute state-to-state kinetic data. This physics based reduced order model ensures that an equilibrium distribution is reached given infinite time, while allowing for non-equilibrium distributions during the relaxation and dissociation processes. A proof-of-concept test case demonstrates the applicability of this model by comparison with state-to-state kinetic data for the N2(Chi (sup 1) Sigma ((sub g), (sup plus sign)) minus N ((exp 4) S (sub u)) system with excellent agreement.

Maximum Entropy↗
NASA NTRS1982

Integrated and spectral energetics studies of the GLAS general circulation model

Integrated and spectral energetics of the Goddard Laboratory for Atmospheric Sciences (GLAS) general circulation model are compared with observations and examined when subdivided beyond hemispheric integrals. In the monthly mean zonal averages, qualitative improvements are found over previous versions of the model in eddy kinetic energy and barotropic conversions although vertical shear above the subtropical jet remains weak. This paper traces this problem to the erroneous growth of kinetic energy above strong jets, shows that this difficulty is common to all orographic general circulation models, and suggests that this problem has significant effects on upper-level long-wave predictability. In particular, it is considered whether hemispherically integrated kinetic energy or low-wavenumber spectral coefficients retain any skill as the predictability limit. In the time-dependent comparisons, the model shows significant skill in predicting the hemispherically integrated eddy kinetic energy out to two weeks for one of two cases.

Tenenbaum, J.↗
NASA NTRS2001

A stochastic differential equation model of diurnal cortisol patterns

Circadian modulation of episodic bursts is recognized as the normal physiological pattern of diurnal variation in plasma cortisol levels. The primary physiological factors underlying these diurnal patterns are the ultradian timing of secretory events, circadian modulation of the amplitude of secretory events, infusion of the hormone from the adrenal gland into the plasma, and clearance of the hormone from the plasma by the liver. Each measured plasma cortisol level has an error arising from the cortisol immunoassay. We demonstrate that all of these three physiological principles can be succinctly summarized in a single stochastic differential equation plus measurement error model and show that physiologically consistent ranges of the model parameters can be determined from published reports. We summarize the model parameters in terms of the multivariate Gaussian probability density and establish the plausibility of the model with a series of simulation studies. Our framework makes possible a sensitivity analysis in which all model parameters are allowed to vary simultaneously. The model offers an approach for simultaneously representing cortisol's ultradian, circadian, and kinetic properties. Our modeling paradigm provides a framework for simulation studies and data analysis that should be readily adaptable to the analysis of other endocrine hormone systems.

NASA Discipline Space Human Factors↗
NASA NTRS1979

Analysis of population mortality kinetics with application to the longevity followup of the Navy's '1,000 aviators'

A simple physiological model of mortality kinetics is used to assess the intuitive concept that the aging rates of populations are proportional to their mortality rates. It is assumed that the vitality of an individual can be expressed as a simple summation of the weighted functional capacities of its organs and homeostatic systems that are indispensable for survival. It is shown that the mortality kinetics of a population can be derived by a linear transformation of the frequency distribution of vitality, assuming a uniform constant rate of decline of the physiological functions. A simple comparison of two populations is not possible when they have different vitality frequency distributions. Analysis of the data using the model suggests that the differences in decline of survivorship with age between the military pilot population, a medically insured population, and the control population can be accounted for by the effect of physical selection on the vitality frequency distribution of the screened populations.

Economos, A. C.↗
NASA NTRS2002

High Temperature Corrosion of Silicon Carbide and Silicon Nitride in Water Vapor

Silicon carbide (SiC) and silicon nitride (Si3N4) are proposed for applications in high temperature combustion environments containing water vapor. Both SiC and Si3N4 react with water vapor to form a silica (SiO2) scale. It is therefore important to understand the durability of SiC, Si3N4 and SiO2 in water vapor. Thermogravimetric analyses, furnace exposures and burner rig results were obtained for these materials in water vapor at temperatures between 1100 and 1450 C and water vapor partial pressures ranging from 0.1 to 3.1 atm. First, the oxidation of SiC and Si3N4 in water vapor is considered. The parabolic kinetic rate law, rate dependence on water vapor partial pressure, and oxidation mechanism are discussed. Second, the volatilization of silica to form Si(OH)4(g) is examined. Mass spectrometric results, the linear kinetic rate law and a volatilization model based on diffusion through a gas boundary layer are discussed. Finally, the combined oxidation and volatilization reactions, which occur when SiC or Si3N4 are exposed in a water vapor-containing environment, are presented. Both experimental evidence and a model for the paralinear kinetic rate law are shown for these simultaneous oxidation and volatilization reactions.

Opila, E. J.↗
NASA NTRS1982

Integrated and spectral energetics of the GLAS general circulation model

Integrated and spectral error energetics of the GLAS General circulation model are compared with observations for periods in January 1975, 1976, and 1977. For two cases the model shows significant skill in predicting integrated energetics quantities out to two weeks, and for all three cases, the integrated monthly mean energetics show qualitative improvements over previous versions of the model in eddy kinetic energy and barotropic conversions. Fundamental difficulties remain with leakage of energy to the stratospheric level, particularly above strong initial jet streams associated in part with regions of steep terrain. The spectral error growth study represents the first comparison of general circulation model spectral energetics predictions with the corresponding observational spectra on a day by day basis. The major conclusion is that eddy kinetics energy can be correct while significant errors occur in the kinetic energy of wavenumber 3. Both the model and observations show evidence of single wavenumber dominance in eddy kinetic energy and the correlation of spectral kinetics and potential energy.

Tenenbaum, J.↗
NASA NTRS1996

Rejoining and misrejoining of radiation-induced chromatin breaks. II. Biophysical Model

A biophysical model for the kinetics of the formation of radiation-induced chromosome aberrations is developed to account for the recent experimental results obtained with a combination of the premature chromosome condensation (PCC) and fluorescence in situ hybridization (FISH) techniques. In this model, we consider the broken ends of DNA double-strand breaks (DSBs) to be reactant and make use of the interaction distance hypothesis. The repair/misrepair process between broken ends is suggested to consist of two steps; the first step represents the two break ends approaching each other, and the second step represents the enzymatic processes leading to DNA end-to-end rejoining. Only the second step is reflected in the kinetics observed in experiments using PCC. The model appears to be able to fit existing data for human cells. It is shown that the kinetics of the formation of chromosome aberrations can be explained by a single rate that characterizes both rejoining and misrejoining of DSBs, suggesting that repair and misrepair share the same mechanism. Fast repair (completed in minutes) in a subset of DSBs is suggested as an explanation of the complete exchanges observed with PCC in human lymphocytes immediately after irradiation. The fast repair component seems to be absent in human fibroblasts.

NASA Center JSC↗