Ignition of single coal particle in a hot furnace under normal- and micro-gravity condition

An experimental study on ignition and combustion of single particles was conducted at normal gravity (1-g) and microgravity (μ-g) for three high volatile coals with initial diameter of 1.5 and 2.0 mm, respectively. The non-intrusive twin-color pyrometry method was used to retrieve the surface temperature of the coal particle through processing the images taken by a color CCD camera. At the same time, a mathematical model considering thermal conduction inside the coal particle was developed to simulate the ignition process.Both experiments and modeling found that ignition occurred homogeneously at the beginning and then heterogeneously for the testing coal particles burning at μ-g. Experimental results confirmed that ignition temperature decreased with increasing volatile content and increasing particle size. However, contradicted to previous studies, this study found that for a given coal with certain particle size, ignition temperature was about 50–80 K lower at μ-g than that at 1-g.The model predictions agreed well with the μ-g experimental data on ignition temperature. The criterion that the temperature gradient in the space away from the particle surface equaled to zero was validated to determine the commence of homogeneous ignition. Thermal conduction inside the particle could have a noticeable effect for determining the ignition temperature. With the consideration of thermal conduction, the critical size for the phase transient from homogeneous to heterogeneous is about 700 μm at ambient temperature 1500 K and oxygen concentration 0.23.     

Characteristics of laminar lifted flames in coflow jets with initial temperature variation

Characteristics of laminar lifted flames of propane highly diluted with nitrogen have been investigated by varying the initial temperature in coflow jets. The result showed that the lifted flame maintained the tribrachial structure up to the initial temperature of 900 K and the liftoff height decreased with initial temperature and dilution ratio. The overall behavior of liftoff heights correlated well with the jet velocity scaled by the stoichiometric laminar burning velocity, emphasizing the importance of the stoichiometric laminar burning velocity on the propagation speed of tribrachial flame. The exponent of the liftoff height with jet velocity in the relation of increased with initial fuel mole fraction, which has been attributed to the differential diffusion between propane and diluent nitrogen. Consequently, nitrogen concentration varied along the stoichiometric contour, which affected the propagation speed. Also, the exponent increased with initial temperature due to the sensitiveness of the propagation speed variation with nitrogen dilution on initial temperature. The liftoff conditions have been observed for the jet velocity even smaller than the stoichiometric laminar burning velocity at relatively low initial temperatures. This can be attributed to the effect of the buoyancy. Liftoff velocities accounting for the relative buoyancy effect were found to have a satisfactory correlation with regardless of initial temperatures and nitrogen dilution.     

Thermal convection and related instabilities in models of crystal growth from the melt on earth and in microgravity: Past history and current status

The paper presents a comparative study of a number of theoretical/experimental/numerical results concerning the dynamics of natural (gravitational), Marangoni and related mixed convection in various geometrical models of widely‐used technologies for the production of single‐crystalline materials (Horizontal and vertical Bridgman growth, Czochralski method, Floating Zone Technique). Emphasis is given to fundamental knowledge provided over the years by landmark analyses as well as to very recent contributions. Such a knowledge is of paramount importance since it is validating new, more complex models, accelerating the current trend towards predictable and reproducible phenomena and finally providing an adequate scientific foundation to industrial processes which are still conducted on a largely empirical basis. A deductive approach is followed with fluid‐dynamic systems of growing complexity being treated as the discussion progresses. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of Schmidt number on the structure and propagation of density currents

The results of a numerical study of two- and three-dimensional Boussinesq density currents are described. They are aimed at exploring the role of the Schmidt number on the structure and dynamics of density driven currents. Two complementary approaches are used, namely a spectral method and a finite-volume interface capturing method. They allow for the first time to describe density currents in the whole range of Schmidt number 1 ≤ Sc ≤ ∞ and Reynolds number 10² ≤ Re ≤ 10⁴. The present results confirm that the Schmidt number only weakly influences the structure and dynamics of density currents provided the Reynolds number of the flow is large, say of O(10⁴) or more. On the contrary low- to moderate-Re density currents are dependant on Sc as the structure of the mixing region and the front velocities are modified by diffusion effects. The scaling of the characteristic density thickness of the interface has been confirmed to behave as (ScRe)^−1/2. Three-dimensional simulations suggest that the patterns of lobes and clefts are independent of Sc. In contrast the Schmidt number is found to affect dramatically (1) the shape of the current head as a depression is observed at high-Sc, (2) the formation of vortex structures generated by Kelvin–Helmholtz instabilities. A criterion is proposed for the stability of the interface along the body of the current based on the estimate of a bulk Richardson number. This criterion, derived for currents of arbitrary density ratio, is in agreement with present computed results as well as available experimental and numerical data.      

Large Eddy Simulation of Turbulent Convective Cavity Flow

The large eddy simulation model with Smagorinsky subgrid-scale model was applied to two-dimensional turbulent convective cavity flow. The Reynolds number is lying from 1×10⁴ to 4×10⁵ and Archimedes number from 0 to 0.4. The simulation results were compared with the k?? model results and experimental results wherever possible. Flow results were in good agreement with experimental data across the mid-planes. Effects of Smagorinsky constant and grid resolution were investigated.     

Effects of buoyancy and thermophysical property variations on the flow of supercritical carbon dioxide

The flow and heat transfer behaviours of fluids at supercritical pressure have been studied using direct numerical simulations (DNS), in which one or more thermal properties are artificially frozen to discern the various physical mechanisms from each other so as to better understand the complex phenomena. Different from previous similar studies on this topic, this study focuses on the axial flow development resulted from the large variations of thermophysical properties. The contribution of the flow inertia has been quantified by analysing the momentum balance for each case studied, which has been found to be significant throughout the entire length of the pipe in cases when buoyancy is considered. The effect of the inertia on momentum in turn impacts on turbulence production, generally delaying flow laminarisation. Such an influence of flow development is non-trivial and cannot be omitted in flow analysis and heat transfer calculations. This suggests that the results of simplified analyses based on a spatially developed flow cannot be directly applied to such flows despite they can be very useful in developing fundamental understanding of the physics. Similarly, this also explains that in some cases, buoyancy parameters based on local flow quantities cannot describe heat transfer deterioration accurately. The effect of variable viscosity alone can cause turbulence reduction by flattening the velocity profile, but it will not turn the velocity profile to an M-shape, which can only be achieved by buoyancy.     

A study on single fuel droplets combustion in vertical direct current electric fields

Both the electric force working on flames and the natural buoyancy are body forces, so there is a possibility to control the natural buoyancy by applying than electric field. It is important to discuss the body force in the flame because it induces the convective flow around flames. In this circumstance, combustion behavior of single droplets in vertical direct current electric fields was investigated. Ethanol, n-octane, and toluene were used as fuels, and the flame shape and the burning rate constants were measured. The distance between electrodes is 50 mm, and the applied voltage ranged between −4 and 6 kV as the bottom electrode is ground. When the direction of the electric field is opposite to the natural buoyancy direction, the applied voltage exists that make the flame symmetrical in the vertical direction, and the burning rate constants have local minima for ethanol and n-octane at the voltage. However, the minimum burning rate constants are larger than those under microgravity. This means that the electric force roughly balances with natural buoyancy, but it does not completely balance with the same. When the direction of the electric field is in the same direction as the natural buoyancy, there exist some experimental results, which cannot be explained by the assumption that electric field induces the body force only through the positive ions. This suggests that the additional body force is induced by the negative ions. The effects of negative charged soot particles on the combustion behavior are also discussed for the sooty flame of toluene.     

PLIF and PIV measurements of the self-preserving structure of steady round buoyant turbulent plumes in crossflow

Measurements of the mean concentration of source fluid and mean velocity fields were obtained for the first time in the self-preserving region of steady round buoyant turbulent plumes in uniform crossflows using Planar-Laser-Induced-Fluorescence (PLIF) and Particle-Image-Velocimetry (PIV), respectively. The experiments involved salt water sources injected into water/ethanol crossflows within a water channel. Matching the index of refraction of the source and ambient fluids was required in order to avoid image distortion and laser intensity nonuniformities. Further experimental methods and procedures are explained in detail. The self-preserving structure properties of the flow were correlated successfully based on the scaling analysis of [Fischer, H.B., List, E.J., Koh, R.C., Imberger, J., Brooks, N.H., 1979. Mixing in Inland and Coastal Waters, Academic Press, New York, pp. 315–389]. The resulting self-preserving structure consisted of two counter-rotating vortices having their axes nearly aligned with the crossflow direction that move away from the source in the streamwise (vertical) direction due to the action of buoyancy. This alignment, was a strong function of the source/crossflow velocity ratio, u₀/v_∞. Finally, the counter-rotating vortex system was responsible for substantial increases in the rate of mixing of the source fluid with the ambient fluid compared to axisymmetric round buoyant turbulent plumes in still environments, e.g., transverse dimensions in the presence of the self-preserving counter-rotating vortex system were 2–3 times larger than the transverse dimensions of self-preserving axisymmetric plumes at similar streamwise distances from the source.     

Pattern formation and mass transfer under stationary solutal Marangoni instability

According to the seminal theory by Sternling and Scriven [1], solutal Marangoni convection during mass transfer of surface-active solutes may occur as either oscillatory or stationary instability. With strong support of Manuel G. Velarde, a combined initiative of experimental works, in particular to mention those of Linde, Wierschem and coworkers, and theory has enabled a classification of dominant wave types of the oscillatory mode and their interactions. In this way a rather comprehensive understanding of the nonlinear evolution of the oscillatory instability could be achieved. A comparably advanced state-of-the-art with respect to the stationary counterpart seemed to be out of reach a short time ago. Recent developments on both the numerical and experimental side, in combination with assessing an extensive number of older experiments, now allow one to draw a more unified picture. By reviewing these works, we show that three main building blocks exist during the nonlinear evolution: roll cells, relaxation oscillations and relaxation oscillations waves. What is frequently called interfacial turbulence results from the interaction between these partly coexisting basic patterns which may additionally occur in different hierarchy levels. The second focus of this review lies on the practical importance of such convection patterns concerning their influence on mass transfer characteristics. Particular attention is paid here to the interaction between Marangoni and buoyancy effects which frequently complicates the pattern formation even more. To shed more light on these dependencies, new simulations regarding the limiting case of stabilizing density stratification and vanishing buoyancy are incorporated.     

Three-dimensional thermocapillary and buoyancy convections and interface shape in horizontal Bridgman crystal growth

Computer simulation is conducted to study three-dimensional (3D) thermocapillary and buoyancy convections and their effects on the growth interface for horizontal Bridgman crystal growth. The free-boundary model is based on a finite volume approximation of continuity, momentum, and energy equations on a collocated grid. Crystal growth of GaAs is used as an example. From calculated results, it is observed that the effect of buoyancy convection on the growth interface is significant. With the thermocapillary effect, the 3D flow structures are not changed much, but its effect on the growth interface is not trivial. Due to the convections, the growth interface is always concave, and its deflection is affected significantly by the growth rate and thermal environment. A simple strategy of interface control is illustrated. Furthermore, slight crucible tilting can also affect the 3D flows leading to an asymmetric growth interface.