|
[Sponsors] |
[swak4Foam] Wave generation with 1.7.x and groovyBC |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 18 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 17 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 18 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 17 ![]() |
When I runned your case, I put the tangentialVelocity in the U file at (0 0 0). Maybe your problem comes from that.
Rgrds |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 18 ![]() |
where do you put tangentialVelocity uniform (0 0 0); ?? On which patch?
Thanks E |
|
![]() |
![]() |
![]() |
![]() |
#6 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 17 ![]() |
In your U file, patch atmosphere, I replace (1 0 0) by (0 0 0)
|
|
![]() |
![]() |
![]() |
![]() |
#7 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 18 ![]() |
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 Last edited by nuovodna; September 9, 2010 at 11:23. Reason: wrong attached file |
|
![]() |
![]() |
![]() |
![]() |
#8 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 17 ![]() |
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.
|
|
![]() |
![]() |
![]() |
![]() |
#9 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 18 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#10 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 17 ![]() |
Sorry the file was to big.
|
|
![]() |
![]() |
![]() |
![]() |
#11 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 17 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#12 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 18 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#13 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 38 ![]() ![]() |
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 ![]() Best regards, Niels |
|
![]() |
![]() |
![]() |
![]() |
#14 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 18 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#15 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 38 ![]() ![]() |
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 |
|
![]() |
![]() |
![]() |
Tags |
groovybc, interfoam, wave, wavebc, waves |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] Change in alpha and U with groovyBC in twoPhaseEulerFoam | dani2702 | OpenFOAM Community Contributions | 0 | November 17, 2016 03:30 |
Wave generation inlet using interFoam | Katzenpaul | OpenFOAM Running, Solving & CFD | 2 | February 6, 2013 09:28 |
[swak4Foam] GroovyBC for 2D wave flume! | Hisham | OpenFOAM Community Contributions | 13 | January 20, 2012 05:04 |
[swak4Foam] groovyBC in 1.7.x ??? | fisch | OpenFOAM Community Contributions | 1 | June 7, 2011 08:58 |
Wrong wave generation | nuovodna | OpenFOAM | 2 | August 9, 2010 10:53 |