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

Wave generation with 1.7.x and groovyBC

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 9, 2010, 08:57
Default Wave generation with 1.7.x and groovyBC
  #1
Senior Member
 
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 8
nuovodna is on a distinguished road
Hi, i setted up a 2d box with an inlet patch that generates a 1st order Stokes wave using groovyBC. The simulation stop to run (giving floating point exception) after 0.55 sec

This is the case

http://db.tt/Jiicy8A
nuovodna is offline   Reply With Quote

Old   September 9, 2010, 10:14
Default
  #2
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 7
keul is on a distinguished road
Hi Emanuele,

I runned your case with some modifications. In the alpha1 file, I changed the inlet condition by a zeroGradient and in the p_rgh file I changed the outlet condition by buyoantPressure. Now your case works fine.

Best regards
keul is offline   Reply With Quote

Old   September 9, 2010, 10:34
Default still crash
  #3
Senior Member
 
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 8
nuovodna is on a distinguished road
Hi mikael, thanks for your fast reply. I made the corrections that you suggested me but it still crashes. Now at t=6.8 sec

Regards
nuovodna is offline   Reply With Quote

Old   September 9, 2010, 10:42
Default
  #4
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 7
keul is on a distinguished road
When I runned your case, I put the tangentialVelocity in the U file at (0 0 0). Maybe your problem comes from that.

Rgrds
keul is offline   Reply With Quote

Old   September 9, 2010, 10:50
Default
  #5
Senior Member
 
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 8
nuovodna is on a distinguished road
where do you put tangentialVelocity uniform (0 0 0); ?? On which patch?

Thanks

E
nuovodna is offline   Reply With Quote

Old   September 9, 2010, 10:55
Default
  #6
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 7
keul is on a distinguished road
In your U file, patch atmosphere, I replace (1 0 0) by (0 0 0)
keul is offline   Reply With Quote

Old   September 9, 2010, 11:14
Default it works but not like i want
  #7
Senior Member
 
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 8
nuovodna is on a distinguished road
Changing the value of tangentialVelocity it doesnt crash anymore. But the free surface prediction is wrong.
This is the graph of computed free surface (red) against analytical solution (green) at t=20 sec
Attached Images
File Type: png contour20sec.png (7.6 KB, 47 views)

Last edited by nuovodna; September 9, 2010 at 11:23. Reason: wrong attached file
nuovodna is offline   Reply With Quote

Old   September 9, 2010, 11:32
Default
  #8
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 7
keul is on a distinguished road
Maybe you could try with the program I use. It generates regular waves with the deep water hypothesis. You could compile the waveBC folder with wmake. I put a floating body tutorial.
keul is offline   Reply With Quote

Old   September 9, 2010, 11:40
Default waveBC condition
  #9
Senior Member
 
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 8
nuovodna is on a distinguished road
I compiled waveBC yesterday but in the tutorial folder alpha1 file is missing. I grabbed it from shipHydrodynamicIG svn. Can you attach an example of waveBC on a case??

Thanks

Regards

E
nuovodna is offline   Reply With Quote

Old   September 9, 2010, 11:40
Default
  #10
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 7
keul is on a distinguished road
Sorry the file was to big.
Attached Files
File Type: gz prog.tar.gz (37.4 KB, 63 views)
keul is offline   Reply With Quote

Old   September 9, 2010, 11:44
Default
  #11
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 7
keul is on a distinguished road
For the alpha1 file, you could use the same file as your case with zeroGradient on all patches except eventually atmosphere with an inletOutlet condition.

Good luck
keul is offline   Reply With Quote

Old   September 9, 2010, 11:44
Default Allrun error
  #12
Senior Member
 
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 8
nuovodna is on a distinguished road
I launched Allrun but the alpha1.org~ file is missing


cp: cannot stat `0/alpha1.org~': No such file or directory

Edit: Ok i ll try. Thanks!

Edit: It works. Now i m going to make some test

Last edited by nuovodna; September 9, 2010 at 11:57. Reason: progress
nuovodna is offline   Reply With Quote

Old   September 9, 2010, 12:06
Default
  #13
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi

First of all, you have defined omega (w) differently in the alpha1.ORG and the U files in your input groovy file. This is most definitely one error in the input, which can result in the strange profiles in your post #7.

However, be aware, that linear wave theory will not yield the analytical solution in shallow water as the linear theory does not fulfill the boundary conditions for the non-linear wave problem at the inlet.

I am utterly surprised that the waveBC still exists out there, it has changed a bit since those years back, however I can still find my old comments in Danish I hope it works for you.

Best regards,

Niels
ngj is online now   Reply With Quote

Old   September 14, 2010, 05:02
Default Wave
  #14
Senior Member
 
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 8
nuovodna is on a distinguished road
Hi Niels, thanks for your answer. I re-setted the case using an water depth >= Wavelength/2 but the results are still wrong (i'm using groovyBC). Do you know another system available in OF to generate waves?
Thanks in advance

Regards

Emanuele
nuovodna is offline   Reply With Quote

Old   September 14, 2010, 05:31
Default
  #15
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,619
Rep Power: 25
ngj will become famous soon enoughngj will become famous soon enough
Hi Emanuele

Yes, I have a different system, however currently I am working on a journal article, and the source will not be released prior to publication/acceptance of the article. I made a presentation of the framework at the Gothenburg workshop, see [1].

I have been working on a specific validation case for a long time without succes, so I will not be able to give a time line for the release.

Good luck,

Niels

[1]: http://web.student.chalmers.se/group...SlidesOFW5.pdf
ngj is online now   Reply With Quote

Reply

Tags
groovybc, interfoam, wave, wavebc, waves

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
Non-normal inflow with groovyBC bjr OpenFOAM 9 January 17, 2011 11:35
OpenFOAM 1.6 and 1.7 with interFoam, groovyBC give different strange results Arnoldinho OpenFOAM 7 December 9, 2010 17:29
Wrong wave generation nuovodna OpenFOAM 2 August 9, 2010 10:53
Scale discrete inlet velocity profile with groovyBC cboss OpenFOAM 1 June 20, 2010 13:02
groovyBC tamsilian OpenFOAM Programming & Development 1 May 27, 2010 03:15


All times are GMT -4. The time now is 05:08.