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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗

NASA NTRS1994

The stabilizing effect of compressibility in turbulent shear flow

Direct numerical simulation of turbulent homogeneous shear flow is performed in order to clarify compressibility effects on the turbulence growth in the flow. The two Mach numbers relevant to homogeneous shear flow are the turbulent Mach number M(t) and the gradient Mach number M(g). Two series of simulations are performed where the initial values of M(g) and M(t) are increased separately. The growth rate of turbulent kinetic energy is observed to decrease in both series of simulations. This 'stabilizing' effect of compressibility on the turbulent energy growth rate is observed to be substantially larger in the DNS series where the initial value of M(g) is changed. A systematic companion of the different DNS cues shows that the compressibility effect of reduced turbulent energy growth rate is primarily due to the reduced level of turbulence production and not due to explicit dilatational effects. The reduced turbulence production is not a mean density effect since the mean density remains constant in compressible homogeneous shear flow. The stabilizing effect of compressibility on the turbulence growth is observed to increase with the gradient Mach number M(g) in the homogeneous shear flow DNS. Estimates of M(g) for the mixing and the boundary layer are obtained. These estimates show that the parameter M(g) becomes much larger in the high-speed mixing layer relative to the high-speed boundary layer even though the mean flow Mach numbers are the same in the two flows. Therefore, the inhibition of turbulent energy production and consequent 'stabilizing' effect of compressibility on the turbulence (over and above that due to the mean density variation) is expected to be larger in the mixing layer relative to the boundary layer in agreement with experimental observations.

Sarkar, S.↗

NASA NTRS1993

Thermal effects in rapid directional solidification - Linear theory

We study the morphological instability of the planar solid/liquid interface for a unidirectionally-solidified dilute binary mixture. We use a model developed by Boettinger et al. (1985, 1986), Aziz (1982), and Jackson et al. (1980), which allows for nonequilibrium effects on the interface through velocity-dependent segregation and attachment kinetics. Two types of instabilities are found in the linear stability analysis: (1) a cellular instability, and (2) an oscillatory instability driven by disequilibrium effects. Merchant and Davis (1990) characterized these instabilities subject to the frozen-temperature approximation (FTA). The present work relaxes the FTA by including the effects of latent heat and the full temperature distribution. Thermal effects slightly postpone the onset of the cellular instability but dramatically postpone the onset of the oscillatory instability; however, the absolute-stability conditions, at which at high speed the cellular and oscillatory instabilities are suppressed, remain unchanged from the FTA.

Huntley, D. A.↗

NASA NTRS2000

Measurements and Modeling of Nitric Oxide Formation in Counterflow, Premixed CH4/O2/N2 Flames

Laser-induced fluorescence (LIF) measurements of NO concentration in a variety of CH4/O2/N2 flames are used to evaluate the chemical kinetics of NO formation. The analysis begins with previous measurements in flat, laminar, premixed CH4/O2/N2 flames stabilized on a water-cooled McKenna burner at pressures ranging from 1 to 14.6 atm, equivalence ratios from 0.5 to 1.6, and volumetric nitrogen/oxygen dilution ratios of 2.2, 3.1 and 3.76. These measured results are compared to predictions to determine the capabilities and limitations of the comprehensive kinetic mechanism developed by the Gas Research Institute (GRI), version 2.11. The model is shown to predict well the qualitative trends of NO formation in lean-premixed flames, while quantitatively underpredicting NO concentration by 30-50%. For rich flames, the model is unable to even qualitatively match the experimental results. These flames were found to be limited by low temperatures and an inability to separate the flame from the burner surface. In response to these limitations, a counterflow burner was designed for use in opposed premixed flame studies. A new LIF calibration technique was developed and applied to obtain quantitative measurements of NO concentration in laminar, counterflow premixed, CH4/O2/N2 flames at pressures ranging from 1 to 5.1 atm, equivalence ratios of 0.6 to 1.5, and an N2/O2 dilution ratio of 3.76. The counterflow premixed flame measurements are combined with measurements in burner-stabilized premixed flames and counterflow diffusion flames to build a comprehensive database for analysis of the GRI kinetic mechanism. Pathways, quantitative reaction path and sensitivity analyses are applied to the GRI mechanism for these flame conditions. The prompt NO mechanism is found to severely underpredict the amount of NO formed in rich premixed and nitrogen-diluted diffusion flames. This underprediction is traced to uncertainties in the CH kinetics as well as in the nitrogen oxidation chemistry. Suggestions are made which significantly improve the predictive capability of the GRI mechanism in near-stoichiometric, rich, premixed flames and in atmospheric-pressure, diffusion flames. However, the modified reaction mechanism is unable to model the formation of NO in ultra-rich, premixed or in high-pressure, nonpremixed flames, thus indicating the need for additional study under these conditions.

Thomsen, D. Douglas↗

NASA NTRS2010

Use of LS-DYNA(Registered TradeMark) to Assess the Energy Absorption Performance of a Shell-Based Kevlar(TradeMark)/Epoxy Composite Honeycomb

The forward and vertical impact stability of a composite honeycomb Deployable Energy Absorber (DEA) was evaluated during a full-scale crash test of an MD-500 helicopter at NASA Langley?s Landing and Impact Research Facility. The lower skin of the helicopter was retrofitted with DEA components to protect the airframe subfloor upon impact and to mitigate loads transmitted to Anthropomorphic Test Device (ATD) occupants. To facilitate the design of the DEA for this test, an analytical study was conducted using LS-DYNA(Registered TradeMark) to evaluate the performance of a shell-based DEA incorporating different angular cell orientations as well as simultaneous vertical and forward impact conditions. By conducting this study, guidance was provided in obtaining an optimum design for the DEA that would dissipate the kinetic energy of the airframe while maintaining forward and vertical impact stability.

Polanco, Michael↗

NASA NTRS1983

Homogeneous nucleation kinetics

Homogeneous nucleation kinetics are rederived in a manner fundamentally similar to the approach of classical nucleation theory with the following modifications and improvements. First, the cluster is a parent phase cluster and does not require energization to the parent state. Second, the thermodynamic potential used to describe phase stability is a continuous function along the pathway of phase decomposition. Third, the kinetics of clustering corresponds directly to the diffusional flux of monomers through the cluster distribution and are formally similar to classical theory with the resulting kinetic equation modified by two terms in the preexponential factor. These terms correct for the influence of a supersaturation dependent clustering within the parent phase and for the influence of an asymmetrical cluster concentration as a function of cluster size at the critical cluster size. Fourth, the supersaturation dependence of the nucleation rate is of the same form as that given by classical nucleation theory. This supersaturation dependence must however be interpreted in terms of a size dependent surface tension. Finally, there are two scaling laws which describe supersaturation to either constant nucleation rate or to the thermodynamically determined physical spinodal.

Rasmussen, D. H.↗

NASA NTRS1983

A kinetic cross-field streaming instability

In a high-beta plasma the so-called modified-two-stream instability, which results from strongly magnetized electrons drifting relative to unmagnetized ions across a homogeneous magnetic field, is misnamed because the mode is highly kinetic, particularly when the relative streaming velocity exceeds the Alfven speed of the plasma. This kinetic cross-field streaming instability is investigated in detail, examining the effect of the electromagnetic terms and the stability boundaries in both low- and high-beta plasmas. An approximate dispersion relation showing the relation of this mode to the whistler is derived and solutions of it are compared with those obtained from the exact dispersion relation. The kinetic mode, unlike the usual modified-two-stream instability, is not stabilized by electromagnetic effects when the relative electron-ion drift speed exceeds the Alfven speed.

Wu, C. S.↗

NASA NTRS1996

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗

NASA NTRS1992

Growth from a hypercooled melt near absolute stability

The stability of a solid-liquid interface in a hypercooled melt is studied, taking into account attachment kinetics, surface energy, and surface energy in the heat balance. There is a basic-state solution with the planar interface moving at constant speed. Linear-stability theory gives a long-wave absolute-stability limit. Near this point a string model is introduced in which a thermal-boundary-layer approximation is used and an evolution equation for the interface is obtained. In a limiting case this interface equation reduces to a Kuromoto-Sivashinsky equation. Comparison with experimental and numerical results are discussed, and a conceptual picture of unconstrained growth for all undercoolings is addressed.

Umantsev, A.↗

NASA NTRS1995

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗

NASA NTRS1983

A linear analysis of the transition curve for the baroclinic annulus

A primitive equation linear model is shown to be capable of quantitatively predicting the neutral curve for a baroclinic annulus. Flows unstable to small nonaxisymmetric perturbations are assumed to lie on the nonaxisymmetric side of a transition curve, while the stable flows are on the axisymmetric side. The stability in the presence of perturbations is tested with linearized hydrostatic primitive equations. Attention is given to the fastest growing eigenmodes, demonstrating that the linear and nonlinear wave have similar structures at a particular point in the nonaxisymmetric regime. All the waves are similar to the Eddy wave except at the extreme to the transition curve, where the waves maintain little structure with height. The eddy kinetic energy is baroclinic with a potential energy suppresed by the large static stability of the basic state. Eddy dissipation is noted to be significant near all boundaries.

Miller, T. L.↗

NASA NTRS1987

A comparative study of rapidly and slowly rotating dynamical regimes in a terrestrial general circulation model

As a preliminary step in the development of a general circulation model for general planetary use, a simplified vesion of thef GISS Model I GCM has been run at various rotation periods to investigate differences between the dynamical regimes of rapidly and slowly rotating planets. To isolate the dynamical processes, the hydrologic cycle is suppressed and the atmosphere is forced with perpetual annual mean solar heating. All other parameters except the rotation period remain fixed at their terrestrial values. Experiments were conducted for rotation periods of 2/3, 1, 2, 4, 8, 16, 64 and 256 days. The results are in qualitative agreement with similar experiments carried out previously with other GCMs and with certain aspects of one Venus GCM simulation. As rotation rate decreases, the energetics shifts from baroclinc to quasi-barotropic when the Rossby radius of deformation reaches planetary scale. The Hadley cell expands poleward and replaces eddies as the primary mode of large-scale heat transport. Associated with this is a poleward shift of the baroclinic zone and jet stream and a reduction of the equator-pole temperature contrast. Midlatitude jet strength peaks at 8 days period, as does the weak positive equatorial zonal wind which occurs at upper levels at all rotation periods. Eddy momentum transport switches from poleward to equatorward at the same period. Tropospheric mean static stability generally increases in the tropics and decreases in midlatitudes as rotation rate decreases, but the global mean static stability is independent of rotation rate. The peak in the eddy kinetic energy spectrum shifts toward lower wavenumbers, reaching wavenumber 1 at a period of 8 days. Implications of these results for the dynamics of Venus and Titan are discussed. Specifically, it is suggested that the extent of low-level convection determines whether the Gierasch mechanism contributes significantly to equatorial superrotation on these planets.

Del Genio, Anthony D.↗

NASA NTRS1987

Prediction of transition on a flat plate under the influence of free-stream turbulence using low-Reynolds-number two-equation turbulence models

The prediction characteristics and capabilities of two popular two-equation low-Reynolds-number turbulence models (Jones, Launder 1972; Lam, Bremhorst 1981) have been evaluated with respect to the prediction of transition of a flat plate under the influence of free-stream turbulence. The sensitivity of the predictions to free-stream turbulence intensity, initial starting location of the calculation, and the assumed initial starting profiles for k and epsilon has been determined and presented. Although both models predict the correct qualitative characteristics of transition, they also exhibit significant quantitative deficiencies with regard to both the predicted location and the length of transition. A modification to the production term in the turbulent kinetic energy equation is proposed which is based on a simple stability criterion and correlated to the free-stream turbulence level. The modification becomes inactive in the fully turbulent regime, but is shown to improve both the qualitative and quantitative characteristics of the transition predictions.

Schmidt, Rodney C.↗

NASA NTRS1993

Structure and length scales of rotational discontinuities

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

Krauss-Varban, D.↗

NASA NTRS1992

Microstructural transitions during containerless processing of undercooled Fe-Ni alloys

In order to identify the possible hierarchy of microstructural pathways and transitions, a systematic evaluation of the microstructural evolution in undercooled Fe-Ni alloys was performed on uniformly processed samples covering seven orders of magnitude in volume. At appropriate undercooling levels, alternate solidification pathways become thermodynamically possible and metastable product structures can result from the operation of competitive solidification kinetics. For thermal history evaluation, a heat flow analysis was applied and tested with large Fe-Ni alloy particles (1 to 3 mm) to assess undercooling potential. Alloy powders (10 to 150 microns), with large liquid undercoolings, were studied under the same composition and processing conditions to evaluate the solidification kinetics and microstructural evolution, including fcc/bcc phase selection and the thermal stability of a retained metastable bcc phase. The identification of microstructural transitions with controlled variations in sample size and composition during containerless solidification processing was used to develop a microstructure map which delineates regimes of structural evolutions and provides a unified analysis of experimental observations in the Fe-Ni system.

Thoma, D. J.↗

NASA NTRS1995

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗

NASA NTRS1985

Future fundamental combustion research for aeropropulsion systems

Physical fluid mechanics, heat transfer, and chemical kinetic processes which occur in the combustion chamber of aeropropulsion systems were investigated. With the component requirements becoming more severe for future engines, the current design methodology needs the new tools to obtain the optimum configuration in a reasonable design and development cycle. Research efforts in the last few years were encouraging but to achieve these benefits research is required into the fundamental aerothermodynamic processes of combustion. It is recommended that research continues in the areas of flame stabilization, combustor aerodynamics, heat transfer, multiphase flow and atomization, turbulent reacting flows, and chemical kinetics. Associated with each of these engineering sciences is the need for research into computational methods to accurately describe and predict these complex physical processes. Research needs in each of these areas are highlighted.

Mularz, E. J.↗

NASA NTRS1985

Future fundamental combustion research for aeropropulsion systems

Mularz, E. J.↗