Comparison of a Full Second-Order Moment Model and an Algebraic Stress Two-Phase Turbulence Model for Simulating Bubble-Liquid Flows in a Bubble Column

A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.

EFFECT OF EMPIRICAL COEFFICIENTS ON SIMULATION IN TWO-SCALE SECOND-ORDER MOMENT PARTICLE-PHASE TURBULENCE MODEL

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. The proposed model is used to simulate gas-particle downer reactor flows. The computational results of both particle volume fraction and mean velocity are in agreement with the experimental results. After analyzing effects of empirical coefficient on prediction results, we can come to a conclusion that, inside the limit range of empirical coefficient, the predictions do not reveal a large sensitivity to the empirical coefficient in the downer reactor, but a relatively great change of the constants has important effect on the prediction.

Numerical simulation of the effect of void fraction and inlet velocity on two-phase turbulence in bubble-liquid flows

There are contradicted opinions on whether bubbles enhance or reduce the liquid turbulence. In this paper, the effect of void fraction and inlet velocity on the bubble-liquid two-phase turbulence of the multiple bubble-liquid jets in a two-dimensional channel is studied by using the two-phase second-order moment turbulence model. The results confirm the phenomena observed in experiments and reported in references that at a low void fraction and low inlet velocities the bubbles enhance the liquid turbulence, whereas at a high void fraction and high inlet velocities the bubbles reduce the liquid turbulence.

NUMERICAL SIMULATION OF METHANE-AIR TURBULENT JET FLAME USING A NEW SECOND-ORDER MOMENT MODEL

A new second-order moment model for turbulent combustion is applied in the simulation of methane-air turbulent jet flame. The predicted results are compared with the experimental results and with those predicted using the well-known EBU-Arrhenius model and the original second-order moment model. The comparison shows the advantage of the new model that it requires almost the same computational storage and time as that of the original second-order moment model, but its modeling results are in better agreement with experiments than those using other models. Hence, the new second-order moment model is promising in modeling turbulent combustion with NOx formation with finite reaction rate for engineering application.

Second-order moment modeling of dispersed two-phase turbulence——Part 1:USM,k-ε-k_p,and non-linear k-ε-k_p two-phase turbulence models

Turbulent dispersed multiphase flows,including gas-particle,gas-droplet and bubble-liquid flows,are widely encountered in various engineering facilities.Modeling of two-phase turbulence,in particular the dispersed phase turbulence,is the key problem in the Eulerian-Eulerian simulation of practical dispersed multiphase flows.Although different models were developed and used,the experimental validation shows that they cannot always give satisfactory prediction results.In this paper the present author give a detailed review of the unified second-order moment (USM),k-k p and nonlinear k-k p two-phase turbulence models,proposed by him.The derivation and closure of these models are described in detail and the experimental validation and application of these models are extensively discussed.

Second-order moment modeling of dispersed two-phase turbulence-Part 2-USM-Θ two-phase turbulence model and USM-SGS two-phase stress model

Dense gas-particle flows are frequently encountered in fluidized beds,riser and downer reactors,pneumatic transport and the near-wall zone of dilute gas-particle flows.Particle-particle collision plays an important role in the behavior of two-phase flows.In this paper a USM-Q two-phase turbulence model for dense gas-particle flows is proposed to account for both two-phase turbulence and inter-particle collision.For two-fluid large-eddy simulation of gas-particle flows,the author proposed a unified second-order moment(USM) two-phase SGS stress model and a two-phase k-kp SGS energy-equation stress model.The proposed models can fully account for the interaction between the gas and particle SGS stresses.

SECOND-ORDER MOMENT MODEL FOR DENSE TWO-PHASE TURBULENT FLOW OF BINGHAM FLUID WITH PARTICLES

The USM-θ model of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θ model in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θ model has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θ model was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale （SGS） stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the cur- rent second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body- fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbu-lent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results.

TURBULENCE IN BROKEN LEE WAVES

In this paper, the waves' breaking in the lee waves is successfully simulated by the atmospheric mesoscale numerical model with a second-order tur- bulent closure. It is further proved that the turbulence in the wave-breaking region plays the role of intense mixing for the average field, which leads to the trapping of upward propagating waves and thus promotes the development of the downslope wind. The turbulent structure in the wave-breaking region is discussed and the fol- lowing conclusions are obtained: (1) In the wave-breaking region, the turbulent heat fluxes transfer from inside to outside and the turbulent momentum fluxes transfer from outside to inside. (2) In the wave-breaking region, the turbulent energy mainly comes from the wind shear and the buoyancy promotes the turbulent development only in part of the region. (3) In the upper part of the wave-breaking region, the turbulent momentum fluxes behave as a counter-gradient transfer. (4) The turbulent mixing in the wave-breaking region is non-local.

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.

Effect of turbulence on NO formation in swirling combustion

Turbulence affects both combustion and NO formation. Fluctuation correlations are ideally used for quantitative analysis. From the instantaneous chemical reaction rate expression,ignoring the third-order correlation terms, the averaged reaction rate will have four terms, including the term of averaged-variable product, a concentration fluctuation correlation term, and temperature-concentration fluctuation correlation term. If the reaction-rate coefficient is denoted as K, the temperature fluctuation would be included in the K fluctuation. In order to quantitatively study the effect of turbulence on NO formation in methane-air swirling combustion, various turbulencechemistry models are tested. The magnitudes of various correlations and their effects on the time-averaged reaction rate are calculated and analyzed, and the simulation results are compared with the experimental measurement data. The results show that among various correlation moments, the correlation between the reaction-rate coefficient K fluctuation with the concentration fluctuation is most important and is a strong nonlinear term.

Modeling Evaporating Droplets in Complex Unsteady Flows

In many applications, a moving fluid carries a suspension of droplets of a second phase which may change in size due to evaporation or condensation. Examples include liquid fuel drops in engines and raindrops or ice-crystals in a thunderstorm. If the number of such particles is very large, and, if further, the flow is inhomogeneous, unsteady or turbulent, it may be practically impossible to explicitly compute all of the fluid and particle degrees of freedom in a numerical simulation of the system. Under such circumstances Lagrangian Particle Tracking (LPT) of a small subset of the particles is used to reduce the computational effort. The purpose of this paper is to compare the LPT with an alternate method that is based on an approximate solution of the conservation equation of particle density in phase space by the method of moments (MOM). Closure is achieved by invoking the assumption that the droplet size distribution is locally lognormal. The resulting coupled transport equations for the local mean and variance of the particle size distribution are then solved in conjunction with the usual equations for the fluid and associated scalar fields. The formalism is applied to the test case of a uniform distribution of droplets placed in a non homogeneous temperature field and stirred with a decaying Taylor vortex. As a benchmark, we perform a direct numerical simulation (DNS) of high resolution that keeps track of all the particles together with the fluid flow.

Development of SOM combustion model for Reynolds-averaged and large-eddy simulation of turbulent combustion and its validation by DNS

The derivation and closure methods of the second-order moment (SOM) combus- tion model are proposed. The application of this model to Reynolds averaged (RANS) and large-eddy simulation (LES) of turbulent swirling diffusion combustion, jet diffusion combustion, and bluff-body stabilized premixed combustion is sum- marized. It is indicated that the SOM model is much better than the eddy-beak-up (EBU) and presumed PDF models widely used in commercial software and engi- neering. The SOM modeling results are close to those obtained using the most accurate but much more complex PDF equation model. Moreover, it can save much more computation time than the PDF equation model. Finally, the SOM model is validated by the direct numerical simulation (DNS) of turbulent reacting channel flows.