Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic model data”

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 343 records · Page 19

Numerical Study of Flow Past a Circular Cylinder Using RANS, Hybrid RANS/LES and PANS Formulations

Two multiscale type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(sub epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. Solutions from these models are compared to RANS results and experimental data for a stationary and rotating cylinder. The parameter f(sub k) varies between zero and one and has the characteristic to be equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation example are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The models are compared to RANS results and experimental data for turbulent separated flows (TS) and laminar separated flows (LS) around stationary and rotating cylinders. For a stationary cylinder, the TS case is accurately simulated using the general two-equation k(sub epsilon) turbulence model (eddy-viscosity model). PAB3D accurately predicts the drag coefficient (CD), lift coefficient (CL) and the Strouhal number (St). The LS case was a challenge for the RANS computation with an eddy-viscosity turbulence model. The RANS/LES and PANS performed well and showed marked improvements over the RANS solution. The modified PANS model was the most accurate. For the rotating cylinder, the spin ratio varied from zero to one, and the PANS results were in good agreement with published experimental data. RANS/LES and PANS capture both temporal and spatial fluctuations and produce large-scale structures that do not occur in the RANS simulation. The current results show promise for the capability of PANS in simulating unsteady and complex flow phenomena.

Elmiligui, Alaa A.↗

Prototype Biology-Based Radiation Risk Module Project

Biological effects of space radiation and risk mitigation are strategic knowledge gaps for the Evolvable Mars Campaign. The current epidemiology-based NASA Space Cancer Risk (NSCR) model contains large uncertainties (HAT #6.5a) due to lack of information on the radiobiology of galactic cosmic rays (GCR) and lack of human data. The use of experimental models that most accurately replicate the response of human tissues is critical for precision in risk projections. Our proposed study will compare DNA damage, histological, and cell kinetic parameters after irradiation in normal 2D human cells versus 3D tissue models, and it will use a multi-scale computational model (CHASTE) to investigate various biological processes that may contribute to carcinogenesis, including radiation-induced cellular signaling pathways. This cross-disciplinary work, with biological validation of an evolvable mathematical computational model, will help reduce uncertainties within NSCR and aid risk mitigation for radiation-induced carcinogenesis.

Terrier, Douglas↗

High Reynolds number calculations using the dynamic subgrid-scale stress model

The dynamic subgrid-scale eddy viscosity model was used in the large-eddy simulation of the turbulent flow in a plane channel for Reynolds numbers based on friction velocity and channel half-width ranging between 200 and 2000, a range including values significantly higher than in previous simulations. The computed wall stress, mean velocity, and Reynolds stress profiles compare very well with experimental and direct simulation data. Comparison of higher moments is also satisfactory. Although the grid in the near-wall region is fairly coarse, the results are quite accurate: the turbulent kinetic energy peaks at y(+) of about 12, and the near-wall behavior of the resolved stresses is captured accurately. The model coefficient is o(0.001) in the buffer layer and beyond, where the cutoff wave numbers are in the decaying region of the spectra; in the near-wall region, the cutoff wave numbers are nearer the energy-containing range, and the resolved turbulent stresses become a constant fraction of the resolved stresses. This feature is responsible for the correct near-wall behavior of the model coefficient.

Piomelli, Ugo↗

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES↗

The ultraviolet dayglow. I - Far UV emissions of N and N2

Far ultraviolet rocket spectra of N I and N2 dayglow emissions have been analyzed by using AE-E photoelectron spectra, laboratory-measured excitation cross sections, and photochemical models of atomic nitrogen. A self-consistent picture of both optically thick and thin emission features is found by using a model in which the principal production mechanism for N I 1200-A and 1493-A photons is photodissociative excitation of N2. The aeronomic data require that 50-70% of the excited 4 P atoms produced dissociatively have velocities within the Doppler core of the ambient nitrogen atoms, contrary to the expectation that those atoms are produced with large excess kinetic energy. The derived atomic nitrogen density has a maximum density of 2.7 x 10 to the 7th per cu cm at 170 km, a value that is within 40% of that from recent models of odd nitrogen photochemistry.

Meier, R. R.↗

Investigation of turbulent transport in an axisymmetric sudden expansion

Simultaneous two-component laser velocimeter measurements were made in the incompressible turbulent flowfield following an axisymmetric sudden expansion. Mean velocities, Reynolds stresses, and triple products were measured and are presented at axial positions ranging from x/H = 0.2-14. A balance of the turbulent kinetic energy in the flow was performed. The production, convection, and diffusion of turbulent kinetic energy were computed directly from the experimental data using central differencing. A specially designed correction lens was employed to correct for optical aberrations introduced by the circular tube. This lens system allowed the accurate simultaneous measurement of axial and radial velocities in the test section. The experimental measurements were compared to predictions generated by a code that employed the k-epsilon turbulence model. Agreement was good for mean axial velocities, turbulent kinetic energy, and turbulent shear stresses. However, the modeled turbulent normal stresses where in poor agreement with the measured values. The modeled diffusion of turbulent kinetic energy was underpredicted in the region between the shear layer and the centerline of the flow giving lower values of turbulent kinetic energy downstream of the potential core than measured.

Gould, Richard D.↗

Kinetic Model Development of the Carbon Vapor Deposition Reactor to Predict Performance Versus Design

To enable complete closure of a Sabatier-based atmosphere revitalization system where oxygen recovery is predicated by carbon dioxide reduction to methane and water, there must be a means of recovering the hydrogen from the byproduct methane. Honeywell’s carbon vapor deposition (CVD) reactor is a potential technology that provides a thermochemical means by which methane decomposes at high temperatures (> 1,000 °C) into hydrogen, solid carbon, and a small percentage of hydrocarbons. This advanced methane pyrolysis reactor utilizes a carbon fiber substrate to capture and immobilize the solid pyrolytic carbon rather than allowing continuation toward undesirable soot. The CVD reactor is a highly dynamic system whose characteristic geometry changes while operating. The reactor model was developed to provide predictive performance assessment capabilities at different operating conditions, geometries, and configurations to aid in technology development. Modeling of the CVD reactor considered the multitude of simultaneous and interdependent physicochemical processes that are entailed, including: (1) the homogeneous gas-phase pyrolysis reactions, (2) the heterogeneous deposition reactions, and (3) the transport of gaseous species in the reactor domains. A 1-D CVD reactor model was developed based on a general packed-bed, axial-dispersion plug flow reactor model with a solid substrate balance to account for carbon deposition. The model includes 59 species and 243 reactions constituting the methane pyrolysis reaction network, carbon deposition, and interconversions of oxygenate species. The reactor model was compared to brassboard test data, where it was shown to be able to capture the pertinent physicochemical phenomena as evidenced by its agreement with short-duration run data (i.e., looking at different temperature and flow rate conditions), as well as trend agreement with the brassboard reactor’s long-duration test data. Recommendations are provided on further model development steps as well as its use cases.

Air Revitalization↗

MGS MAG/ER Data Analysis Using a Time and Magnetic Field Dependent Electron Transport Model

The goal of that project was to examine certain details about the dayside electron environment at Mars as seen by the Mars Global Surveyor (MGS) magnetometer/electron reflectometer (MAG/ER) instrument. Specifically, we stated that we would use the Khazanov and Liemohn (K&L) kinetic electron transport model to analyze features in the observations. This code includes a non-uniform magnetic field and time-dependence in the result (different from most other models of this type). It was originally developed for electron motion along field lines in the Earth's magnetosphere (between conjugate ionospheres), and is thus quite appropriate for application to the Mars magnetic field scenario. Numerous code developments were implemented and the Mars version of the K&L model is fully operational. Initial results from this code have focused on the examination of MGS MAG/ER observations in the crustal field region when it is on the dayside. After several presentations at scientific meetings, this study culminated in a JGR publication last year.

Liemohn, Michael W.↗

Computational simulation of corotational vortex merger using 0, 1, and 2 equation turbulence models

A computational study has been conducted to gain a better understanding of the role of turbulence in the merger of corotating vortices and in the decay of a single vortex. The flow field is simulated by using an inviscid flow model and zero-, one-, and two-equation turbulent-flow models. An implicit finite-difference procedure is employed to integrate the unsteady, two-dimensional, incompressible flow equations in a cross-plane. In the zero- and one-equation models, a spatially varying mixing-length model that simulates the streamline curvature effect is used to define an eddy viscosity. The computational merging distance criterion is in good agreement with the experimental one. In the case of the two-equation model, the turbulent kinetic energy and a modified rate of dissipation equation incorporating a curvature correction term are solved. A comparative study of the computational results corresponding to the various flow models and the available experimental data is presented.

Raj, P.↗

Mass spectrometric measurements of the neutral gas composition of the thermosphere and exosphere of Venus

The neutral gas composition and density in the thermosphere of Venus is being measured with a quadrupole mass spectrometer on the Pioneer Venus orbiter. Data are obtained near periapsis once per day approximately 150-250 km above the surface. The principal gases in the thermosphere are CO2, CO, N2, O, N, and He. Atomic oxygen is the major constituent above 155 km on the dayside and also on the nightside up to 180 km when helium becomes the major constituent. The average values of CO2, CO, N2, O, and N remain nearly constant during day and night, but an abrupt change occurs across the terminator from a high dayside value to a low nightside value. The helium density varies in the opposite way, and a distinct bulge was observed at night near the morning terminator. The data have been used as the basis of an empirical model. Large orbit to orbit variations in densities were also observed on the nightside, suggesting perhaps strong turbulent motion in the atmosphere below. Kinetic temperatures inferred from scale heights are approximately 285 K on the dayside and 110 K at night. The average global temperature obtained from the model is 199 K.

Niemann, H. B.↗

Diffusion processes in Al2O3 scales - Void growth, grain growth, and scale growth

The internal microstructure and growth kinetics of Al2O3 scales on Ni-15Cr-13Al (wt percent) are investigated by TEM and analyzed in relation to models of diffusivity. Polished arc-melted specimens were oxidized in 1-atm air at 1100 C for 0.1, 1.0, and 20 hours and ion-thinned for TEM at 100 kV. The frequency distribution of void size and grain size is determined for different oxidation times and scale depths. The kinetics of microvoid growth and of grain and scale growth are plotted and related via simplified models to lattice and grain-boundary oxygen diffusivity, respectively. Good agreement is found between model predictions and data obtained by Oishi and Kingery (1960) on oxygen diffusion in bulk Al2O3. The further implications and limitations of these findings are discssed.

Smialek, J. L.↗

Hydrazine vapor detonations

The detonation velocity and cell widths for hydrazine decomposition were measured over a wide range of temperatures and pressures. The detonation velocity in pure hydrazine was within 5 percent of the calculated C-J velocity. The detonation cell width measurements were interpreted using the Zeldovich-Doering-von Neumann model with a detailed reaction mechanism for hydrazine decomposition. Excellent agreement with experimental data for pure hydrazine was obtained using the empirical relation that detonation cell width was equal to 29 times the kinetically calculated reaction zone length.

Pedley, M. D.↗

Validation of a Computational Fluid Dynamics (CFD) Code for Supersonic Axisymmetric Base Flow

The ability to accurately and efficiently calculate the flow structure in the base region of bodies of revolution in supersonic flight is a significant step in CFD code validation for applications ranging from base heating for rockets to drag for protectives. The FDNS code is used to compute such a flow and the results are compared to benchmark quality experimental data. Flowfield calculations are presented for a cylindrical afterbody at M = 2.46 and angle of attack a = O. Grid independent solutions are compared to mean velocity profiles in the separated wake area and downstream of the reattachment point. Additionally, quantities such as turbulent kinetic energy and shear layer growth rates are compared to the data. Finally, the computed base pressures are compared to the measured values. An effort is made to elucidate the role of turbulence models in the flowfield predictions. The level of turbulent eddy viscosity, and its origin, are used to contrast the various turbulence models and compare the results to the experimental data.

Tucker, P. Kevin↗

Kinetic Roughening and Energetics of Tetragonal Lysozyme Crystal Growth

Lysozyme crystal growth rates over 5 orders of magnitude in range can be described using a layer-by-layer model where growth occurs by 2D nucleation on the crystal surface. Based upon the 2D nucleation model of layer growth, the effective barrier for growth was determined to be gamma = 1.3 plus or minus 0.3 x 10(exp -13) erg/molecule, corresponding to a barrier of 3.2 plus or minus 0.7 k(sub B)T, at 22 C. For solution supersaturation, In c/c(sub eq) greater than or equal to 1.9 plus or minus 0.2, the nucleation model would not predict or consistently estimate the highest observable crystal growth rates. As such, a kinetic roughening hypothesis where crystal growth occurs by a continuous mode was implemented for all growth rate data obtained above In c(sub r)/c(sub eq) greater than or equal to 2. That is, independent of the solution conditions that vary with either buffer pH, temperature or precipitant concentration, crystal growth occurs by the continuous addition of molecules anywhere on the crystal surface, above a roughening solution supersaturation. The energy barrier, E(sub c), for the continuous growth process is determined as 6.1 plus or minus 0.4 x 10(exp -13) erg/molecule or 15 plus or minus 1 k(sub B)T at 22 C.

Gorti, Sridhar↗

Catalyzed Combustion of Bipropellants for Micro-Spacecraft Propulsion

This paper addresses the need to understand the physics and chemistry involved in propellant combustion processes in micro-scale combustors for propulsion systems on micro-spacecraft. These spacecraft are planned to have a mass less than 50 kilograms with attitude control estimated to be in the 10 milli-Newton thrust class. These combustors are anticipated to be manufactured using Micro Electrical Mechanical Systems (MEMS) technology and are expected to have diameters approaching the quenching diameter of the propellants. Combustors of this size are expected to benefit significantly from surface catalysis processes. Miniature flame tube apparatus is chosen for this study because microtubes can be easily fabricated from known catalyst materials and their simplicity in geometry can be used in fundamental simulations for validation purposes. Experimentally, we investigated the role of catalytically active surfaces within 0.4 and 0.8 mm internal diameter microtubes, with special emphases on ignition processes in fuel rich gaseous hydrogen and gaseous oxygen. Flame thickness and reaction zone thickness calculations predict that the diameters of our test apparatus are below the quenching diameter of the propellants in sub-atmospheric tests. Temperature and pressure rise in resistively heated platinum and palladium microtubes was used as an indication of exothermic reactions. Specific data on mass flow versus preheat temperature required to achieve ignition are presented. With a plug flow model, the experimental conditions were simulated with detailed gas-phase chemistry, thermodynamic properties, and surface kinetics. Computational results generally support the experimental findings, but suggest an experimental mapping of the exit temperature and composition is needed.

Schneider, Steven J.↗

Advanced Computational Modeling of Vapor Deposition in a High-pressure Reactor

In search of novel approaches to produce new materials for electro-optic technologies, advances have been achieved in the development of computer models for vapor deposition reactors in space. Numerical simulations are invaluable tools for costly and difficult processes, such as those experiments designed for high pressures and microgravity conditions. Indium nitride is a candidate compound for high-speed laser and photo diodes for optical communication system, as well as for semiconductor lasers operating into the blue and ultraviolet regions. But InN and other nitride compounds exhibit large thermal decomposition at its optimum growth temperature. In addition, epitaxy at lower temperatures and subatmospheric pressures incorporates indium droplets into the InN films. However, surface stabilization data indicate that InN could be grown at 900 K in high nitrogen pressures, and microgravity could provide laminar flow conditions. Numerical models for chemical vapor deposition have been developed, coupling complex chemical kinetics with fluid dynamic properties.

Cardelino, Beatriz H.↗

Advanced Computational Modeling of Vapor Deposition in a High-Pressure Reactor

In search of novel approaches to produce new materials for electro-optic technologies, advances have been achieved in the development of computer models for vapor deposition reactors in space. Numerical simulations are invaluable tools for costly and difficult processes, such as those experiments designed for high pressures and microgravity conditions. Indium nitride is a candidate compound for high-speed laser and photo diodes for optical communication system, as well as for semiconductor lasers operating into the blue and ultraviolet regions. But InN and other nitride compounds exhibit large thermal decomposition at its optimum growth temperature. In addition, epitaxy at lower temperatures and subatmospheric pressures incorporates indium droplets into the InN films. However, surface stabilization data indicate that InN could be grown at 900 K in high nitrogen pressures, and microgravity could provide laminar flow conditions. Numerical models for chemical vapor deposition have been developed, coupling complex chemical kinetics with fluid dynamic properties.

Cardelino, Beatriz H.↗

Comparison of DSMC and CFD Solutions of Fire II Including Radiative Heating

The ability to compute rarefied, ionized hypersonic flows is becoming more important as missions such as Earth reentry, landing high mass payloads on Mars, and the exploration of the outer planets and their satellites are being considered. These flows may also contain significant radiative heating. To prepare for these missions, NASA is developing the capability to simulate rarefied, ionized flows and to then calculate the resulting radiative heating to the vehicle's surface. In this study, the DSMC codes DAC and DS2V are used to obtain charge-neutral ionization solutions. NASA s direct simulation Monte Carlo code DAC is currently being updated to include the ability to simulate charge-neutral ionized flows, take advantage of the recently introduced Quantum-Kinetic chemistry model, and to include electronic energy levels as an additional internal energy mode. The Fire II flight test is used in this study to assess these new capabilities. The 1634 second data point was chosen for comparisons to be made in order to include comparisons to computational fluid dynamics solutions. The Knudsen number at this point in time is such that the DSMC simulations are still tractable and the CFD computations are at the edge of what is considered valid. It is shown that there can be quite a bit of variability in the vibrational temperature inferred from DSMC solutions and that, from how radiative heating is computed, the electronic temperature is much better suited for radiative calculations. To include the radiative portion of heating, the flow-field solutions are post-processed by the non-equilibrium radiation code HARA. Acceptable agreement between CFD and DSMC flow field solutions is demonstrated and the progress of the updates to DAC, along with an appropriate radiative heating solution, are discussed. In addition, future plans to generate more high fidelity radiative heat transfer solutions are discussed.

Liechty, Derek S.↗