Two Phase Darcy flow with gravity -> problem
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:
- 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 |
All times are GMT -4. The time now is 08:47. |