Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic stabilization”

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 361 records · Page 20

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗

Modelling of photodegradation in solar cell modules of substrate and superstrate design made with ethylene-vinyl acetate as pottant material

A computer model which simulates, in principle, the chemical changes in the photooxidation of hydrocarbons using as input data a set of elementary reactions, corresponding kinetic rate data and appropriate initial conditions was developed. The Model was refined and exploited to examine more closely the photooxidation and photostabilization of a hydrocarbon polymer. The results lead to the following observations. (1) The time to failure, tau sub f (chosen as the level of 5% C-H bond oxidation which is within the range anticipated for marked change in mechanical properties) varies as the inverse square root of the light intensity. However, tau sub f is almost unaffected by both the photoinitiator type and concentration. (2) The time to failure decreases with the rate of abstraction of C-H by peroxy radicals but increases with the rate of bimolecular radical termination controlled by diffusion. (3) Of the various stabilization mechanisms considered, the trapping of peroxy radicals is distinctly the most effective, although the concommitant decomposition of hydroperoxide is also desirable.

Somersall, A. C.↗

On the balance of stresses in the plasma sheet.

The stress resulting from magnetic tension on the neutral sheet must, in a steady state, be balanced by any one or a combination of (1) a pressure gradient in the direction along the axis of the tail, (2) a similar gradient of plasma flow kinetic energy, and (3) the tension resulting from a pressure anisotropy within the plasma sheet. Stress balance in the first two cases requires that the ratios h/LX and BZ/BX be of the same order of magnitude, where h is the half-thickness of the neutral sheet, LX is the length scale for variations along the axis of the tail, and BZ and BX are the magnetic field components in the plasma sheet just outside the neutral sheet. The second case requires, in addition, that the plasma flow speed within the neutral sheet be of the order of or larger than the Alfven speed outside the neutral sheet. Stress balance in the third case requires that just outside the neutral sheet the plasma pressure obey the marginal firehose stability condition.

Rich, F. J.↗

Proton affinity determinations and proton-bound dimer structure indications in C2 to C15, (alpha),(omega)-alkyldiamines

The 'kinetic method' was used to determine the proton affinity (PA) of a,coalkyldiamines from collision induced dissociation (CID) studies of protonated heterodimers. These PA values were consistently lower than those reported in the proton affinity scale. The apparent discrepancy was rationalized in terms of differences in the conformation of the protonated diamine monomers. The minimum energy species, formed by equilibrium proton transfer processes, have a cyclic conformation and the ion charge is shared by both amino-groups which are bridged by the proton. On the other hand, the species formed through dissociation of protonated dimers have a linear structure and the charge is localized on one of the amino-groups. Thus, the difference in the PA values obtained by both methods is a measure of the additional stability acquired by the protonated diamines through cyclization and charge delocalization. The major collision dissociation pathway of the protonated diamine monomers involved elimination of an ammonia moiety. Other reactions observed included loss of the second amino-group and several other bond cleavages. CID of the protonated dimers involved primarily formation of a protonated monomer through cleavage of the weaker hydrogen bond and subsequently loss of ammonia at higher collision energies. As observed from the CID studies, doubly charged ions were also formed from the diamines under conditions of the electrospray ionization.

Karpas, Z.↗

Microgravity Flammability of PMMA Rods in Concurrent Flow

Microgravity experiments burning cast PMMA cylindrical rods in axial flow have been conducted aboard the International Space Station in the Microgravity Science Glovebox (MSG) facility using the Burning and Suppression of Solids (BASS) flow duct, as part of the BASS-II experiment. Twenty-four concurrent-flow tests were performed, focusing on finding flammability limits as a function of oxygen and flow speed. The oxygen was varied by using gaseous nitrogen to vitiate the working volume of the MSG. The speed of the flow parallel to the rod was varied using a fan at the entrance to the duct. Both blowoff and quenching limits were obtained at several oxygen concentrations. Each experiment ignited the rod at the initially hemispherical stagnation tip of the rod, and allowed the flame to develop and heat the rod at a sufficient flow to sustain burning. For blowoff limit tests, the astronaut quickly turned up the flow to obtain extinction. Complementary 5.18-second Zero Gravity Facility drop tests were conducted to compare blowoff limits in short and long duration microgravity. For quenching tests, the flow was incrementally turned down and the flame allowed to stabilize at the new flow condition for at least the solid-phase response time before changing it again. Quenching was observed when the flow became sufficiently weak that the flame could no longer provide adequate heat flux to compensate for the heat losses (conduction into the rod and radiation). A surface energy balance is presented that shows the surface radiative loss exceeds the conductive loss into the rod near the limit. The flammability boundary is shown to represent a critical Damkohler number, expressed in terms of the reaction rate divided by the stretch rate. For the blowoff branch, the boundary exhibits a linear dependence on oxygen concentration and stretch rate, indicating that the temperature at blowoff must be fairly constant. For the quenching branch, the dominance of the exponential nature of the Arrhenius kinetics reaction rate indicates that the temperature is critical.

flammability↗

Alternative Fuel Research in Fischer-Tropsch Synthesis

NASA Glenn Research Center has recently constructed an Alternative Fuels Laboratory which is solely being used to perform Fischer-Tropsch (F-T) reactor studies, novel catalyst development and thermal stability experiments. Facility systems have demonstrated reliability and consistency for continuous and safe operations in Fischer-Tropsch synthesis. The purpose of this test facility is to conduct bench scale Fischer-Tropsch (F-T) catalyst screening experiments while focusing on reducing energy inputs, reducing CO2 emissions and increasing product yields within the F-T process. Fischer-Tropsch synthesis is considered a gas to liquid process which reacts syn-gas (a gaseous mixture of hydrogen and carbon monoxide), over the surface of a catalyst material which is then converted into liquids of various hydrocarbon chain length and product distributions1. These hydrocarbons can then be further processed into higher quality liquid fuels such as gasoline and diesel. The experiments performed in this laboratory will enable the investigation of F-T reaction kinetics to focus on newly formulated catalysts, improved process conditions and enhanced catalyst activation methods. Currently the facility has the capability of performing three simultaneous reactor screening tests, along with a fourth fixed-bed reactor used solely for cobalt catalyst activation.

Surgenor, Angela D.↗

Aerodynamic levitation of laser-heated solids in gas jets

The aerodynamic levitation technique is developed for studies of high-temperature material properties and gas/condensed-phase reaction kinetics. Stable levitation is demonstrated in a supersonic jet from a 0.081 cm nozzle with 0.03-0.20 g 0.24-0.47 cm diameter solid spheres at a height between 0.7-2.0 cm above the nozzle and ambient pressures between 1.1-18 Torr. A model of supersonic jet levitation is developed which accurately predicts the values of height vs pressure over the full range of conditions investigated. It is found that the efficiency with which jet momentum is converted into levitation force decreases with the jet/specimen diameter ratio and the jet Reynolds number, and the rate of jet spreading with distance from the nozzle is found to agree with that measured by pitot tube traverses of the jet. In addition, laser heating is shown to reduce the jet momentum required for levitation at a given height and to increase levitation stability. Measurements of sphere levitation in subsonic gas jets show that the required jet momentum flow rate exceeds the specimen weight by about 2/the specimen drag coefficient at its terminal free-fall speed under ambient conditions.

Nordine, P. C.↗

Similarity theory of the buoyantly interactive planetary boundary layer with entrainment

A similarity model is developed for the vertical profiles of turbulent flow variables in an entraining turbulent boundary layer of arbitrary buoyant stability. In the general formulation the vertical profiles, internal rotation of the velocity vector, discontinuities or jumps at a capping inversion and bulk aerodynamic coefficients of the boundary layer are given by solutions to a system of ordinary differential equations in the similarity variable. To close the system, a formulation for buoyantly interactive eddy diffusivity in the boundary layer is introduced which recovers Monin-Obukhov similarity near the surface and incorporates a hypothesis accounting for the observed variation of mixing length throughout the boundary layer. The model is tested in simplified versions which depend only on roughness, surface buoyancy, and Coriolis effects by comparison with planetary-boundary-layer wind- and temperature-profile observations, measurements of flat-plate boundary layers in a thermally stratified wind tunnel and observations of profiles of terms in the turbulent kinetic-energy budget of convective planetary boundary layers. On balance, the simplified model reproduced the trend of these various observations and experiments reasonably well, suggesting that the full similarity formulation be pursued further.

Hoffert, M. I.↗

Equiparatition of energy for turbulent astrophysical fluids: Accounting for the unseen energy in molecular clouds

Molecular clouds are observed to be partially supported by turbulent pressure. The kinetic energy of the turbulence is directly measurable, but the potential energy, which consists of magnetic, thermal, and gravitational potential energy, is largly unseen. We have extended previous results on equipartition between kinetic and potential energy to show that it is likely to be a very good approximation in molecular clouds. We have used two separate approaches to demonstrate this result: For small-amplitude perturbations of a static equilibrium, we have used the energy principle analysis of Bernstein et al. (1958); this derivation applies to perturbations of arbitary wavelength. To treat perturbations of a nonstatic equilibrium, we have used the Lagrangian analysis of Dewar (1970); this analysis applies only to short-wavelength perturbations. Both analysis assume conservation of energy. Wave damping has only a small effect on equipartition if the wave frequency is small compared to the neutral-ion collision frequency; for the particular case we considered, radiative losses have no effect on equipartition. These results are then incorporated in a simple way into analyses of cloud equilibrium and global stability. We discuss the effect of Alfvenic turbulence on the Jeans mass and show that it has little effect on the magnetic critical mass.

Zweibel, Ellen G.↗

Low NO(x) potential of gas turbine engines

The purpose is to correlate emission levels of gas turbine engines. The predictions of NO(x) emissions are based on a review of the literature of previous low NO(x) combustor programs and analytical chemical kinetic calculations. Concepts included in the literature review consisted of lean-premixed-prevaporized (LPP), rich burn/quick quench/lean burn (RQL), and direct injection. The NO(x) emissions were found to be an exponential function of adiabatic combustion temperature over a wide range of inlet temperatures, pressures and (lean) fuel-air ratios. A simple correlation of NO(x) formation with time was not found. The LPP and direct injection (using gaseous fuels) concepts have the lowest NO(x) emissions of the three concepts. The RQL data has higher values of NO(x) than the LPP concept, probably due to the stoichiometric temperatures and NO(x) production that occur during the quench step. Improvements in the quick quench step could reduce the NO(x) emissions to the LPP levels. The low NO(x) potential of LPP is offset by the operational disadvantages of its narrow stability limits and its susceptibility to autoignition/flashback. The Rich-Burn/Quick-Quench/Lean-Burn (RQL) and the direct injection concepts have the advantage of wider stability limits comparable to conventional combustors.

Tacina, Robert R.↗

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗

The Oscillations and Stability of Differentially Rotating Neutron Stars

Studies of the oscillations and stability of neutron stars are motivated by the fact that vibrating neutron stars are a promising source of gravitational waves. One important factor is the influence of differential rotation, which is likely to arise in a neutron star at times, such as the immediate aftermath of the supernova, when we expect strong vibrations. I will discuss two phenomena unique to differentially rotating systems: dynamical shear instabilities, and the existence of a co-rotation band (a frequency band in which mode pattern speed matches the local angular velocity). Using a simple model, we have found dynamical shear instabilities that arise where modes cross into the co-rotation band, if the degree of differential rotation exceeds a certain threshold. We are currently investigating whether this mechanism operates in more realistic stellar models, and whether it is responsible for the dynamical instabilities occurring at low ratios of kinetic to potential energy that have been observed by several authors. I will present the latest results of these studies. Another topic of investigation is the nature of oscillations within the co-rotation band. The band gives rise to a continuous spectrum whose collective physical perturbation exhibits complicated temporal behaviour. I will also discuss the existence of modes within the continuous spectrum that appear physically indistinguishable from the discrete modes outside the band, despite the apparently singular nature of their eigenfunctions.

Watts, A.↗

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused convergence problems in some RANS simulations, was also captured in LES / RANS simulations, which were able to accommodate its effects accurately.

Fulton, Jesse A.↗

Orion Landing Simulation Eight Soil Model Comparison

LS-DYNA finite element simulations of a rigid Orion Crew Module (CM) were used to investigate the CM impact behavior on eight different soil models. Ten different landing conditions, characterized by the combination of CM vertical and horizontal velocity, hang angle, and roll angle were simulated on the eight different soils. The CM center of gravity accelerations, pitch angle, kinetic energy, and soil contact forces were the outputs of interest. The simulation results are presented, with comparisons of the CM behavior on the different soils. The soils analyzed in this study can be roughly categorized as soft, medium, or hard, according to the CM accelerations that occur when landing on them. The soft group is comprised of the Carson Sink Wet soil and the Kennedy Space Center (KSC) Low Density Dry Sand. The medium group includes Carson Sink Dry, the KSC High Density In-Situ Moisture Sand and High Density Flooded Sand, and Cuddeback B. The hard soils are Cuddeback A and the Gantry Unwashed Sand. The softer soils were found to produce lower peak accelerations, have more stable pitch behavior, and to be less sensitive to the landing conditions. This investigation found that the Cuddeback A soil produced the highest peak accelerations and worst stability conditions, and that the best landing performance was achieved on the KSC Low Density Dry Sand.

Mark, Stephen D.↗

The development of gamma-gamma-prime lamellar structures in a nickel-base superalloy during elevated temperature mechanical testing

The kinetics of the formation and subsequent development of the directional coarsening of the gamma-prime precipitate in model Ni-Al-Mo-Ta superalloy single crystals are examined during tensile creep under various stress levels at 982 and 1038 C. Special attention is given to the gamma and gamma-prime relation to creep time and strain in order to trace the changing gamma-gamma-prime morphology. Directional coarsening of gamma-prime is found to begin during primary creep and its rate is shown to increase with an increase in temperature or stress level. The length of gamma-prime thickness increased linearly with time up to a plateau reached after the onset of steady state creep. The raft thickness, equal to the gamma-prime size, remained constant at this initial value up through the onset of the tertiary creep. The interlaminar spacing indicates the stability of directionally coarsened structure.

Mackay, R. A.↗

Some minimization calculations concerning the effect of horizontal eddy momentum fluxes on the equilibrium zonal mean motions

A simplified linear two-level zonally averaged atmospheric model is used to determine the minimum zonal kinetic energy (ZKE) plus the zonal available potential energy (ZAPE) of equilibrium states for any distribution of horizontal eddy momentum fluxes. The model assumes fixed thermal forcing except for coefficients of the longwave radiative response to temperature. Some distinctive features of the model solutions are described. It is found that: the states of minimum ZKE or minimum ZKE + ZAPE contain a jet stream; the structure of the zonal winds is dependent on thermal forcing; the horizontal eddy momentum fluxes which minimize ZKE or ZKE + ZAPE can be up-gradient everywhere with respect to the upper level zonal mean zonal winds, or down-gradient depending on the thermal forcing; ZKE minimization results in a jet stream at 54 deg latitude independent of external parameters. The minimization of ZKE + ZAPE results in a jet stream at 35 deg latitude for a planet of the earth's radius and similar static stability.

Schneider, E. K.↗

Convective flow effects on protein crystal growth

The long-term stability of the interferometric setup for the monitoring of protein morphologies has been improved. Growth or dissolution of a crystal on a 100 A scale can now be clearly distinguished from dimensional changes occurring within the optical path of the interferometer. This capability of simultaneously monitoring the local interfacial displacement at several widely-spaced positions on the crystal surface with high local depth resolution, has already yielded novel results. We found with lysozyme that (1) the normal growth rate is oscillatory, and (2) the mean growth step density is greater at the periphery of a facet than in its center. The repartitioning of Na(+) and Cl(-) ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed to interpret the large body of data in a unified way. The results strongly suggests that (1) the ion to lysozyme ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter on the order of 10 microns. The computational model for diffusive-convective transport in protein crystallization (see the First Report) has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies. These results show that some elements of a moving boundary problem must be incorporated into the model in order to obtain a more realistic description. Our experimental setup for light scattering investigations of aggregation and nucleation in protein solutions has been extensively tested. Scattering intensity measurements with a true Rayleigh scatterer produced systematically increased forward scattering, indicating problems with glare. These have been resolved. Preliminary measurements with supersaturated lysozyme solutions revealed that the scatterers grow with time. Work has begun on a computer program for the unified evaluation of simultaneously obtained, multi-angle static and dynamic light scattering data.

Rosenberger, Franz↗