CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Solving pressure in anisotropic porous media

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 3, 2015, 10:50
Default Solving pressure in anisotropic porous media
  #1
New Member
 
Join Date: Apr 2013
Posts: 17
Rep Power: 4
Nathanael3 is on a distinguished road
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'
Nathanael3 is offline   Reply With Quote

Old   May 5, 2015, 15:23
Default
  #2
New Member
 
Join Date: Apr 2013
Posts: 17
Rep Power: 4
Nathanael3 is on a distinguished road
Bump !

Nathanael.
Nathanael3 is offline   Reply With Quote

Old   May 9, 2015, 04:51
Default
  #3
New Member
 
Join Date: Apr 2013
Posts: 17
Rep Power: 4
Nathanael3 is on a distinguished road
No one has a clue ?

Nath'
Nathanael3 is offline   Reply With Quote

Old   May 12, 2015, 14:14
Default
  #4
New Member
 
Join Date: Apr 2013
Posts: 17
Rep Power: 4
Nathanael3 is on a distinguished road
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'
Nathanael3 is offline   Reply With Quote

Old   May 13, 2015, 02:08
Default
  #5
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 367
Rep Power: 10
arjun is on a distinguished road
Quote:
Originally Posted by Nathanael3 View Post
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'

It is very difficult to tell what the exact problem is.

Fluxes could be related to problem.
arjun is offline   Reply With Quote

Old   May 13, 2015, 08:06
Default
  #6
New Member
 
Join Date: Apr 2013
Posts: 17
Rep Power: 4
Nathanael3 is on a distinguished road
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();
With Poro = porosity, M = gas molar mass, R = ideal gas constant, p = pressure, Kappa = permability tensor, Rho_g = gas density, Mu = viscosity, MassSource = pressure source term.

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();
Thank for your help,
Nath'
Nathanael3 is offline   Reply With Quote

Old   May 14, 2015, 02:07
Default
  #7
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 367
Rep Power: 10
arjun is on a distinguished road
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 (v-p)
arjun is offline   Reply With Quote

Old   May 16, 2015, 07:42
Default
  #8
New Member
 
Join Date: Apr 2013
Posts: 17
Rep Power: 4
Nathanael3 is on a distinguished road
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.
Nathanael3 is offline   Reply With Quote

Old   July 1, 2015, 04:14
Default
  #9
New Member
 
Join Date: Apr 2013
Posts: 17
Rep Power: 4
Nathanael3 is on a distinguished road
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.
Nathanael3 is offline   Reply With Quote

Old   July 1, 2015, 16:44
Default
  #10
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 367
Rep Power: 10
arjun is on a distinguished road
Quote:
Originally Posted by Nathanael3 View Post
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.
If you have access to other software do try it with them too. For example if you have starccm then try this with and without delta-v dissipation.

Also it would be interesting to see how fluent does this.
My guess is that both of these softwares would do it.
arjun is offline   Reply With Quote

Old   July 1, 2015, 23:59
Default Flood Modelling
  #11
New Member
 
Ahmad Ali
Join Date: Jun 2015
Posts: 11
Rep Power: 2
alimalik is on a distinguished road
Hello everyone,

I need some guidance..
Which software is best for river & flood modelling ?

And which software is good for hydrogeology?

Regards,

Ali
alimalik is offline   Reply With Quote

Reply

Tags
anisotropy, media, porous

Thread Tools
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 On
Pingbacks are On
Refbacks are On


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 k-epsilon 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


All times are GMT -4. The time now is 12:17.