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

strange pressure behaviour with symmetricPlane boudary condition - interFoam

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

Like Tree4Likes
  • 4 Post By duongquaphim

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 16, 2011, 11:46
Default strange pressure behaviour with symmetricPlane boudary condition - interFoam
  #1
Member
 
Duong A. Hoang
Join Date: Apr 2009
Location: Delft, Netherlands
Posts: 93
Rep Power: 16
duongquaphim is on a distinguished road
Send a message via Yahoo to duongquaphim
Dear Foamer,

Recently, I tested symmetricPlane boundary condition with interFoam solver and found quite strange pressure behaviour. My problem is rather simple:

- A 2D straight channel with dimension of 310x15 micrometer
- Uniform inlet velocity of 0.003m/s and zerogradient pressure
- Uniform outlet pressure of 0Pa and zerogradient velocity
- The top is symmetricPlane
- The bottom is wall
- Only one phase stay in the channel

I ran the case and the results (velocity and pressure) look fine until 1.5ms. After that, there is an increase of pressure and velocity at both inlet and outlet (as in the attachments). Then those high pressure and velocity propagate into the channel and the results look even much strange. The case is also included if anyone want to take a look.

And I should mention here that I had the above problem with two-phase flow simulations (droplet in a channel) so I decided to test again with only one phase (this case) and still got the same weird results.

It really took me quite sometime but still not understand the reason of that behaviour. It would be very appreciative to hear comments from you on that.

Thanks,

Duong
Attached Images
File Type: jpg 1.jpg (17.0 KB, 270 views)
File Type: jpg 2.jpg (19.2 KB, 236 views)
Attached Files
File Type: zip symmetric_cfd.zip (13.2 KB, 54 views)
duongquaphim is offline   Reply With Quote

Old   December 1, 2011, 09:49
Default
  #2
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16
Andrea_85 is on a distinguished road
Hi Duong,

i faced the same problem as you with symmetric Bc and interFoam. Some strange behaviours near the symmetry wall in the pressure and velocity solution that i can not understand. Did you solved it?


best
andrea
Andrea_85 is offline   Reply With Quote

Old   December 1, 2011, 10:02
Default
  #3
Member
 
Rob
Join Date: Sep 2011
Posts: 55
Rep Power: 14
robbirobocop is on a distinguished road
I faced the same problem in a 0.5 million cell case.
Changing "symmetryPlane" to "zeroGradient" fixed the strange velocity behaviour at the symmetry patch. So maybe check that one out!
robbirobocop is offline   Reply With Quote

Old   December 1, 2011, 11:19
Default
  #4
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16
Andrea_85 is on a distinguished road
Hi again,

What you mean with "zeroGradient"? on the velocity? What you used for pressure?

In my case i have an inlet and an outlet (top and bottom of a rectangular domain) and the flow is from top to bottom, so doesn't matter what happens on the right and left side, but i prefer not to use "no flow". I solved the issue using cyclic condition but i'm curious about why symmetry plane doesn't work.

best

andrea
Andrea_85 is offline   Reply With Quote

Old   December 2, 2011, 04:04
Default
  #5
Member
 
Rob
Join Date: Sep 2011
Posts: 55
Rep Power: 14
robbirobocop is on a distinguished road
I used zeroGradient for everything (k, eps, alpha1, U, p_rgh).

I think symmetryPlane might not work if your "patch" is not meshed well. But that's just my opinion and might as well be wrong.

In duongquaphim's case the grid is probably not the problem.

I was trying out his case but since it was build with OF 1.6. I was too lazy to actually change fvSolution and so forth to get it working with OF 2.0.

So maybe Andrea you could check your case with the "zeroGradient" BC and tell me if it works better with it?
robbirobocop is offline   Reply With Quote

Old   December 2, 2011, 06:27
Default
  #6
Member
 
Duong A. Hoang
Join Date: Apr 2009
Location: Delft, Netherlands
Posts: 93
Rep Power: 16
duongquaphim is on a distinguished road
Send a message via Yahoo to duongquaphim
Hi,

I am struggling with this problem for more than two months without any solutions. I already reported a bug at http://openfoam.com/bugs/ with job ID 0000297 but still didn't get respond yet. It might be good if you guys can also join and report that problem.

@ Andrea: I did not resolve it yet. Tried different way but still get the same errors.

@ Rob: I dont think zeroGradient will be a good alternative for symmetryPlane since there is a difference between them. In symmetryPlane, the normal velocity of that plane should be 0 and zero gradient for other variables. So using zeroGradient for all of variables is not a proper way I think.

I tried to mimic symmetryPlane bc by applying bc like

- U: mixed boundary condition: zero normal fixed value, zero tangential gradient
- p and alpha1: zerogradient

but still get the same behaviour: after a while there are some pressure accumulation close to the outlet and then propagate back into the domain. That creates a strange pressure and velocity distribution.

I am looking forward to your comments on this.

Regards,

Duong
duongquaphim is offline   Reply With Quote

Old   December 2, 2011, 09:31
Default
  #7
Senior Member
 
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 16
Andrea_85 is on a distinguished road
Hi Duong,
i agree on the fact that zeroGradient and symmetryPlane are not the same BC, so it depends on what you want to have on that patch.
in my case the symmetry plane is close to another patch (a wall) and so maybe this closeness could create problems. The strange think is that the problem appears only using a fine mesh and not if i use a coarse mesh. Have you tried changing the discretization?

best

andrea
Andrea_85 is offline   Reply With Quote

Old   December 2, 2011, 09:43
Default
  #8
Member
 
Duong A. Hoang
Join Date: Apr 2009
Location: Delft, Netherlands
Posts: 93
Rep Power: 16
duongquaphim is on a distinguished road
Send a message via Yahoo to duongquaphim
Hi Andrea,

That is true. In my case, I like to simulate flow at the T-junction. And Reynolds number is small, I have laminar flow. Therefore, I will have a symmetry plane which split T-junction into two symmetric part. That's why I am interested in symmetryPlane bc. And I can't use zeroGradient in my case.

I've tried a couple of discretization schemes but it is not work. for one phase flow simulation in a straight channel, using Momentum predictor (in FvSolution) helps but if I applied that to a T-junction, that predictor creates another problem.

As far as I found, in the straight channel case, there is a change in direction of velocity at the outlet: small velocity in normal direction appears at the intersecting cell between my symmetryPlane and the outlet. But at the symmetryPlane patch, normal velocity is 0. I think due to that, this velocity increases and increases and also propagates back to the internal domain. That leads to the accumulation of velocity and pressure.

I don't know if it is really the reason for this problem. I haven't found a solution yet . Also I didn't notice on the mesh-dependent part you mentioned. I will give it to try.

Cheers,

Duong
duongquaphim is offline   Reply With Quote

Old   July 26, 2012, 11:58
Default
  #9
Member
 
Duong A. Hoang
Join Date: Apr 2009
Location: Delft, Netherlands
Posts: 93
Rep Power: 16
duongquaphim is on a distinguished road
Send a message via Yahoo to duongquaphim
Hi,

For all who are concerning on this BC behavior, I got a solution:

1) Use momentumPredictor (set in fvSolution)
2) Increase number of corrections in PISO (nCorrectors in fvSolution) to around 6-8. You can play with that number by looking at the initial residual of p in the log file. If the residual decreases to 10^(-5), it is enough I think.

By doing that I got the much better results. Also according to Hrvoje's comments, this problem usually appears in microscale cases.

Cheers,

Duong
pcaron, styleworker, Yage and 1 others like this.
duongquaphim is offline   Reply With Quote

Old   August 20, 2013, 04:10
Default
  #10
Member
 
X Meng
Join Date: Jun 2012
Location: Scotland
Posts: 89
Rep Power: 13
mxylondon is on a distinguished road
I am struggling on this issue now, will let you guys know the result.
mxylondon is offline   Reply With Quote

Old   August 20, 2013, 15:00
Default
  #11
Member
 
X Meng
Join Date: Jun 2012
Location: Scotland
Posts: 89
Rep Power: 13
mxylondon is on a distinguished road
Oh, yes, Duong's method is great! I have tried it with good result.
mxylondon is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
private boundary condition in controlDict strange behaviour markc OpenFOAM Programming & Development 12 November 23, 2011 11:18
Strange behaviour because of contact angle (interfoam) Kim123 OpenFOAM Running, Solving & CFD 0 January 12, 2011 11:16
interFoam: strange pressure with high fluid viscosity? ckroener OpenFOAM 3 April 7, 2010 13:51
How to give pressure boundary condition in SIMPLER algorithm shyamprasad Main CFD Forum 0 March 17, 2009 00:40
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00


All times are GMT -4. The time now is 23:45.