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

OpenFOAM 1.6 and 1.7 with interFoam, groovyBC give different strange results

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 11, 2010, 10:28
Unhappy OpenFOAM 1.6 and 1.7 with interFoam, groovyBC give different strange results
  #1
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 8
Arnoldinho is on a distinguished road
Hi all,

I have recently updated to OpenFOAM 1.7.0, after not having done simulations in Openfoam for a longer time.

I have a problem with a test case of a wave tank using interFoam and groovyBC on the inlet (left side) to generate waves in the flume. On 1.6 (or maybe it was 1.5), everything went fine, waves are running into the channel (left picture below). On 1.7.0, the same case with corrected p > p_rgh, gamma > alpha1 gives strange results and of course aborts after some time with a 'floating point exception'. Waves might be generated directly on the left inlet boundary, but the water is running up the front wall and down the back wall.

For net generation, Salome was used. The case files can be found at http://130.75.108.10/~material/pdf_downloads/20_wellentank_fine_complete.tar.gz

Pictures are attached. Left: OpenFoam 1.7, right, OpenFoam 1.6, for a similar time step.


Have you got any ideas of whats wrong here?

Arne


Attached Images
File Type: jpg Bildschirmfoto-1.jpg (23.5 KB, 203 views)

Last edited by Arnoldinho; August 11, 2010 at 12:05. Reason: New picture
Arnoldinho is offline   Reply With Quote

Old   August 11, 2010, 12:38
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Arne,

I'm not familiar with groovyBC (yet) but I've got a feeling that it might be an issue with the compiler!

What gcc version(s) did you use to build those OpenFOAM versions and respective groovyBCs?

If you are wandering why I ask... this is the reason: OpenFOAM 1.6.x, 1.7.0 and 1.7.x are not fully prepared to work with gcc-4.5.x

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   August 11, 2010, 13:11
Default
  #3
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 8
Arnoldinho is on a distinguished road
Hi Bruno,

thanks for your hint. I'm running OpenFOAM on a 64bit Ubuntu 10.04.

A gcc -v gives me: gcc version 4.4.3 (Ubuntu 4.4.3-4ubuntu5). For the OpenFOAM 'installation', I did it the way described in http://www.openfoam.com/download/ubuntu.php.

Btw: Could anywone who is using groovyBC (1.6) and OpenFoam 1.7 test my case?

Arne
Arnoldinho is offline   Reply With Quote

Old   August 12, 2010, 04:53
Default
  #4
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,912
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by Arnoldinho View Post
Hi all,

I have recently updated to OpenFOAM 1.7.0, after not having done simulations in Openfoam for a longer time.

I have a problem with a test case of a wave tank using interFoam and groovyBC on the inlet (left side) to generate waves in the flume. On 1.6 (or maybe it was 1.5), everything went fine, waves are running into the channel (left picture below). On 1.7.0, the same case with corrected p > p_rgh, gamma > alpha1 gives strange results and of course aborts after some time with a 'floating point exception'. Waves might be generated directly on the left inlet boundary, but the water is running up the front wall and down the back wall.

For net generation, Salome was used. The case files can be found at http://130.75.108.10/~material/pdf_downloads/20_wellentank_fine_complete.tar.gz

Pictures are attached. Left: OpenFoam 1.7, right, OpenFoam 1.6, for a similar time step.


Have you got any ideas of whats wrong here?

Arne


No offense. But have you checked the direction of gravity?
gschaider is offline   Reply With Quote

Old   August 12, 2010, 06:14
Default Check groovyBC
  #5
mks
New Member
 
Join Date: Nov 2009
Posts: 10
Rep Power: 7
mks is on a distinguished road
The direction of gravity may by wrong, check the orientation of your mesh and correct (or change the g-direction in constant/g). Check also the GroovyBC in 0/U, 0/alpha1. It also depends on the axes orientation(

valueExpression "(pos().y<=A*cos(-w*ti.......
.
.
variables "l=10;A=0.5;g=vector(0,-9.81,0);k=2*.....
)
mks is offline   Reply With Quote

Old   August 12, 2010, 09:47
Default
  #6
Senior Member
 
Arne Stahlmann
Join Date: Nov 2009
Location: Hanover, Germany
Posts: 209
Rep Power: 8
Arnoldinho is on a distinguished road
Thanks mks,

it was the orientation of 'g' in the file. I shouldn't copy files from a tutorial case without looking if they fit for my case...

You know, always this problem with the wood and the trees...


Another question: What kind out boundary condition would you suggest for the outlet of the flume? The waves shall 'run out' of the flume without influencing the next incoming waves - so some kind of a spongle layer. Is there already one implemented?

Arne
Arnoldinho is offline   Reply With Quote

Old   August 12, 2010, 11:14
Default
  #7
mks
New Member
 
Join Date: Nov 2009
Posts: 10
Rep Power: 7
mks is on a distinguished road
well... if there is something ready to use, i have not found it... it is still an issue for me how to let the waves go "through" the wall. I just made the flume tank long enough:-)
mks is offline   Reply With Quote

Old   December 9, 2010, 17:29
Default
  #8
New Member
 
Join Date: Dec 2010
Posts: 4
Rep Power: 6
abhi25_itbhu is on a distinguished road
Hey all,

I am trying to simulate a wave tank in OpenFoam 1.7.1.
I am able to get groovyBC working on my system. I followed the example of 2D stokes wave available on the WiKi pages. I really loved this boundary condition as it is relatively easier to understand and implement.
Problem is the reflection from the outlet boundary. Currently I am using zeroGradient boundary condition on the outlet boundary.
I was wondering if anybody came up with a solution to transmit the waves.
Has anybody tried the waveTransmissive boundary condition at the outlet boundary? Is it pertinent to this situation? If yes, I would like to know how to implement it.

I also read that the inclusion of sponge layer near the outlet boundary may solve the problem. Is anybody familiar with it?

Thanks a lot in advance.

Abhi
abhi25_itbhu 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
OpenFOAM 1.7 installation on Redhat linux maxims OpenFOAM Installation 2 November 30, 2012 05:29
OpenFOAM 1.7 - openSUSE 11.3 - gcc 4.5.0 alberto OpenFOAM 12 July 28, 2010 11:59
OpenFoam 1.7 & Ubuntu 9.1 bendel_boy OpenFOAM 1 July 23, 2010 06:27
OpenFOAM 1.5 vs. 1.6 vs. 1.7 skarnani OpenFOAM 2 July 7, 2010 17:46
Compatibility from OF 1.6 to OF 1.7? Chrisi1984 OpenFOAM Installation 2 July 7, 2010 03:07


All times are GMT -4. The time now is 21:07.