|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Elwardi Fadeli
Join Date: Dec 2016
Location: Boumerdes, Algeria
Posts: 41
Rep Power: 9 ![]() |
Hi all,
I'm trying to solve the basic incompressible, immiscible two-phase flow that can be modeled with the following equations in a coupled manner: SwEqn: (water equation) ![]() pEqn: (oil equation) ![]() For now, what works for me is to solve: 1- SwEqn explicitly for Sw: All coeffs ![]() ![]() 2- The equation (SwEqn+pEqn) implicitly for p (ddt(Sw) eliminated from pressure equation only if the flow is incompressible ![]() ![]() Of course, the explicit calculation of Sw sets some constraints on deltaT and introduces instabilities when attempting to model ![]() My first attempts to couple these two equations (using foam-extend-4.0) through the terms: 1- pInSw: ![]() 2- SwInp: ![]() did not prove to be successful ![]() The general layout of the coupled solver looks like this: 1- Create the block matrix for the working variable pSw. 2- Assemble pEqn and SwEqn without the coupling terms 3- Calculate coupling terms separately. 4- Add equations and coupling terms to the block matrix 5- Solve it and retrieve solution. Till now, the testing case works fine with the segregated solver but gives weird results with the coupled version. I'm currently investigating the boundary conditions but i don't think there is a problem there. Any help/orientation to reference material is greatly appreciated |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18 ![]() |
Hi ELwardi,
You can have a look at the porousMultiphaseFoam package https://github.com/phorgue/porousMultiphaseFoam that solves the multiphase Darcy equations with the IMPES algorithm. The corresponding paper is there : https://www.sciencedirect.com/scienc...10465514003403 Cheers, Cyprien
__________________
www.cypriensoulaine.com/openfoam |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
Member
Elwardi Fadeli
Join Date: Dec 2016
Location: Boumerdes, Algeria
Posts: 41
Rep Power: 9 ![]() |
Quote:
Thanks for your reply, I have already seen the awesome toolkit ![]() A very simple case where IMPES methods fall short is when applying a semi-implicit well rate (BHP-controled well, you may improve deltaT calculation for this if you want). |
||
![]() |
![]() |
![]() |
Tags |
block matrix, block-coupling, coupled, foam ext 4.0 |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Issues on the simulation of high-speed compressible flow within turbomachinery | dowlee | OpenFOAM Running, Solving & CFD | 11 | August 6, 2021 06:40 |
Solving non-linear coupled equations using blockmatrix solver (OpenFOAM-3.1ext) | Rolanzo | OpenFOAM Running, Solving & CFD | 0 | April 12, 2018 20:49 |
Advice regarding algorithm for solution of discretised equations in compressible flow | lorbekl | Main CFD Forum | 10 | January 22, 2018 01:59 |
DPM coupling with multiphaseInterFoam for gas-oil-water-particle 4phase flow | particlecfd | OpenFOAM Programming & Development | 0 | September 9, 2016 13:09 |
Flow meter Design | CD adapco Group Marketing | Siemens | 3 | June 21, 2011 08:33 |