Effects of gas and particle emissions on wall radiative heat flux in oxy-fuel combustion

Oxy-fuel combustion exhibits combustion and heat transfer characteristics different from air-fuel combustion due to high concentrations of CO₂ and H₂O. This study evaluated the effect of gas and particle emissions on radiative heat transfer in oxy-fuel combustion of coal. For a hexahedral furnace, prescribed gas compositions based on combustion calculation were used to simplify the combustion reactions. The values of radiative heat fluxes (q_rad) were compared for different combustion modes, flue gas recirculation (FGR) methods, particle concentrations, furnace sizes and O₂ concentrations in the oxidizer. The radiation was calculated by the discrete ordinate method with gaseous emission predicted by the weighted sum of gray gases models (WSGGMs). The results showed that employing an optimized WSGGM is essential for the accurate prediction of q_rad in oxy-fuel combustion for gaseous fuels. The conventional WSGGM showed large errors for larger furnace volumes or under dry FGR conditions. With higher particle concentrations such as in pulverized coal combustion, however, q_rad was dominated by emission of particles. The effect of gas emissivity was not critical in the furnace with a mean beam length of 8.3m. Oxy-fuel combustion with wet FGR had higher q_rad than dry FGR. The O₂ concentration in the oxidizer was a key parameter for oxy-fuel combustion since increasing its value linearly increased q_rad.     

Study on the radiative transfer through non-gray gas mixtures within an irregular 3-D enclosure by using the modified weighted sum of gray gas method

In this study the modified weighted sum of gray gases model (WSGGM) using the gray gas regrouping technique and the discrete ordinates interpolation method (DOIM) are applied to analyze the radiative transfer within an irregular 3-D enclosure filled with non-gray gas mixture of CO₂, H₂O and N₂. The computer code developed in this study is successfully applied for solving the non-gray gas radiation within a 3-D rectangular enclosure and the gray gas radiation within an irregular 3-D enclosure by showing fairly good agreements with the existing results. In this paper the radiative transfer within an irregular 3-D enclosure filled with non-isothermal non-gray gases with uniform mixtures of CO₂, H₂O and N₂ is studied to demonstrate the applicability of the modified sum of gray gases model for irregular systems and to examine the effect of the concentrations of CO₂ and H₂O on the radiative transfer within modern combustors. Results show that the wall heat fluxes and the radiative heat source terms are increased as the concentrations of CO₂ and H₂O are increased. Results also show that the radiative fluxes caused by the mixture gases with high concentrations CO₂ and H₂O which can be observed in oxy-fuel combustion systems can reach up to nearly twice of those found in ordinary air-fuel combustion systems.     

Numerical analysis of the flow and heat transfer characteristics for forced convection-radiation in entrance region of an internally finned tubes

The flow and heat transfer characteristics of combined forced convection and radiation in the entrance region of internally finned tubes are investigated numerically in this paper. The uniform flow is considered for an inlet flow condition. A three dimensional parabolic problem is solved by a marching-type procedure envolving a series of two dimensional elliptic problems in the cross-stream plane. The SIMPLER-algorithm and Raithby's pressure-velocity coupling method are employed to analyze the flow and heat transfer characteristics. For the calculation of radiative heat transfer, the P₁-approximation and the weighted sum of gray gases method (WSGGM) are used. The effects of fin height, number of fins, optical thickness, reference temperature, and Planck number on the flow and heat transfer characteristics are examined. It was found that the effect of fin-height on the heat transfer characteristic is more dominant than that of number of fins. The present results show that the optimal non-dimensional fin height and number of fins are 0.4 and 16, respectively.     

Narrow band radiative solutions within a cubical enclosure filled with real gas mixtures

Radiative transfer by nongray gas mixtures with nonuniform concentration and temperature profiles is studied by using the statistical narrow-band model and the ray-tracing method with the sufficiently accurate T₆₀ quadrature set. Transmittances through the nonhomogeneous gas mixtures are calculated by using the Curtis-Godson approximation. Three different cases with different temperature and concentration profiles are considered to obtain benchmark solutions for the radiative transfer by nongray gas mixtures. The solutions obtained from this study are verified and found to be very well matched with the previous solutions for uniform gas mixtures. The results presented in this paper can be used as benchmark solutions in developing various solution methods for radiative transfer by nongray gas mixtures.     

Combustion of a single droplet in the presence of an oscillating flow

The two-dimensional, unsteady, laminar conservation equations for mass, momentum, energy and species transport in the gas phase are solved numerically in spherical coordinates. This is to study the heat and the mass transfer, and the combustion around a single spherical droplet. The droplet mass and momentum equations are also solved simultaneously with the gas phase equations in order to investigate the effects of droplet entrainment in the oscillating flow with and without a steady velocity. The numerical solution for a single droplet combustion gives the droplet diameter variation as well as the gas phase velocity, temperature and species concentrations as a function of time. The effects of frequency, amplitude of oscillating flow, velocity ratio of oscillating flow amplitude to the steady velocity, ambient temperature and initial droplet diameter on the droplet combustion are also investigated. The droplet burning history is not governed by thed ²-law in the presence of oscillating flow, unlike to the case under quiescent ambient conditions.     

Automated instrumentation for high-temperature Seebeck coefficient measurements

The fabrication of automated instrumentation for high-temperature Seebeck coefficient measurements is presented. K-type thermocouples were used to measure the average temperature of the sample and the Seebeck voltage. The temperature dependence of the Seebeck coefficients of the thermocouple and its negative leg were obtained by the integration method. A steady-state-based differential technique was used for Seebeck coefficient measurement. The use of limited components and a thin heater simplified the sample holder design and minimized heat loss. The power supplied to the heater determined the temperature difference across the sample and the measurement was carried out by achieving a steady state. A LabVIEW-based program was constructed to automate the measurements. The complete setup was fabricated using commonly available materials. This instrument is standardized for materials with a wide range of Seebeck coefficients and temperature differences. High temperature measurements for iron, constantan, bismuth, and Bi_0.36Sb_1.45Te₃ were carried out and the results were in good agreement with standard values.     

Experimental and numerical studies on the thermal analysis of the plate in indirectly fired continuous heat treatment furnace

Experimental and numerical studies were performed by considering convective and radiative heat transfer to predict the transient thermal behavior of a plate in an indirectly fired continuous heat treatment furnace. The temperature profiles in the plate were determined by solving the transient one-dimensional heat conduction equation in conjunction with appropriate boundary conditions by using a time marching scheme. The results obtained from the transient analysis were substantiated by comparing with experimental results. Additionally, parametric investigations were performed to examine how the thermal behavior of the plate is affected by plate and refractory emissivities, charging temperature and residence time of the plate, gas temperature of the work and drive sides of the heat treatment furnace, and plate thickness. This paper was recommended for publication in revised form by Associate Editor Ohchae Kwon Young-Deuk Kim is a graduate student at Hanyang University in Seoul, Korea. He earned his B.S. in Mechanical Engineering from Korea Maritime University in 2002 and his M.S. in mechanical engineering from Hanyang university in 2004. His current research areas are modeling of automotive aftertreatment catalysts, optimal design of thermal systems, and phase change modeling with free surface flow. Deok-Hong Kang is a senior researcher at the RIST (Research Institute of Industrial Science and Technology) in Pohang, Korea. He earned his B.S. and M.S. in mechanical engineering from Hanyang University in 1989 and 1993, respectively, and his Ph.D. in mechanical engineering from POSTECH in 2004. His current research areas are mathematical modeling for combustion control, furnace optimization control system, and energy saving engineering in all kinds of furnaces. Woo-Seung Kim is a professor in mechanical engineering at Hanyang University in Ansan, Korea. He earned his B.S. in Mechanical Engineering in 1981 from Hanyang University and his M.S. and Ph.D. in mechanical engineering from North Carolina State University in 1986 and 1989, respectively. His current research areas are modeling of automotive aftertreatment systems, inverse heat transfer problems, optimal design of thermal systems, and phase change heat transfer problems with free surface flow.     

An experimental study on laminar CH<Subscript>4</Subscript>/O<Subscript>2</Subscript>/N<Subscript>2</Subscript> premixed flames under an electric field

This work investigates the electric field effect on gas temperature, radiative heat flux and flame speed of premixed CH₄/O₂/N₂ flames in order to gain a better insight into the mechanism of controlling the combustion process by electrophysical means. Experiments were performed on laminar Bunsen flames (Re<2200) of lean to rich mixture composition (φ =0.8–1.2) with slight oxygen enrichment (Ω=0.21-0.30). The Schlieren flame angle technique was used to determine the flame speed, and thermocouple measurements at the post flame gas were conducted. The radiative heat flux was measured by using a heat flux meter. At high field strengths, coincident with the appearance and enhancement of flame surface curvatures, an apparent change in flame speed and gas temperature was observed. However, the application of an electric field had no significant effect on flame speed and temperature when the flame geometry was unaltered. This was supported by radiative heat flux showing negligible electric field effects. The modification in flame temperature and flame speed under electric field was attributed to the field-induced flame stretch due to the body forces produced by the ionic winds. This additional flame stretch, coupled with the influence of non-unity Lewis number, accounts for such changes. This reinforces the idea that the action of an electric field on flames with a geometry that remains practically undeformed produces very minimal effect on flame speed, temperature and radiative heat flux. A possible mechanism of combustion control by the application of flame stretch using electric field was introduced.     

Development of the delamination evaluation parameters (I) — the delamination aspect ratio and the delamination shape factors —

Although the previous researches evaluated the fatigue behavior of A1/GFRP laminates using the traditional fracture mechanism, their researches were not sufficient to do it : the damage zone of A1/GFRP laminates was occurred at the delamination zone instead of the crack-metallic damages. Thus, previous researches were not applicable to the fatigue behavior of Al/GFRP laminates. The major purpose of this study was to evaluate delamination behavior using the relationship between crack length(a) and delamination width(b) in A1/GFRP laminate. The details of investigation were as follows : 1) Relationship between the crack length(a) and the delamination width(b), 2) Relationship between the delamination aspect ratio(b/a) and the delamination area rate ( (A_d)_n/ (A _d ) a_ll) , 3) The effect of delamination aspect ratio(b/a) on the delamination shape factor (f _s ) and the delamination growth rate (dA _d /da). As results, it was known that the delamination aspect ratio(b/a) was decreased and the delamination area rate ( (A _d ) _n / (A _d ) _All ) was increased as the normalized crack size(a/W) was increased. And, the delamination shape factors (f _s ) of the ellipse-II (f _s3 ) was greater than of the ellipse-1 (f _s2 ) but that of the triangle (f_s1) was less than of the ellipse-I (f _s2 ).     

Determining cutting speeds for the CO2 laser machining of decorative ceramic tile

This paper presents a comparison of theoretically predicted optimum cutting speeds for decorative ceramic tile with experimentally derived data. Four well-established theoretical analyses are considered and applied to the laser cutting of ceramic tile, i.e. Rosenthal's moving point heat-source model, and the heat-balance approaches of Powell, Steen and Chryssolouris. The theoretical results are subsequently compared and contrasted with actual cutting data taken from an existing laser machining database. Empirical models developed by the author are described which have been successfully used to predict cutting speeds for various thicknesses of ceramic tile.Notation A absorptivity - a thermal diffusivity (m²/s) - C specific heat (J/kgK) - d cutting depth (mm) - E _cut specific cutting energy (J/kg) - k thermal conductivity (W/mK) - J laser beam intensity (W/ m²) - L latent heat of vaporisation (J/kg) - l length of cut (mm) - n coordinate normal to cutting front - P laser power (W) - P _b laser power not interacting with the cutting front (W) - q heat input (J/s) - R radial distance (mm) - r beam radius (mm) - s substrate thickness (mm) - S _crit critical substrate thickness (mm) - T temperature (°C) - T _o ambient temperature (°C) - T _p peak temperature (°C) - T _s temperature at top surface (°C) - t time (s) - V cutting speed (mm/min) - V _opt optimum cutting speed (mm/min) - w kerf width (mm) - X, Y, Z coordinate location - x, y, z coordinate distance (mm) - conductive loss function - radiative loss function - convective loss function - angle between -coordinate andx-coordinate (rad) - coordinate parallel to bottom surface - angle of inclination of control surface w. r. t.X-axis (rad) - coupling coefficient - translated coordinate distance (mm) - density (kg/m³) - angle of inclination of control surface w.r.t.Y-axis (rad)     

柴油机缸内辐射换热的多维数值模拟研究

对柴油机来说,辐射换热极为重要,在缸内总传热量中占有非常大的比重,直接关系到发动机热效率及因传热引起的各种热负荷、热强度问题,同时,辐射换热对燃烧系统的研究也十分重要,辐射热流量会深刻影响内燃机的燃烧性能,对发动机的各种燃烧产物的形成产生至关重要的影响。为此,利用离散传递法实现柴油机缸内辐射换热的多维数值模拟,通过多维模拟计算同时考察燃烧室部件表面发射率及喷雾提前角对柴油机缸内辐射换热的影响。结果表明：活塞的辐射热流量峰值高于缸盖的辐射热流;缸盖的辐射热流量的最大值并不在中心位置处,而是随时间变化;随着壁面辐射率的增加,缸内向燃烧室部件辐射换热量逐渐增大;喷雾提前角直接影响所有燃烧室部件表面的辐射热流密度。     

Tribochemical Reactions and Lubricating Effects of Fluorinated Methanes for Al2O3 Ceramic

It has been found that CF₃CH₂F (HFC-134a) gas is an effective lubricant for several ceramics because of the formation of fluorine-containing tribochemical products. To understand the influence of the molecular structure of fluorine-containing gases on the lubricating characteristics, the lubricating effects and tribochemical reactions of some fluorinated methanes for Al₂O₃ ceramic were studied. X-ray Photoelectron Spectroscopy (XPS) and Time-of-Flight Secondary Ion Mass Spectroscopy (TOF-SIMS) were used to identify the chemical structure of tribochemical products. It was found that the friction and wear of Al₂O₃ ceramic were dependent on the molecular structure of the reacting fluorocarbon gases. The lowest friction was obtained in CH₂F₂. CHF₃ showed the best anti-wear effect. The results of surface analyses indicate that greater amounts of tribochemical products are produced when Al₂O₃ rubs in CH₂F₂ than in CF₄ and CHF₃ gases. C–C/C–H/carbon and AlF₃, which are mainly formed on the sliding surface in CH₂F₂, are responsible for the low friction. The mechanism of tribochemical reactions of the different environmental molecules on the nascent surface of Al₂O₃ was also discussed.     

A steady-state combustion modelling of composite solid propellants

By depicting the transfer of heat and combustion reaction to take place within thin gas layers close to the propellant surface burning in a steady-state fashion, a mathematical equation has been deduced to describe the burning rate of solid propellant as a function of initial grain temperature and chamber pressure. It has been also assumed that chemical reaction could take place in premixing-diffusing zone but were carried out mainly in the reaction-flame zone. All these phenomena taken place in each zone of combustion have been assumed to be steady-state. In the present investigation, the equation, γ=k·(1/R(T _i +C))ⁿ. exp (-E _a/R(T _i +C)(P/z) is being presented and it is compared with experimental data. The proposed model has been tested and evaluated vis-a-vis strand burner data for three different propellants based on CTPB, and it has been found that the deviation of the computed burning rates from the measured rates ranged up to 2%.     

The effect of gasses on the viscosity of dimethyl ether

Dimethyl ether (DME) has been recognised as a clean substitute for diesel oil as it does not form soot during combustion. DME has a vapour pressure of 6 bar at 25 °C; so pressurisation is necessary to keep DME liquid at ambient temperature. Inert gases are good candidates as pressurising media, but their effect on DME viscosity is unknown.Argon (Ar), nitrogen (N₂), carbon dioxide (CO₂), hydrogen (H₂) and propane (C₃H₈) have been investigated at pressure levels of 12–15 bar. A Cannon-Manning semi-micro capillary glass viscometer, size 25, enclosed in a cylindrical pressure container, of glass, submerged completely in a constant temperature bath, has been used. A distinct reduction of efflux times was found only for the gas, CO₂. The reduction in efflux time was about 9%.The kinematic viscosity of pure DME was determined to be: 0.188±0.001 cSt, 25 °C. A previously reported viscosity of pure DME has been corrected for the surface tension effect. Viscosity determination was initially based on a direct comparison of efflux times of DME with that of distilled water. The calculation gave a revised viscosity of 0.186±0.002 cSt, 25 °C, consistent with the above experimental result.     

Sliding wear characteristics of plasma-sprayed Al2O3 and Cr2O3 coatings against copper alloy under severe conditions

Al₂O₃ and Cr₂O₃ coatings were deposited by atmospheric plasma spraying and their tribological properties dry sliding against copper alloy were evaluated using a block-on-ring configuration at room temperature. It was found that the wear resistance of Al₂O₃ coating was superior to that of the Cr₂O₃ coating under the conditions used in the present study. This mainly attributed to its better thermal conductivity of Al₂O₃ coating, which was considered to effectively facilitate the dissipation of tribological heat and alleviate the reduction of hardness due to the accumulated tribological heat. As for the Al₂O₃ coating, the wear mechanism was plastic deformation along with some micro-abrasion and fatigue-induced brittle fracture, while the failure of Cr₂O₃ coating was predominantly the crack propagation-induced detachment of transferred films and splats spallation.     

Some stereological and statistical consequences derived from Cartier's formula

Let X and Y be two random variables with distribution functions F_x and F_Y, respectively. If Cartier's formula E[Y|X]=X holds, then F_x is necessarily less dispersed than F_Y. In this paper, the main consequences of Cartier's formula are derived. Some examples where such a formula holds are also given, namely in stereology (comparison of various estimators of Minkowsky's functionals) and in geostatistics (the change of support problem).     

Growth and properties analysis of metal-organic chemical vapor deposited MgxZn1-xO films on C-Al2O3 substrates

Mg_xZn_1−xO (0 〈 x ⩽ 0.12) thin films with the wurtzite structure have been successfully grown on c-Al₂O₃ substrates by metal-organic chemical vapor deposition (MOCVD). X-ray photoelectron spectroscopy (XPS), X-ray diffraction (XRD), photoluminescence (PL) spectrometry, and transmission measurements are performed to study the characteristics of Mg _x Zn_1−x O thin films. Results show that with increasing Mg content, the diffraction peak of Mg _x Zn_1−x O thin films shifts towards a higher diffraction angle (the biggest shift is 0.22°), and the full width at half maximum (FWHM) of the diffraction peak is broadened. Meanwhile, a blue-shift occurs at the near-band-edge (NBE) emission peak and the largest blue-shift of the band gap of the Mg _x Zn_1−x O films is 113 meV with Mg content x50.12. Therefore, the energy band gap of the Mg _x Zn_1−x O films is determined by Mg content in the thin films and the energy band gap increases with an increase of Mg content.     

Geometrical effects on spray characteristics of air-pressurized swirl flows

In an effort to obtain the significant features associated with the ALR and length/diameter ratio of the final discharge orifice in swirling flows, experimental observations using a 3-D PDPA system were carried out. Profiles of SMD distributions depending on l _o /d _o , correlation between SMD and turbulence intensities in terms of l _o /d _o and correlations between droplet size and turbulence components were quantitatively analyzed. As discussed in a previous literature, an axisymmetric swirl angle of 30° was selected for this investigation because of its strong turbulence levels in the flow-field and finer droplet disintegrations. Three ALRs of 0.093, 0.106, and 0.122 as well as the length/diameter ratio of 0.15, 0.45, and 0.60 were chosen as parameters. Due to the complex interactions in swirling flows under these variables, this experimental observation will be of fundamental importance to the understanding of geometrical effects on spray trajectory. From the observations, it is indicated that increasing the ALR causes the spray development to be more dependent on number density and volume flux. The results indicated that the SMD decreases discernibly with smaller l _o /d _o , substantiating the fact that turbulence intensities are inversely proportional to the SMD. But, l _o /d _o is quite proportional to the SMD.     

Tribo-physical and tribo-chemical aspects of WC-based cermet/Ti3SiC2 tribo-pair at elevated temperatures

The tribological behaviors of tungsten carbide (WC) based cermet/Ti₃SiC₂ tribo-pair at elevated temperatures were investigated. Lead oxide (PbO) was added as a solid lubricant. The tribo-physical and tribo-chemical changes on sliding surfaces were studied in detail. The results indicated that adhesive and abrasive wear due to removal of metallic binder and pullout of grains were the dominant wear mechanisms at room temperature. At high temperature, tribo-physical changes (i.e. mechanical mixture or sintering) and tribo-chemical reactions including complex reaction and oxidations were induced by frictional heat combined with high environmental temperature. As a result PbWO₄ was formed as a reaction product and acted as a solid lubricant. PbWO₄ and tribo-oxides along with the physically changed layer on the sliding surfaces were favorable to reduce wear of both materials. At high temperature, the wear mechanism varied from adhesive and abrasive wear at room temperature to lubrication by tribo-layer containing PbWO₄, tribo-oxides, and sinters at high temperatures.     

Combined convection and radiation in a tube with circumferential fins and circular disks

Combined convection and radiation heat transfer in a circular tube with circumferential fins and circular disks is investigated for various operating conditions. Using a finite volume technique for steady laminar flow, the governing equations are solved in order to study the flow and temperature fields. TheP- 1 approximation and the weighted sum of gray gases model (WSGGM) are used for solving the radiation transport equation. The results show that the total Nusselt number of combined convection and radiation is higher than that of pure convection. If the temperatures of the combustion gas and the wall in a tube are high, radiation becomes dominant. Therefore, it is necessary to evaluate the effect of radiation on the total heat transfer. Key Words: Convection, Radiation, Nongray Radiation,P- 1 Approximation, Weighted Sum of Gray Gases Model