A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structu...A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations(LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3.More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.展开更多
Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling facto...Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.展开更多
By comparing the energy spectrum and total kinetic energy, the effects of numerical errors （which arise from aliasing and discretization errors）, subgrid-scale （SGS） models, and their interactions on direct numeri...By comparing the energy spectrum and total kinetic energy, the effects of numerical errors （which arise from aliasing and discretization errors）, subgrid-scale （SGS） models, and their interactions on direct numerical simulation （DNS） and large eddy simulation （LES） are investigated, The decaying isotropic turbulence is chosen as the test case. To simulate complex geometries, both the spectral method and Pade compact difference schemes are studied. The truncated Navier-Stokes （TNS） equation model with Pade discrete filter is adopted as the SGS model. It is found that the discretization error plays a key role in DNS. Low order difference schemes may be unsuitable. However, for LES, it is found that the SGS model can represent the effect of small scales to large scales and dump the numerical errors. Therefore, reasonable results can also be obtained with a low order discretization scheme.展开更多
Moving-particle semi-implicit(MPS) method is a new mesh-free numerical method based on Lagrangian particle. In this paper, MPS method is applied to the study on numerical wave tank. For the purpose of simulating numer...Moving-particle semi-implicit(MPS) method is a new mesh-free numerical method based on Lagrangian particle. In this paper, MPS method is applied to the study on numerical wave tank. For the purpose of simulating numerical wave, we combine the MPS method with large eddy simulation(LES) which can simulate the turbulence in the flow. The intense pressure fluctuation is a significant shortcoming in MPS method. So, we improve the original MPS method by using a new pressure Poisson equation to ease the pressure fluctuation. Divergencefree condition representing fluid incompressible is used to calculate pressure smoothly. Then, area-time average technique is used to deal with the calculation. With these improvements, the modified MPS-LES method is applied to the simulation of numerical wave. As a contrast, we also use the original MPS-LES method to simulate the wave in a numerical wave tank. The result shows that the new method is better than the original MPS-LES method.展开更多
Large eddy simulation (LES) explicitly calculates the large-scale vortex field and parameterizes the small-scale vortices.In this study,LES and κ-ε models were developed for a specific geometrical configuration of b...Large eddy simulation (LES) explicitly calculates the large-scale vortex field and parameterizes the small-scale vortices.In this study,LES and κ-ε models were developed for a specific geometrical configuration of backward-facing step (BFS).The simulation results were validated with particle image velocimetry (PIV) measurements and direct numerical simulation (DNS).This LES simulation was carried out with a Reynolds number of 9000 in a pressurized water tunnel with an expansion ratio of 2.00.The results indicate that the LES model can reveal largescale vortex motion although with a larger grid-cell size.However,the LES model tends to overestimate the top wall separation and the Reynolds stress components for the BFS flow simulation without a sufficiently fine grid.Overall,LES is a potential tool for simulating separated flow controlled by large-scale vortices.展开更多
The numerical simulation of modern aero-engine combustion chamber needs accurate description of the interaction between turbulence and chemical reaction mechanism. The Large Eddy Simulation(LES) method with the Transp...The numerical simulation of modern aero-engine combustion chamber needs accurate description of the interaction between turbulence and chemical reaction mechanism. The Large Eddy Simulation(LES) method with the Transported Probability Density Function(TPDF) turbulence combustion model is promising in engineering applications. In flame region, the impact of chemical reaction should be considered in TPDF molecular mixing model. Based on pioneer research, three new TPDF turbulence-chemistry dual time scale molecular mixing models were proposed tentatively by adding the chemistry time scale in molecular mixing model for nonpremixed flame. The Aero-Engine Combustor Simulation Code(AECSC) which is based on LES-TPDF method was combined with the three new models. Then the Sandia laboratory's methane-air jet flames: Flame D and Flame E were simulated. Transient simulation results show that all the three new models can predict the instantaneous combustion flow pattern of the jet flames. Furthermore,the average scalar statistical results were compared with the experimental data. The simulation result of the new TPDF arithmetic mean modification model is the closest to the experimental data:the average error in Flame D is 7.6% and 6.6% in Flame E. The extinction and re-ignition phenomena of the jet flames especially Flame E were captured. The turbulence time scale and the chemistry time scale are in different order in the whole flow field. The dual time scale TPDF combustion model has ability to deal with both the turbulence effect and the chemistry reaction effect, as well as their interaction more accurately for nonpremixed flames.展开更多
Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics h...Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.展开更多
Large Eddy Simulations(LES) in conjunction with the Flamelet Progress Variable(FPV) approach have been performed to investigate the flame and large-scale flow structures in the bluff-body stabilized non-premixed flame...Large Eddy Simulations(LES) in conjunction with the Flamelet Progress Variable(FPV) approach have been performed to investigate the flame and large-scale flow structures in the bluff-body stabilized non-premixed flames, HM1 and HM3. The validity of the numerical methods is first verified by comparing the predicted velocity and composition fields with experimental measurements. Then the evolution of the flame and large-scale flow structures is analyzed when the flames approach blow-off. The analysis of instantaneous and statistical data indicates that there exists a shift of the control mechanism in the recirculation zone in the two flames. In the recirculation zone, HM1 flame is mainly controlled by the mixing effect and ignition mainly occurs in the outer shear layer. In HM3 flame, both the chemical reactions and mixing are important in the recirculation zone. The Proper Orthogonal Decomposition(POD) results show that the fluctuations in the outer shear layer are more intense in HM1, while the flow structures are more obvious in the outer vortex structure in HM3, due to the different control mechanism in the recirculation zone.It further shows that the flow structures in HM1 spread larger in the intense mixing zone due to higher temperature and less extinction.展开更多
A non-linear turbulent hydrodynamic model was proposed by Zhan and Li (1993) using the large eddy simulation approach and directional eddy viscosity method. The model is simplified in this paper and is applied to the ...A non-linear turbulent hydrodynamic model was proposed by Zhan and Li (1993) using the large eddy simulation approach and directional eddy viscosity method. The model is simplified in this paper and is applied to the problem of a side discharge into open-channel flow,where a recirculation zone develops downstream of the discharge. The numerical results are compared with experimental results (Strazisar and Prahl (1973) and Mikhail, Chu, and Savage (1975)) and previous numerical results obtained using a k-model (McGuirk and Rodi, 1978). The agreement is good.展开更多
By the Volume of Fluid (VOF) multiphase flow model two-dimensional gravity currents with three phases including air are numerically simulated in this article. The necessity of consideration of turbulence effect for hi...By the Volume of Fluid (VOF) multiphase flow model two-dimensional gravity currents with three phases including air are numerically simulated in this article. The necessity of consideration of turbulence effect for high Reynolds numbers is demonstrated quantitatively by LES (the Large Eddy Simulation) turbulence model. The gravity currents are simulated for h not equal H as well as h = H, where h is the depth of the gravity current before the release and H is the depth of the intruded fluid. Uprising of swell occurs when a current flows horizontally into another lighter one for h not equal H. The problems under what condition the uprising of swell occurs and how long it takes are considered in this article. All the simulated results are in reasonable agreement with the experimental results available.展开更多
Bunsen burner is a typical geometry for investigating the turbulence-flame interaction.In most experimental studies,only turbulence intensity u′and integral scale l0 are used to characterize the turbulent flow field,...Bunsen burner is a typical geometry for investigating the turbulence-flame interaction.In most experimental studies,only turbulence intensity u′and integral scale l0 are used to characterize the turbulent flow field,regardless of the perforation geometry of perforated plates.However,since the geometry influences the developing process and vortex broken,the plate geometry has to be considered when discussing the flame-turbulence interaction.In order to investigate conditions at the same l0 and u′using different geometries,large eddy simulation of CH_(4)/air flames with dynamic TF combustion model was performed.The model validation shows good agreement between Large Eddy Simulation(LES)and experimental results.In the non-reacting flows,the Vortex Stretching of circular-perforated plate condition is always larger than that of slot-perforated plate condition,which comes from the stresses in the flow fields to stretch the vorticity vector.In reacting flows,at the root of the flame,the Vortex Stretching plays a major role,and the total vorticity here of circular-perforated plate condition is still larger(53.8%and 300%larger than that of the slot-perforated plate at x/D=0 and x/D=2.5,respectively).More small-scale vortex in circular-perforated plate condition can affect and wrinkle the flame front to increase the Probability Density Function(PDF)at large curvatures.The 3D curvature distributions of both cases bias to negative values.The negative trend of curvatures at the instant flame front results from the Dilatation term.Also,the value of the Vortex Stretching and the Dilatation at the flame front of circular-perforated plate condition is obviously larger.展开更多
A computational model of three-dimensional, time-dependent flame spread in microgravity environment is presented. The solid is assumed to be a thermally-thin, pyrolysing cellulosic sheet. The gas phase model includes ...A computational model of three-dimensional, time-dependent flame spread in microgravity environment is presented. The solid is assumed to be a thermally-thin, pyrolysing cellulosic sheet. The gas phase model includes the full Navier-Stokes equations with density and pressure variations and six-flus model of radiation heat transfer. The solid phase model consists of continuity and energy eqllations whose solution provides boundary conditions for the gas phase equations. In the numerical procedure, the gas-and solid-phase equations are solved sepaxately and iteratively at each time step. Predictions have been made of flame spread in slow forced flow under gravitational acceleration normal to fuel surface and flame spread in a quiescent environment in an enclosed chamber under gravitational acceleration parallel to fuel surface. Numerical simulations show that, under microgravity, slow-flow conditions, flame spread process is highly unsteady with the upstream flame spreads faster than the downstream flame after a period of ignition. It has also been shown that the level of microgravity has a significant effect on the name spread process.展开更多
At present,there is a shortage of experimental and simulation studies on fire spread in medium-and large-scale compartments while the existing models of the fire spread are limited for typical engineering applications...At present,there is a shortage of experimental and simulation studies on fire spread in medium-and large-scale compartments while the existing models of the fire spread are limited for typical engineering applications.This paper proposes a new model for large-scale fire spread on medium density fibreboard(MDF)panels.Validating the model with single burning item(SBI)experiments,it is found that the numerical simulation closely predicts the experimental heat release rate(HRR)with some error near the peak.The predicted heat flux and distance of lateral flame spread are consistent with the experiments and an existing model.The effects of kinetic properties and heat of combustion are identified through a sensitivity analysis.The decrease of activation energy and increase of pre-exponential factor make the MDF easier to pyrolyze and the increase of heat of combustion enhances the flame temperature and thus provide more heat feedback to the sample surface.The low activation energy(71.9 kJ/mol)and high heat of combustion(46.5 MJ/kg)of the model ensure the occurrence of flame spread.Furthermore,the model was validated using medium-scale compartment fire experiments and the results showed that the model can accurately predict the HRR after flashover(the error is within 7%).While the burner is ignited,the predictions of in-compartment gas temperature and heat flux are more accurate.However,when the burner is extinguished,the modelled in-compartment gas temperature is lower than the experimental values,resulting in a lower heat flux prediction.The model leads to easier flame spread;therefore,the modelled flame spreads faster in the compartment compared to the experiment,and thus the HRR increases more rapidly.展开更多
To improve the NO modelling in turbulent flames,the flamelet/progress variable(FPV)model is extended by introducing NO mass fraction into the progress variable and incorporating an additional NO transport equation.Two...To improve the NO modelling in turbulent flames,the flamelet/progress variable(FPV)model is extended by introducing NO mass fraction into the progress variable and incorporating an additional NO transport equation.Two sets of flamelet databases are tabulated with progress variables based on major species and NO mass fraction,respectively.The former is used for the acquisition of the main thermochemical variables,while the latter is employed for NO modelling.Moreover,an additional transport equation is solved to obtain the NO mass fraction,with the source term corrected using the scale similarity method.Model assessments are first conducted on laminar counterflow diffusion flames to identify lookup-related errors and assess the suitability of progress variable definitions.The results show that the progress variables based on major species and NO could correctly describe the main thermochemical quantities and NO-related variables,respectively.Subsequently,the model is applied to the large eddy simulation(LES)of Sandia flames.The results indicate that the extended FPV model improves the NO prediction,with a mean error for NO prediction at 55%,significantly lower than those of existing FPV models(130%and 385%).The LES with the extended FPV model quantitatively captures NO suppression in the mid-range of Reynolds numbers from 22400(Flame D)to 33600(Flame E),but underestimates the NO suppression at higher Reynolds numbers from 33600 to 44800(Flame F).This underprediction is primarily attributed to the underestimation of local extinction levels in flames with high Reynolds numbers.展开更多
Turbulent flow in a 3-D blade passage of a Francis hydro turbine was simulated with the Large Eddy Simulation （LES） to investigate the spatial and temporal distributions of the turbulence when strongly distorted wak...Turbulent flow in a 3-D blade passage of a Francis hydro turbine was simulated with the Large Eddy Simulation （LES） to investigate the spatial and temporal distributions of the turbulence when strongly distorted wakes in the inflow sweep over the passage, In a suitable consideration of the energy exchanging mechanism between the large and small scales in the complicated passage with a strong 3-D curvature, one-coefficient dynamic Sub-Grid-Scale （SGS） stress model was used in this article. The simulations show that the strong wakes in the inflow lead to a flow separation at the leading zone of the passage, and to form a primary vortex in the span-wise direction. The primary span-wise vortex evolves and splits into smaller vortex pairs due to the constraint of no-slip wall condition, which triggers losing stability of the flow in the passage. The computed pressures on the pressure and suction sides agree with the measured data for a working test turbine model.展开更多
This paper presents a numerical study on the hydrodynamic behaviours of a round buoyant jet under the effect of JONSWAP random waves. A three-dimensional large eddy simulation （LES） model is developed to simulate th...This paper presents a numerical study on the hydrodynamic behaviours of a round buoyant jet under the effect of JONSWAP random waves. A three-dimensional large eddy simulation （LES） model is developed to simulate the buoyant jet in a stagnant ambient and JONSWAP random waves. By comparison of velocity and concentration fields, it is found that the buoyant jet exhibits faster decay of centerline velocity, wider lateral spreading and larger initial dilution under the wave effect, indicating that wave dynamics improves the jet entrainment and mixing in the near field, and subsequently mitigate the jet impacts in the far field. The effect of buoyancy force on the jet behaviours in the random waves is also numerically investigated. The results show that the wave effect on the jet entrainment and mixing is considerably weakened under the existence of buoyancy force, resulting in a slower decay rate of centerline velocity and a narrower jet width for the jet with initial buoyancy.展开更多
An algebraic modelling o?large eddy simulation (LES) has been proposed and discussed in detail. Its basic idea and application to numerical simulation of turbulence in a straight channel and a curved channel are prese...An algebraic modelling o?large eddy simulation (LES) has been proposed and discussed in detail. Its basic idea and application to numerical simulation of turbulence in a straight channel and a curved channel are presented. Primary computational results show the feasibility of this modelling. It is our hope that this modelling) if further improved, will be useful for the numerical simulation of more complicated turbulence.展开更多
The turbulent large eddy simulation (LES) technique and the finite element method (FEM) of computational fluid dynamics (CFD) are used to predict the three-dimensional flow field in a vector flow clean-room under em...The turbulent large eddy simulation (LES) technique and the finite element method (FEM) of computational fluid dynamics (CFD) are used to predict the three-dimensional flow field in a vector flow clean-room under empty state and static state conditions. The partly expanded Taylor-Galerkin (TG) discretization scheme is combined with implicit stream-upwind diffusion in the finite element formulation of the basic equations with Gauss filtering. The vortex viscosity subgrid model is used in the numerical simulation. The numerical results agree well with the available experimental data, showing that the LES method can more accurately predict the size and location of large eddies in clean-rooms than the standard k-ε two equation model.展开更多
Flame stabilization is the key to extending scramjets to hypersonic speeds;accordingly,this topic has attracted much attention in theoretical research and engineering design.This study performed large eddy simulations...Flame stabilization is the key to extending scramjets to hypersonic speeds;accordingly,this topic has attracted much attention in theoretical research and engineering design.This study performed large eddy simulations(LESs)of lifted hydrogen jet combustion in a stepped-wall combustor,focusing on the flame stabilization mechanisms,especially for the autoignition effect.An assumed probability density function(PDF)approach was used to close the subgrid chemical reaction source.The reliability of the solver was confirmed by comparing the LES results with experimental data and published simulated results.The hydrogen jet and the incoming stream were first mixed by entraining large-scale vortices in the shear layer,and stable combustion in the near-wall region was achieved downstream of the flame induction region.The autoignition cascade is a transition of fuel-rich flame to stoichiometric ratio flame that plays a role in forming the flame base,which subsequently causes downstream flame stabilization.Three cases with different jet total temperatures are compared,and the results show that the increase in the total temperature reduces the lift-off distance of the flame.In the highest total temperature case,an excessively large scalar dissipation rate inhibits the autoignition cascade,resulting in a fuel-rich low-temperature flame.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.91441117 and 51576182)
文摘A sub-grid scale(SGS) combustion model, which combines the artificial thickened flame(ATF) model with the flamelet generated manifold(FGM) tabulation method, is proposed. Based on the analysis of laminar flame structures, two self-contained flame sensors are used to track the diffusion and reaction processes with different spatial scales in the flame front, respectively. The dynamic formulation for the proposed SGS combustion model is also performed. Large eddy simulations(LESs) of Bunsen flame F3 are used to evaluate the different SGS combustion models. The results show that the proposed SGS model has the ability in predicting the distributions of temperature and velocity reasonably, while the predictions for the distributions of some species need further improvement. The snapshots of instantaneous normalized progress variables reveal that the flame is more remarkably and severely wrinkled at the flame tip for flame F3.More satisfactory results obtained by the dynamic model indicate that it can preserve the premixed flame propagation characteristics better.
基金Project supported by the National Natural Science Foundation of China(Nos.91441117 and51576182)
文摘Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.
基金Project supported by the National Natural Science Foundation of China (No.10502029)the Scientific Research Foundation for Returned Overseas Chinese Scholars of Ministry of Education of China
文摘By comparing the energy spectrum and total kinetic energy, the effects of numerical errors （which arise from aliasing and discretization errors）, subgrid-scale （SGS） models, and their interactions on direct numerical simulation （DNS） and large eddy simulation （LES） are investigated, The decaying isotropic turbulence is chosen as the test case. To simulate complex geometries, both the spectral method and Pade compact difference schemes are studied. The truncated Navier-Stokes （TNS） equation model with Pade discrete filter is adopted as the SGS model. It is found that the discretization error plays a key role in DNS. Low order difference schemes may be unsuitable. However, for LES, it is found that the SGS model can represent the effect of small scales to large scales and dump the numerical errors. Therefore, reasonable results can also be obtained with a low order discretization scheme.
基金the National Natural Science Foundation of China(Nos.50979059 and 11272213)
文摘Moving-particle semi-implicit(MPS) method is a new mesh-free numerical method based on Lagrangian particle. In this paper, MPS method is applied to the study on numerical wave tank. For the purpose of simulating numerical wave, we combine the MPS method with large eddy simulation(LES) which can simulate the turbulence in the flow. The intense pressure fluctuation is a significant shortcoming in MPS method. So, we improve the original MPS method by using a new pressure Poisson equation to ease the pressure fluctuation. Divergencefree condition representing fluid incompressible is used to calculate pressure smoothly. Then, area-time average technique is used to deal with the calculation. With these improvements, the modified MPS-LES method is applied to the simulation of numerical wave. As a contrast, we also use the original MPS-LES method to simulate the wave in a numerical wave tank. The result shows that the new method is better than the original MPS-LES method.
基金supported by the National Natural Science Foundation of China(Grant No.51379128)
文摘Large eddy simulation (LES) explicitly calculates the large-scale vortex field and parameterizes the small-scale vortices.In this study,LES and κ-ε models were developed for a specific geometrical configuration of backward-facing step (BFS).The simulation results were validated with particle image velocimetry (PIV) measurements and direct numerical simulation (DNS).This LES simulation was carried out with a Reynolds number of 9000 in a pressurized water tunnel with an expansion ratio of 2.00.The results indicate that the LES model can reveal largescale vortex motion although with a larger grid-cell size.However,the LES model tends to overestimate the top wall separation and the Reynolds stress components for the BFS flow simulation without a sufficiently fine grid.Overall,LES is a potential tool for simulating separated flow controlled by large-scale vortices.
基金co-supported by the National Key R&D Program of China(Nos.2017YFB0202400 and 2017YFB0202402)the National Natural Science Foundation of China(No.91741125)the Project of Newton International Fellowship Alumnus from Royal Society(No.AL120003)
文摘The numerical simulation of modern aero-engine combustion chamber needs accurate description of the interaction between turbulence and chemical reaction mechanism. The Large Eddy Simulation(LES) method with the Transported Probability Density Function(TPDF) turbulence combustion model is promising in engineering applications. In flame region, the impact of chemical reaction should be considered in TPDF molecular mixing model. Based on pioneer research, three new TPDF turbulence-chemistry dual time scale molecular mixing models were proposed tentatively by adding the chemistry time scale in molecular mixing model for nonpremixed flame. The Aero-Engine Combustor Simulation Code(AECSC) which is based on LES-TPDF method was combined with the three new models. Then the Sandia laboratory's methane-air jet flames: Flame D and Flame E were simulated. Transient simulation results show that all the three new models can predict the instantaneous combustion flow pattern of the jet flames. Furthermore,the average scalar statistical results were compared with the experimental data. The simulation result of the new TPDF arithmetic mean modification model is the closest to the experimental data:the average error in Flame D is 7.6% and 6.6% in Flame E. The extinction and re-ignition phenomena of the jet flames especially Flame E were captured. The turbulence time scale and the chemistry time scale are in different order in the whole flow field. The dual time scale TPDF combustion model has ability to deal with both the turbulence effect and the chemistry reaction effect, as well as their interaction more accurately for nonpremixed flames.
文摘Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range. By decomposing the subgrid energy transfer and nonlinear interaction into ‘forward’ and ‘backward’ groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.
基金supported by the National Natural Science Foundation of China(Nos.91441202 and 51476087)
文摘Large Eddy Simulations(LES) in conjunction with the Flamelet Progress Variable(FPV) approach have been performed to investigate the flame and large-scale flow structures in the bluff-body stabilized non-premixed flames, HM1 and HM3. The validity of the numerical methods is first verified by comparing the predicted velocity and composition fields with experimental measurements. Then the evolution of the flame and large-scale flow structures is analyzed when the flames approach blow-off. The analysis of instantaneous and statistical data indicates that there exists a shift of the control mechanism in the recirculation zone in the two flames. In the recirculation zone, HM1 flame is mainly controlled by the mixing effect and ignition mainly occurs in the outer shear layer. In HM3 flame, both the chemical reactions and mixing are important in the recirculation zone. The Proper Orthogonal Decomposition(POD) results show that the fluctuations in the outer shear layer are more intense in HM1, while the flow structures are more obvious in the outer vortex structure in HM3, due to the different control mechanism in the recirculation zone.It further shows that the flow structures in HM1 spread larger in the intense mixing zone due to higher temperature and less extinction.
基金jointly supported by research grants from the Hong Kong Polytechnic Universitythe Hong Kong Research Grants Council.
文摘A non-linear turbulent hydrodynamic model was proposed by Zhan and Li (1993) using the large eddy simulation approach and directional eddy viscosity method. The model is simplified in this paper and is applied to the problem of a side discharge into open-channel flow,where a recirculation zone develops downstream of the discharge. The numerical results are compared with experimental results (Strazisar and Prahl (1973) and Mikhail, Chu, and Savage (1975)) and previous numerical results obtained using a k-model (McGuirk and Rodi, 1978). The agreement is good.
基金This paper was supported by the National Natural Science Foundation of China (Grant No.19972061)
文摘By the Volume of Fluid (VOF) multiphase flow model two-dimensional gravity currents with three phases including air are numerically simulated in this article. The necessity of consideration of turbulence effect for high Reynolds numbers is demonstrated quantitatively by LES (the Large Eddy Simulation) turbulence model. The gravity currents are simulated for h not equal H as well as h = H, where h is the depth of the gravity current before the release and H is the depth of the intruded fluid. Uprising of swell occurs when a current flows horizontally into another lighter one for h not equal H. The problems under what condition the uprising of swell occurs and how long it takes are considered in this article. All the simulated results are in reasonable agreement with the experimental results available.
基金supported by National Science and Technology Major Project(J2019-III-0014-0058)Natural Science Foundation of Science and Technology Department of Shaanxi Province(2022JQ-712)Scientific Research Program of Shaanxi Provincial Education Department(21JK0642)。
文摘Bunsen burner is a typical geometry for investigating the turbulence-flame interaction.In most experimental studies,only turbulence intensity u′and integral scale l0 are used to characterize the turbulent flow field,regardless of the perforation geometry of perforated plates.However,since the geometry influences the developing process and vortex broken,the plate geometry has to be considered when discussing the flame-turbulence interaction.In order to investigate conditions at the same l0 and u′using different geometries,large eddy simulation of CH_(4)/air flames with dynamic TF combustion model was performed.The model validation shows good agreement between Large Eddy Simulation(LES)and experimental results.In the non-reacting flows,the Vortex Stretching of circular-perforated plate condition is always larger than that of slot-perforated plate condition,which comes from the stresses in the flow fields to stretch the vorticity vector.In reacting flows,at the root of the flame,the Vortex Stretching plays a major role,and the total vorticity here of circular-perforated plate condition is still larger(53.8%and 300%larger than that of the slot-perforated plate at x/D=0 and x/D=2.5,respectively).More small-scale vortex in circular-perforated plate condition can affect and wrinkle the flame front to increase the Probability Density Function(PDF)at large curvatures.The 3D curvature distributions of both cases bias to negative values.The negative trend of curvatures at the instant flame front results from the Dilatation term.Also,the value of the Vortex Stretching and the Dilatation at the flame front of circular-perforated plate condition is obviously larger.
文摘A computational model of three-dimensional, time-dependent flame spread in microgravity environment is presented. The solid is assumed to be a thermally-thin, pyrolysing cellulosic sheet. The gas phase model includes the full Navier-Stokes equations with density and pressure variations and six-flus model of radiation heat transfer. The solid phase model consists of continuity and energy eqllations whose solution provides boundary conditions for the gas phase equations. In the numerical procedure, the gas-and solid-phase equations are solved sepaxately and iteratively at each time step. Predictions have been made of flame spread in slow forced flow under gravitational acceleration normal to fuel surface and flame spread in a quiescent environment in an enclosed chamber under gravitational acceleration parallel to fuel surface. Numerical simulations show that, under microgravity, slow-flow conditions, flame spread process is highly unsteady with the upstream flame spreads faster than the downstream flame after a period of ignition. It has also been shown that the level of microgravity has a significant effect on the name spread process.
基金This work was supported by the National Key Research&Development(R&D)Plan of China under(No.2020YFC 1522800)the National Natural Science Foundation of China(NSFC)(No.51876148)the Science and Technology Project of State Grid Anhui Electric Corporation of China(No.52120518001S).
文摘At present,there is a shortage of experimental and simulation studies on fire spread in medium-and large-scale compartments while the existing models of the fire spread are limited for typical engineering applications.This paper proposes a new model for large-scale fire spread on medium density fibreboard(MDF)panels.Validating the model with single burning item(SBI)experiments,it is found that the numerical simulation closely predicts the experimental heat release rate(HRR)with some error near the peak.The predicted heat flux and distance of lateral flame spread are consistent with the experiments and an existing model.The effects of kinetic properties and heat of combustion are identified through a sensitivity analysis.The decrease of activation energy and increase of pre-exponential factor make the MDF easier to pyrolyze and the increase of heat of combustion enhances the flame temperature and thus provide more heat feedback to the sample surface.The low activation energy(71.9 kJ/mol)and high heat of combustion(46.5 MJ/kg)of the model ensure the occurrence of flame spread.Furthermore,the model was validated using medium-scale compartment fire experiments and the results showed that the model can accurately predict the HRR after flashover(the error is within 7%).While the burner is ignited,the predictions of in-compartment gas temperature and heat flux are more accurate.However,when the burner is extinguished,the modelled in-compartment gas temperature is lower than the experimental values,resulting in a lower heat flux prediction.The model leads to easier flame spread;therefore,the modelled flame spreads faster in the compartment compared to the experiment,and thus the HRR increases more rapidly.
基金supported by the National Natural Science Foundation of China(51906075)the National Key Research and Development Program of China(2019YFE0100100).
文摘To improve the NO modelling in turbulent flames,the flamelet/progress variable(FPV)model is extended by introducing NO mass fraction into the progress variable and incorporating an additional NO transport equation.Two sets of flamelet databases are tabulated with progress variables based on major species and NO mass fraction,respectively.The former is used for the acquisition of the main thermochemical variables,while the latter is employed for NO modelling.Moreover,an additional transport equation is solved to obtain the NO mass fraction,with the source term corrected using the scale similarity method.Model assessments are first conducted on laminar counterflow diffusion flames to identify lookup-related errors and assess the suitability of progress variable definitions.The results show that the progress variables based on major species and NO could correctly describe the main thermochemical quantities and NO-related variables,respectively.Subsequently,the model is applied to the large eddy simulation(LES)of Sandia flames.The results indicate that the extended FPV model improves the NO prediction,with a mean error for NO prediction at 55%,significantly lower than those of existing FPV models(130%and 385%).The LES with the extended FPV model quantitatively captures NO suppression in the mid-range of Reynolds numbers from 22400(Flame D)to 33600(Flame E),but underestimates the NO suppression at higher Reynolds numbers from 33600 to 44800(Flame F).This underprediction is primarily attributed to the underestimation of local extinction levels in flames with high Reynolds numbers.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 50579025 and 90210005 ).
文摘Turbulent flow in a 3-D blade passage of a Francis hydro turbine was simulated with the Large Eddy Simulation （LES） to investigate the spatial and temporal distributions of the turbulence when strongly distorted wakes in the inflow sweep over the passage, In a suitable consideration of the energy exchanging mechanism between the large and small scales in the complicated passage with a strong 3-D curvature, one-coefficient dynamic Sub-Grid-Scale （SGS） stress model was used in this article. The simulations show that the strong wakes in the inflow lead to a flow separation at the leading zone of the passage, and to form a primary vortex in the span-wise direction. The primary span-wise vortex evolves and splits into smaller vortex pairs due to the constraint of no-slip wall condition, which triggers losing stability of the flow in the passage. The computed pressures on the pressure and suction sides agree with the measured data for a working test turbine model.
基金supported by the National Key Basic Research Program of the Ministry of Science and Technology of China(Grant No.2010CB429001)the Special Fund of State Key Laboratory of China(Grant No.2011585812)+2 种基金the Fundamental Research Funds for the Central Universities(Grant No.2011B05614)the 111 Project of the Ministry of Educationthe State Administration of Foreign Experts Affairs,China(Grant No.B12032)
文摘This paper presents a numerical study on the hydrodynamic behaviours of a round buoyant jet under the effect of JONSWAP random waves. A three-dimensional large eddy simulation （LES） model is developed to simulate the buoyant jet in a stagnant ambient and JONSWAP random waves. By comparison of velocity and concentration fields, it is found that the buoyant jet exhibits faster decay of centerline velocity, wider lateral spreading and larger initial dilution under the wave effect, indicating that wave dynamics improves the jet entrainment and mixing in the near field, and subsequently mitigate the jet impacts in the far field. The effect of buoyancy force on the jet behaviours in the random waves is also numerically investigated. The results show that the wave effect on the jet entrainment and mixing is considerably weakened under the existence of buoyancy force, resulting in a slower decay rate of centerline velocity and a narrower jet width for the jet with initial buoyancy.
基金Project supported by the National Natural Science Foundation of China and Scientific Fund of Tsinghua University
文摘An algebraic modelling o?large eddy simulation (LES) has been proposed and discussed in detail. Its basic idea and application to numerical simulation of turbulence in a straight channel and a curved channel are presented. Primary computational results show the feasibility of this modelling. It is our hope that this modelling) if further improved, will be useful for the numerical simulation of more complicated turbulence.
文摘The turbulent large eddy simulation (LES) technique and the finite element method (FEM) of computational fluid dynamics (CFD) are used to predict the three-dimensional flow field in a vector flow clean-room under empty state and static state conditions. The partly expanded Taylor-Galerkin (TG) discretization scheme is combined with implicit stream-upwind diffusion in the finite element formulation of the basic equations with Gauss filtering. The vortex viscosity subgrid model is used in the numerical simulation. The numerical results agree well with the available experimental data, showing that the LES method can more accurately predict the size and location of large eddies in clean-rooms than the standard k-ε two equation model.
基金National Natural Science Foundation of China(Nos.91741205 and 11522222)。
文摘Flame stabilization is the key to extending scramjets to hypersonic speeds;accordingly,this topic has attracted much attention in theoretical research and engineering design.This study performed large eddy simulations(LESs)of lifted hydrogen jet combustion in a stepped-wall combustor,focusing on the flame stabilization mechanisms,especially for the autoignition effect.An assumed probability density function(PDF)approach was used to close the subgrid chemical reaction source.The reliability of the solver was confirmed by comparing the LES results with experimental data and published simulated results.The hydrogen jet and the incoming stream were first mixed by entraining large-scale vortices in the shear layer,and stable combustion in the near-wall region was achieved downstream of the flame induction region.The autoignition cascade is a transition of fuel-rich flame to stoichiometric ratio flame that plays a role in forming the flame base,which subsequently causes downstream flame stabilization.Three cases with different jet total temperatures are compared,and the results show that the increase in the total temperature reduces the lift-off distance of the flame.In the highest total temperature case,an excessively large scalar dissipation rate inhibits the autoignition cascade,resulting in a fuel-rich low-temperature flame.
