1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: Please help me in solving these

  1. Dec 3, 2007 #1
    Also ,see the attachments for clarity...

    1. The problem statement, all variables and given/known data

    I would like to find a numerical solution for the 3 equations using the conditions.

    ‘w’ refers to water, g to gas (CO2) and ‘o’ to oil
    P=Pressure; u=Velocity; s=Saturation (wetting phase);
    L=distance between CO2 injection and oil production;
    t=time; x=horizontal distance in X-direction; µ=viscosity; ρ=density; ø=porosity
    F=fractional flow of water=relative water mobility/sum of relative mobilities

    F is a function of s and F is ratio of relative premeability and viscosity.

    All the variables are known except s, u and P.

    Initial and boundary conditions
    (x is distance and t is time,P is pressure,s is saturation):
    At t=0, s=1
    At x=0. s=s_wi
    At x=0, P=P_1
    At x=L, P=P_2

    2. Relevant equations

    Mass conservation laws for water and for CO2:
    \phi\frac{\partial s} {\partial t} + \frac{\partial u}{\partial x} \( F(s) \) =0 [/tex]

    \phi \frac{\partial \rho} {\partial t}(P)(1-s)+\frac{\partial \rho}{\partial x}P(u)

    The Darcy Law for both phases, water and gas is

    u = -k (\frac{s(k_(rw))}{\mu_w}+\frac{s(k_(rg))}{\mu_g})
    (\frac{\partial P}{\partial x}) [/tex]

    3. The attempt at a solution

    Consider finite steps,
    (\Delta x) and (\Delta t).

    [tex] {s_i} ^ k =s (i \Delta x,k \Delta t) [/tex] and the same for P and u.

    Then (for the explicit method), we can write approximately using discretization as

    [tex]\frac{\partial s}{\partial t}= \frac{({s_i}^(k+1)-{s_i}^(k))}{\Delta t}[/tex]


    \frac{\partial u}{\partial x}F(s)=\frac{uF_(i+1)^k-uF_(i)^k}{\Delta x} [/tex]

    On substitution in (1), we get an equation for s at (i,k+1) .

    Now , i tried to do the same for the other 2 equations but could not separate the
    variables u and p.Also did not know how to use the initial and boundary conditions.

    But i think the procedure could be like:

    The solution at the layer k=0 (t=0) is known from initial conditions.

    Assume that the solution at layer k has been calculated. In order to find the solution at the layer k+1,

    1) Find the values of saturation s_i^k+1, for each i, from Eq. (s);
    2) Find the values of rho_i^k+1= rho(P_i^k+1) from Eq. (r);
    3) Re-calculate P_i^k+1 based on the known values of rho_i^k+1;
    4) Find the values of u_i^k+1 from Eq. (u).

    Thank you...

    Attached Files:

  2. jcsd
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?
Draft saved Draft deleted