Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic methods”

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 199 records · Page 11

Computer-aided mechanogenesis of skeletal muscle organs from single cells in vitro

Complex mechanical forces generated in the growing embryo play an important role in organogenesis. Computerized application of similar forces to differentiating skeletal muscle myoblasts in vitro generate three dimensional artificial muscle organs. These organs contain parallel networks of long unbranched myofibers organized into fascicle-like structures. Tendon development is initiated and the muscles are capable of performing directed, functional work. Kinetically engineered organs provide a new method for studying the growth and development of normal and diseased skeletal muscle.

Vanderburgh, Herman H.↗

Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model

The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.

Lipatov, Alexander S.↗

Microstructural Characterization and Modeling of SLM Superalloy 718

Superalloy 718 is an excellent candidate for selective laser melting (SLM) fabrication due to a combination of excellent mechanical properties and workability. Predicting and validating the microstructure of SLM-fabricated Superalloy 718 after potential post heat-treatment paths is an important step towards producing components comparable to those made using conventional methods. At present, obtaining accurate volume fraction and size measurements of gamma-double-prime, gamma-prime and delta precipitates has been challenging due to their size, low volume fractions, and similar chemistries. A technique combining high resolution distortion corrected SEM imaging and with x-ray energy dispersive spectroscopy has been developed to accurately and independently measure the size and volume fractions of the three precipitates. These results were further validated using x-ray diffraction and phase extraction methods and compared to the precipitation kinetics predicted by PANDAT and JMatPro. Discrepancies are discussed in context of materials properties, model assumptions, sampling, and experimental errors.

SEM↗

Comparison of Climatological Planetary Boundary Layer Depth Estimates Using the GEOS-5 AGCM

Planetary boundary layer (PBL) processes, including those influencing the PBL depth, control many aspects of weather and climate and accurate models of these processes are important for forecasting changes in the future. However, evaluation of model estimates of PBL depth are difficult because no consensus on PBL depth definition currently exists and various methods for estimating this parameter can give results that differ by hundreds of meters or more. In order to facilitate comparisons between the Goddard Earth Observation System (GEOS-5) and other modeling and observational systems, seven PBL depth estimation methods are used to produce PBL depth climatologies and are evaluated and compared here. All seven methods evaluate the same atmosphere so all differences are related solely to the definition chosen. These methods depend on the scalar diffusivity, bulk and local Richardson numbers, and the diagnosed horizontal turbulent kinetic energy (TKE). Results are aggregated by climate class in order to allow broad generalizations. The various PBL depth estimations give similar midday results with some exceptions. One method based on horizontal turbulent kinetic energy produces deeper PBL depths in the winter associated with winter storms. In warm, moist conditions, the method based on a bulk Richardson number gives results that are shallower than those given by the methods based on the scalar diffusivity. The impact of turbulence driven by radiative cooling at cloud top is most significant during the evening transition and along several regions across the oceans and methods sensitive to this cooling produce deeper PBL depths where it is most active. Additionally, Richardson number-based methods collapse better at night than methods that depend on the scalar diffusivity. This feature potentially affects tracer transport.

Planetary Boundary Layer Depth↗

Cell Kinetic and Histomorphometric Analysis of Microgravitational Osteopenia: PARE.03B

Previous methods of identifying cells undergoing DNA synthesis (S-phase) utilized H-3 thymidine (3HT) autoradiography. 5-Bromo-2'-deoxyuridine (BrdU) immunohistochemistry is a nonradioactive alternative method. This experiment compared the two methods using the nuclear volume model for osteoblast histogenesis in two different embedding media. Twenty Sprague-Dawley rats were used, with half receiving 3HT (1 micro Ci/g) and the other half BrdU (50 microgram/g). Condyies were embedded (one side in paraffin, the other in plastic) and S-phase nuclei were identified using either autoradiography or immunohistochemistry. The fractional distribution of preosteoblast cell types and the percentage of labeled cells (within each cell fraction and label index) were calculated and expressed as mean q standard error. Chi-Square analysis showed only a minor difference in the fractional distribution of cell types. However, there were significant differences (p less than 0.05) by ANOVA, in the nuclear labeling of specific cell types. With the exception of the less-differentiated A+A'cells, more BrdU label was consistently detected in paraffin than in plastic-embedded sections. In general, more nuclei were labeled with 3H-thymidine than with BrdU in both types of embedding media. Labeling index data (labeled cells/total cells sampled x 100) indicated that BrdU in paraffin, but not plastic gave the same results as 3HT in either embedding method. Thus, we conclude that the two labeling methods do not yield the same results for the nuclear volume model and that embedding media is an important factor whenusing BrdU. As a result of this work, 3HT was chosen for used in the PARE.03 flight experiments.

Roberts, W. Eugene↗

Cell Kinetic and Histomorphometric Analysis of Microgravitational Osteopenia: PARE.03B

Previous methods of identifying cells undergoing DNA synthesis (S-phase) utilized 3H-thymidine (3HT) autoradiography. 5-Bromo-2'-deoxyuridine (BrdU) immunohistochemistry is a nonradioactive alternative method. This experiment compared the two methods using the nuclear volume model for osteoblast histogenesis in two different embedding media. Twenty Sprague-Dawley rats were used, with half receiving 3HT (1 micro-Ci/g) and the other half BrdU (50 micro-g/g). Condyles were embedded (one side in paraffin, the other in plastic) and S-phase nuclei were identified using either autoradiography or immunohistochemistry. The fractional distribution of preosteoblast cell types and the percentage of labeled cells (within each cell fraction and label index) were calculated and expressed as mean +/- standard error. Chi-Square analysis showed only a minor difference in the fractional distribution of cell types. However, there were,significant differences (p less than 0.05) by ANOVA, in the nuclear labeling of specific cell types. With the exception of the less-differentiated A+A' cells, more BrdU label was consistently detected in paraffin than in plastic-embedded sections. In general, more nuclei were labeled with 3H-thymidine than with BrdU in both types of embedding media (Fig 2.). Labeling index data (labeled cells/total cells sampled x 100) indicated that BrdU in paraffin, but not plastic gave the same results as 3HT in either embedding method. Thus, we conclude that the two labeling methods do not yield the same results.

Roberts, W. Eugene↗

Consistent transport coefficients in astrophysics

A consistent theory for dealing with transport phenomena in stellar atmospheres starting with the kinetic equations and introducing three cases (LTE, partial LTE, and non-LTE) was developed. The consistent hydrodynamical equations were presented for partial-LTE, the transport coefficients defined, and a method shown to calculate them. The method is based on the numerical solution of kinetic equations considering Landau, Boltzmann, and Focker-Planck collision terms. Finally a set of results for the transport coefficients derived for a partially ionized hydrogen gas with radiation was shown, considering ionization and recombination as well as elastic collisions. The results obtained imply major changes is some types of theoretical model calculations and can resolve some important current problems concerning energy and mass balance in the solar atmosphere. It is shown that energy balance in the lower solar transition region can be fully explained by means of radiation losses and conductive flux.

Fontenla, Juan M.↗

Structure and length scales of rotational discontinuities

Results are presented of a series of 1D hybrid (kinetic ions and fluid electrons) simulations, using a nondynamic method to form rotational discontinuities (RDs) at an angle of 60 deg between the normal direction and the upstream magnetic field. Ion kinetic effects are found to be important with regard to the stability of the current layer, its thickness and scaling, and its dependence on initial conditions. A detailed examination of the downstream wave train and gradient scale size showed that RDs with alpha1 = +/-180 deg are stable within the framework of the simulation. With proper initialization, electron-sense rotations of -270 deg were found that do not disintegrate immediately but rather appear to consist of a superposition of a solitary wave with -360-deg rotation and a 90-deg ion-sense rotation.

Krauss-Varban, D.↗

Electrode kinetics of oxygen reduction - A theoretical and experimental analysis of the rotating ring-disc electrode method

In order to calculate most of the rate constants for the intermediate formation of H2O2 in the electroreduction of O2 to H2O, the theoretical treatments of the rotating ring-disc electrode method by Damjanovic et al. (1966, 1967), Bagotskii et al. (1968, 1969), and Wroblowa et al. (1976) are modified. Rotating ring-disc electrode experimental data obtained for O2 reduction in Pt in 0.55 M H2SO4 are used to illustrate the calculations of rate constants according to the above theoretical treatments. A simple reaction model as proposed by the first author is consistent with the experimental data. The results indicate that O2 (97 percent) reduces to H2O in a direct four-electron transfer reaction. The adsorption of O2 is probably the rate-determining step in the potential region more negative than 0.5 V vs. reversible hydrogen electrode.

Hsueh, K.-L.↗

Measuring Epoxy-Curing Kinetics

Key reaction parameters estimated from single run. Single DSC Curve used to construct linearized plot from which kinetic reaction parameters are obtained: vertical axis intercept proportional to activation energy, while slope equals reaction order. DSC methods show promise for rapid screening and estimation of cure kinetic parameters of epoxy resins.

Cizmecioglu, M.↗

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

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗

Group-kinetic theory of turbulence

The two phases are governed by two coupled systems of Navier-Stokes equations. The couplings are nonlinear. These equations describe the microdynamical state of turbulence, and are transformed into a master equation. By scaling, a kinetic hierarchy is generated in the form of groups, representing the spectral evolution, the diffusivity and the relaxation. The loss of memory in formulating the relaxation yields the closure. The network of sub-distributions that participates in the relaxation is simulated by a self-consistent porous medium, so that the average effect on the diffusivity is to make it approach equilibrium. The kinetic equation of turbulence is derived. The method of moments reverts it to the continuum. The equation of spectral evolution is obtained and the transport properties are calculated. In inertia turbulence, the Kolmogoroff law for weak coupling and the spectrum for the strong coupling are found. As the fluid analog, the nonlinear Schrodinger equation has a driving force in the form of emission of solitons by velocity fluctuations, and is used to describe the microdynamical state of turbulence. In order for the emission together with the modulation to participate in the transport processes, the non-homogeneous Schrodinger equation is transformed into a homogeneous master equation. By group-scaling, the master equation is decomposed into a system of transport equations, replacing the Bogoliubov system of equations of many-particle distributions. It is in the relaxation that the memory is lost when the ensemble of higher-order distributions is simulated by an effective porous medium. The closure is thus found. The kinetic equation is derived and transformed into the equation of spectral flow.

Tchen, C. M.↗

Shock wave structure using nonlinear model Boltzmann equations

The structure of a strong plane shock wave in a monatomic rarefied perfect gas is one of the simplest problems able to be posed in kinetic theory, and one of the hardest to solve. Its simplicity lies in the absence of solid boundaries, geometrical complications, or internal molecular energy. Its difficulty arises from the great departure of the gas from equilibrium within the shock, which invalidates many of the techniques used successfully elsewhere in kinetic theory. In addition to this theoretical challenge, the modern development of ballistics and hypersonic flight has helped to stimulate extensive theoretical and experimental interest in the shock problem. The experimenters in turn have encountered great difficulties on account of the very small physical dimensions of shocks. In fact, until very recently indeed, any close comparisons of theoretical and experimental shock structure results have been rather unprofitable due to the inadequacies of both theory and experiment. During the last few years this situation has been appreciably improved by development of the Monte Carlo method. This allows idealized 'experiments' to be performed on large computers instead of in wind tunnels, using a known intermolecular force law. The most developed of these methods has been shown to be equivalent theoretically to the Boltzmann equation and to give results which agree extremely closely with measurements of high accuracy. Thus Monte Carlo results not only form the soundest basis for our present theoretical knowledge of shock wave structure, but, for purposes of developing other theories, can also be considered a very valuable experimental resource. However, such results remain very expensive to obtain. In this thesis we develop more economical kinetic theory methods for the approximate prediction of shock structure, and compare our results with those of the Monte Carlo method.

Segal, Ben Maurice↗

Particle/Continuum Hybrid Simulation in a Parallel Computing Environment

The objective of this study was to modify an existing parallel particle code based on the direct simulation Monte Carlo (DSMC) method to include a Navier-Stokes (NS) calculation so that a hybrid solution could be developed. In carrying out this work, it was determined that the following five issues had to be addressed before extensive program development of a three dimensional capability was pursued: (1) find a set of one-sided kinetic fluxes that are fully compatible with the DSMC method, (2) develop a finite volume scheme to make use of these one-sided kinetic fluxes, (3) make use of the one-sided kinetic fluxes together with DSMC type boundary conditions at a material surface so that velocity slip and temperature slip arise naturally for near-continuum conditions, (4) find a suitable sampling scheme so that the values of the one-sided fluxes predicted by the NS solution at an interface between the two domains can be converted into the correct distribution of particles to be introduced into the DSMC domain, (5) carry out a suitable number of tests to confirm that the developed concepts are valid, individually and in concert for a hybrid scheme.

Baganoff, Donald↗

Probing Sensitivity of Discharge Characteristics to Model Selection using Uncertainty Quantification in an aprotic Li-Oxygen Battery

Currently, there are several models in the literature, such as kinetic models, microstructural models, and mass transport models that describe a Li-air battery's discharge behavior. Many of these models are calibrated and tested at low current densities and cannot be easily transferred to high current densities. Even at low current densities, there is no quantitative method for a researcher to choose a reaction kinetic model such as classical Butler-Volmer and its derivatives, and modified Marcus-Hush-Chidsey, a resistance model for lithium peroxide such as electron transport via tunneling or linear resistivity, a surface coverage model (lithium peroxide growth) such as partial coverage or full coverage, and mass transport model (discussed in Ref. [1]). Also, it is time-consuming to test different models at high current density (1C) due to a lack of well-tested models and well-calibrated model parameters. For this presentation, we will develop an analytical model, which acts as a surrogate model for a sophisticated finite element model to predict discharge time and discharge voltage. Next, we use an uncertainty quantifying technique called reduced-order stochastic optimization [2, 3] to determine the uncertainty in model parameters for rate kinetics, lithium peroxide resistivity, and parasitic resistance. Finally, a finite element simulation is performed to determine the error introduced by the surrogate model and its influence on the uncertainty in the model parameters.

M Mehta↗

Composition/Structure/Dynamics of comet and planetary satellite atmospheres

This research program addresses two cases of tenuous planetary atmospheres: comets and Io. The comet atmospheric research seeks to analyze a set of spatial profiles of CN in comet Halley taken in a 7.4-day period in April 1986; to apply a new dust coma model to various observations; and to analyze observations of the inner hydrogen coma, which can be optically thick to the resonance scattering of Lyman-alpha radiation, with the newly developed approach that combines a spherical radiative transfer model with our Monte Carlo H coma model. The Io research seeks to understand the atmospheric escape from Io with a hybrid-kinetic model for neutral gases and plasma given methods and algorithms developed for the study of neutral gas cometary atmospheres and the earth's polar wind and plasmasphere. Progress is reported on cometary Hydrogen Lyman-alpha studies; time-series analysis of cometary spatial profiles; model analysis of the dust comae of comets; and a global kinetic atmospheric model of Io.

Combi, Michael R.↗

An approximate analysis of the diffusing flow in a self-controlled heat pipe.

Constant-density two-dimensional axisymmetric equations are presented for the diffusing flow of a class of self-controlled heat pipes. The analysis is restricted to the vapor space. Condensation of the vapor is related to its mass fraction at the wall by the gas kinetic formula. The Karman-Pohlhausen integral method is applied to obtain approximate solutions. Solutions are presented for a water heat pipe with neon control gas.

Somogyi, D.↗