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

Pressure in Channel with simpleFoam and cyclic

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 2, 2014, 19:30
Default Pressure in Channel with simpleFoam and cyclic
  #1
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 7
ArathoN is on a distinguished road
I've done some simulations over a channel using RANS models first with pimpleFoam however because of some problems i switched to simpleFoam (here the thread i created for pimple link)with cyclic Boundaries. I noticed:
  1. it's hard to lower the residuals of pressure and Uy. Tightening the fvsolution tolerances and lowering aspect ratio may help. For Re=15000 this approach worked and the residual decrease to 1e-6 however for Re=64000 i couldn't accomplish this no matter what scheme or IC i've put.
  2. the profile of the Uy and p residuals present an oscillatory behavior they systematically have peaks.
  3. The pressure field is weird (see pictures LS stand for launderSharmaKE and SSY for kOmegaSST) i don't know why it is like that
  4. tHe pressure gradient at the first step is always high and Ubar is alway low after some time they will stabilize.

Does any of you have an explanation about this p field? I'm really lost here.

For the info on the flow i'll link the other thread link
Attached Images
File Type: jpg p_front_sst.jpg (14.0 KB, 24 views)
File Type: jpg p_frontview_LS.jpg (13.3 KB, 22 views)
File Type: jpg p_inlet_sst.jpg (12.6 KB, 18 views)
File Type: png linear.png (7.8 KB, 24 views)
ArathoN is offline   Reply With Quote

Old   November 4, 2014, 08:23
Default
  #2
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,097
Rep Power: 16
RodriguezFatz will become famous soon enough
Hey,

I see the same for a periodic pipe, LaunderSharmaKE
U_no_limiter.jpg
p_no_limiter.jpg
This is a cut through the (short) pipe. Upper profile is "U", flow goes in "x"-direction. X-direction is very poorly resolved by the mesh (see the lower pressure profile). Here you also see the wiggles.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   November 4, 2014, 18:40
Default
  #3
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 7
ArathoN is on a distinguished road
I refined the mesh along x and y direction but the pressure still had that same behavior. I searched more on the forum and in the net but didn't find anything useful. Maybe The pressure gradient source is bugged with simpleFoam, I'll check my case with pimplefoam and report back.
ArathoN is offline   Reply With Quote

Old   November 7, 2014, 11:08
Default
  #4
Member
 
davide basso
Join Date: Jan 2012
Posts: 47
Rep Power: 5
rolloblues is on a distinguished road
This is what my residuals and middle plane profiles of velocity and pressure look like (simpleFoam, Re=44000 approx)

The inlet patch is at the bottom and it is mapped with the top outlet patch.
I didn't use cyclic boundary conditions for inlet and outlet since it would have killed the pressure gradient between in and out.

If you did maybe that's the reason for those strange behaviours. Infact especially in ArathoN's there is barely any pressure gradient at all
Attached Images
File Type: png mapSSQ_residuals.png (13.7 KB, 14 views)
File Type: jpg mapSSQ_middlePlane.jpg (15.5 KB, 19 views)
rolloblues is offline   Reply With Quote

Old   November 14, 2014, 13:41
Default
  #5
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 7
ArathoN is on a distinguished road
Can you upload your /0 folder, I tried to use mappedBC but I had very bad residual. The mapped variable oscillates around 0.01-0.1 almost like your U. In my case I mapped all the field variables except for nut and I mapped to the nearest cell of the outlet (I've put the offset a Lil bit away from the inlet position, I've read that it can cause problems if you put it at the exact position of the outlet). I'll try to do the simulation again and I'll report, for now the u-profile is quite good.
ArathoN is offline   Reply With Quote

Old   November 14, 2014, 14:04
Default
  #6
Member
 
davide basso
Join Date: Jan 2012
Posts: 47
Rep Power: 5
rolloblues is on a distinguished road
There you go
Attached Files
File Type: zip ssdMap0.zip (2.9 KB, 10 views)
rolloblues is offline   Reply With Quote

Old   November 14, 2014, 14:23
Default
  #7
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 7
ArathoN is on a distinguished road
Quote:
Originally Posted by rolloblues View Post
There you go
Interesting you set average values to 0. Is this right? I usually set that flag to false. What interpolation did you use in block mesh dictionary for mappedPatch.
ArathoN is offline   Reply With Quote

Old   November 14, 2014, 15:00
Default
  #8
Member
 
davide basso
Join Date: Jan 2012
Posts: 47
Rep Power: 5
rolloblues is on a distinguished road
I can't give you reference about the zero value for the average. To be honest since the thing worked I didn't pay much attention to it.

Regarding block mesh dictionary I used:

sampleRegion region0;
sampleMode nearestPatchFace;
samplePatch top_end;
rolloblues is offline   Reply With Quote

Old   November 14, 2014, 18:10
Default
  #9
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 7
ArathoN is on a distinguished road
Quote:
Originally Posted by rolloblues View Post
I can't give you reference about the zero value for the average. To be honest since the thing worked I didn't pay much attention to it.

Regarding block mesh dictionary I used:

sampleRegion region0;
sampleMode nearestPatchFace;
samplePatch top_end;
Ah interesting didn't know you can directly specify the patch from which to sample. I usually used the offset option. Thanks I'll try tonight with mappedPatch and I'll report tomorrow.
ArathoN is offline   Reply With Quote

Old   November 21, 2014, 23:25
Default
  #10
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 7
ArathoN is on a distinguished road
Quote:
Originally Posted by rolloblues View Post
I can't give you reference about the zero value for the average. To be honest since the thing worked I didn't pay much attention to it.

Regarding block mesh dictionary I used:

sampleRegion region0;
sampleMode nearestPatchFace;
samplePatch top_end;
I confirm with mapped patch the pressure field is perfect however the residuals will have an oscillatory behaviour and they will not converge to 1e-6. I compared the U-profile of the mappedpatch and cyclic cases and they are identical.

So in the end the pressure field behavior with cyclic is a problem of how the fvOption works with the solver simpleFoam. In my opinion it is needed a pressure correction loop to reassess the continuity.

P. S ah because I made the same error in blockmeshdict it is needed to set the offset value for mappedPatch otherwise it will not work. Moreover if setaverage is 0 there is no averaging.
ArathoN 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
fully developed channel with cyclic using simplefoam ArathoN OpenFOAM Running, Solving & CFD 4 February 21, 2014 10:53
ribbed channel / simpleFoam / boundary conditions beeo OpenFOAM Pre-Processing 20 July 17, 2013 08:39
Cyclic boundary conditions concerning channel with constant flow rate pascool OpenFOAM Pre-Processing 1 November 16, 2012 14:04
How to setup cyclic BCs in simpleFOAM hhuang84 OpenFOAM Running, Solving & CFD 12 July 16, 2012 07:01
multiple cyclic fans and outlets irregular ( checkerboard ) pressure on patch faces ajl42 OpenFOAM 0 November 30, 2011 09:05


All times are GMT -4. The time now is 01:16.