Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Kinetic inverse 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 19 records

On the analysis of para-ammonia observations

The intensities and optical depths of the (1, 1), (2, 2), and (2, 1) inversion transitions of ammonia can be calculated quite accurately without solving the equations of statistical equilibrium. A two-temperature partition function suffices. The excitation of the K-ladders can be approximated by using a temperature obtained from a two-level model with the (2, 1) and (1, 1) levels. Distribution of populations between the ladders is described with the kinetic temperature. This enables one to compute the (1, 1) and (2, 1) inversion transition excitation temperatures and optical depths. To compute the (2, 2) brightness temperatures, the fractional population of the (2, 2) doublet is computed from the population of the (1, 1) doublet using the 'true rotation temperature,' which is calculated using a three-level model with the (2, 1), (2, 2), and (1, 1) levels. In spite of some iterative steps, the calculation is quite fast.

Kuiper, T. B. H.↗

Inverse energy cascades in three-dimensional turbulence

Fully three-dimensional magnetohydrodynamic (MHD) turbulence at large kinetic and low magnetic Reynolds numbers is considered in the presence of a strong uniform magnetic field. It is shown by numerical simulation of a model of MHD that the energy inverse cascades to longer length scales when the interaction parameter is large. While the steady-state dynamics of the driven problem is three-dimensional in character, the behavior has resemblance to two-dimensional hydrodynamics. These results have implications in turbulence theory, MHD power generator, planetary dynamos, and fusion reactor blanket design.

Hossain, Murshed↗

Estimation of Lower-Body Kinetics from Loading Profile and Kinematics Alone, Without Measured Ground Reaction Forces

Biomechanical models of human motion can estimate kinetic outcomes, such as joint moments, joint forces and muscle forces. Typically, one performs an inverse dynamics (ID) analysis to compute joint moments from joint angles and measured external forces. Sometimes it is impractical to measure ground reaction forces and moments (GRF&M). We devised an empirical method for performing ID analysis of resistance exercises without measured GRF&M. The method solves the multibody dynamics equations of motion with four key assumptions about the GRF&M that reduce the number of unknowns. The assumptions are 1) negligible ground reaction moments, 2) fixed lateral/medial location of the center of pressure (COP), 3) equal fore/aft location of the COP between the feet, and 4) constant angle of the GRF vector relative to the vertical axis in the frontal plane. We used evaluation trials from a spaceflight countermeasure resistance training device to test this approach. Four participants performed squat and deadlift exercises at various loads. We compared results from traditional ID analysis to results without measured GRF&M using our method. We found that joint moment trajectories in the sagittal plane were qualitatively similar in shape between the two methods, and the amount of root mean squared error (RMSE), measured by difference in joint moment impulse, was typically under 10 percent. Non-sagittal joint moment trajectories, which are much lower in overall magnitude, were not qualitatively similar in shape between the two methods. Non-sagittal moments displayed much higher RMSE, with typical values well over 50 percent. These findings were further supported by validation metrics (Sprague and Geers' P and M metrics, Pearson's r correlation coefficient). Based on these findings, we concluded that useful kinetic results are obtained from ID analysis of squat and deadlift exercises, even when GRF&M are not measured, as long as the outcomes of interest lie in the sagittal plane.

computational modeling↗

Correlation effects in the theory of combined Doppler and pressure broadening. I - Classical theory

An investigation is conducted of the combined effects of radiator-perturber collisions and radiator translational motion in the context of foreign gas broadening of optical transitions in neutral radiators. Questions concerning the speed-dependent collision frequency are considered and aspects of general theory are explored, taking into account the correlation function, the ensemble average, and the kinetic equation formalism. An elementary solution is discussed along with a one-perturber approximation, inverse power law model calculations, and a comparison with the Voigt profile.

Ward, J.↗

Vibrational kinetics in CO electric discharge lasers - Modeling and experiments

A model of CO laser vibrational kinetics is developed, and predicted vibrational distributions are compared with measurements. The experimental distributions were obtained at various flow locations in a transverse CW discharge in supersonic (M = 3) flow. Good qualitative agreement is obtained in the comparisons, including the prediction of a total inversion at low discharge current densities. The major area of discrepancy is an observed loss in vibrational energy downstream of the discharge which is not predicted by the model. This discrepancy may be due to three-dimensional effects in the experiment which are not included in the model. Possible kinetic effects which may contribute to vibrational energy loss are also examined.

Stanton, A. C.↗

Conversion of magnetic field energy into kinetic energy in the solar wind

The outflow of the solar magnetic field energy (the radial component of the Poynting vector) per steradian is inversely proportional to the solar wind velocity. It is a decreasing function of the heliocentric distance. When the magnetic field effect is included in the one-fluid model of the solar wind, the transformation of magnetic field energy into kinetic energy during the expansion process increases the solar wind velocity at 1 AU by 17 percent.

Whang, Y. C.↗

Measurement techniques for analysis of fission fragment excited gases

Spectroscopic analysis of fission fragment excited He, Ar, Xe, N2, Ne, Ar-N2, and Ne-N2 have been conducted. Boltzmann plot analysis of He, Ar and Xe have indicated a nonequilibrium, recombining plasma, and population inversions have been found in these gases. The observed radiating species in helium have been adequately described by a simple kinetic model. A more extensive model for argon, nitrogen and Ar-N2 mixtures was developed which adequately describes the energy flow in the system and compares favorably with experimental measurements. The kinetic processes involved in these systems are discussed.

Schneider, R. T.↗

Maxwell's second- and third-order equations of transfer for non-Maxwellian gases

Condensed algebraic forms for Maxwell's second- and third-order equations of transfer are developed for the case of molecules described by either elastic hard spheres, inverse-power potentials, or by Bird's variable hard-sphere model. These hardly reduced, yet exact, equations provide a new point of origin, when using the moment method, in seeking approximate solutions in the kinetic theory of gases for molecular models that are physically more realistic than that provided by the Maxwell model. An important by-product of the analysis when using these second- and third-order relations is that a clear mathematical connection develops between Bird's variable hard-sphere model and that for the inverse-power potential.

Baganoff, D.↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

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↗

Universal Responses of Cyclic-Oxidation Models Studied

Oxidation is an important degradation process for materials operating in the high-temperature air or oxygen environments typical of jet turbine or rocket engines. Reaction of the combustion gases with the component material forms surface layer scales during these oxidative exposures. Typically, the instantaneous rate of reaction is inversely proportional to the existing scale thickness, giving rise to parabolic kinetics. However, more realistic applications entail periodic startup and shutdown. Some scale spallation may occur upon cooling, resulting in loss of the protective diffusion barrier provided by a fully intact scale. Upon reheating, the component will experience accelerated oxidation due to this spallation. Cyclic-oxidation testing has, therefore, been a mainstay of characterization and performance ranking for high-temperature materials. Models simulate this process by calculating how a scale spalls upon cooling and regrows upon heating (refs. 1 to 3). Recently released NASA software (COSP for Windows) allows researchers to specify a uniform layer or discrete segments of spallation (ref. 4). Families of model curves exhibit consistent regularity and trends with input parameters, and characteristic features have been empirically described in terms of these parameters. Although much insight has been gained from experimental and model curves, no equation has been derived that can describe this behavior explicitly as functions of the key oxidation parameters.

Smialek, James L.↗

Group-kinetic theory and modeling of atmospheric turbulence

A group kinetic method is developed for analyzing eddy transport properties and relaxation to equilibrium. The purpose is to derive the spectral structure of turbulence in incompressible and compressible media. Of particular interest are: direct and inverse cascade, boundary layer turbulence, Rossby wave turbulence, two phase turbulence; compressible turbulence, and soliton turbulence. Soliton turbulence can be found in large scale turbulence, turbulence connected with surface gravity waves and nonlinear propagation of acoustical and optical waves. By letting the pressure gradient represent the elementary interaction among fluid elements and by raising the Navier-Stokes equation to higher dimensionality, the master equation was obtained for the description of the microdynamical state of turbulence.

Tchen, C. M.↗

Relativistic particles and gamma-ray in quasars and active galactic nuclei

A model for a class of quasars and active galactic nuclei is described in which a shock around a massive black hole randomizes the infall kinetic energy of spherically accreting matter producing a nonthermal spectrum of high energy protons. These protons may be responsible for the secondary production (via tau + or - decay) of the radio emitting high energy electrons and also of high energy gamma rays (via pi decay and inverse Compton interactions of the electrons). The correlation between radio and gamma ray emission implied by the model is in good agreement with observations of 3C273. Observation of the flux of high energy neutrinos from quasars may provide a test for the model.

Protheroe, R. J.↗

Relativistic particles and gamma-rays in quasars and active galactic nuclei

A model for a class of quasars and active galactic nuclei is described in which a shock around a massive black hole randomizes the infall kinetic energy of spherically accreting matter producing a nonthermal spectrum of high energy protons. These protons may be responsible for the secondary production (via tau + or - decay) of the radio emitting high energy electrons and also of high energy gamma rays (via Pi decay and inverse Compton interactions of the electrons). The correlation between radio and gamma ray emission implied by the model is in good agreement with observations of 3C273. Observation of the flux of high energy neutrinos from quasars may provide a test for the model.

Protheroe, R. J.↗

Implementation of a kappa-epsilon turbulence model to RPLUS3D code

The RPLUS3D code has been developed at the NASA Lewis Research Center to support the National Aerospace Plane (NASP) project. The code has the ability to solve three dimensional flowfields with finite rate combustion of hydrogen and air. The combustion process of the hydrogen-air system are simulated by an 18 reaction path, 8 species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though a fully implicit numerical algorithm is used. A k-epsilon turbulence model has recently been incorporated into the code. Initial validations have been conducted for a flow over a flat plate. Results of the validation studies are shown. Some difficulties in implementing the k-epsilon equations to the code are also discussed.

Chitsomboon, Tawit↗

Turbulence modeling in supersonic combusting flows

To support the National Aerospace Plane project, the RPLUS3D CFD code has been developed at NASA Lewis. The code has the ability to solve three-dimensional flowfields with finite rate combustion of hydrogen and air. The combustion processes of the hydrogen-air system are simulated by an 18-reaction path, 8-species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though it uses a fully implicit numerical algorithm. A k-epsilon (two equation) turbulence model is incorporated into the RPLUS3D code.

Chitsomboon, Tawit↗

On the origin of relativistic particles and gamma-rays in quasars

A model for a class of quasars and active galactic nuclei is proposed in which a shock around a massive black hole randomizes the infall kinetic energy of spherically accreting matter producing a nonthermal spectrum of high energy protons. It is suggested that these protons are responsible for the secondary production (via proton decay) of the radio emitting high energy electrons and also of high energy gamma-rays (via proton decay and inverse Compton interactions of the electrons). The correlation between radio and gamma-ray emission implied by the model is in good agreement with observations of 3C273. Observation of the flux of high energy neutrinos from quasars may provide a test for the model.

Protheroe, R. J.↗

Numerical modeling of D-mappings with applications to chemical kinetics

Numerical modeling of D-mappings was studied and applied to solving nonlinear stiff systems. These mappings were locally linearized for convergence analysis, and some applications were made to chemical kinetics. The technique avoids using multistep implicit codes that require inversion of Jacobian matrices, but depends on the Jacobians for its convergence analysis.

Dey, S. K.↗