(a)                                                                           (b)                                                                           (c)

(a)  Oil (fully non-wetting fluid, red) is displaced by water (fully wetting fluid, blue) in a porous media (solid is white).

(b)  water (fully wetting fluid, blue) is displaced by Oil (fully non-wetting fluid, red) in a porous media (solid is white).

(c)  Oil and Water are neutrally-wetting ( Water displace Oil )

 

Shan-Chen multi-component multiphase are used.  Density of Oil and Water are equal to 2 mu/lu^3. (mu means  mass unit; lu means lattice unit )

Totally 20,000 ts (time steps) are performed in the simulation.

 

A Non-Wetting fluid displacing a Wetting fluid

                                              

                case 1                                                       case 2                                                                   case3

 

                                               

               case 4                                                       case 5                                                                    case 6

 Figure . 6 cases of different pressure boundary for Blue Non-wetting fluid is displacing the Red wetting fluid. The pressure (density ) are listed in the following Table. The domain is 100*100, pressure inlet/ pressure outlet boundary conditions are applied. All the Figures illustrate the final states of the displacement (20,000 time steps )

  ρ1 is the non-wetting fluid density,  ρ2 is the wetting fluid density.

case

ρ1 (Blue area)

ρ2 (Blue area)

ρ2 (Red area)

ρ1 (Red area)

 dp

R

1/R

1

2.032

0.0092

1.969

0.01

0.02123

15.94517

0.06271

2

2.0366

0.0098

1.961

0.011

0.02561

13.625

0.07339

3

2.012

0.01

1.992

0.01

0.00657

51.32325

0.01948

4

2.0165

0.0094

1.985

0.01

0.01075

31.13605

0.03212

5

2.056

0.012

1.9422

0.011

0.03661

9.00319

0.11107

6

2.06

0.01

1.937

0.01

0.04038

8.03056

0.12452

 

The following Figure is a drawing from the above data. , please notice in our simulations:

 

 The cohesive force acting on the σth component is defined as (Martys and Chen, 1996),

 

 


 

Figure : The Now-Wetting fluid (red) was initialized as a cubic area in the center area, due to cohesion force

the non-wetting fluid would form a sphere finally as the animation shown above.   

 

Figure : 3D Test of Laplace’s Law