|
[Sponsors] |
[swak4Foam] Wave generation with 1.7.x and groovyBC |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 9, 2010, 09:57 |
Wave generation with 1.7.x and groovyBC
|
#1 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
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 |
|
September 9, 2010, 11:14 |
|
#2 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 16 |
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 |
|
September 9, 2010, 11:34 |
still crash
|
#3 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
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 |
|
September 9, 2010, 11:42 |
|
#4 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 16 |
When I runned your case, I put the tangentialVelocity in the U file at (0 0 0). Maybe your problem comes from that.
Rgrds |
|
September 9, 2010, 11:50 |
|
#5 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
where do you put tangentialVelocity uniform (0 0 0); ?? On which patch?
Thanks E |
|
September 9, 2010, 11:55 |
|
#6 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 16 |
In your U file, patch atmosphere, I replace (1 0 0) by (0 0 0)
|
|
September 9, 2010, 12:14 |
it works but not like i want
|
#7 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
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 12:23. Reason: wrong attached file |
|
September 9, 2010, 12:32 |
|
#8 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 16 |
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.
|
|
September 9, 2010, 12:40 |
waveBC condition
|
#9 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
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 |
|
September 9, 2010, 12:40 |
|
#10 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 16 |
Sorry the file was to big.
|
|
September 9, 2010, 12:44 |
|
#11 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 16 |
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 |
|
September 9, 2010, 12:44 |
Allrun error
|
#12 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
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 12:57. Reason: progress |
|
September 9, 2010, 13:06 |
|
#13 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 37 |
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 |
|
September 14, 2010, 06:02 |
Wave
|
#14 |
Senior Member
Emanuele
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
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 |
|
September 14, 2010, 06:31 |
|
#15 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 37 |
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 |
|
|
Similar Threads | ||||
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 04:30 |
Wave generation inlet using interFoam | Katzenpaul | OpenFOAM Running, Solving & CFD | 2 | February 6, 2013 10:28 |
[swak4Foam] GroovyBC for 2D wave flume! | Hisham | OpenFOAM Community Contributions | 13 | January 20, 2012 06:04 |
[swak4Foam] groovyBC in 1.7.x ??? | fisch | OpenFOAM Community Contributions | 1 | June 7, 2011 09:58 |
Wrong wave generation | nuovodna | OpenFOAM | 2 | August 9, 2010 11:53 |