
[Sponsors] 
May 3, 2015, 10:50 
Solving pressure in anisotropic porous media

#1 
New Member
Join Date: Apr 2013
Posts: 17
Rep Power: 4 
Hi,
I am trying to solve a pressure equation in a 2D axi reacting porous medium. The permeability of the medium is highly anisotropic: the ratio between the directions is 1/5000. The solvers runs "without problem". Yet, the computed pressure field is problematic. Artefacts appear from time to time. Then, the velocity field is wrong and species transport gets wrong too. Here is a picture of the problem. The artefact is on the third column. I tried my luck with numerical schemes. Yet, it does not seem to change anything. I would like to as you for some advices. Note: the case is transient. Thanks, Nath' 

May 5, 2015, 15:23 

#2 
New Member
Join Date: Apr 2013
Posts: 17
Rep Power: 4 
Bump !
Nathanael. 

May 9, 2015, 04:51 

#3 
New Member
Join Date: Apr 2013
Posts: 17
Rep Power: 4 
No one has a clue ?
Nath' 

May 12, 2015, 14:14 

#4 
New Member
Join Date: Apr 2013
Posts: 17
Rep Power: 4 
Hi,
I am circling around the problem. In case someone is interested. It seems to have something to do with the source term in the pressure equation. Nath' 

May 13, 2015, 02:08 

#5 
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 367
Rep Power: 10 

May 13, 2015, 08:06 

#6 
New Member
Join Date: Apr 2013
Posts: 17
Rep Power: 4 
Hi,
Thank you for your answer. I have to admit that I am not familiar with flux correction. Could you lend me a hand on this problem, please ? If it helps, here is how I solve the pressure equation : Code:
fvScalarMatrix peqn ( fvm::ddt(Poro * M / R / T, p)  fvm::laplacian(Kappa * Rho_g / Mu, p) == MassSource ); peqn.relax(); peqn.solve(); p.relax(); Relaxation factors are set to 1. Then I compute the velocity using Darcy's law : Code:
U =  Kappa & fvc::grad(p) / Mu; phiF = linearInterpolate(U) & mesh.Sf(); Nath' 

May 14, 2015, 02:07 

#7 
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 367
Rep Power: 10 
I am not very familiar with openFOAM, in fact never used it.
If you assume that fluxes could be the problem (very likely actually), then the only way it could happen is that the Rhie Chow part adds some spurious flux that it shall not be adding. In this scenario the origin would be anisotropy in Ap variable for momentum that is used in Rhie Chow flux. I understand that in case of porous flow this Ap var could vary too much and thus some cells could see too much spurious flux coming or going out of them. One way to handle that problem is to limit the amount of Rhie Chow flux allowed. Say for example shall be no more than physical flux (computed by interpolating velocity on face center) Note that in many simulations, Rhie Chow flux is more than physical flux and is necessary for coupling of equations (vp) 

May 16, 2015, 07:42 

#8 
New Member
Join Date: Apr 2013
Posts: 17
Rep Power: 4 
Hi,
Alright. Thank you for the advice. I am currently looking into it. I will come back to you once I am done. Thanks again, Nathanaël. 

July 1, 2015, 04:14 

#9 
New Member
Join Date: Apr 2013
Posts: 17
Rep Power: 4 
Hi dear Foamers,
After a review of literature and with the help of arjun, I have been able to conclude that this problem is a state of the art problem. There is currently no way to solve it. Nathanaël. 

July 1, 2015, 16:44 

#10  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 367
Rep Power: 10 
Quote:
Also it would be interesting to see how fluent does this. My guess is that both of these softwares would do it. 

July 1, 2015, 23:59 
Flood Modelling

#11 
New Member
Ahmad Ali
Join Date: Jun 2015
Posts: 11
Rep Power: 2 
Hello everyone,
I need some guidance.. Which software is best for river & flood modelling ? And which software is good for hydrogeology? Regards, Ali 

Tags 
anisotropy, media, porous 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Wrong fluctuation of pressure in transient simulation  caitao  OpenFOAM Running, Solving & CFD  2  March 5, 2015 22:33 
How to write k and epsilon before the abnormal end  xiuying  OpenFOAM Running, Solving & CFD  8  August 27, 2013 15:33 
pisoFoam with kepsilon turb blows up  Some questions  Heroic  OpenFOAM Running, Solving & CFD  26  December 17, 2012 04:34 
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3  bookie56  OpenFOAM Installation  8  August 13, 2011 04:03 
Orifice Plate with a fully developed flow  Problems with convergence  jonmec  OpenFOAM Running, Solving & CFD  3  July 28, 2011 05:24 