Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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 55 records · Page 3

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↗

The Atacama Cosmology Telescope: Cosmological Parameters from the 2008 Power Spectrum

We present cosmological parameters derived from the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz and 218 GHz over 296 deg(exp 2) with the Atacama Cosmology Telescope (ACT) during its 2008 season. ACT measures fluctuations at scales 500 < l < 10,000. We fit a model for the lensed CMB, Sunyaev-Zel'dovich (SZ), and foreground contribution to the 148 GHz and 218 GHz power spectra, including thermal and kinetic SZ, Poisson power from radio and infrared point sources, and clustered power from infrared point sources. At l = 3000, about half the power at 148 GHz comes from primary CMB after masking bright radio sources. The power from thermal and kinetic SZ is estimated to be Beta(sub 3000) is identical to 6.8 +/- 2.9 mu K (exp 2), where Beta (sub l) is identical to l(l + 1) C(sub l)/2pi. The IR Poisson power at 148 GHz is Bewta(sub 3000) 7.8 +/- 0.7 muK(exp 2) (C(sub l) = 5.5 +/- 0.5 nK(exp 2)), and a clustered IR component is required with Beta (sub 3000) = 4.6 +/- 0.9 muK(exp 2), assuming an analytic model for its power spectrum shape. At 218 GHz only about 15% of the power, approximately 27 mu K(exp 2), is CMB anisotropy at l = 3000. The remaining 85% is attributed to IR sources (approximately 50% Poisson and 35% clustered), with spectral index alpha = 3.69 +/- 0.14 for flux scaling as S(nu) varies as nu(sup alpha). We estimate primary cosmological parameters from the less contaminated 148 GHz spectrum, marginalizing over SZ and source power. The ACDM cosmological model is a good fit to the data (chi square/dof = 29/46), and ACDM parameters estimated from ACT+Wilkinson Microwave Anisotropy Probe (WMAP) are consistent with the seven-year WMAP limits, with scale invariant n(sub s) = 1 excluded at 99.7% confidence level (CL) (3 sigma). A model with no CMB lensing is disfavored at 2.8 sigma. By measuring the third to seventh acoustic peaks, and probing the Silk damping regime, the ACT data improve limits on cosmological parameters that affect the small-scale CMB power. The ACT data combined with WMAP give a 6 sigma detection of primordial helium, with Y(sub p) = 0.313 +/- 0.044, and a 4 sigma detection of relativistic species, assumed to be neutrinos, with N(sub eff) = 5.3 +/- 1.3 (4.6 +/- 0.8 with BAO+H(sub 0) data). From the CMB alone the running of the spectral index is constrained to be d(sub s) / d ln k = -0,034 +/- 0,018, the limit on the tensor-to-scalar ratio is r < 0,25 (95% CL), and the possible contribution of Nambu cosmic strings to the power spectrum is constrained to string tension G(sub mu) < 1.6 x 10(exp -7) (95% CL),

Dunkley, J.↗

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗

Hydrothermal organic synthesis experiments

Ways in which heat is useful in organic synthesis experiments are described, and experiments on the hydrothermal destruction and synthesis of organic compounds are discussed. It is pointed out that, if heat can overcome kinetic barriers to the formation of metastable states from reduced or oxidized starting materials, abiotic synthesis under hydrothermal conditions is a distinct possibility. However, carefully controlled experiments which replicate the descriptive variables of natural hydrothermal systems have not yet been conducted with the aim of testing the hypothesis of hydrothermal organic systems.

Shock, Everett L.↗

The Effect of Gravity on the Combustion Synthesis of Porous Ceramics and Metal Matrix Composites

Combustion synthesis (self propagating, high temperature synthesis-SHS) is a novel technique that is capable of producing many advanced materials. The ignition temperature (Tig) of such combustion synthesis reactions is often coincident with that of the lowest melting point reactant. The resultant liquid metal wets and spreads around the other solid reactant particles of higher melting points, thereby improving the reactant contact and kinetics, followed by formation of the required compounds. This ignition initiates a combustion propagating wave whose narrow reaction front rapidly travels through the reactants. Since this process is highly exothermic, the heat released by combustion often melts the reactant particles ahead of the combustion front and ignites the adjacent reactant layer, resulting in a self-sustaining reaction. Whenever a fluid phase (liquid or gas) is generated by the reaction system, gravity-driven phenomena can occur. Such phenomena include convective flows of fluid by conventional or unstable convection and settling of the higher density phases. A combustion process is often associated with various kinds of fluid flow. For instance, if the SHS reaction is carried out under inert or reactive gas atmospheres, or a volatile, e.g., B2O3, is deliberately introduced as a reactant, convective flows of the gas will occur due to a temperature gradient existing in the atmosphere when a combustion wave is initiated. The increased gas flow will produce a porous (or expanded) SHS product. Owing to the highly exothermic nature of many SHS reactions, liquid phase(s) can also form before, at, or after the combustion front. The huge temperature gradient at the combustion front can induce convective flows (conventional or unstable) of the liquid phase. Each of these types of convective fluid flow can change the combustion behavior of the synthesizing reaction, and, therefore, the resultant product microstructure. In addition, when two or more phases of different density are produced at or ahead of the propagating combustion front settling of the higher density phase will occur resulting in a non-uniform product microstructure and properties.

Moore, J. J.↗

Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods

A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.

DeBonis, James R.↗

Characterization of Chevron Nozzle Performance

The flow fields from a set of 50.8 mm diameter circular nozzles were characterized using both streamwise and cross-stream 3-component Stereo Particle Image Velocimetry (SPIV). The test matrix of nozzles included 15 different chevron nozzle designs and one baseline, circularly symmetric nozzle. Nozzle jet operating conditions ranged from acoustic Mach numbers of 0.9 to 1.5, with static temperature ratios ranging between 0.84 and 2.7. Detailed surveys of the single jet flows were performed to capture three-dimensional features of the turbulent exhaust jet evolution. Cross-flow planar measurements were obtained at twelve axial locations, ranging from 0.1 to 20 nozzle diameters downstream of the nozzle exit planes. Streamwise measurements, along the jet centerlines, were obtained at ten partially overlapping downstream locations, providing complete axial surveys over a region extending beyond 20 nozzle diameters downstream of the nozzle exit planes. In both optical configurations, the measurement planes were sized to completely capture the fully turbulent jet shear layer growth. The measured three-dimensional mean and turbulent velocity fields, along with computed second order statistics including axial vorticity and turbulent kinetic energy, were evaluated for all test points. Well-defined streamwise vortex structures in the jet shear layers were measured and reported.

Jet Noise↗

Test particle study of minor ions in solar wind turbulence

We perform a parameter study of the temporal evolution of a test particle distribution function in MHD turbulence. The turbulent fields are calculated using a pseudo-spectral method and periodic boundary conditions on a regular grid of 180(exp 3) points, appropriate for incompressible, homogeneous and isotropic turbulence. Initially, the kinetic and the magnetic energy are equal on the average. Both, deterministic and random initial conditions are used, in the former case with zeros of the magnetic field located at grid points, in the latter case located by interpolation between grid points. The evolution of the minor ion distribution function is studied in detail as these turbulent fields evolve, developing strong current and vorticity sheets. Using the full collisionless equation of motion for the test particles, the efficiency of nonlinear interactions can be studied. The results are compared to theoretical predictions and are then discussed in connection with the observations of the dynamical properties of solar wind minor ions derived from in situ observations.

Zurbuchen, Th.↗

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

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

Kim, S.-W.↗

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

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

Kim, S.-W.↗

Thermophysical Property Measurements of Molten Semiconductors in 1-g and Reduced-g Condition

Understanding and controlling the formation kinetics of varieties of crystal imperfections such as point defects, non uniform distribution of doping atoms, and impurity atoms in growing crystals are very important. Theoretical (numerical) modeling of the crystal growth process is an essential step to achieving these objectives. In order to obtain reliable modeling results, input parameters, i.e. various thermophysical parameters, must be accurate. The importance of accurate thermophysical properties of semiconductors in crystal growth cannot be overly emphasized. The total hemispherical emissivity, for instance, has a dramatic impact on the thermal environment. It determines the radiative emission from the surface of the melt which determines to a large extent the profile of the solidified crystal. In order to understand the convection and the turbulence in a melt, viscosity becomes an important parameter. The liquid surface tension determines the shape of the liquid-atmosphere interface near the solid-liquid-atmosphere triple point. Currently used values for these parameters are rather inaccurate, and this program intends to provide more reliable measurements of these thermophysical properties. Thus, the objective of this program is in the accurate measurements of various thermophysical properties which can be reliably used in the modeling of various crystal growth processes. In this program, thermophysical properties of molten semiconductors, such as Si, Ge, Si-Ge, and InSb will be measured as a function of temperature using the High Temperature Electrostatic Levitator at JPL. Each material will be doped by different kinds of impurities at various doping levels. Thermophysical properties which will be measured include: density, thermal expansion coefficient, surface tension, viscosity, specific heat, hemispherical total emissivity, and perhaps electrical and thermal conductivities. Many molten semiconductors are chemically reactive with crucibles. As a result, these dispersed impurities in the melts tend to substantially modify the properties of pure semiconductors. Sample levitation done in a vacuum clearly helps maintain the sample purity. However, in the 1-g environment, all gravity caused effects such as convection, sedimentation and buoyancy are still present in the sample. In addition, large forces needed to levitate a sample in the presence of the gravity can cause additional flows in the melt. The use of the High Temperature Electrostatic Levitator (HTESL) for the present research is a recent development and little is known about the flows induced by the electrostatic forces. In this ground base program, we will define the limits of HTESL technology as various thermophysical properties of molten semiconductors are measured.

Rhim, Won-Kyu↗

Artificial Gravity as a Bone Loss Countermeasure in Simulated Weightlessness

The impact of microgravity on the human body is a significant concern for space travelers. We report here initial results from a pilot study designed to explore the utility of artificial gravity (AG) as a countermeasure to the effects of microgravity, specifically to bone loss. After an initial phase of adaptation and testing, 15 male subjects underwent 21 days of 6 head-down bed rest to simulate the deconditioning associated with space flight. Eight of the subjects underwent 1 h of centrifugation (AG, 1 gz at the heart, 2.5 gz at the feet) each day for 21 days, while 7 of the subjects served as untreated controls (CN). Blood and urine were collected before, during, and after bed rest for bone marker determinations. At this point, preliminary data are available on the first 8 subjects (6 AG, and 2 CN). Comparing the last week of bed rest to before bed rest, urinary excretion of the bone resorption marker n-telopeptide increased 95 plus or minus 59% (mean plus or minus SD) in CN but only 32 plus or minus 26% in the AG group. Similar results were found for another resorption marker, helical peptide (increased 57 plus or minus 0% and 35 plus or minus 13% in CN and AG respectively). Bone-specific alkaline phosphatase, a bone formation marker, did not change during bed rest. At this point, sample analyses are continuing, including calcium tracer kinetic studies. These initial data demonstrate the potential effectiveness of short-radius, intermittent AG as a countermeasure to the bone deconditioning that occurs during bed rest.

Smith, S. M.↗

Small Scale Response and Modeling of Periodically Forced Turbulence

The response of the small scales of isotropic turbulence to periodic large scale forcing is studied using two-point closures. The frequency response of the turbulent kinetic energy and dissipation rate, and the phase shifts between production, energy and dissipation are determined as functions of Reynolds number. It is observed that the amplitude and phase of the dissipation exhibit nontrivial frequency and Reynolds number dependence that reveals a filtering effect of the energy cascade. Perturbation analysis is applied to understand this behavior which is shown to depend on distant interactions between widely separated scales of motion. Finally, the extent to which finite dimensional models (standard two-equation models and various generalizations) can reproduce the observed behavior is discussed.

Bos, Wouter↗

Detailed Chemical Kinetic Modeling of Hydrazine Decomposition

The purpose of this research project is to develop and validate a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. Hydrazine is used extensively in aerospace propulsion, and although liquid hydrazine is not considered detonable, many fuel handling systems create multiphase mixtures of fuels and fuel vapors during their operation. Therefore, a thorough knowledge of the decomposition chemistry of hydrazine under a variety of conditions can be of value in assessing potential operational hazards in hydrazine fuel systems. To gain such knowledge, a reasonable starting point is the development and validation of a detailed chemical kinetic mechanism for gas-phase hydrazine decomposition. A reasonably complete mechanism was published in 1996, however, many of the elementary steps included had outdated rate expressions and a thorough investigation of the behavior of the mechanism under a variety of conditions was not presented. The current work has included substantial revision of the previously published mechanism, along with a more extensive examination of the decomposition behavior of hydrazine. An attempt to validate the mechanism against the limited experimental data available has been made and was moderately successful. Further computational and experimental research into the chemistry of this fuel needs to be completed.

Meagher, Nancy E.↗

Effects of Gas-Phase Radiation and Detailed Kinetics on the Burning and Extinction of a Solid Fuel

This is the first attempt to analyze both radiation and detailed kinetics on the burning and extinction of a solid fuel in a stagnation-point diffusion flame. We present a detailed and comparatively accurate computational model of a solid fuel flame along with a quantitative study of the kinetics mechanism, radiation interactions, and the extinction limits of the flame. A detailed kinetics model for the burning of solid trioxane (a trimer of formaldehyde) is coupled with a narrowband radiation model, with carbon dioxide, carbon monoxide, and water vapor as the gas-phase participating media. The solution of the solid trioxane diffusion flame over the flammable regime is presented in some detail, as this is the first solution of a heterogeneous trioxane flame. We identify high-temperature and low-temperature reaction paths for the heterogeneous trioxane flame. We then compare the adiabatic solution to solutions that include Surface radiation only and gas-phase and surface radiation using a black surface model. The analysis includes discussion of detailed flame chemistry over the flammable regime and, in particular, at the low stretch extinction limit. We emphasize the low stretch regime of the radiatively participating flame, since this is the region representative of microgravity flames. When only surface radiation is included, two extinction limits exist (the blow-off limit, and the low stretch radiative limit), and the burning rate and maximum flame temperatures are lower, as expected. With the inclusion of surface and gas-phase radiation, results show that, while flame temperatures are lower, the burning rate of the trioxane diffusion flame may actually increase at low stretch rate due to radiative feedback from the flame to the surface.

Rhatigan, Jennifer L.↗

Critical analysis of nitramine decomposition data: Activation energies and frequency factors for HMX and RDX decomposition

A summary of a literature review on thermal decomposition of HMX and RDX is presented. The decomposition apparently fits first order kinetics. Recommended values for Arrhenius parameters for HMX and RDX decomposition in the gaseous and liquid phases and for decomposition of RDX in solution in TNT are given. The apparent importance of autocatalysis is pointed out, as are some possible complications that may be encountered in interpreting extending or extrapolating kinetic data for these compounds from measurements carried out below their melting points to the higher temperatures and pressure characteristic of combustion.

Schroeder, M. A.↗

Analysis of two-equation turbulence models for recirculating flows

The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.

Thangam, S.↗

A parametric analysis of radiative structure in aerobrake shock layers

A broad-spectrum version of the NEQAIR code was modified to account for self-absorption and applied to AFE flowfields calculated by the LAURA code with a variety of kinetic models. The resulting radiative fluxes were obtained in a decoupled fashion from the flowfield solver along the vehicle's stagnation streamline. The radiative flux obtained was broken down by causative process to study the radiative structure of the AFE's flowfield for the various kinetic models. In addition, the radiative fluxes for several points on a typical AFE trajectory were analyzed to examine how the radiative structure changes as the vehicle completes its aeropass. Only two radiative processes dominated the stagnation radiative flux, and the flow field conditions near the wal were found to exert considerable influence over the radiative flux to the wall.

Greendyke, Robert B.↗