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

interFoam 2.1.x gives wrong results on low Froude numbers

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 11, 2012, 08:55
Exclamation interFoam 2.1.x gives wrong results on low Froude numbers
  #1
New Member
 
Ivan
Join Date: Nov 2011
Posts: 8
Rep Power: 5
lt.quibbler is on a distinguished road
Hello, dear foamers!

I've faced a very interesting problem using interFoam 2.1.x for the prediction of the resistance for Kriso Container Ship (KCS).

On relatively high Froude numbers ( 0.22 - 0.28 ) solver gives very good results compared to the experiment, but when I diminish the Froude number the results get worse.

For example, the coefficient for the total resistance on Froude number 0.1 is twice higher
than it should be. And it is so due to high values of pressure force, which is in this case 70-100 percent of the viscous force and oscillating.

When I tried to use pimpleFoam ( for 1-phase problem with symmetry plane) on the same Froude number for the same ship I got correct result. The error for resistance is about 10 percent.

Has anybody faced similar problems?
I just want to know if it is my fault or it is a well-known problem.

Thanks in advance!

Last edited by lt.quibbler; June 11, 2012 at 12:03.
lt.quibbler is offline   Reply With Quote

Old   June 12, 2012, 10:47
Default interFoam 2.1.x. shows strange results
  #2
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
...well after others reproduced the same bug, I sent a bug report concerning interFoam, in this case only affecting the outlet boundary condition but maybe your case is affected too. See the bug report here: http://www.openfoam.org/mantisbt/view.php?id=554
vonboett is offline   Reply With Quote

Old   June 12, 2012, 13:16
Default
  #3
New Member
 
Ivan
Join Date: Nov 2011
Posts: 8
Rep Power: 5
lt.quibbler is on a distinguished road
Albrecht, I have taken a look on the bug. What BC's for p_rgh do you use?
lt.quibbler is offline   Reply With Quote

Old   June 13, 2012, 05:17
Default
  #4
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
Hi lt.quibbler,

I have tried with several settings, in the case shown I use for the outlet:

type outletInlet;
outletValue uniform 0;

and for the atmosphere (top of domain) that reflects the phase aswell I use:

type totalPressure;
p0 uniform 0;
U U;
phi phi;
rho rho;
psi none;
gamma 1;
value uniform 0;

I would be happy if you could tell me if you can reproduce the bug.

Thanks,

Albrecht
vonboett is offline   Reply With Quote

Old   June 13, 2012, 07:23
Default
  #5
New Member
 
Ivan
Join Date: Nov 2011
Posts: 8
Rep Power: 5
lt.quibbler is on a distinguished road
Hi, Albrecht!

Unfortunately, I cant reproduce this bug even using outletInlet. Perhaps you could send me your case? I would take a look if the bug can be reproduced for the same case on my version of FOAM.
lt.quibbler is offline   Reply With Quote

Old   June 13, 2012, 09:46
Default
  #6
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
Hi Ivan,

Thank you for testin!. The case worked on OF 1.7 but produces no outlet if used in OF 2.0 or 2.1.x. Tha attached case (only the basic files) is the working one with x- coordinates smaller zero. To reproduce the bug shift the vertices in blockMesh into positive the positive quadrant.
Attached Files
File Type: gz flumeOutletTestXSmallerZero.gz (7.9 KB, 13 views)
vonboett is offline   Reply With Quote

Old   June 13, 2012, 15:13
Default
  #7
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
Hello Ivan,

meanwhile I reproduced the bug on a Mac OS x 10.6.5 intel, using OpenFOAM 1.7.x.
I attatched the corresponding case folder with positive x coordinates. Again, the bug vanishes when the vertices are moved by -7 m in x direction to the negative quadrant. To create the mesh and run the case, create the polyMesh folder and move the blockMeshDict file in it, had to remove the mesh files due to attatchment size.

Could you post me your case files that dont show the bug, eg. blockMeshDict, alpha1, p_rgh and U as well as svSolution and svSceme ?

Thanks,

Albrecht
Attached Files
File Type: zip flumeOutletTestX+1,7,x.zip (25.1 KB, 7 views)
vonboett is offline   Reply With Quote

Old   June 13, 2012, 15:24
Default
  #8
New Member
 
Ivan
Join Date: Nov 2011
Posts: 8
Rep Power: 5
lt.quibbler is on a distinguished road
Hi Albrecht!!

I think I know what is your problem. You mirrored the mesh, but you forgot to change g from ( 4.905 0 -8.495 ) to ( -4.905 0 -8.495 ) ( the same reflection as for the mesh ). In the case if u dont change the gravity force direction your gravity force makes the fluid go to opposite directon from the outlet.

At least it worked for me.

Please, check if it is the case for u.
Regards
lt.quibbler is offline   Reply With Quote

Old   June 13, 2012, 16:16
Default
  #9
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
... no I didn't mirror the mesh, I yust translated it. However, using for p_rgh zeroGradient at the outlet works, at least in OF 1.7.x, I will check this on OF 2.1.x tomorrow, thought I have already tried this.
But the reflection on the atmosphere gives me still a headache.
vonboett is offline   Reply With Quote

Old   June 13, 2012, 16:28
Default
  #10
New Member
 
Ivan
Join Date: Nov 2011
Posts: 8
Rep Power: 5
lt.quibbler is on a distinguished road
Perhaps, I did the wrong manipulations in blockMeshDict =) Could you, please, send me the version of the case in which I wont need to change anything? Just to run and see the bug
lt.quibbler is offline   Reply With Quote

Old   June 13, 2012, 16:36
Default
  #11
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
Hi Ivan thanks for your time! I cant attatch the grid its too big, yust create the polyMesh folder and move the blockMeshDict from Constants to Constant/polyMesh and then call blockMesh before calling interFoam.
What version of OpenFOAM do you use? This case flumeOutletTestX+1,7 is modified to fit to for OF 1.7.x.

Best wishes,

Albrecht
Attached Files
File Type: zip flumeOutletTestX+1,7,x.zip (25.1 KB, 11 views)
vonboett is offline   Reply With Quote

Old   June 13, 2012, 17:11
Default
  #12
New Member
 
Ivan
Join Date: Nov 2011
Posts: 8
Rep Power: 5
lt.quibbler is on a distinguished road
I finally reproduced the bug on interFoam 2.1.x

Unfortunately, I cant try another BC's right now.

I think that something is wrong with the direction of the normal to the boundary...

Best regards,
Ivan
lt.quibbler is offline   Reply With Quote

Old   June 14, 2012, 08:06
Default
  #13
Senior Member
 
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 177
Rep Power: 6
vonboett is on a distinguished road
Hi Ivan,

thanks for confirming. Henry does not see this as a bug, when using zeroGradient at the outlet for p_rgh the strange behavior vanishes. He is right, it works aswell on OF 2.1.x when using zeroGradient at the outlet for p_rgh. My problem is I have to explain this to my PhD supervisor, and the only guessed explanation I can think of is that it has to do with reintroducing the plus/minus sign after performing calculations with abs(p_rgh) at the cell faces somwhere deep inside the code, but guessing is not explaining.

Anyway I am happy that I can continue with my simulations, but if someone ever finds an explanation I'd be happy to hear from it.

Best wishes,
Albrecht
vonboett 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
Unexpected results from interFoam mconroy OpenFOAM Running, Solving & CFD 2 March 31, 2011 11:05
BuoyantBoussinesqSimpleFoam and axial-symmetric results wrong mass flow Thomas Baumann OpenFOAM 6 December 21, 2009 11:31
why dynamic mesh method give wrong results? weiyang1980 Main CFD Forum 0 September 22, 2009 21:06
'integrate variables' filter on a polyhedral mesh gives wrong results... jbf OpenFOAM Paraview & paraFoam 0 September 4, 2009 04:08
Are the VOF results of interFoam grid independent asaha OpenFOAM Running, Solving & CFD 1 June 26, 2009 04:27


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