CFD Online Logo CFD Online URL
Home > Wiki > Two phase flow

Two phase flow

From CFD-Wiki

(Difference between revisions)
Jump to: navigation, search
(Eddy Interaction Model)
Line 67: Line 67:

Revision as of 13:42, 25 December 2007



article in progress

Importance of two phase flow in industrial configurations

Two phase flow phenomena occur in various industrial application in all fluid mechanics application fiels. Aerospace, automotive, nuclear applications, etc. In all this domain, prediction of two phase behaviour is important. Prediction of liquid spray in an internal combustion engine should enable us to have a better control on combustion process and then to reduce pollutant emissions. Controlling water - steam equilibrium in a coller system enable to prevent from industrial accident, de-icing/anti-icing of aircraft on the ground etc. Any other examples can be quoted here.

Overview of the different available approach

Two main family can be distinguished to model two phase flow, depending of the two phase configuration approach. In case of dispersed configuration a lagrangian approach is suitable. Such an approach consists in following dropplets (or bubbles) during then movement. This is done by applying external force on the particle and solving acceleration, then velocity and finally position. On the other hand, two phase flow can be solve with an eulerien approach. As in all eulerian framework, this approach consists in considering inlet and outlet flux in a given volume. In such an eulerian approach, two family can be distinguished : Mixture model and Two fluids model, those two approach will be detailed in corresponding section bellow.

Lagrangian dispersed two-phase flow modelling

The main goal in Lagrangien approaches is to statiscally particle history in given flow fields. The starting point in Lagrangien approach is the fundamental law of dynamics:

 m_p\frac{dV_p}{dt} =  F(t)

 \frac{dx_p}{dt} = V_p

Where  {F(t)} stands for the resulting force on the particle.

Resolution of these equations requires the knowledge of the instantaneous velocity of the fluid at particle position. The problem is thus to track fluid particles along the discrete particle trajectory. A fluid particle instantaneously owns the velocity of the surrounding fluid and the simulation of its trajectory relies on a quite simple equation such as:  {x_i(t+}{\delta}{t)}={x_i(t)}+{u_i}{\delta}{t}

The instantaneous fluid velocity  {u_i} is decomposed into a mean part which is known (from turbulence model prediction) and a fluctuating part  {u_i}^{'} . So generating fluctuating part of the fluid velocity is the core of the problem. Generation of the fluid particle velocity fluctuations is based on a Gaussian PDF for the fluid velocities, but different random schemes can then be used to guess the fluid velocity, which are characterized by the underlying fluid Lagrangian correlation function   {R_{fl}}.


The first approach was proposed by Gosman and Ioannides (1981) and has been called “eddy life time”. Each fluid velocity fluctuation is kept constant on a time step which is equal to the Lagrangian integral time scale   \tau_l . The resulting Lagrangian correlation function is linearly decreasing from 1 to 0 in a time delay equal to  2\tau_l  :


A first extension has been proposed by Ormancey & Martinon (1984) and is commonly used in the Lagrangian community. In that scheme, a Poisson distribution of the time interval is introduced: each fluid velocity fluctuation is kept constant until a random number (uniformly distributed between 0 and 1) is smaller than  {\delta}{t} / {\tau_l} . The resulting fluid Lagrangian correlation function is exponentially decreasing from 1 to 0:


Later, Desjonqueres (1987) proposed a process which can handle any given correlation function. This can be done through a correlation matrix, as described by Berlemont et al (1990) or Boughattas et al. (2006). Moreover, Desjonqueres (1987) introduced the Frenkiel family of correlation function (1948) which are defined by:

 {R_{fl}(\tau)}=exp[\frac{-\tau}{(m^2+1)\tau_l}]  cos[\frac{m\tau}{(m^2+1)\tau_l}]

Where m is a loop parameter, giving the exponential decrease for m=0. Let us note that this function presents negative loops and m=1 is the value used by Picart (1981) during confrontations of simulations with the experimental data.

Eddy Interaction Model

Historically, the first approach which has been developed and widely used in engineering calculations is the Eddy Interaction Model, first described by Gosman and Ioannides (1981). In the Eddy Interaction Model, the discrete particle is assumed to interact with a succession of eddies. Each eddy is characterized by a velocity (fluctuating), a time scale (lifetime) and a length scale (size). The fluid fluctuating velocity is randomly sampled with a Gaussian PDF of mean and variance that are deduced from turbulence model, as previously described.

Eddy Interation Model (Graham & James 1996)

The eddy life time and the eddy length scale can be estimated from the local turbulence properties and defined by:

 {l_{e}}={c_\mu^{3/4}} \frac{k^{3/2}}{\epsilon}=0.164\frac{k^{3/2}}{\epsilon}


The time  {t_{c}} for a particle to cross the eddy is calculated from the particle velocity (at the beginning of the time step) and the length scale  {l_{e}}.

 {t_{c}}=-\tau_{p}ln[ 1- \frac{l_{e}}{\tau_{p} {\left|\vec u_f-\vec u_p \right|}}]

The particle is assumed to interact with the eddy for a time which is the minimum of the eddy life time and the eddy transit time {t_{c}}.


During that interaction the fluid fluctuating velocity is kept constant and the discrete particle is moved with respect to its equation of motion. Then a new fluctuating fluid velocity is sampled and the process is repeated.

unfortunately, the Eddy Interaction Model still cannot handle anisotropy or inertia effects properly.

Time Correlated Model

Time Corrolated Model (Burry & Bergeles 1993)
Time Corrolated Model (Berlemont et al. 1991)

The Correlated-Time Model, called also two-step methods is similar to the Eddy Interaction Model in a sense that the focus is also to solve the particle trajectory, the difference lies in how the fluctuating fluid velocity is determined along the particle trajectory. this model takes into account the effect of the anisotropy of turbulence on the particulate phase.

The Correlated-Time Model is based on the simultaneous realization of a fluid trajectory and a particle trajectory. The fluid velocity is transferred from the fluid position to the particle position with respect to Eulerian correlation, and the process is then repeated (Burry & Bergeles 1993).

That approach has been extended by Berlemont et al. 1990. he follow the two particle, monitoring their relative separation between time step, and only when the distance beteen them exceeds a certain length scale {L_{D}},than the next fluid trajectory started from the present particle location.

This scheme is imaging the crossing trajectory effects in a very physical way, but it also allows incorporating the continuity effect through the Eulerian correlation.




Where {Cste} ranges between 0.2 and 0.6

One Step Model

A third approach who mixed the Lagrangian and Eulerian step of the previous method, can also be defined to give a simpler stochastic scheme.

Mixture model for two phase flow

Basics of the mixture model

MAC approach

VOF method

Eulerian Two fluids approach

Basics of the two fluids approach

Interfacial exchange closures

Turbulence modelling in such a context



  • Berlemont A., Desjonqueres PH., Guesbet G. (1990), "Particle Lagrangian simulation in turbulent Flows", Int. J. Multiphase Flow, 16-1, pp 19-34.
  • Boughattas N., Gazzah M. H., Said R. (2006), "Effects of a co-flow on particles or droplets dispersion and on droplets vaporization in turbulent air flow", ICAMEM2006, Hammamet, Tunisia.
  • Boughattas N., Gazzah M. H., Said R. (2007), "Lagrangian prediction of particulate two-phase flow", Fifth Mediterranean Combustion Symposium, Monastir, Tunisia.
  • Burry D., Bergeles G. (1993), "Dispersion of particles in anisotropic turbulence", Int. J. Multiphase Flow, 19-4, pp 651-664.
  • Desjonqueres PH. (1987), "Modélisation Lagrangienne du comportement de particules discrètes en écoulement turbulent", Thèse de doctorat, Faculté des Sciences de L’université de Rouen, France.
  • Graham D. I., James P. W. (1996), "Turbulent dispersion of particles using eddy interaction models", Int. J. Multiphase Flow, 22-1, pp 157-175.
  • Frenkiel F.N. (1948), "Etude statistique de la turbulence - fonctions spectrales et coefficients de corrélation", Rapport technique. ONERA n° 34.
  • Gosman A. D., Ionnides I. E. (1981), "Aspects of computer simulation of liqued fuelled combustors", AIAA aerospace sciences meeting, paper 81-0323, St.louis,MO.
  • Ormancey A., Martinon A. (1984), "Prediction of particle dispersion in turbulent flows", Physico-Chemical Hydrodynamics, 5-314, pp 229-244.
  • Picart A., Berlement A., Guesbet G. (1986), "Modelling and predicting turbulence fields and the dispersion of discrete particles transported by turbulent Flows", Int. J. Heat and Mass Transfer, 12-2, pp 237-261.
My wiki