CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Two Phase Darcy flow with gravity -> problem

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 17, 2010, 08:58
Default Two Phase Darcy flow with gravity -> problem
  #1
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18
Cyp is on a distinguished road
Hi everybody!

I am currently working on a Darcy two phase flow solver with IMPES method (IMplicit Pressure Explicit Saturation). I want to simulate a two phase flow in porous media in 3D. It is why I use permeability tensor fields. These permeabilities depend on the liquid saturation

Here is the snippet I programmed in OpenFOAM:
Code:
        M_beta  = K*pow(S,2)/mu_beta;
        M_gamma = K*pow((scalar(1.0)-S),2)/mu_gamma;

        phi_beta  = linearInterpolate( (rho_beta*M_beta & g))  & mesh.Sf();
        phi_gamma = linearInterpolate( (rho_gamma*M_gamma & g)) & mesh.Sf();


        fvScalarMatrix pEqn
        (
            fvm::laplacian((M_beta + M_gamma),p)
         == fvc::div(phi_beta) + fvc::div(phi_gamma)
        );    
        pEqn.relax();
        pEqn.solve();

        fvScalarMatrix SEqn
        (
           fvm::ddt(eps,S) - fvm::laplacian(diffS,S)
        == (fvc::laplacian(M_beta,p)) -  fvc::div(phi_beta)
        );
        SEqn.relax();
        SEqn.solve();

        U_beta = -(M_beta&fvc::grad(p))+((rho_beta*M_beta)&g);
        U_gamma = -(M_gamma&fvc::grad(p))+((rho_gamma*M_gamma)&g);
where
- S and p are volScalarFields
- M_beta and M_gamma are volTensorFields.
- K is a dimensionedTensor

The compilation is ok.

As a validation case, I simulate a the flow through a 1D vertical column: I inject a liquid at the top of column initially filled up with gas. the simulation is correct on a case without gravity (g= (0 0 0)). However, the code diverges when I tried to simulate a flow with gravity (g = (0 -9.81 0)).

Is there something wrong with the gravity in my code ?

Best Regards,
Cyp
Cyp is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
two phase stratified flow. convergence problem moas CFX 0 September 7, 2008 16:32
error in two phase flow problem mvee FLUENT 0 March 28, 2007 06:17
More Probs - Gas-Solid 2 Phase flow Paresh Jain CFX 4 March 27, 2004 18:10
Problem on boundry of two phase flow youngan CFX 2 June 30, 2003 03:11
Two phase flow models Atholl Main CFD Forum 2 May 7, 2002 04:49


All times are GMT -4. The time now is 06:22.