Page tree

You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 22 Next »

(1) \frac{\partial (\rho_A \, \phi) }{\partial t} + \nabla (\rho_A \, {\bf u}_A) = 0, \quad A = \{ O, G, W \}
(2) \int_V \, \frac{\partial (\rho_A \, \phi) }{\partial t} \, dV = - \int_V \, \nabla (\rho_A \, {\bf u}_A) \, dV = - \int_{\partial V} \, \rho \, {\bf u}_A \, d {\bf A}
(3) V \cdot \delta (\rho_A \, \phi) = \delta \, m_A \Rightarrow V \cdot \delta \left( \phi \, \sum_\alpha \rho_{A,\alpha} \, s_\alpha \right) = \mathring{\rho}_A \cdot \delta \, q_A, \quad \alpha = \{ o, g, w \}
(4) \Rightarrow \delta \left( \phi \, \sum_\alpha \frac{\rho_{A,\alpha}}{\mathring{\rho}_A} \, s_\alpha \right) = V^{-1} \cdot \delta \, q_A \Rightarrow \delta \left( \phi \, \sum_\alpha \frac{\mathring{V}_{A,\alpha}}{V_\alpha} \, s_\alpha \right) = V^{-1} \cdot \delta \, q_A
(5) \delta \left( \phi \, \sum_\alpha \frac{\mathring{V}_{A,\alpha}}{V_\alpha} \, s_\alpha \right) = V^{-1} \cdot \delta \, q_A


For the MBO fluid:

(6) \xi_{A, \alpha} = \frac{\mathring{V}_{A,\alpha}}{V_\alpha}
(7) \xi_{O, o} = \frac{\mathring{V}_{O, o}}{V_o} = \frac{1}{B_o}



(8) \xi_{O, g} = \frac{\mathring{V}_{O, g}}{V_g} = \frac{\mathring{V}_{O, g}}{\mathring{V}_{G, g}} \cdot \frac{\mathring{V}_{G, g}}{V_g} = \frac{R_s }{B_g}
(9) \xi_{G, o} = \frac{\mathring{V}_{G, o}}{V_o} = \frac{\mathring{V}_{G, o}}{\mathring{V}_{O, o}} \cdot \frac{\mathring{V}_{O, o}}{V_o} =\frac{R_s }{B_o}
(10) \xi_{G, g} = \frac{\mathring{V}_{G, g}}{V_g} = \frac{1}{B_g}
(11) \xi_{W, w} = \frac{\mathring{V}_{W, w}}{V_w} = \frac{1}{B_w}




Next step is to write the equations explicitly for MBO fluid.


(12) \delta \left[ \phi \cdot \left( \xi_{O,o} \, s_o + \xi_{O,g} \, s_g \right) \right] = V^{-1} \, \delta q_O
(13) \delta \left[ \phi \cdot \left( \xi_{G,o} \, s_o + \xi_{G,g} \, s_g \right) \right] = V^{-1} \, \delta q_G
(14) \delta \left[ \phi \cdot \xi_{W,w} \, s_w \right] = V^{-1} \, \delta q_W

(15) \phi \cdot \left[ \frac{1}{B_o} \, s_o + \frac{R_v}{B_g} \, s_g \right] = V^{-1} \, \delta q_O + \phi_i \cdot \left[ \frac{1}{B_{oi}} \, s_o + \frac{R_{vi}}{B_{gi}} \, s_{gi} \right]
(16) \phi \cdot \left[ \frac{R_s}{B_o} \, s_o + \frac{1}{B_g} \, s_g \right] = V^{-1} \, \delta q_G + \phi_i \cdot \left[ \frac{R_{si}}{B_{oi}} \, s_o + \frac{1}{B_{gi}} \, s_{gi} \right]
(17) \phi \cdot \frac{1}{B_w} \, s_w = V^{-1} \, \delta q_W + \phi_i \cdot \frac{1}{B_{wi}} \, s_{wi}

(18) \phi_n \cdot \left[ \frac{1}{B_o} \, s_o + \frac{R_v}{B_g} \, s_g \right] = V_e^{-1} \, \delta q_O + \left[ \frac{1}{B_{oi}} \, s_o + \frac{R_{vi}}{B_{gi}} \, s_{gi} \right]
(19) \phi_n \cdot \left[ \frac{R_s}{B_o} \, s_o + \frac{1}{B_g} \, s_g \right] = V_e^{-1} \, \delta q_G + \left[ \frac{R_{si}}{B_{oi}} \, s_o + \frac{1}{B_{gi}} \, s_{gi} \right]
(20) \phi_n \cdot \frac{1}{B_w} \, s_w = V_e^{-1} \, \delta q_W + \frac{1}{B_{wi}} \, s_{wi}

With new definitions:

(21) \frac{1}{B_o} \, s_o + \frac{R_v}{B_g} \, s_g = G_O/\phi_n
(22) G_O = V_e^{-1} \, \delta q_O + \left[ \frac{1}{B_{oi}} \, s_o + \frac{R_{vi}}{B_{gi}} \, s_{gi} \right]
(23) \frac{R_s}{B_o} \, s_o + \frac{1}{B_g} \, s_g = G_G/\phi_n
(24) G_G = V_e^{-1} \, \delta q_G + \left[ \frac{R_{si}}{B_{oi}} \, s_o + \frac{1}{B_{gi}} \, s_{gi} \right]
(25) \frac{1}{B_w} \, s_w = G_W/\phi_n
(26) G_W = V_e^{-1} \, \delta q_W + \frac{1}{B_{wi}} \, s_{wi}


The equations can be finally explicitly express sturations:

(27) s_o = \frac{B_o \, (G_o - R_v \, G_G)}{\phi_n \, (1- R_s \, R_v)}
(28) s_g = \frac{B_g \, (G_G - R_s \, G_O)}{\phi_n \, (1- R_s \, R_v)}
(29) s_w = \frac{B_w \, G_W}{\phi_n}


Now summing up and taking into account that  s_o + s_g + s_w = 1 one arrives to a single equation:

(30) \frac{B_o \, (G_o - R_v \, G_G)}{\phi_n \, (1- R_s \, R_v)} + \frac{B_g \, (G_G - R_s \, G_O)}{\phi_n \, (1- R_s \, R_v)} + \frac{B_w \, G_W}{\phi_n} =1
(31) B_o \, (G_o - R_v \, G_G) + B_g \, (G_G - R_s \, G_O)  + B_w \, G_W \, (1- R_s \, R_v)  = \phi_n \, (1- R_s \, R_v)
(32) (B_o - R_s \, B_g) \, G_O +(B_g - R_V \, B_o) \, G_G + (B_w \, G_W - \phi_n )\, (1- R_s \, R_v) = 0



  • No labels