Laminar burning velocity and explosion index of LPG-air and propane-air mixtures

The determination of burning velocity is very important for the calculations used in hazardous waste explosion protection and fuel tank venting, which has a direct impact on environmental protection. The scope of the present study encompass an extensive study to map the variations of the laminar burning velocity and the explosion index of LPG-air and propane-air mixtures over wide ranges of equivalence ratio (Φ = 0.7-2.2) and initial temperature (T_i = 295-400 K) and pressure (P_i = 50-400 kPa). For this purpose a cylindrical combustion bomb was developed. The reliability and accuracy of the built up facility together with the calculation algorithm are confirmed by comparing the values of the laminar burning velocity obtained for a standard fuel (propane at normal pressure normal temperature conditions, NPT) with those available in the literature. The burning velocity was determined using different models depending on the pressure history (P-t) of the central ignition combustion process at the minimum ignition energy.The data obtained for the laminar burning velocity is correlated to S_L = S_L0(T/T₀)^α(P/P₀)^β where S_L0 is the burning velocity at NPT, α and β are the temperature and pressure exponents respectively. The value of β is observed to slightly vary with the equivalence ratio for both fuels. However, propane exhibits higher pressure dependency than that of LPG. The maximum laminar burning velocity found for propane is nearly 455 mm/s at Φ = 1.1, while that for LPG is nearly 432 mm/s at 4.5% fuel percent (Φ ≈ 1.5). The maximum explosion index, commonly called the “explosion severity parameter”, is calculated from the determined laminar burning velocity and is found to be 93 bar m/s for propane, and nearly 88 bar m/s for LPG.     

Combustion characteristics of droplets composed of light cycle oil and diesel light oil in a hot-air chamber☆

Development of gas turbines fueled with light cycle oil (LCO) and oil mixture of LCO and diesel light oil (LO) requires an understanding of the droplet burning and vaporization characteristics of those oils. The present study is devoted to comparing the burning characteristics of isolated fuel droplets composed of an LCO and an LO. The tests were conducted in an atmospheric hot-air chamber preset at 1173 K, and the examined LCO had a lower cetane number but higher volatility and aromatics content compared to LO. It was demonstrated that the burning of the LCO droplet was sootier, while that of the LO droplet was more disruptive. At the tested temperature, coke formation was indistinct for both the oils, whereas slightly higher ignition delay time was shown for the LO droplet. The microexplosive burning more or less complicated the time-series droplet size d, an explicit burning rate constant, however, was still definable according to the d²-law to show the overall regression speed of the droplet surface area d² with burning time t. The rate constant exhibited little difference for smaller LCO and LO droplets but was greater for LO when the droplet was larger. The rate constant also gradually increased with increasing the initial droplet diameter d₀, which caused the relative size d/d₀ to be unified (normalized) into a single curve by a burning time t/d₀ⁿ (1.0<n<2.0). Analysis revealed that this unification resulted from the respective overlaps of the unsteady and quasi-steady burning phases for differently sized droplets. Further, it was clarified that the unification and analysis are generally valid to isolated liquid fuel droplet burning in hot ambiences.     

The effect of temperature on the adiabatic laminar burning velocities of CH4−air and H2−air flames

The effect of temperature on the adiabatic laminar burning velocities of CH₄ + air and H₂ + air flames was analyzed. Available measurements were interpreted using correlation S_L = S_L0 (T/T₀)^α. Particular attention was paid to the variation of the power exponent α with equivalence ratio at fixed (atmospheric) pressure. Experimental data and proposed empirical expressions for α as a function of equivalence ratio were summarized. They were compared with predictions of detailed kinetic models in methane + air and hydrogen + air flames. Unexpected non-monotonic behavior of α was found in rich methane + air flames. Modeling results are further examined using sensitivity analysis to elucidate the reason of particular dependences of the power exponent α on equivalence ratio.     

A simplified correlation for fixed bed pressure drop

An experimental study was conducted on the pressure drop characteristics of a variety of vertical packed beds in turbulent flow of air. The materials of different particle diameter, D_p, with a range of sphericity Φ, 0.55 ≤ Φ ≤ 1.00 were used in random loose packing to produce beds of different lengths, L, with a range of porosity, ε, 0.36 ≤ ε ≤ 0.56. In the covered test cases the cross-sectional velocity distribution at the exit plane of the packed beds and the pressure drop ΔP_Bed were measured in a particle Reynolds number range of Re_p, 675 ≤ Re_p ≤ 7772. The particular emphasis of the study was given to determine the influence of ε, Φ, D_p, L, Re_p on ΔP_Bed. In this respect the measurements of ΔP_Bed were compared with the well-known Ergun's Equation and the data were expressed in terms of correlations through introduced dimensionless parameters of pressure coefficient, ΔP^? and exit Reynolds number Re_exit. The proposed correlations of ΔP^? = ΔP^?(εRe_pD_p / L) and Re_exit = Re_exit(Re_pD_p / L) are found to be appropriate for the determination of ΔP_Bed and mean exit velocity, U, respectively with an acceptable fit of experimental data in an error margin less than ± 20%. The methodology is presented in this paper as an alternative approach to the available literature on packed beds.     

Two-dimensional unsteady flow of power-law fluids over a cylinder

The unsteady flow of incompressible power-law fluids over an unconfined circular cylinder in cross-flow arrangement has been studied numerically. The two-dimensional (2-D) field equations have been solved using a finite volume method based solver (FLUENT 6.3). In particular, the effects of the power-law index (0.4?n?1.8) and Reynolds number (40?Re?140) on the detailed kinematics of the flow (streamline, surface pressure and vorticity patterns) and on the macroscopic parameters (drag and lift coefficients, Strouhal number) are presented in detail. The periodic vortex shedding and the evolution of detailed kinematics with time are also presented to provide insights into the nature of flow. The two-dimensional flow transits from steady to unsteady behaviour at a critical value of the Reynolds number Re∼(40-50) and the von-Karman vortex street is observed beyond the critical Reynolds number (Re). Obviously, both the lift coefficient and Strouhal number values are zero for the steady flow, but their values increase with the increasing Reynolds number (Re) in the unsteady flow regime. For highly shear-thickening fluids (n=1.8), the flow becomes unsteady at Re=40 and unsteadiness in the flow appears at Re=50 for all values of power-law index (n). As expected, the evolution of the kinematics and vortex shedding show a complex dependence on the flow parameters near the transition in the flow. For a fixed value of the Reynolds number (Re), the drag coefficient increases and lift coefficient decreases with increasing value of the power-law index (n). For a fixed value of the power-law index (n), the drag coefficient gradually increases with the Reynolds number (Re). Similar to the drag coefficient, lift coefficient also shows a complex dependence on the power-law index (n) near the transition zone. The value of the Strouhal number (St) decreases with the increasing value of the power-law index (n) at a fixed value of the Reynolds number (Re).     

Effects of initial temperature on the structure of laminar CH4-air premixed flames

The growing popularity of natural gas as a eco-friendly fuel, is of paramount motivation of present investigation. In the present paper, the effect of initial temperature on the flame structure have been investigated in which laminar one-dimensional planar propagating flames of CH₄/air mixtures is simulated numerically using detailed chemical kinetic scheme and realistic transport models. The burning velocities are fundamentally important in developing models to predict progress of combustion. Hence, the burning velocities as a function of initial temperature of unburnt gas have been computed for stoichiometric mixture. The present predictions of burning velocities are compared with reported experimental data of Stone et al. [Combust. Flame. 114 (1998) 546], Hill and Huang [Combust. Sci. Technol. 60 (1980) 7] and Rallis and Garforth [Combust. Flame 31 (1978) 53]. The present prediction lies within the scatter of experimental data. A correlation in the form of S_u/S_u,0=(T_u/T_u,0)^1.575 has been developed to describe the dependence of initial temperature on the burning velocity for stoichiometric mixture. The structures of flame are investigated in details for initial temperature of 300 and 600 K which clearly indicate that detailed chemical kinetics are essential for prediction of the effects of initial temperature on the burning velocities. The present study will help in designing and developing the regenerative combustion systems.     

Electrochemical incineration of p-cresol and o-cresol in the filter-press-type FM01-LC electrochemical cell using BDD electrodes in sulfate media at pH 0

This work shows a comparative study of the incineration of 2-mM p-cresol and o-cresol in 1 M-H₂SO₄ in aqueous media. Microelectrolysis studies indicated that both the p-cresol and o-cresol oxidation were carried out via hydroxyl radicals (OH) formed by water oxidation in the boron-doped diamonds (BDD)-H₂O-H₂SO₄-p-cresol and o-cresol interface. In both cases, the potential and current density ranges, where great amounts of OH are formed, were between 2.3 V ≤ E ≤ 2.75 V versus SHE and J = 10 mA cm⁻². Electrolyses in an undivided FM01-LC reactor were performed at different Reynolds values 27,129 ≤ Re ≤ 42,631, and at J = 10 mA cm⁻². For p-cresol and o-cresol, the rate of degradation was slow, however it increases slightly as a function of the Re, indicating that the oxidation involves a complex pathway; current efficiency also rises as a function of the Re. For p-cresol, the mineralization at Re = 42,631 reached 90%, with 71% current efficiency and an energy consumption of 7.84 kWh m⁻³; whereas o-cresol was mineralized to 84%, with 67% current efficiency and an energy consumption of 6.56 kWh m⁻³. The results obtained in this work demonstrated that o-cresol is more recalcitrant than p-cresol.     

Influence of droplet size on the release of atomic sodium from a burning black liquor droplet in a flat flame

The influence of initial droplet size on the release of atomic sodium from black liquor solids (BLS) during each stage of black liquor combustion has been assessed using a planar laser-induced fluorescence (PLIF) technique. Three different initial diameters of black liquor droplets, 1.3, 1.7 and 2.2 mm were burned in a flat flame at equivalence ratios of 0.8, 0.9 and 1.25. The temporal release of the atomic sodium under fuel rich conditions was found to be different from that under fuel-lean conditions, especially during the smelt coalescence stage. For each stage of black liquor combustion, the measured release rate of atomic sodium increases with decreasing d_i. This implies that significant release of atomic sodium could occur during the in-flight combustion of small droplets, which are known to be generated in recovery boilers from either the carryover or the ejecta.     

Drag of non-spherical solid particles of regular and irregular shape

The drag of a non-spherical particle was reviewed and investigated for a variety of shapes (regular and irregular) and particle Reynolds numbers (Re_p). Point-force models for the trajectory-averaged drag were discussed for both the Stokes regime (Re_p ? 1) and Newton regime (Re_p ? 1 and sub-critical with approximately constant drag coefficient) for a particular particle shape. While exact solutions were often available for the Stokes regime, the Newton regime depended on: aspect ratio for spheroidal particles, surface area ratio for other regularly-shaped particles, and min-med-max area for irregularly shaped particles. The combination of the Stokes and Newton regimes were well integrated using a general method by Ganser (developed for isometric shapes and disks). In particular, a modified Clift-Gauvin expression was developed for particles with approximately cylindrical cross-sections relative to the flow, e.g. rods, prolate spheroids, and oblate spheroids with near-unity aspect ratios. However, particles with non-circular cross-sections exhibited a weaker dependence on Reynolds number, which is attributed to the more rapid transition to flow separation and turbulent boundary layer conditions. Their drag coefficient behavior was better represented by a modified Dallavalle drag model, by again integrating the Stokes and Newton regimes. This paper first discusses spherical particle drag and classification of particle shapes, followed by the main body which discusses drag in Stokes and Newton regimes and then combines these results for the intermediate regimes.     

Modeling and design of transdermal drug delivery patches containing an external heating device

Process modeling and design concepts were implemented to aid in the manufacturing of heat-enhanced transdermal drug-delivery systems. The simulated prototype consists of a corticosterone-loaded polymer patch applied to the skin and connected to a heating device in which an exothermic reaction occurs. To achieve a desired transdermal flux of 1.2 × 10⁻⁵ mg/cm² h, this contribution focuses on the influences of the (1) initial reaction rate (−r_A0), (2) mass of filler material in the device (m), (3) initial concentration (C₀) of medicament in the patch and (4) overall heat transfer coefficient (U). A regression technique yielded the following results: −r_A0 = 3.000 × 10⁻² kg/m³ s, m = 1.251 × 10⁻⁸ kg, U = 6.124 × 10 J/m² K s and C₀ = 1.966 × 10⁻¹ kg/m³. When m was fixed at 12.5 g, the optimum design required the following specifications: r_A0 = 2.765 × 10⁻² kg/m³ s, U = 1.402 × 10³ J/m² K s and C₀ = 1.941 × 10⁻¹ kg/m³. The priority (S_i) of the input factors (i) in reaching the target delivery rate is: SC₀>S−rA0>Sm>SU.     

Thermoelectric properties of layered Ca3.95RE0.05Mn3O10 compounds (RE = Ce,Nd, Sm,Eu, Gd,Dy)

In this work, polycrystalline samples of the substituted n = 3 Ruddlesden-Popper Ca_4−xRE_xMn₃O₁₀ phase were prepared by solid-state reaction (RE, rare earth = Ce, Nd, Sm, Eu, Gd, Dy). Single phased samples were synthesized for sintering times larger than 150 h at 1350 °C. Complete thermoelectric characterizations were performed from 5 to 390 K, in terms of electrical resistivity (ρ), Seebeck coefficient (S) and thermal conductivity (κ). As expected, the substitution of Ca by different rare earth elements leads to a significant modification of the thermoelectric properties. With substitution level as low as 1.25 at.%, a remarkable decrease of the electrical resistivity is observed. The influence of this cationic substitution on the thermal conductivity (κ), Seebeck coefficient (S), and the figure of merit ZT is also discussed. In this study, the best one reaches 5.8 × 10⁻³ at 300 K for the Ca_3.95Eu_0.05Mn₃O₁₀ composition, a value 6 times higher than the ZT exhibited by the beginning Ca₄Mn₃O₁₀ sample.     

Calcination and sintering characteristics of limestone under O2/CO2 combustion atmosphere

The O₂/CO₂ coal combustion technology is an innovative combustion technology that can control CO₂, SO₂ and NO_x emissions simultaneously. Calcination and sintering characteristics of limestone under O₂/CO₂ atmosphere were investigated in this paper. The pore size, the specific pore volume and the specific surface area of CaO calcined were measured by N₂ adsorption method. The grain size of CaO calcined was determined by XRD analysis. The specific pore volume and the specific surface area of CaO calcined in O₂/CO₂ atmosphere are less than that of CaO calcined in air at the same temperature. And the pore diameter of CaO calcined in O₂/CO₂ atmosphere is larger than that in air. The specific pore volume and the specific surface area of CaO calcined in O₂/CO₂ atmosphere increase initially with temperature, and then decline as temperature exceeds 1000 °C. The peaks of the specific pore volume and the specific surface area appear at 1000 °C. The specific surface area decreases with increase in the grain size of CaO calcined. The correlations of the grain size with the specific surface area and the specific pore volume can be expressed as L = 744.67 + 464.64 lg(1 / S) and L = − 608.5 + 1342.42 lg(1 / ε), respectively. Sintering has influence on the pore structure of CaO calcined by means of influencing the grain size of CaO.     

One-step fabrication of multihollow polystyrene particles from miniemulsion system with nonionic surfactant

This paper presented a new facile approach to fabricate polystyrene (PS) multihollow particles from miniemulsion polymerization. The surfactant used in this miniemulsion system was OP-10, a nonionic surfactant of nonyl phenyl polyoxyethylene with an average of 10 ethylene oxide units per molecule (hydrophilic-lipophilic balance, HLB = 13.9). Due to the partition characteristic of OP-10 in miniemulsion, W/O/W structured monomer droplets could be formed after ultrasonic dispersion. Under irradiation by γ-ray, styrene polymerized through a droplet nucleation mechanism, which was the feature of miniemulsion polymerization. The formation of multihollow structure was affected by the content of OP-10 (W_OP) and pre-added PS (W_PS). It was found that when W_OP was above 2 wt% and W_PS was in the region of 2-10 wt% based on the monomer, multihollow PS particles could be obtained. The molecular weight of the PS latex particles synthesized was determined to be 10⁵ g/mol by GPC.     

Determination of the elastic constants of portlandite by Brillouin spectroscopy

The single crystal elastic constants C_ij and the shear and adiabatic bulk modulus of a natural portlandite (Ca(OH)₂) crystal were determined by Brillouin spectroscopy at ambient conditions. The elastic constants, expressed in GPa, are: C₁₁ = 102.0(± 2.0), C₁₂ = 32.1(± 1.0), C₁₃ = 8.4(± 0.4), C₁₄ = 4.5(± 0.2), C₃₃ = 33.6(± 0.7), C₄₄ = 12.0(± 0.3), C₆₆ = (C₁₁-C₁₂)/2 = 35.0(± 1.1), where the numbers in parentheses are 1σ standard deviations. The Reuss bounds of the adiabatic bulk and shear moduli are K_0S = 26.0(± 0.3) GPa and G₀ = 17.5(± 0.4) GPa, respectively, while the Voigt bounds of these moduli are K_0S = 37.3(± 0.4) GPa and G₀ = 24.4(± 0.3) GPa. The Reuss and Voigt bounds for the aggregate Young's modulus are 42.8(± 1.0) GPa and 60.0(± 0.8) GPa respectively, while the aggregate Poisson's ratio is equal to 0.23(± 0.01). Portlandite exhibits both large compressional elastic anisotropy with C₁₁/C₃₃ = 3.03(± 0.09) equivalent to that of the isostructural hydroxide brucite (Mg(OH)₂), and large shear anisotropy with C₆₆/C₄₄ = 2.92(± 0.12) which is 11% larger than brucite. The comparison between the bulk modulus of portlandite and that of lime (CaO) confirms a systematic linear relationship between the bulk moduli of brucite-type simple hydroxides and the corresponding NaCl-type oxides.     

The Effect of Ambient Conditions on the Burning Rate of Gel Fuel Droplets

An experimental study was conducted to determine the dependence of the burning rate coefficient of gel fuel droplets on the pressure at different ambient oxygen nitrogen mixtures. Experiments were conducted using a pressure chamber, in which the droplet was suspended and the combustion process was video‐photographed by a high‐speed digital video camera. The tests were conducted at pressures between 0.1–4 MPa at different ambient oxygen nitrogen compositions (air, 40 % O₂ – 60 % N₂, and 60 % O₂ – 40 % N₂). The fuel was a compound of 95 % kerosene and 5 % gellant. At sub‐critical pressure conditions, the burning rate coefficient was found to increase with increasing ambient oxygen mass fraction. At supercritical conditions, no dependence of the burning rate coefficient on the ambient mixture was found. The results indicate that the burning rate coefficient depends on the oxygen partial pressure, at least at low pressures.     

Single droplet drying: Transition from the effective diffusion model to a modified receding interface model

Drying experiments on single droplets of aqueous amorphous polymer solution show morphological changes towards the end of drying that result in an under-prediction of the drying rate using an effective diffusion based model. Alternately, other researchers argue that the receding interface model more accurately reflects the physics of drying by predicting a fixed droplet radius once a specified surface condition is reached, usually the saturation concentration. However, this surface condition is not adequate for many skin forming materials. The conditions at which droplet radial contraction ceases will be determined by the balance between internal moisture loss causing a collapsing pressure and the mechanical strength of the surface skin. Because measurements and prediction of surface stress are difficult, it is proposed that they are related to the state of the polymer solution at the surface which is defined by the proximity of the surface temperature to its glass transition temperature, (T − T_g). In this work, an effective diffusion model is used to predict ideal shrinkage until a critical temperature difference or (T − T_g)_crit is reached where the surface of the droplet becomes fixed and the skin grows towards the droplet centre, that is, as a receding interface. For maltodextrin DE5, a (T − T_g)_crit of 20 °C was found to provide an accurate prediction of the drying rate. While these results show (T − T_g)_crit is indicative of mechanical stress development, it points to a need for further understanding of mechanical stress development in skin forming polymers during drying.     

A numerical study of particle wall-deposition in a turbulent square duct flow

The deposition of dense solid particles in a downward, fully developed turbulent square duct flow at Re_τ = 360, based on the mean friction velocity and the duct width, is studied using large eddy simulations of the fluid flow. The fluid and the particulate phases are treated using Eulerian and Lagrangian approaches, respectively. A finite-volume based, second-order accurate fractional step scheme is used to integrate the incompressible form of the unsteady, three-dimensional, filtered Navier-Stokes equations on an 80 × 80 × 128 grid. A dynamic subgrid kinetic energy model is used to account for the unresolved scales. The Lagrangian particle equation of motion includes the drag, lift, and gravity forces and is integrated using the fourth-order accurate Runge-Kutta scheme. Two values of particle to fluid density ratio (ρ_p/ρ_f = 1000 and 8900) and five values of dimensionless particle diameter (d_p/δ × 10⁶ = 100, 250, 500, 1000 and 2000, δ is the duct width) are studied. Two particle number densities, consisting of 10⁵ and 1.5 × 10⁶ particles initially in the domain, are examined.Variations in the probability distribution function (PDF) of the particle deposition location with dimensionless particle response time, i.e. Stokes number, are presented. The deposition is seen to occur with greater probability near the center of the duct walls, than at the corners. The average streamwise and wall-normal deposition velocities of the particles increase with Stokes number, with their maxima occurring near the center of the duct wall. The computed deposition rates are compared to previously reported results for a circular pipe flow. It is observed that the deposition rates in a square duct are greater than those in a pipe flow, especially for the low Stokes number particles. Also, wall-deposition of the low Stokes number particles increases significantly by including the subgrid velocity fluctuations in computing the fluid forces on the particles. Two-way coupling and, to a greater extent, four-way coupling are seen to increase the deposition rates.     

Synthesis and electrical conductivity of La-Sr-X-Mg-O (X = Ti, Zr, Al) perovskite solid solution

Perovskite solid solutions of (La_0.6Sr_0.4)(X_1−yMg_y)O_3−δ (X = Ti, Zr, Al) were prepared by a coprecipitation method using corresponding aqueous solutions and ammonium carbonate solution. The freeze-dried powders were sintered in air at 1000-1500 °C for 1-36 h. Single phase solid solutions were produced in the compositions of (La_0.6Sr_0.4)(Zr_0.6Mg_0.4)O_3−δ and (La_0.6Sr_0.4)(Al_0.9Mg_0.1)O_3−δ where (3 − δ) < 3. For the compositions of X = Ti and Zr for y = 0.1 where (3 − δ) > 3, two phases including perovskite solid solution were produced at 1400-1500 °C. The stability of perovskite solid solution was closely related to the fraction of lattice oxygen atom (3 − δ). A relatively high conductivity was measured for (La_0.6Sr_0.4)(Al_0.9Mg_0.1)O_3−δ (σ = 4.15 × 10⁻⁴ S/cm at 600 °C, activation energy 113.4 kJ/mol). The influence of fraction of oxide ion vacancy on the activation energy was small for δ = 0.1-0.3 of perovskite solid solution.     

Preparation of poly(ethylene terephthalate) nonwoven fabric from endless microfibers obtained by CO2 laser-thinning method

Poly(ethylene terephthalate) (PET) nonwoven fabric was prepared from microfibers obtained by using a carbon dioxide laser-thinning method. The PET nonwoven fabric obtained was made of the endless mircofibers with a uniform diameter without droplets. The fiber diameter can be varied by controlling airflow rate into the air jet and supplying speed of an original fiber into a laser-irradiating point. The fiber diameter decreased, and its birefringence increased as the airflow rate increased and the supplying speed decreased. When the microfiber prepared by irradiating the laser operated at a power density of 4.8 W cm⁻² to the original fiber supplied at S_s = 0.15 m min⁻¹ was dragged at an airflow rate of 30 L min⁻¹, the thinnest microfiber with a diameter of 3.6 μm was obtained.     

Variations in oxygenated blend composition to meet energy and combustion characteristics very similar to the diesel fuel

By the method of data collation, research into changes in life histories (ignition delay plus time of combustion) of the compounded fuel droplets (diesel fuel-biodiesel fuel (RME)-bioethanol), as well as diesel engine D-144 brake specific fuel consumption rates was performed and obtained results being compared to diesel fuel by an analogous manner.An optimum composition of the multi-component blend B30 + 7.5E demonstrating specific fuel consumption rates and droplet combustion characteristics very similar to diesel fuel was derived. In comparison to B30, a newly derived combustible blend demonstrated fairly improved emissions of exhaust gases. For low load mode: smoke opacity (−10%), NO_X (−2%), CO (−20%), and HC (−12.5%). For average load mode: smoke opacity (−10%), NO_X (−2%), CO (−22%), and HC (−14.5%). For high load mode: smoke opacity (−18%), NO_X (−2%), CO (−22%), and HC (−18%).