Advances in Large-eddy Simulation of Two-phase Combustion (I) LES of Spray Combustion

Spray combustion is widely used in power, transportation, chemical and metallurgical, iron and steel making, aeronautical and astronautical engineering. In recent years, large-eddy simulation (LES) becomes more and more attractive, because it can give the instantaneous flow and flame structures, and may give more accurate statistical results than the Reynolds averaged Navier-Stokes (RANS) modeling. In this paper, the present status of the studies on LES of spray combustion is reviewed, and the future research needs are discussed.

Advances in LES of Two-phase Combustion(Ⅱ) LES of Complex Gas-Particle Flows and Coal Combustion

Large-eddy simulation(LES) is under its rapid development and is recognized as a possible second generation of CFD methods used in engineering.Large-eddy simulation of two-phase flows and combustion is particularly important for engineering applications.Some investigators,including the present authors,give their review on LES of spray combustion in gas-turbine combustors and internal combustion engines.However,up to now only a few papers are related to the state-of-the-art on LES of gas-particle flows and combustion.In this paper a review of the advances in LES of complex gas-particle flows and coal combustion is presented.Different sub-grid scale(SGS) stress models and combustion models are described,some of the main results are summarized,and some research needs are discussed.

A Review on SOM-LES of Turbulent Two-Phase Combustion

Turbulent two-phase combustion is widely encountered in spray and pulverized-coal combustors,and large-eddy simulation(LES)becomes a powerful CFD method for its simulation,because LES can give unsteady flame structures and more reasonable statistical results than Reynolds-averaged modeling.Present combustion models in LES either lack of generality or are computationally too expensive.A statistical moment model based on the idea of turbulence modeling called“second-order moment(SOM)combustion model”was developed by the present authors for LES of two-phase combustion.In this paper,a review is given on our published research results for SOM-LES of two-phase combustion,including the description of the SOM-LES model,its application,validation of statistical results by experiments,as well as the phenomena obtained by instantaneous results.

Heat transfer model of two-phase flow across tube bundle in submerged combustion vaporizer

In order to optimize the design of the submerged combustion vaporizer(SCV), an experimental apparatus was set up to investigate the heat transfer character outside the tube bundle in SCV. Several experiments were conducted using water and CO_2 as the heat transfer media in the tubes, respectively. The results indicated that hot air flux, the initial liquid level height and the tube pitch ratio had great influence on the heat transfer coefficient outside the tube bundle(ho). Finally, the air flux associated factor β and height associated factor γ were introduced to propose a new hocorrelation. After verified by experiments using cold water, high pressure CO_2 and liquid N_2 as heat transfer media, respectively, it was found that the biggest deviation between the predicted and the experimental values was less than 25%.

NUMERICAL MODELING OF 3-D TURBULENT TWO-PHASE FLOW AND COAL COMBUSTION IN A PULVERIZED-COAL COMBUSTOR

In the present paper, a multifluid model of two-phase flows with pulverized-coal combustion, based on a continuum-trajectory model with reacting particle phase, is developed and employed to simulate the 3-D turbulent two-phase hows and combustion in a new type of pulverized-coal combustor with one primary-air jet placed along the wall of the combustor. The results show that: (1) this continuum-trajectory model with reacting particle phase can be used in practical engineering to qualitatively predict the flame stability, concentrations of gas species, possibilities of slag formation and soot deposition, etc.; (2) large recirculation zones can be created in the combustor, which is favorable to the ignition and flame stabilization.

PREDICTIONS OF 3-D STRONGLY SWIRLING GAS-SOLID TWO-PHASE FLOW WITH GAS COMBUSTION

ＰＲＥＤＩＣＴＩＯＮＳＯＦ３┐ＤＳＴＲＯＮＧＬＹＳＷＩＲＬＩＮＧＧＡＳ┐ＳＯＬＩＤＴＷＯ┐ＰＨＡＳＥＦＬＯＷＷＩＴＨＧＡＳＣＯＭＢＵＳＴＩＯＮＷａｎｇＺｈｅｎｙｕ（王振宇）（ＳｈａｎｇｈａｉＷｕｊｉｎｇＴｈｅｒｍａｌＰｏｗｅｒＰｌａｎｔ）ＨｕａｎＢｏ...

NUMERICAL MODELING OF SUSPENSION TWO-PHASE FLOWS AND THEIR COMBUSTION CURRENT STATUS AND DEVELOPMENT

The background, history and development of numerical modeling of suspension two-phase flows and their combustion are reviewed. Different models, including single-fluid model, particle trajectory model, multifluid model and continuum-trajectory model, their numerical procedure and application are discussed. The advantage and drawback of these models are analyzed. Finally, the prospect of development is described and the contributions made by the present author are also incorporated.

LARGE-EDDY SIMULATION OF TWO-PHASE REACTING FLOW IN MODEL COMBUSTOR

The gas-droplet two-phase reacting flow in a model combustor with the V-gutter flame holder is studied by an Eulerian-Lagrangian large-eddy simulation (LES) approach. The k-equation subgrid-scale model is used to simulate the subgrid eddy viscosity, and the eddy-break-up (EBU) combustion subgrid-scale model is used to determine the chemical reaction rate. A two-step turbulent combustion subgrid-scale model is employed for calculating carbon monoxide CO concentration, and the NO subgrid-scale pollutant formation model for the evaluation of the rate of NO formation. The heat flux model is applied to the prediction of radiant heat transfer. The gas phase is solved with the SIMPLE algorithm and a hybrid scheme in the staggered grid system. The liquid phase equations are solved in a Lagrangian frame in reference of the particle-source-in-cell (PSIC) algorithm. From simulation results, the exchange of mass, moment and energy between gas and particle fields for the reacting flow in the afterburner with a V-gutter flame holder can be obtained. By the comparison of experimental and simulation results, profile temperature and pollutant of the outlet are quite in agreement with experimental data. Results show that the LES approach for predicting the two-phase instantaneous reacting flow and pollutant emissions in the afterburner is feasible.

Unraveling the role of dual Ti/Mg metals on the ignition and combustion behavior of HTPB-boron-based fuel

Metal additives play an essential role in explosive and propellant formulations. Boron(B) is widely used in propellant applications owing to its high energetic content. The addition of B to explosives and propellants increases their energy density, making them more efficient and powerful. Nevertheless, B forms oxide layers on its surface during combustion, slowing down the combustion rate and reducing rocket motor efficiency. To overcome this issue, other metal additives such as aluminum(Al), magnesium(Mg),and titanium(Ti) are revealed to be effective in boosting the combustion rate of propellants. These additives may improve the combustion rate and therefore enhance the rocket motor’s performance. The present study focused on preparing and investigating the ignition and combustion behavior of pure hydroxyl-terminated polybutadiene(HTPB)-B fuel supplemented with nano-titanium and nanomagnesium. The burn rates of HTPB-B fuel samples were evaluated on the opposed flow burner(OFB)under a gaseous oxygen oxidizer, for which the mass flux ranges from 22 kg/(m^(2)·s) to 86 kg/(m^(2)·s). The addition of Ti and Mg exhibited higher regression rates, which were attributed to the improved oxidation reaction of B due to the synergetic metal combustion effect. The possible combustion/oxidation reaction mechanism of B-Mg and B-Ti by heating the fuel samples at 900℃ and 1100℃ was also examined in a Nabertherm burnout furnace under an oxygen atmosphere. The post-combustion products were collected and further subjected to X-ray diffraction(XRD) and field emission scanning electron microscopy(FE-SEM) analyses to inspect the combustion behavior of B-Ti and B-Mg. It has been observed that the B oxide layer at the interface between B-Ti(B-Mg) is removed at lower temperatures, hence facilitating oxygen transfer from the surroundings to the core B. Additionally, Ti and Mg decreased the ignition delay time of B, which improved its combustion performance.

Fabrication and characterization of multi-scale coated boron powders with improved combustion performance:A brief review

Boron has high mass and volume calorific values,but it is difficult to ignite and has low combustion efficiency.This literature review summarizes the strategies that are used to solve the above-mentioned problems,which include coatings of boron by using fluoride compounds,energetic composites,metal fuels,and metal oxides.Coating techniques include recrystallization,dual-solvent,phase transfer,electrospinning,etc.As one of the effective coating agents,the fluorine compounds can react with the oxide shell of boron powder.In comparison,the energetic composites can effectively improve the flame temperature of boron powder and enhance the evaporation efficiency of oxide film as a condensed product.Metals and metal oxides would react with boron powder to form metal borides with a lower ignition point,which could reduce its ignition temperature.

The regulation of ferrocene-based catalysts on heat transfer in highpressure combustion of ammonium perchlorate/hydroxyl-terminated polybutadiene/aluminum composite propellants

The regulation of the burning rate pressure exponent for the ammonium perchlorate/hydroxylterminated polybutadiene/aluminum(AP/HTPB/Al)composite propellants under high pressures is a crucial step for its application in high-pressure solid rocket motors.In this work,the combustion characteristics of AP/HTPB/Al composite propellants containing ferrocene-based catalysts were investigated,including the burning rate,thermal behavior,the local heat transfer,and temperature profile in the range of 7-28 MPa.The results showed that the exponent breaks were still observed in the propellants after the addition of positive catalysts(Ce-Fc-MOF),the burning rate inhibitor((Ferrocenylmethyl)trimethylammonium bromide,Fc Br)and the mixture of Fc Br/catocene(GFP).However,the characteristic pressure has increased,and the exponent decreased from 1.14 to 0.66,0.55,and 0.48 when the addition of Ce-FcMOF,Fc Br and Fc Br/GFP in the propellants.In addition,the temperature in the first decomposition stage was increased by 7.50℃ and 11.40℃ for the AP/Fc Br mixture and the AP/Fc Br/GFP mixture,respectively,compared to the pure AP.On the other hand,the temperature in the second decomposition stage decreased by 48.30℃ and 81.70℃ for AP/Fc Br and AP/Fc Br/GFP mixtures,respectively.It was also found that Fc Br might generate ammonia to cover the AP surface.In this case,a reaction between the methyl in Fc Br and perchloric acid caused more ammonia to appear at the AP surface,resulting in the suppression of ammonia desorption.In addition,the coarse AP particles on the quenched surface were of a concave shape relative to the binder matrix under low and high pressures when the catalysts were added.In the process,the decline at the AP/HTPB interface was only exhibited in the propellant with the addition of Ce-Fc-MOF.The ratio of the gas-phase temperature gradient of the propellants containing catalysts was reduced significantly below and above the characteristic pressure,rather than 3.6 times of the difference in the blank propellant.Overall,the obtained results demonstrated that the pressure exponent could be effectively regulated and controlled by adjusting the propellant local heat and mass transfer under high and low pressures.

Micro-aluminum powder with bi-or tri-component alloy coating as a promising catalyst:Boosting pyrolysis and combustion of ammonium perchlorate

A novel design of micro-aluminum(μAl)powder coated with bi-/tri-component alloy layer,such as:Ni-P and Ni-P-Cu(namely,Al@Ni-P,Al@Ni-P-Cu,respectively),as combustion catalysts,were introduced to release its huge energy inside Al-core and promote rapid pyrolysis of ammonium perchlorate(AP)at a lower temperature in aluminized propellants.The microstructure of Al@Ni-P-Cu demonstrates that a three-layer Ni-P-Cu shell,with the thickness of~100 nm,is uniformly supported byμAl carrier(fuel unit),which has an amorphous surface with a thickness of~2.3 nm(catalytic unit).The peak temperature of AP with the addition of Al@Ni-P-Cu(3.5%)could significantly drop to 316.2℃ at high-temperature thermal decomposition,reduced by 124.3℃,in comparison to that of pure AP with 440.5℃.It illustrated that the introduction of Al@Ni-P-Cu could weaken or even eliminate the obstacle of AP pyrolysis due to its reduction of activation energy with 118.28 kJ/mol.The laser ignition results showed that the ignition delay time of Al@Ni-P-Cu/AP mixture with 78 ms in air is shorter than that of Al@Ni-P/AP(118 ms),decreased by 33.90%.Those astonishing breakthroughs were attributed to the synergistic effects of adequate active sites on amorphous surface and oxidation exothermic reactions(7597.7 J/g)of Al@Ni-P-Cu,resulting in accelerated mass and/or heat transfer rate to catalyze AP pyrolysis and combustion.Moreover,it is believed to provide an alternative Al-based combustion catalyst for propellant designer,to promote the development the propellants toward a higher energy.

A reduced combustion mechanism of ammonia/diesel optimized with multi-objective genetic algorithm

For the deep understanding on combustion of ammonia/diesel,this study develops a reduced mechanism of ammonia/diesel with 227 species and 937 reactions.The sub-mechanism on ammonia/interactions of N-based and C-based species(N—C)/NOx is optimized using the Non-dominated Sorting Genetic Algorithm II(NSGA-II)with 200 generations.The optimized mechanism(named as 937b)is validated against combustion characteristics of ammonia/methane(which is used to examine the accuracy of N—C interactions)and ammonia/diesel blends.The ignition delay times(IDTs),the laminar flame speeds and most of key intermediate species during the combustion of ammonia/methane blends can be accurately simulated by 937b under a wide range of conditions.As for ammonia/diesel blends with various diesel energy fractions,reasonable predictions on the IDTs under pressures from 1.0 MPa to5.0 MPa as well as the laminar flame speeds are also achieved by 937b.In particular,with regard to the IDT simulations of ammonia/diesel blends,937b makes progress in both aspects of overall accuracy and computational efficiency,compared to a detailed ammonia/diesel mechanism.Further kinetic analysis reveals that the reaction pathway of ammonia during the combustion of ammonia/diesel blend mainly differs in the tendencies of oxygen additions to NH_2 and NH with different equivalence ratios.

Evaluation of frictional pressure drop correlations for air-water and air-oil two-phase flow in pipeline-riser system

Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.

Effect and Mechanism of Rare Earth Hydrotalcite Inhibiting Coal Spontaneous Combustion

A hydrotalcite(layered double hydroxide, LDH) inhibitor which is suitable for the whole process of coal spontaneous combustion and a LDH inhibitor containing rare earth lanthanum elements were prepared. The inhibition effect and mechanism were analyzed by scanning electron microscopy(SEM),X-ray diffraction(XRD), thermal performance analysis, in-situ diffuse reflectance infrared spectroscopy and temperature-programmed experiment. The results have shown that the inhibitor containing lanthanum can play a good inhibitory role in every stage of coal oxidation. During the slow oxidation of coal samples, the inhibitor containing lanthanum ions can slow down the oxidation process of coal and increase the initial temperature of coal spontaneous combustion. At the same time, because the hydroxyl groups in LDHs are connected with-COO-groups on the coal surface through hydrogen bonds, the stability of coal is improved. With the increase of temperature, LDHs can remove interlayer water molecules and reduce the surface temperature of coal. CO release rate of coal samples decreases significantly after adding inhibitor containing lanthanum element, and the maximum inhibition rate of the inhibitor is 58.1%.

Evaluating Ignition and Combustion Performance with Al-Metal- Organic Frameworks and Nano-Aluminum in HTPB Fuel

Incorporating aluminum metal-organic frameworks(Al-MOFs)as energetic additives for solid fuels presents a promising avenue for enhancing combustion performance.This study explores the potential benefits of Al-MOF(MIL-53(Al))energetic additive on the combustion performance of hydroxyl-terminated polybutadiene(HTPB)fuel.The HTPB-MOF fuel samples were manufactured using the vacuum-casting technique,followed by a comprehensive evaluation of their ignition and combustion properties using an opposed flow burner(OFB)setup utilizing gaseous oxygen as an oxidizer.To gauge the effectiveness of Al-MOFs as fuel additives,their impact is compared with that of nano-aluminum(nAl),another traditional additive in HTPB fuel.The results indicate that the addition of 15%(mass fraction)nAl into HTPB resulted in the shortest ignition delay time(136 ms),demonstrating improved ignition performance compared to pure HTPB(273 ms).The incorporation of Al-MOF in HTPB also reduced ignition delay times to 227 ms and 189 ms,respectively.Moreover,under high oxidizer mass flux conditions(79—81 kg/(m^(2)s)),HTPB fuel with 15%nAl exhibited a substantial 83.2%increase in regression rate compared to the baseline HTPB fuel,highlighting the positive influence of nAl on combustion behavior.In contrast,HTPB-MOF with a 15%Al-MOF additive showed a 32.7%increase in regression rate compared to pure HTPB.These results suggest that HTPB-nAl outperforms HTPB-MOF in terms of regression rates,indicating a more vigorous and rapid burning behavior.

Investigation of oxy-fuel combustion for methane and acid gas in a diffusion flame

Co-combustion of methane(CH4)and acid gas(AG)is required to sustain the temperature in Claus reaction furnace.In this study,oxy-fuel combustion of methane and acid gas has been experimentally studied in a diffusion flame.Three equivalence ratios(ER=1.0,1.5,2.0)and CH_(4)-addition ratios(CH_(4)/AG=0.3,0.5,0.7)were examined and the flame was interpreted by analyzing the distributions of the temperature and species concentration along central axial.CH_(4)-AG diffusion flame could be classified into three sections namely initial reaction,oxidation and complex reaction sections.Competitive oxidation of CH_(4)and H_(2)S was noted in the first section wherein H_(2)S was preferred and both were mainly proceeding decomposition and partial oxidation.SO_(2)was formed at oxidation section together with obvious presence of H2 and CO.However,H2 and CO were inclined to be sustained under fuel rich condition in the complex reaction section.Reducing ER and increasing CH4/AG contributed to higher temperature,H_(2)S and CH_(4)oxidation and CO_(2)reactivity.Hence a growing trend for CH_(4)and AG to convert into H_(2),CO and SO_(2)could be witnessed.And this factor enhanced the generation of CS2 and COS in the flame inner core by interactions of CH4 and CO_(2)with sulfur species.COS was formed through the interactions of CO and CO_(2)with sulfur species.The CS_(2)production directly relied on reaction of CH_(4)with sulfur species.The concentration of COS was greater than CS_(2)since CS_(2)was probably inhibited due to the presence of H_(2).COS and CS_(2)could be consumed by further oxidation or other complex reactions.

Comparison of Combustion Characteristics of Tars Produced with Tobacco Stem Biomass Gasification

In order to study the combustion characteristics of tar in biomass gasifier inner wall and gasification gas,“tobacco stem semi-tar inside furnace”,“tobacco stem tar inside furnace”and“tobacco stem tar out-of-furnace”were subjected to thermogravimetric experiments,and the combustion characteristics and kinetic characteristics were analyzed.The result shows that“tobacco stem semi-tar inside furnace”has the highest value and“tobacco stem tar out-of-furnace”is has the lowest value on ignition characteristics,combustion characteristics and combustible stability;“tobacco stem semi-tar inside furnace”has the lowest value and“tobacco stem tar outside furnace”has the highest value on burnout characteristics;“tobacco stem tar outside furnace”has the highest value and“tobacco stem tar inside furnace”has the lowest value on integrated combustion characteristics.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Numerical Simulations of Wave Impact Forces on the Open-Type Sea Access Road Using A Two-Phase SPH Model

A numerical study based on a two-dimensional two-phase SPH(Smoothed Particle Hydrodynamics)model to analyze the action of water waves on open-type sea access roads is presented.The study is a continuation of the analyses presented by Chen et al.(2022),in which the sea access roads are semi-immersed.In this new configuration,the sea access roads are placed above the still water level,therefore the presence of the air phase becomes a relevant issue in the determination of the wave forces acting on the structures.Indeed,the comparison of wave forces on the open-type sea access roads obtained from the single and two-phase SPH models with the experimental results shows that the latter are in much better agreement.So in the numerical simulations,a two-phaseδ-SPH model is adopted to investigate the dynamical problems.Based on the numerical results,the maximum horizontal and uplifting wave forces acting on the sea access roads are analyzed by considering different wave conditions and geometries of the structures.In particular,the presence of the girder is analyzed and the differences in the wave forces due to the air cushion effects which are created below the structure are highlighted.