A turbulence model for inclined, bubbly flows

A new turbulence model for the flow of a two phase (liquid-liquid) flow in an inclined pipe is presented. An eddy viscosity is used to model the effects of shear induced turbulence and bubble induced turbulence. The cross-pipe momentum transport arising from the buoyant rise of bubbles across the axial flow is also modelled. Numerical simulations have been carried out in both one and two dimensions. One and two dimensional numerical simulations are presented.On leave from the University of Leeds, Leeds LS2 9JT, U.K.     

超声速流动中非线性EASM湍流模式应用研究

针对超声速复杂流动区域精确模拟的需要,发展了基于k-ω可压缩修正形式的非线性显式代数雷诺应力模式（EASM）,提高了该模式对超声速复杂流动的数值模拟精度。通过对二维超声速凹槽和三维双椭球的数值计算表明,与SA和SST常规线性涡黏性湍流模式比较,非线性的EASM模式对大分离以及剪切层流动结构的刻画能力更精细,对剪切层再附区的压力及摩擦系数分布模拟更加精确;EASM模式能够准确地模拟二次激波引起的压强和热流分布情况。     

Comparative assessment of SAS and DES turbulence modeling for massively separated flows

Numerical studies of the flow past a circular cylinder at Reynolds number 1.4 × 10~5 and NACA0021airfoil at the angle of attack 60?have been carried out by scale-adaptive simulation (SAS) and detached eddy simulation (DES), in comparison with the existing experimental data. The new version of the model developed by Egorov and Menter is assessed, and advantages and disadvantages of the SAS simulation are analyzed in detail to provide guidance for industrial application in the future. Moreover, the mechanism of the scale-adaptive characteristics in separated regions is discussed, which is obscure in previous analyses. It is concluded that: the mean flow properties satisfactorily agree with the experimental results for the SAS simulation, although the prediction of the second order turbulent statistics in the near wake region is just reasonable. The SAS model can produce a larger magnitude of the turbulent kinetic energy in the recirculation bubble, and, consequently, a smaller recirculation region and a more rapid recovery of the mean velocity outside the recirculation region than the DES approach with the same grid resolution. The vortex shedding is slightly less irregular with the SAS model than with the DES approach,probably due to the higher dissipation of the SAS simulation under the condition of the coarse mesh.     

A general one-equation turbulence model for free shear and wall-bounded flows

The purpose of thiswork is to introduce a complete and general one-equation model capable of correctly predicting a wide class of fundamental turbulent flows like boundary layer, wake, jet, and vortical flows. The starting point is the mature and validated two-equation k−ω turbulence model of Wilcox. The newly derived one-equation model has several advantages and yields better predictions than the Spalart-Allmaras model for jet and vortical flows while retaining the same efficiency and quality of the results for near-wall turbulent flows without using a wall distance. The derivation and validation of the new model using findings computed by the Spalart-Allmaras and the k−ω models are presented and discussed for several free shear and wall-bounded flows.     

A two-equation turbulence model for jet flows laden with vaporizing droplets

A two-equation turbulence model for steady incompressible two-phase flows including phase change has been recently developed by Mostafa & Elghobashi (1984). This model is tested for the flow of a turbulent axisymmetric gaseous jet laden with evaporating liquid droplets. To avoid the problem of density fluctuations of the carrier phase at this stage, only isothermal flow is considered and vaporization is assumed to be due to the vapor concentration gradient. The continuous size distribution of the droplets is approximated by finite size groups. Each group is considered as a continuous phase interpenetrating and interacting with the carrier phase. Two test cases have been predicted by the model. The first is for a Freon-11 spray issuing from a round nozzle, where experimental data are available at distances equal to or greater than 170 nozzle diameters. Good agreement between the data and the predictions was achieved. The second is for a methanol spray where no experiments are available yet and the predictions consider the flow region close to the nozzle ($\text{z/D < 40}$). The results of the methanol spray include distributions of the mean velocity, volume fractions of the different phases, concentration of the evaporated material in the carrier phase, turbulence intensity and shear stress of the carrier phase, droplet diameter distribution, and the jet spreading rate. In this case the results are analyzed based on a qualitative comparison with the corresponding single phase jet flow.     

A USM-Θ two-phase turbulence model for simulating dense gas-particle flows

A second-order moment two-phase turbulence model for simulating dense gas-particle flows (USM- model), combining the unified second-order moment two-phase turbulence model for dilute gas-particle flows with the kinetic theory of particle collision, is proposed. The interaction between gas and particle turbulence is simulated using the transport equation of two-phase velocity correlation with a two-time-scale dissipation closure. The proposed model is applied to simulate dense gas-particle flows in a horizontal channel and a downer. Simulation results and their comparison with experimental results show that the model accounting for both anisotropic particle turbulence and particle-particle collision is obviously better than models accounting for only particle turbulence or only particle-particle collision. The USM- model is also better than the k--kp- model and the k--kp-p- model in that the first model can simulate the redistribution of anisotropic particle Reynolds stress components due to inter-particle collision, whereas the second and third models cannot.The project supported by the Special Funds for Major State Basic Research of China (G-1999-0222-08), the National Natural Science Foundation of China (50376004), and Ph.D. Program Foundation, Ministry of Education of China (20030007028)     

A filtered mass density function approach for modeling separated two-phase flows for LES I: Mathematical formulation

The overall objective of this study is to develop a full velocity-scalar filtered mass density function (FMDF) formulation for large eddy simulation (LES) of a separated two-phase flow. Required in the development of the two-phase FMDF transport equation are the local instantaneous equations of motion for a two-phase flow previously derived by Kataoka. In Kataoka’s development, phase interaction terms are cast in terms of a Dirac delta distribution on the phase interface. For this reason, it is difficult to close these coupling terms in the instantaneous formulation and this difficulty is propagated into the phase-coupling terms in the FMDF transport equation. To address this point a new derivation of the local instantaneous equations for a separated two-phase flow is given. The equations are shown to be consistent with the formulation given by Kataoka, and in the development, a direct link between the conditionally surface-filtered coupling terms, arising in the FMDF formulation, and LES phase-coupling terms is established. Clarification of conditions under which conditionally filtered interphase conversion terms in the marginal FMDF transport equations may be disregarded in a separated continuum-dispersed phase flow is discussed. Modeling approaches and solutions procedures to solve the two-phase FMDF transport equation via Monte-Carlo methods are outlined.     

A two-scale second-order moment two-phase turbulence model for simulating dense gas-particle flows

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the single-scale second-order moment two-phase turbulence model superposed with a particle collision model (USM-θ model) in most regions. The project supported by the Special Funds for Major State Basic Research, China (G-1999-0222-08), and the Postdoctoral Science Foundation (2004036239) The English text was polished by Keren Wang     

Thermal contact conductance for cylindrical and spherical contacts

A prediction methodology based on Monte-Carlo simulation model, developed for flat conforming surfaces in contact, is modified and extended to predict contact conductance between curvilinear surfaces like cylinders and spheres. Experiments are also conducted in vacuum for the measurement of contact conductance between stainless steel and aluminium cylindrical contacts and stainless steel spherical contacts over a range of contact pressures. The contact conductance between cylindrical and spherical bodies is, in general, about an order of magnitude lower than for flat surfaces in contact. Increase of surface roughness and decrease in contact pressure lowers the contact conductance. However, the influence of these parameters is larger than those obtained for flat surfaces. The prediction for different parametric conditions agree closely with those measured in the experiments.     

A stability condition for turbulence model: From EMMS model to EMMS-based turbulence model

The closure problem of turbulence is still a challenging issue in turbulence modeling. In this work, a stability condition is used to close turbulence. Specifically, we regard single-phase flow as a mixture of turbulent and non-turbulent fluids, separating the structure of turbulence. Subsequently, according to the picture of the turbulent eddy cascade, the energy contained in turbulent flow is decomposed into different parts and then quantified. A turbulence stability condition, similar to the principle of the energy-minimization multi-scale (EMMS) model for gas–solid systems, is formulated to close the dynamic constraint equations of turbulence, allowing the inhomogeneous structural parameters of turbulence to be optimized. We name this model as the “EMMS-based turbulence model”, and use it to construct the corresponding turbulent viscosity coefficient. To validate the EMMS-based turbulence model, it is used to simulate two classical benchmark problems, lid-driven cavity flow and turbulent flow with forced convection in an empty room. The numerical results show that the EMMS-based turbulence model improves the accuracy of turbulence modeling due to it considers the principle of compromise in competition between viscosity and inertia.     

A novel two-indenter based micro-pump for lab-on-a-chip application: modeling and characterizing flows for a non-Newtonian fluid

Meccanica - Inspired by the feeding mechanisms of a nematode, a novel two-indenter (2I) micro-pump is analyzed theoretically for transport and mixing of a non-Newtonian fluid for the purpose of...     

Two-phase turbulence models for simulating dense gas-particle flows

The two-fluid model is widely adopted in simulations of dense gas-particle flows in engineering facili- ties. Present two-phase turbulence models for two-fluid modeling are isotropic. However, turbulence in actual gas-particle flows is not isotropic. Moreover, in these models the two-phase velocity correlation is closed using dimensional analysis, leading to discrepancies between the numerical results, theoretical analysis and experiments. To rectify this problem, some two-phase turbulence models were proposed by the authors and are applied to simulate dense gas-particle flows in downers, risers, and horizontal channels; Experimental results validate the simulation results. Among these models the USM-O and the two-scale USM models are shown to give a better account of both anisotropic particle turbulence and particle-particle collision using the transport equation model for the two-phase velocity correlation.     

Current capabilities for modelling turbulence in industrial flows

This paper summarizes some of the more widely applied or promising schemes for computing heat and momentum transport in industrially relevant flows. Such flows typically involve complex flow domains, severe pressure gradients and regions of flow separation and reattachment. Models tuned by reference to equilibrium, simple shear flows cannot in general be relied upon to predict accurately the effects of these complexities on the transport processes. The main conclusions drawn are that second-moment closure offers a far more reliable basis for computing non-equilibrium turbulent flows than eddy-viscosity schemes, especially in flows with very complex strain fields or those substantially affected by external force fields. Moreover, where significant variations in shear stress occur across the near-wall viscosity-affected sublayer, the usual practice of employing wall functions needs to be replaced by at least a two-equation model in order to capture, even qualitatively, the consequent effects on wall heat transfer or skin friction coefficient.     

Parametric data-based turbulence modelling for vortex dominated flows

Computational fluid dynamics simulations employing eddy-viscosity turbulence models remain the baseline numerical tool in the aerospace industry, mainly due to their numerical stability and computational efficiency. However, many industrially relevant cases require a level of accuracy that is not routinely achieved by global turbulence models. The simulation of leading-edge vortices shed at low aspect ratio wings is one such class of flows that remains a challenge for turbulence modelling. A local approach is proposed in which a parametrised eddy-viscosity turbulence model is calibrated using experimental results of configurations and flow conditions similar to the one being analysed. In this paper, the Spalart–Allmaras one-equation model is enhanced with additional source terms, which are exclusively active in the vortex field. An automatic optimisation procedure with experimental data as reference is then applied. The resulting optimised model improves the eddy viscosity distribution for a limited but relevant range of configurations and flow conditions.     

Two-dimensional viscoelastic flows: experimentation and modeling

Extensive experimental data on the birefringence in converging and diverging flows of a polymeric melt have been obtained. The birefringence and pressure drop measurements were carried out in working cells of planar geometry having different contraction angles and contraction ratios. For investigation of diverging or abrupt expansion flow, the direction of flow in the cells was reversed. The theoretical predictions are based upon the Leonov constitutive equation and a finite element scheme with streamwise integration.In contrast to Newtonian and second-order fluids, viscoelastic fluids at high shear rates show significant differences in pressure drop and birefringence (i.e. stresses) in converging and diverging flows. For a constant flow rate, the pressure drop is higher and the birefringence smaller in diverging flows than in converging flows. This difference increases with increasing flow rate. Further, for the same contraction ratio but different contraction angles, the birefringence maximum increases considerably with contraction angle. In addition, an increase in contraction ratio has the same effect.The viscoelastic constitutive equation of Leonov has been shown to describe all the above viscoelastic effects observed in the experiments. In general, a reasonable agreement between theory and experiment has been obtained, which shows the usefulness of the Leonov model in describing actual flows.     

Particle motion and turbulence in dense two-phase flows

Measurements of particle mean and r.m.s. velocity were obtained by laser-Doppler anemometry in a descending solid-liquid turbulent flow in a vertical pipe with volumetric concentrations of suspended spherical particles of 270 μm mean diameter in the range 0.1–14%. Similar measurements were obtained in the flow downstream of an axisymmetric baffle of 50% area blockage placed in the pipe with volumetric concentrations of 310 μm particles up to 8% and of 665 μm particles up to 2%. In order to enable measurements in high particle concentrations without blockage of the laser beams the refractive index of the particles was matched to that of the carrier fluid.

The results show that the particle mean velocity profiles become more uniform and the particle r.m.s. velocity decreases with increasing concentration in both flow cases. The particle mean velocity in the pipe flow also decreases with concentration and the relative velocity, the difference between the particle velocity and the fluid velocity in single-phase flow, decreases with increasing Reynolds number. The length of the recirculation region downstream of the baffle was shorter than in single-phase flow by 11 and 24% for particle concentrations of 4 and 8%, respectively. The particle mean velocities were hardly affected by size for concentrations up fo 2%, but the r.m.s. velocities were lower with the larger particles.     

A combustion problem in cylindrical and spherical symmetry

Summary We study a system of partial differential equations describing the behaviour of a perfect, viscous, polytropic, compressible, chemically reactive gas in a bounded container, under assumptions of cylindrical symmetry. The global existence in the time of a classical solution is proved by some a priori estimates. One extension at the spherical case is given.

---

Sommario Si studia un sistema di equazioni a derivate parziali che descrive il comportamento di un gas perfetto, viscoso, politropico, comprimibile che reagisce chimicamente in un recipiente limitato, con simmetria cilindrica. Si prova I'esistenza globale nel tempo di una soluzione classica per mezzo di stime a priori. Si estende poi questo procedimento al caso delta simmetria sferica.

     

Spectral transport model for turbulence

The spectrum of inhomogeneous turbulence is modeled by an approach that is not limited to regimes of large Reynolds numbers or small mean-flow strain rates. In its simplest form and applied to incompressible flow, the model depends on five phenomenological constants defining the strength of turbulence coupling to mean flow, turbulence transport in physical and wave-number space, and mixing of stress-tensor components. The implications for homogeneous isotropic turbulence are investigated in detail and found to correspond well to the conclusions from more fundamental theories. Under appropriate limiting conditions, a turbulent system described by the model will relax over time into a state of approximate spectral equilibrium permitting a reduction to a one-point model for the system that is substantially like the familiar K- model. This yields preliminary estimates of the present model's parameters and points to the way to improved modeling of flows beyond the applicability of the K- method.