A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm ECP, 2011-2012
DERIVATION OF CONTINUUM MECHANICS BALANCE EQUATIONS 1. First principles (4) Leibniz rule and Gauss theorem. On material and arbitrary control volumes. 2. Local instantaneous balance equations (single-phase). The closure issue (I) Fixed volume with an interface (discontinuity surface). 3. Local instantaneous balance for each phase and the interface (jump conditions). Space averaging: 1D balance equations. Time averaging: 3D local balance equations (Reynolds style). Composite averaging: two-fluid model. 4. The closure issue (II) Two-phase flow balance equations 1/39
MATHEMATICAL TOOLS 5 K L L L K K Displacement velocity of a surface, S: ( ) M v S t u,v Depends on the choice of parameters. Implicit equation: f(x, y, z, t) 0 inside V f(x, y, z, t + t) = f(x 0, y 0, z 0, t), BN O J, J BN O J + f(m 0 ) M + f t t + Geometrical displacement velocity (intrinsic, scalar): v S n = lim t 0 M t f = t f Two-phase flow balance equations 2/39
LEIBNIZ RULE L 5 L 5 3D-extension of the derivation of integrals theorem: d f f dv = dt V (t) V (t) t dv + fv S n ds S(t) 8 J 5 J Differential geometry theorem, S arbitrary. n points outwardly (always). Use: commutes time derivative and space integration. Material control volumes arbitrary volumes. Two-phase flow balance equations 3/39
GAUSS THEOREM 5 8 Divergence is the flux per unit volume: 1 B lim n B ds ɛ 0 V ɛ Divergence theorem, Gauss-Otstrodradski (Green) : B dv = n B ds V (t) S S(t) Differential geometry theorem, S and V arbitrary, n et on the same side. n, points outwards. B, arbitrary tensor. Use: some particular volume integrals surface integrals. Two-phase flow balance equations 4/39
MATERIAL VOLUMES-ARBITRARY VOLUMES Let V m (t), limited by S m (t) be a material volume : v Sm n = v n. d f f dv = dt t dv + fv n ds V m (t) V m (t) S m (t) Consider V (t) which coincides with V m (t) at t. d f f dv = dt t dv + V (t) V (t) S(t) fv S n ds Identity: for all V (t) which coincides with V m (t) at time t, d dt V m (t) f dv = d dt V (t) f dv + S(t) f(v v S ) n ds Two-phase flow balance equations 5/39
A SIMPLE EXAMPLE: MASS BALANCE Principle: the mass of a material volume is constant. d ρ dv = 0 dt V m (t) Use the identity with f = ρ, d ρ dv dt V (t) }{{} Mass of V, m + S(t) ρ(v v S ) n ds }{{} Net mass flux leaving S, M = 0 The time variation of the mass of V, m, equals the net incoming mass rate, M. dm dt + M = 0, dm dt = M First principles can be formulated on material or arbitrary volumes. statements are equivalent. Both Two-phase flow balance equations 6/39
MASS BALANCE The time variation of the mass equals the net mass flow rate entering in the volume V ( V ). d ρ dv = ρ(v v S ) n ds. (1) dt Particular cases, For a fixed volume, v S n = 0, For a material volume, v S n = v n V S Two-phase flow balance equations 7/39
SPECIES BALANCE The time variation of the mass of component α equals (i) the net incoming mass rate of α and (ii) the production in the volume V ( V ). d ρ α dv = ρ α (v α v S ) n dv + r α dv dt V Add all equations for α gives the mixture mass balance. α r α = 0. Chemicals redistribution, no overall net mass production. S V Two-phase flow balance equations 8/39
LINEAR MOMENTUM BALANCE The time variation of the linear momentum equals the sum of (i) the incoming momentum flux, (ii) the applied forces ( V ). d ρv dv = ρv(v v S ) n ds + n T ds + ρg dv (2) dt V T: stress tensor, contact forces. g: volume forces. NB: vector equation. S S V Two-phase flow balance equations 9/39
ANGULAR MOMENTUM BALANCE The time variation of the momentum moment equals the sum of (i) the net incoming flux of moment of momentum and (ii) the applied torques ( V ). d ρr v dv = ρr v(v v S ) n ds + r (n T) ds + r ρg dv dt V S S V (3) When torques results only of applied forces (non polar fluids). Take two get the third. The stress tensor is symmetrical. The linear momentum balance is satisfied. The angular momentum balance is satisfied. Two-phase flow balance equations 10/39
TOTAL ENERGY BALANCE Equivalent to the first principle of thermodynamics: the time variation of the total energy of a closed system equals the sum of (i) the thermal power added and (ii) the power of external forces applied to the system. The time variation of the total energy (internal and mechanical) equals the sum of (i) the incoming total energy flux, (ii) the mechanical power of the applied forces and (iii) the thermal power given to the system.( V ). d dt V ρ (u + 12 ) v2 dv = + ρ (u + 12 ) v2 (v v S ) n ds S (n T) v ds + ρg v dv S V S q n ds + q : volume heat sources (Joulean heating, radiation absorption, etc.) NOT of thermodynamical origin, heat of reaction, phase transition of any order... The process is arbitrary: reversible or not. V q dv (4) Two-phase flow balance equations 11/39
ENTROPY BALANCE AND SECOND PRINCIPLE The time variation of the entropy of a closed and isolated system is non negative. The time variation of the entropy equals (i) the net inflow of entropy, (ii) the entropy given to the system in a reversible manner, (iii) the entropy sources ( V ). d dt V ρs dv = S ρs(v v S ) n ds σ 0. The second principle is only σ 0. When reversible, σ = 0. S n j s ds + V q T dv + V σ dv, (5) Two-phase flow balance equations 12/39
GENERALIZED BALANCE EQUATION Balance equations have similar forms, d ρψ dv = n ρ(v v S )ψ ds dt V S S n j ψ ds + V φ ψ dv. Balance ψ j ψ φ ψ Mass 1 Species α ω α j α r α L. momentum v T ρg A. momentum r v T R ( ) r ρg Total energy u + 1 2 v2 q T v ρg v + q Entropy s j s σ + q T (*)R, R ij = ɛ ijk r k Two-phase flow balance equations 13/39
PRIMARY LOCAL EQUATIONS Leibniz rule, V ρψ t dv = S Gauss theorem, V D f, [ ρψ t Instantaneous local primary balances, V n ρvψ ds S n j ψ ds + + (ρvψ) + j ψ φ ψ ] V dv = 0 φ ψ dv. ρψ t = (ρvψ) }{{} Convection j ψ }{{} Diffusion +φ ψ }{{} Source Balance on a fixed and infinitesimal volume, strictly equivalent to first principles. Two-phase flow balance equations 14/39
TOTAL FLUX FORM Total flux form (Bird et al., 2007), stationary flows. ρψ t = j t ψ + φ ψ Balance total flux convective flux diffusive flux j t ψ ρψv j ψ Mass n = ρv Species n α = ρω α v j α Momentum φ = ρvv T Total energy e = ρv ( u + 1 2 v2) q T v Entropy j t s = ρsv j s NB: Some authors may use different sign conventions for fluxes. Don t pick up an equation from a text without care... Two-phase flow balance equations 15/39
CONVECTIVE FORM Combine with the mass balance, ρ t = (ρv) Expand products in the balance equation, ρψ t = (ρvψ) j ψ + φ ψ ρ ψ t + ψ ρ t = ψ (ρv) ρv ψ j ψ + φ ψ Definition of the convective derivative: Df Dt = f t + v f ρ Dψ Dt = j ψ + φ ψ Balance on a material volume (infinitesimal). Only diffusive fluxes. Practical form to derive secondary equations. Two-phase flow balance equations 16/39
SUMMARY OF CONTINUUM MECHANICS EQUATIONS For a pure fluid, on an arbitrary control volume, Mass balance (1) Linear momentum balance (2) Angular momentum balance (3) Total energy balance (4) Entropy inequality (5) Local primary balance equations, (1) Mass balance (6) (2) Momentum balance (7) (3) Stress tensor symmetry (4) Total energy balance (8) (5) Entropy inequality (9) Two-phase flow balance equations 17/39
CONTINUUM MECHANICS EQUATIONS Secondary balance equations, for a pure fluid, Mechanical energy balance, (10) v momentum balance. Internal energy balance (11), total energy balance (8)-(10). Enthalpy balance (12). (11), h u + p/ρ Entropy balance (13), (11), du = T ds pdv (Gibbs). Comparing to entropy inequality (9), provides j s and σ. Two-phase flow balance equations 18/39
THE CLOSURE ISSUE (I) In balance equations, Local variables, v, v α, p, u, etc. Unknown fluxes, j α, T, q, j s. NB: T = pi + V Unknown sources, r α, σ. First principles cannot provide expressions for fluxes. The CME are not closed. An extended interpretation of the second principle, provides entropy sources. For a pure fluid, T σ = q T + V : v. provides the thermodynamic equilibrium conditions, σ = 0, provides constraints on possible closure to ensure return to equilibrium. Linearity assumption, transport properties, T = µ( v + v ) + (ζ 2 µ) v, q = κ T, µ, ζ, κ 0 3 Transport properties must be measured or modeled beyond CME scope. Two-phase flow balance equations 19/39
TWO-PHASE LOCAL BALANCE EQUATIONS ) J 8 J ) E J 8 J ) J d dt For V 1 (t) (not fixed), Leibniz rule: d ρ 1 ρ 1 dv = dt V 1 V 1 t dv + Example: mass balance, V = V 1 V 2, A = A 1 A 2 fixed. Interfaces,surface of discontinuity. d ρ dv = ρv n ds, V dt V Split contributions from V 1 and V 2 : ρ dv + d ρ dv = ρv n ds ρv n ds V 1 dt V 2 A 1 A 2 A i (t) A ρ 1 v Ai n 1 da Gauss theorem: A 1 ρv 1 n 1 ds = V 1 (ρ 1 v 1 ) ds A i ρ 1 v 1 n da Apply the same procedure for V 2, sum up all contributions, Two-phase flow balance equations 20/39
TWO-PHASE MASS BALANCE Collect integral terms wrt dimension, V, 2 k=1 V k ( ) ρk t + (ρ kv k ) dv A i (ρ 1 (v 1 v i ) + ρ 2 (v 2 v i )) da = 0 Local mass balance, k = 1, 2, for all points in V k (PDE), ρ k t + (ρ kv k ) = 0 For all points of the interface, jump condition, ρ 1 (v 1 v i ) n 1 + ρ 2 (v 2 v i ) n 2 = 0 }{{}}{{} ṁ 1 ṁ 2 Jump condition is the mass balance of the interface, ṁ k = ρ k (v k v i ) n k. Two-phase flow balance equations 21/39
LOCAL BALANCE EQUATIONS Use the generalized local balance equation, same procedure, V 2 k=1 V k ( ρk ψ k A i t + (ρ k ψ k v k ) + (j ψk ) φ k ) dv 2 (ṁ k ψ k + n k j ψk + φ i ) da = 0 k=1 At every points of each phase, ρ k ψ k t + (ρ k ψ k v k ) + (j ψk ) φ k = 0 At every points of the interface, jump condition, balance of the interface. 2 (ṁ k ψ k + n k j ψk + φ i ) = 0 k=1 φ i : entropy source at the interface. Two-phase flow balance equations 22/39
JUMP CONDITIONS Mass balance, ρ 1 (v 1 v i ) n 1 + ρ 2 (v 2 v i ) n 2 = 0 ṁ 1 + ṁ 2 = 0 No phase change: ṁ k = 0, ṁ 1 = ṁ 2 = 0 Assumption: no slip at the interface (φ i = 0), (v 1 v i ) n 1 = 0, (v 2 v i ) n 2 = 0 (v 1 v 2 ) n 1 = 0 v 1 = v 2 Two-phase flow balance equations 23/39
JUMP CONDITIONS (CT D) Momentum balance, ṁ 1 v 1 + ṁ 2 v 2 n 1 T 1 n 2 T 2 = 0 When no viscosity, T = pi + V, v = v t + v n, v n = n(v n), ṁ 1 (v1 n v2 n ) + (p 1 p 2 )n 1 = 0 v1 t = v2 t General case, ṁ 1 (v 1 v 2 ) + (p 1 p 2 )n 1 n 1 (V 1 V 2 ) = 0 Two-phase flow balance equations 24/39
JUMP CONDITIONS (CT D) Particular case: 1D flow, v k (x) interface : ṁ 1 (v 1 v 2 ) n 1 + (p 1 p 2 ) n 1 (V 1 V 2 ) n 1 = 0 1D incompressible flow, dv k dx = 0 V k = 0, ṁ 1 (v 1 v 2 ) n 1 + (p 1 p 2 ) = 0 From the mass balance, definition: ṁ k = ρ k (v k v i ) n k, Results, pressure jump, recoil force, p 1 p 2 ρ 1 ρ 2 whatever ṁ 1. ṁ 1 ( 1 ρ 1 1 ρ 2 ) = (v 1 v 2 ) n 1 p 1 p 2 = ρ 1 ρ 2 ρ 1 ρ 2 ṁ 2 1. Two-phase flow balance equations 25/39
MOMENTUM BALANCE AND SURFACE TENSION ) J 8 J + J ) E J 8 J ) J Momentum balance on a fixed volume. Forces: = σn dl + n k T k ds + ρ k F k dv. C(t) A 1 A 2 V 1 V 2 Théorème de Gauss Aris (1962), Delhaye (1974) : σn dl = ( S σ nσ S n) ds C(t) A i (t) S : surface gradient, S : surface divergence. Momentum balance interface: ṁ 1 v 1 + ṁ 2 v 2 n 1 T 1 n 2 T 2 = S σ + nσ S n S σ: Marangoni force, nσ S n: capillary pressure, Laplace pressure jump. nσ S n = 2Hn H: mean curvature of the surface. Circular cylinder: 1/R, sphere: 2/R. Two-phase flow balance equations 26/39
EXAMPLE: 2D INTERFACES Momentum jump at the interface, ṁ 1 v 1 + ṁ 2 v 2 n 1 T 1 n 2 T 2 + dσ dl τ σ R n = 0 For a non viscous fluid, T = pi + V,no phase change, n 1 (p 1 p 2 ) + dσ dl τ σ R n = 0 Laplace relation, : (p 1 p 2 ) = σ R n 1 n Inconsistency, // : µ k = 0 dσ dl = 0 Marangoni effect for viscous fluids, σ(t ), σ(c), (n 1 V 1 + n 2 V 2 ) τ + dσ dl = 0 Be careful with the parameter selection. Pressure is always higher in the concavity side (balloon). Two-phase flow balance equations 27/39
JUMP CONDITIONS (CT D) Total energy balance: ṁ 1 (u 1 + 1 ( 2 v2 1 )+ṁ 2 u 2 + 1 ) 2 v2 2 +q 1 n 1 +q 2 n 2 n 1 T 1 v 1 n 2 T 2 v 2 = 0 Enthalpy form, 3 common assumptions, phase change is the dominant effect, variation of mechanical energy can be neglected, the effect of pressure and viscous stress jump can be neglected (no surface tension), ṁ 1 h 1 + ṁ 2 h 2 + q 1 n 1 + q 2 n 2 = 0 More on the derivation, see Delhaye (1974, 2008). Thermodynamic equilibrium condition at the interface: v1 t = v2, t T 1 = T 2, g 1 g 2 = 1 ( 1 2ṁ2 1 ρ 2 1 ) ( n2 V 2 n 2 2 ρ 2 n ) 1 V 1 n 1 1 ρ 2 ρ 1 Two-phase flow balance equations 28/39
USE OF LOCAL EQUATIONS First principles balance on arbitrary control volumes Local phase equations, Jump conditions at the interface (see also the Rankine-Hugoniot eqs). Flows with simple interface configuration Stability of a liquid film, Growth/collapse of a vapor bubble (nucleate boiling, cavitation). More general problems, Tremendously large number of interfaces, non-equilibrium. Large scale fluctuations, intermittency, engineers seek for mean values. Space averaging of local equations (area-averaged): 1D models Time averaging of local equations: CMFD (3D codes) Space and time averaging, composite averaging: two-fluid 1D model, 1D codes, system codes. Two-phase flow balance equations 29/39
AREA-AVERAGED BALANCE EQUATIONS Area-averaging operator: < f k > 2 = 1 A k A k f k da How to get a balance equation for a mean value? Average the local balance on A k. Example, mass balance, + ρ k t + (ρ kv k ) = 0 Integrate on A k, ρ k t da + (ρ k v k ) da = 0 A k A k + ) + Limiting forms of the Leibniz rule and Gauss theorems, ρ k da + + ρ k w k da + = 0 t A } k z A {{}} k {{} A k <ρ k > 2 A k <ρ k w k > 2 Two-phase flow balance equations 30/39
MATHEMATICAL TOOLS Limiting form of the Leibniz rule, f k f k da = t A k t da + dl f k v i n k C k n k n kc A k + Limiting form of the Gauss theorem, BdA = n z BdA + A k z A k n k B dl C k n k n kc First useful identity, B = n z A k z = n k n z C k dl n k n kc (1) + ) + Second useful identity, B = pi A k p da = z A k pn z da + C k pn k dl n k n kc (2) Two-phase flow balance equations 31/39
AREA-AVERAGED BALANCE EQUATIONS (CT D) Integrate on A k, A k ρ k t da + Leibniz rule and Gauss theorem, A k (ρ k v k ) da = 0 t A k < ρ k > 2 + z A k < ρ k w k > 2 = C k ṁ k dl n k n kc Γ k : production rate of phase k [kg/s/m] per unit length of pipe. dl Γ k = ṁ k C k n k n kc No phase change: ṁ k = 0 Γ k = 0. Mass balance of the interface, ṁ 1 + ṁ 2 = 0 Γ 1 + Γ 2 0. The area-averaged mass balance is not closed. Two-phase flow balance equations 32/39
AREA-AVERAGED BALANCE EQUATIONS Based on the general form of the local balance equations., t A k < ρ k ψ k > 2 + z A k < n z ρ k v k ψ k > 2 + z A k < n z j ψk > 2 A k < φ k > 2 dl dl = (ṁ k ψ k + n k j ψk ) n k j ψk C i n k n kc C k n k n kc δa k = C i C k, C k wall fraction wetted by k, C i, interface. + ) + + E + + + E Two-phase flow balance equations 33/39
MOMENTUM BALANCE Note on the momentum balance, vector equation, t A k < ρ k v k > 2 + z A k < ρ k w k v k > 2 z A k < n z T k > 2 A k < ρ k g k > 2 dl dl = (ṁ k v k n k T k ) + n k T k C i n k n kc C k n k n kc Projection on n z : right dot product, w k = v k n z, stress tensor decomposition, t A k < ρ k w k > 2 + z A k < ρ k wk 2 > 2 + z A k < p k > 2 z A k < n z V k n z > 2 dl dl A k < ρ k g z > 2 = (ṁ k w k n k T k n z ) + n k T k n z C i n k n kc C k n k n kc Identity (1), assume < p k > 2 = p C, dl p k n k n z = p C n k n z C k C i n k n kc C k C i dl A k =< p k > 2 n k n kc z Other choices are possible, introduce an excess pressure, p i, p C =< p k > 2 +p i Two-phase flow balance equations 34/39
MOMENTUM BALANCE (CT D) Momentum balance, single-pressure, < p k > 2 = p C, t A k < ρ k w k > 2 + z A k < ρ k wk 2 > 2 +A k z < p k > 2 z A k < n z V k n z > 2 dl dl A k < ρ k g z > 2 = (ṁ k w k n k V k n z ) + n k V k n z C i n k n kc C k n k n kc Transfers are dominant in the radial direction, quasi fully developed flows, z A k < n z V k n z > 2 ( ) w ν t 0 z z Give an example of a situation where this term might not be neglected. Closures are required: interactions at the interface and the wall (wall friction). Two-phase flow balance equations 35/39
B TIME AVERAGED BALANCE EQUATIONS 6 Conditional time-averaging, = 1 f k dt = T k f X T X kf k dt k [T k ] T X = X kf k kdt X k Plain time-average, J J J J J J J 6 f = 1 fdt J 6 T T 6 How to get a balance equation for a mean value? Average the local balance on [T k ]. (ρ k v k ) dt = 0 [T k ] ρ k t dt + Limiting forms of the Leibniz rule and the Gauss, theorem, ρ k dt + + ρ k v k dt + = 0 t [T k ] } {{ } T k ρ X k [T k ] [T k ] } {{ } T k ρ k v X k Two-phase flow balance equations 36/39
MATHEMATICAL TOOLS Limiting form of the Leibniz rule (derivation of an integral wrt upper limit), f k [T k ] t dt = f k dt v i n k f k t [T k ] v i n k disc. [T ] }{{} ±1 Limiting form of the Gauss theorem, B k dt = [T k ] [T k ] B k dt + disc. [T ] n k B k v i n k Time-averaged mass balance, production rate of phase k, interfacial interactions are homogenized [kg/s/m 3 ], α k ρ k X t + ( α k ρ k v k X ) = 1 T disc. [T ] ṁ k v i n k Time-averaged balance equations, starting point of the Reynolds decomposition, T? α k ρ k ψ k X t ( ) ( ) X X X 1 + α k ρ k v k ψ k + α k j ψk α k φ k = T disc. [T ] ṁ k ψ k + n k j ψk v i n k Two-phase flow balance equations 37/39
COMPOSITE AVERAGES: THE TWO-FLUID MODEL Example: mass balance, space-averaged and time averaged, t A k < ρ k > 2 + z A dl k < ρ k w k > 2 = ṁ k C k n k n kc Time averaged and space averaged, t A< α kρ k X > 2 + z A< α kρ k w k X > 2 = A< 1 T disc. [T ] ṁ k v i n k > 2 LHS are identical, the RHS should also. Proof, identity on interaction terms, Local specific interfacial area, γ = 1 1 T v i n k disc. [T ] Possible closure of interaction terms: interfacial area mean flux, 1 T disc. [T ] ṁ k v i n k = ṁki T disc. [T ] 1 v i n k = γṁ ki Two-phase flow balance equations 38/39
REFERENCES Aris, R. 1962. Vectors, tensors, and the basic equations of fluid mechanics. Prentice-Hall. Bird, R. B., Stewart, W. E., & Lightfoot, E. N. 2007. Transport phenomena. Revised second edn. John Wiley & Sons. Delhaye, J. M. 1974. Jump conditions and entropy sources in two-phase systems, local instant formulation. Int. J. Multiphase Flow, 1, 359 409. Delhaye, J.-M. 2008. Thermohydraulique des réacteurs nucléaires. Collection génie atomique. EDP Sciences. Two-phase flow balance equations 39/39