Prediction of turbulent diffusion flames with a four-equation turbulence model

For turbulent flame prediction a four-equation turbulence model consisting of the turbulent kinetic energy, the dissipation rate and the density-velocity correlations using unweighted statistics is suggested. The fluctuating concentrations are taken into account by using Spalding's (1971) equation for the variance of the mixture fraction f and by defining the probability density function of f as beta-function. The thermodynamic model for the combustion of H₂ with air assumes infinitely fast reaction in a global single step. Preliminary results are shown for the calculation of Kent and Bilger (1972) turbulent H₂-air diffusion flame. For this flame with fluctuating density the inclusion of the velocity-density correlations appears to be essential, because they have a strong effect especially on the development of the turbulent viscosity and the mean velocity and mass fraction profiles. Consequently the mean concentrations are in good agreement with the experimental results.     

An experimental study of flame turbulence

An experimental study of the turbulence generated by an enclosed, premixed, propane-air flame has been carried out in a combustion chamber of 25 × 20 cm cross section. Care was taken to reduce any effects of the axial pressure gradient. By suitable changes in the grid geometry, the turbulence intensity and scale of the approach flow were varied independently. The results of these experiments show that a strong link exists between the mechanisms of turbulent flame propagation and flame-generated turbulence. Thus three distinct regions may be identified, each having different characteristics in regard to the effects of turbulence scale on flame-generated turbulence. For each region, a physical mechanism for flame-generated turbulence is proposed. In particular, it is observed that over a wide range of intensity and scale of the approach turbulence, (a) the relative turbulence intensity in the flame zone varies in the range 1–2 times the relative turbulence intensity of the cold flow, (b) in the region of intermediate turbulence levels ( ) the flame-generated turbulence intensity reaches a minimum value which is equal to the approach stream turbulence intensity, (c) the flame-generated turbulence intensity reaches a maximum value when the rate of production of turbulent vorticity is equal to about half its rate of viscous dissipation.     

Population inversion and gain distribution in supersonic mixed flow systems

Experimental results are reported for small signal gain distribution across a cavity of a mixed flow gasdynamic laser system at different turbulent supersonic mixing regimes. It is shown that the temperature range of the GDL generation regime can be extended up to 7000°K, and gain coefficients as high as 3.5 m⁻¹ be attained in a “double-freeze” supersonic gas flow. Basic advantages are discussed as well as the opportunity to obtain higher efficiencies in thermally pumped laser systems.     

An approach to the autoignition of a turbulent mixture

This paper considers the turbulent homogeneous mixing of two reactants undergoing a one step, second order, irreversible, exothermic chemical reaction with a rate constant of the Arrhenius type. A statistically stationary turbulent velocity field is assumed given and unaffected by mass or heat production due to the chemical reaction. Relative density fluctuations are neglected. A Hopf-like functional formalism is presented, with application to both statistically inhomogeneous and statistically homogeneous flows. Single and double point probability density function differential equations are derived from those functional equations. The limit of very large activation energies is considered; a low degree of statistical correlation between temperature and concentration fields during the ignition period is hypothesized. After making use of the homogeneity assumption a closure problem is still present due to the nonlocalness of the molecular diffusion term. The problem is rendered closed by assuming a Gaussian conditional expected value for the temperature at a point given the temperature at a neighboring point. The closure is seen to preserve very important mathematical and physical properties. A linear first order hyperbolic differential equation with variable coefficients for the probability density function of the temperature field is obtained. A second Damköhler number based on Taylor's microscale turns out to be an important controlling parameter. A numerical integration for different values of the second Damköhler number and the initial stochastic parameters is carried out. The mixture is seen to evolve towards an eventual thermal runaway, the detailed behavior however being different for different systems. Some peculiarities during the ignition period evolution are uncovered.     

Chemi-acoustic instability structure in irreversibly reacting systems

A basic understanding of the structure of the interactions between chemical and acoustic instabilities and the effects derived therefrom is sought for a medium undergoing one-step irreversible chemical reaction. Detailed examination of the acoustic-chemical system after the completion of the reaction shows the distinct presence of the chemical, acoustic, and mixed modes of instability. The chemical mode appears as a stationary yet spatially inhomogeneous entropy distribution, even though the medium initially (before reaction) is homogeneous throughout. The acoustic mode appears as a composite of both right- and left-travelling pressure or velocity waves, even though the initial acoustic wave is only right-travelling. The left-travelling wave is generated due to partial reflection or scattering as the right-travelling wave propagates in a spatially inhomogenous medium during the chemical reaction and is sustained even after the reaction is completed. The mixed modes of density and temperature fluctuations apparently retain to some degree the characteristics of both the acoustic and the chemical modes and may, under certain conditions, be dominated by either one of the modes of behavior. This is determined largely by the parameter Ω, the ratio of chemical to acoustic time scales.During reaction, it is found that the observed complicated behavior can be interpreted fruitfully in terms of the mode concept developed for the post-reaction behavior. It is observed that the temporal development of the physical fluctuating variables is dependent upon spatial position (the effect being stronger as Ω is reduced). Further, it is found that at specific values of Ω, the following effects are maximized: (a) acoustic amplification; (b) wave reflection; (c) reaction evolution enhancement.The energy contained in the fluctuations is found to be composed of an acoustic and a chemical part. The latter dominates during reaction. However, after reaction, the chemical part exactly balances any change in the acoustic part which occurs during reaction, thus resulting in no net change in the fluctuation energy. The chemical part seems to represent the loss (or gain) of mean thermal energy which was diverted by chemi-acoustic interactions into the increased (or decreased) acoustic energy of the system.     

Detonation engine fed by acetylene–oxygen mixture

The advantages of a constant volume combustion cycle as compared to constant pressure combustion in terms of thermodynamic efficiency has focused the search for advanced propulsion on detonation engines. Detonation of acetylene mixed with oxygen in various proportions is studied using mathematical modeling. Simplified kinetics of acetylene burning includes 11 reactions with 9 components. Deflagration to detonation transition (DDT) is obtained in a cylindrical tube with a section of obstacles modeling a Shchelkin spiral; the DDT takes place in this section for a wide range of initial mixture compositions. A modified ka-omega turbulence model is used to simulate flame acceleration in the Shchelkin spiral section of the system. The results of numerical simulations were compared with experiments, which had been performed in the same size detonation chamber and turbulent spiral ring section, and with theoretical data on the Chapman–Jouguet detonation parameters.     

Two-phase flow effect on hybrid rocket combustion

This study numerically explores the aerodynamic and combustion processes in a hybrid rocket combustor, under a two-phase turbulent flow environment, considering the evaporation, combustion and drag of droplet and droplet ignition criterion. The predictions of temperature, reaction mode, reactant mass fraction, velocity, oxidizer consumption, fuel regression and droplet number distribution enhance understanding of the two-phase combustion aerodynamics inside the combustor. A parametric study of the inlet spray pattern, including spray cone angle, spray injection velocity and droplet size, is performed to improve the operation of reactant mixing and higher fuel regression rate. Analytical results indicate that both the oxidizer consumption and the fuel regression increase with increasing spray cone angle and spray injection velocity in the practical range of operation. However, for stoichiometric operation, the superior spray cone angle is within 20–60°, and spray injection velocity within 20–40 m/s, under a volume-mean droplet radius of 50 μm. The power dependence of solid-fuel regression on total mass flux is found to decrease with rising of droplet mean size.     

Life and intelligence in the universe from nanobiological principles: A survey and budget of concepts and perspectives

The new-born bioscience called Nanobiology has tackled the problems of the possibility of existence of extraterrestrial life and intelligence and of biosystem distribution in the Universe, as such questions actually belong to the realm of Theoretical Biology. The central, and yet unanswered points of such science have been reinvestigated by attempting knowledge and control of the hard-to-determine nanoscale-level classical and quantum interactions, which would supposedly give mechanistic, definite answers, both informationally and energetically, to the vexing questions put by biosystems to science: is the “living state” a physically definible concept, and how to define it? Are nanoscale kinetics or even detailed mechanics involved in the origin of life? What about intelligence, consciousness and their nanophysical roots? Are “life” and “intelligence” engineerable properties, or is any Artificial Intelligence program bound to mere metaphors? Self-organization, studied at the thermodynamic and the hydrodynamic level, showed the possibility of chemical evolution from amino acids, probably of cometary and/or meteoritic origin, up to spatiotemporal organization, autopoiesis and biological evolution, but didn't explain the origins of life. Questioning the uniqueness of the earthly evolutionary chemistry is cardinal for the ETI dilemma, as from a budgetary appraisal of perspectives in bionanoscale chaotic undecidable dynamics, quantum gravity and quantum vacuum, both “living state” and “intelligence” look like nonlocal, spacetime-linked cosmic phenomena.     

空气涡轮火箭发动机内外涵气流掺混研究

通过无化学反应、均匀进气条件下肼单组元空气涡轮火箭发动机混流燃烧室内流场的数值计算,得到了流向涡与正交涡系产生、衰减演变过程及其对内外涵气流掺混效率的影响规律。结果表明,大尺度阵列二次环流诱导形成的流向涡对内外涵气流掺混起主导作用,大波瓣穿透率的斜切波瓣混流器的综合性能较优。结合热试车结果,分析了包括波瓣混流器在内的两类掺混方案的强化掺混效率。分析表明,非均匀进气条件对小尺寸空气涡轮火箭发动机掺混燃烧效率影响很大。     

富氢/富氧燃气温度对同轴直流气-气喷嘴性能的影响

通过求解使用k-ε湍流模型的Navier-Stokes方程组对采用同轴直流气-气单喷嘴燃烧室的燃烧流场进行数值模拟,对比分析了富氢/富氧燃气推进剂与常温氢气/氧气推进剂条件下的燃烧流场、燃烧室室壁和喷注面板处的燃气温度,研究了富氢/富氧燃气温度变化对燃烧流场和燃烧室热载的影响。数值结果表明:富氢/富氧燃气气-气喷嘴的燃烧性能较好,但热载较高;富氢/富氧燃气温度一定范围内提高对燃烧性能影响不明显,而热载增加。     

Nonisotropic turbulence in swirling flows

Recent experimental, inverse and prediction studies have disputed isotropy assumptions for turbulent swirling flows. To cater for the simulation (and hence solution) of these flows, simple turbulence models may be modified or more advanced models developed. This paper reviews recent work in this area, clarifies the entent of nonisotropy and recommends suitable turbulence models to effect closure of the governing Reynolds equations.     

Aquatic modules for bioregenerative life support systems based on the C.E.B.A.S. biotechnology

Most concepts for bioregenerative life support systems are based on edible higher land plants which create some problems with growth and seed generation under space conditions. Animal protein production is mostly neglected because of the tremendous waste management problems with tetrapods under reduced weightlessness. Therefore, the “Closed Equilibrated Biological Aquatic System” (C.E.B.A.S.) was developed which represents an artificial aquatic ecosystem containing aquatic organisms which are adpated at all to “near weightlessness conditions” (fishes Xiphophorus helleri, water snails Biomphalaria glabrata, ammonia oxidizing bacteria and the rootless non-gravitropic edible water plant Ceratophyllum demersum). Basically the C.E.B.A.S. consists of 4 subsystems: a ZOOLOGICASL COMPONENT (animal aquarium), a BOTANICAL COMPONENT (aquatic plant bioreactor), a MICROBIAL COMPONENT (bacteria filter) and an ELECTRONICAL COMPONENT (data acquisition and control unit). Superficially, the function principle appears simple: the plants convert light energy into chemical energy via photosynthesis thus producing biomass and oxygen. The animals and microorganisms use the oxygen for respiration and produce the carbon dioxide which is essential for plant photosynthesis. The ammonia ions excreted by the animals are converted by the bacteria to nitrite and then to nitrate ions which serve as a nitrogen source for the plants. Other essential ions derive from biological degradation of animal waste products and dead organic matter. The C.E.B.A.S. exists in 2 basic versions: the original C.E.B.A.S. with a volume of 150 liters and a self-sustaining standing time of more than 13 month and the so-called C.E.B.A.S. MINI MODULE with a volume of about 8.5 liters. In the latter there is no closed food loop by reasons of available space so that animal food has to be provided via an automated feeder. This device was flown already successfully on the STS-89 and STS-90 spaceshuttle missions and the working hypothesis was verified that aquatic organisms are nearly not affected at all by space conditions, i . e. that the plants exhibited biomass production rates identical to the ground controls and that as well the reproductive, and the immune system as the the embryonic and ontogenic development of the animals remained undisturbed. Currently the C.E.B.A.S. MINI MODLULE is prepared for a third spaceshuttle fligt (STS-107) in spring 2001. Based on the results of the space experiments a series of prototypes of aquatic food production modules for the implementation into BLSS were developed. This paper describes the scientific disposition of the STS-107 experiments and of open and closed aquaculture systems based on another aquatic plant species, the Lemnacean Wolffia arrhiza which is cultured as a vegetable in Southeastern Asia. This plant can be grown in suspension culture and several special bioreactors were developed for this purpose. W. arrhiza reproduces mainly vegetatively by buds but also sexually from time to time and is therefore especially suitable for genetic engineering, too. Therefore it was used, in addition, to optimize the C.E.B.A.S. MINI MODULE to allow experiments with a duration of 4 month in the International Space Station the basic principle of which will be explained. In the context of aquaculture systems for BLSS the continuous replacement of removed fish biomass is an essential demand. Although fish reproduction seems not to be affected in the short-term space experiments with the C.E.B.A.S. MIMI MODULE a functional and reliable hatchery for the production of siblings under reduced weightlessness is connected with some serious problems. Therefore an automated “reproduction module” for the herbivorous fish Tilapia rendalli was developed as a laboratory prototype. It is concluded that aquatic modules of different degrees of complexity can optimize the productivity of BLSS based on higher land plants and that they offer an unique opportunity for the production of animal protein in lunar or planetary bases.     

Effect of enclosure shape on natural convection velocities in microgravity

A numerical analysis was performed to compare natural convection velocities in two-dimensional enclosures of various shape. The following shapes were investigated: circle, square, horizontal and upright 2 × 1 aspect ratio rectangles, horizontal and upright half-circles, diamond (square oriented with diagonal vertical) and triangle (equilateral and horizontal base). In all cases, the length scale in the various dimensionless parameters, such as Rayleigh number, is defined as the diameter of the equal area circle. Natural convection velocities were calculated for Rayleigh numbers of 100 and 500 with the temperature difference taken to be across (a) the maximum horizontal dimension, (b) the median horizontal line (line through centroid) and (c) the horizontal distance such that the temperature gradient is the same for shapes of equal area. A Rayleigh number of 1000 is within the “low Rayleigh number” range for agreement with first order theory for circular enclosures. A Rayleigh number of 5000 is slightly out of this range. For the class of shapes including the square, upright half-circle and upright rectangle, the computed velocities were found to agree very closely with that of the equal area circle when the temperature difference is taken to be across the maximum horizontal dimension [condition (a)]. The velocities for the horizontal rectangle and half circle were found to be approximately one-half that of the equal area circle for the same condition. Better overall agreement among all shapes was obtained by setting the temperature difference across a distance such that the temperature gradients were equal for shapes of equal area.     

Influences of H2O mass fraction and chemical kinetics mechanism on the turbulent diffusion combustion of H2–O2 in supersonic flows

Hydrogen is one of the most promising fuels for the airbreathing hypersonic propulsion system, and it attracts an increasing attention of the researchers worldwide. In this study, a typical hydrogen-fueled supersonic combustor was investigated numerically, and the predicted results were compared with the available experimental data in the open literature. Two different chemical reaction mechanisms were employed to evaluate their effects on the combustion of H₂–O₂, namely the two-step and the seven-step mechanisms, and the vitiation effect was analyzed by varying the H₂O mass fraction. The obtained results show that the predicted mole fraction profiles for different components show very good agreement with the available experimental data under the supersonic mixing and combustion conditions, and the chemical reaction mechanism has only a slight impact on the overall performance of the turbulent diffusion combustion. The simple mechanism of H₂–O₂ can be employed to evaluate the performance of the combustor in order to reduce the computational cost. The H₂O flow vitiation makes a great difference to the combustion of H₂–O₂, and there is an optimal H₂O mass fraction existing to enhance the intensity of the turbulent combustion. In the range considered in this paper, its optimal value is 0.15. The initiated location of the reaction appears far away from the bottom wall with the increase of the H₂O mass fraction, and the H₂O flow vitiation quickens the transition from subsonic to supersonic mode at the exit of the combustor.     

Perturbation solutions for variable energy blast waves

The trajectory of and the flow field behind blast waves with time varying energy input is determined. Freeman's (1968) Lagrangean coordinate formulation is modified to include both the geometric factor, α, for plane, cylindrical and spherical shocks and also non-integer values of β, the energy input parameter, in a single computational algorithm. Numerical problems associated with vanishing density at the inner mass boundary or “piston face” are then examined and solved. Second order perturbation solutions about the solution for an infinite strength shock are then obtained in Sakurai's (1965) inverse shock Mach number expansion parameter for 0 β < α + 1. Tables and graphs of significant numerical coefficients are presented for comparison to, and extension of, results of other authors. Graphs of typical shock trajectories and flow field density, pressure and velocity variations are also presented and discussed.     

高过载条件下固体火箭发动机绝热层失效研究

根据化学反应／两相流耦合,建立了铝粒子燃烧模型,通过化学反应速率模型的模化湍流燃烧,对含铝推进剂固体火箭发动机在高过载条件下的内流场进行数值研究。结果表明,过载条件下燃烧室局部异常铝滴积聚及剧烈的化学放热反应是导致绝热层异常烧蚀的主要原因。     

The influence of pressure on combustion intensity

The influence of pressure in the range of 3–15 kgf/cm² on combustion intensity is studied experimentally in a medium-sized rig. The apparatus is described and temperature measurements by different thermocouple techniques are discussed; gas composition was monitored by gas chromatography. The experiments were performed at different cross sections allowing to map temperatures and gas concentrations. The results show that temperature and carbon dioxide concentration increase more rapidly as the pressure is raised. Carbon monoxide appears as an intermediate and is concentrated near the combustion axis. The combustion zone becomes shorter with increasing pressure and the combustion intensity increases correspondingly.