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

porousSimpleFoam: oscillating velocity in the porous zone

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

Like Tree6Likes
  • 2 Post By gschaider
  • 4 Post By gschaider

Reply
 
LinkBack Thread Tools Display Modes
Old   December 3, 2010, 16:49
Default porousSimpleFoam: oscillating velocity in the porous zone
  #1
New Member
 
Sergei D.
Join Date: Mar 2009
Posts: 4
Rep Power: 8
Se9a is on a distinguished road
Hello!
I try to use a porousSimpleFoam solver to simulate a flow in a 2D channel from two equal parts, a free flow part and a porous part (wall's type set in the 'slip' for simplicity):



Base case is a angleDictImplicit.

My results are a good p, but strange oscillating U.
On the next sctrutured mesh:


I get this solution:


Velocity field is:



One can see what velocity is good but slightly oscillated on porous boundary. So, question is:

Are such oscillations caused by a stepwise porous source? If yes, how can I specify smooth porous source?

But more bad situation have place on non-structured mesh.
Mesh is


Solution is:




One can see very oscillating velocity in the porous region.
So, why this happens? Non-orthogonal correction not helped.

Thanks.
Se9a is offline   Reply With Quote

Old   January 3, 2011, 14:28
Default
  #2
New Member
 
Marc-Florian Uth
Join Date: Jan 2010
Posts: 10
Rep Power: 7
Marc10 is on a distinguished road
Hi!

I have exactly the same problem. Did you find a solution for that?

kind regards,

marc
Marc10 is offline   Reply With Quote

Old   February 10, 2011, 08:44
Default
  #3
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
This problem has been reported to http://www.openfoam.com/mantisbt/view.php?id=134. A fix is also described there
cheng1988sjtu and Bahram like this.
gschaider is offline   Reply With Quote

Old   January 27, 2014, 09:27
Default
  #4
New Member
 
Faraj
Join Date: Feb 2010
Posts: 22
Rep Power: 7
Filankes is on a distinguished road
This is not a problem.


Just change under-relaxation factor to 0.0001, and you will not have this oscillations.

OpenFOAM software is the best software among all CFD, and being unexperienced in CFD, does not meam that OpenFOAM has bugs, that should be reported like that...
Filankes is offline   Reply With Quote

Old   January 28, 2014, 06:11
Default
  #5
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by Filankes View Post
This is not a problem.


Just change under-relaxation factor to 0.0001, and you will not have this oscillations.
This is absolute brilliant: under-relaxation is like alcohol: if you take enough you see no problems anymore. Getting things done will take forever, but who cares?
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   January 29, 2014, 19:02
Default
  #6
New Member
 
Faraj
Join Date: Feb 2010
Posts: 22
Rep Power: 7
Filankes is on a distinguished road
Quote:
Originally Posted by gschaider View Post
This is absolute brilliant: under-relaxation is like alcohol: if you take enough you see no problems anymore. Getting things done will take forever, but who cares?
Hi Bernhard,

You are right it takes a lot of time of calculations.

In industrial application it is not good idea to use it.

However, if you use 10 000 - 17 000 cells (hex) for the case proposed in this topic, and make underrelaxation = 0.0001, it will take 30 minutes to have an Exelent solution.

Unfortunately, tetrahedral cells were used, and seems to me under-relaxation 0.0001 will take a day or two for this simulation.

I feel like it is not industrial application, since it is tetrahedral, and 1-2 days is possible to give to get better velocity profile)) corrected pressure loss))

se9a
btw, I forgot to mension - I always use 7 cells in flow direction of porous zone + 0.0001 underrelaxation.

other regions can have very coarse mesh, since they are of no interest I guess
Filankes is offline   Reply With Quote

Old   August 8, 2014, 03:26
Default
  #7
New Member
 
Bahram's Avatar
 
Bahram Haddadi
Join Date: Feb 2014
Location: Vienna, Austria
Posts: 17
Rep Power: 3
Bahram is on a distinguished road
Thanks Bernhard

It worked like magic
Bahram is offline   Reply With Quote

Old   September 23, 2014, 04:18
Default
  #8
New Member
 
Bahram's Avatar
 
Bahram Haddadi
Join Date: Feb 2014
Location: Vienna, Austria
Posts: 17
Rep Power: 3
Bahram is on a distinguished road
Quote:
Originally Posted by gschaider View Post
This problem has been reported to http://www.openfoam.com/mantisbt/view.php?id=134. A fix is also described there
Dear Bernhard

Thanks for the help. I tried it, somehow it works, but as it is also mentioned in there it causes a huge conservation error. Do you have any idea where does this error comes from? any help appreciated!

Best regards
Bahram
Bahram is offline   Reply With Quote

Old   September 23, 2014, 14:01
Default
  #9
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by Bahram View Post
Dear Bernhard

Thanks for the help. I tried it, somehow it works, but as it is also mentioned in there it causes a huge conservation error. Do you have any idea where does this error comes from? any help appreciated!

Best regards
Bahram
Have a look at the link http://www.openfoam.com/mantisbt/view.php?id=134

It is not my claim and I never checked how big HUGE is
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   September 24, 2014, 03:10
Default
  #10
New Member
 
Bahram's Avatar
 
Bahram Haddadi
Join Date: Feb 2014
Location: Vienna, Austria
Posts: 17
Rep Power: 3
Bahram is on a distinguished road
Thanks a lot for quick respond, I have already checked that page but as I said before that idea causes some errors and there is not any clue in there to solve this fluctuating velocity without loosing or gaining some mass!!!
Bahram is offline   Reply With Quote

Reply

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
Heat source in porous zone anger OpenFOAM Running, Solving & CFD 11 December 16, 2013 10:49
Problem in IMPORT of ICEM input file in FLUENT csvirume FLUENT 2 September 9, 2009 01:08
Oscillating velocity and porous zones alberto OpenFOAM Running, Solving & CFD 3 August 18, 2008 15:25
Temperature drop after porous zone MN FLUENT 0 December 10, 2003 13:28
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11


All times are GMT -4. The time now is 11:36.