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

Problems with wave BC without a relaxation zone

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

Like Tree2Likes
  • 1 Post By Hisham
  • 1 Post By Hisham

Reply
 
LinkBack Thread Tools Display Modes
Old   October 7, 2013, 07:28
Default Problems with wave BC without a relaxation zone
  #1
Senior Member
 
Hisham's Avatar
 
Hisham El Safti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 248
Blog Entries: 10
Rep Power: 9
Hisham is on a distinguished road
Dear Foamers,

I have been trying to model waves in a 3-D tank using the boundary conditions from waves2Foam (i.e. waveAlpha and waveVelocity). Because I do not need to absorb waves at the inlet, I have not configured a relaxation zone at the inlet (left). The use of the BCs only is OK for 2D flumes. However, when I run a 3D geometry (with a symmetry plane) the case crashes from an artificial high velocity near the water surface (air phase) and very close to the BC. The error is gone if a small relaxation zone is used at the inlet near the BC (1 m. long).

The error also happens for:
- Wall BC instead of the symmetry BC and
- Use of codedBC for wave BC instead of waves2Foam BCs (not a waves2Foam problem)

I have tried to neglect the air phase convection as well as relax the air velocities at each time step (set to zero). This seems to take the case for more steps further but still fails in the end.
The mesh is structured (check mesh OK). I also tried different schemes/solvers (but not all).

I am using OF-2.1.1 and thinking of trying a different version but have doubts if the version is the problem!

I would appreciate help with this problem!

Best regards,
Hisham
alia likes this.
Hisham is offline   Reply With Quote

Old   October 17, 2013, 14:39
Default
  #2
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,685
Rep Power: 26
ngj will become famous soon enoughngj will become famous soon enough
Hi Hisham,

I have first seen your post now.

I have had the same problems many time, and I am sorry to tell you that I do not have a solution for it. Therefore, I have not developed/tested for this type of use in waves2Foam for many years. The only suggestion I have is to look in the literature for absorbing boundary conditions with reflection compensation. As this works without relaxation zones, a solution might be available.

If you find a solution, please tell me, such that I can improve waves2Foam.

Kind regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   October 17, 2013, 15:06
Default
  #3
Senior Member
 
Hisham's Avatar
 
Hisham El Safti
Join Date: Apr 2011
Location: Braunschweig, Germany
Posts: 248
Blog Entries: 10
Rep Power: 9
Hisham is on a distinguished road
Quote:
Originally Posted by ngj View Post
Hi Hisham,

I have first seen your post now.

I have had the same problems many time, and I am sorry to tell you that I do not have a solution for it. Therefore, I have not developed/tested for this type of use in waves2Foam for many years. The only suggestion I have is to look in the literature for absorbing boundary conditions with reflection compensation. As this works without relaxation zones, a solution might be available.

If you find a solution, please tell me, such that I can improve waves2Foam.

Kind regards,

Niels

Hi Niels,

I really appreciate your response. I have doubts that an absorbing boundary can solve that because I tried to increase the domain length to avoid reflection, but the case crashes too soon for any reflection to take place (at least that is what I experience or understand/think you mean by reflection). I can totally understand your idea about reflection because of solvers that use active absorption and are reported to work! Nevertheless!!

I wonder if it will work if I were to use two relaxation zones on both sides like the 3DTank tutorial (but for a small length). I will try the tutorial without the inlet zone.

As for my problem I am developing a "hack" for it (while admiring beauty of waves2Foam ) as a new waveTheory class that I will send to you after it's published. However, it is still a relaxation approach.

Best regards,
Hisham
alia likes this.
Hisham is offline   Reply With Quote

Old   October 17, 2013, 18:57
Default
  #4
ngj
Senior Member
 
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Deltares, Delft, The Netherlands
Posts: 1,685
Rep Power: 26
ngj will become famous soon enoughngj will become famous soon enough
Cool, looking forward to your contribution.

I hope that your project are still on track

Kind regards,

Niels
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.
ngj is offline   Reply With Quote

Old   January 21, 2015, 09:23
Default
  #5
Senior Member
 
Join Date: Jul 2011
Posts: 105
Rep Power: 6
haze_1986 is on a distinguished road
Hi Hisham, I am having a similar problem. Using groovyBC or IHFoam to generate waves (just Stokes I into a 3D harbour), the simulations did not crash but sampling the inlet point I am never able to obtain the specified Stokes I waves. I have not tried waves2foam but do you think inlet relaxation zones could be the answer to my question?
haze_1986 is offline   Reply With Quote

Reply

Tags
waves, waves2foam

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
Cht tutorial in 15 braennstroem OpenFOAM Running, Solving & CFD 197 June 10, 2015 03:02
Mesh conversion problem (fluent3DMeshToFoam) Aadhavan OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 1 December 12, 2012 11:38
Reference on analytical solutions to wave/ Euler problems Hooman Main CFD Forum 0 August 8, 2012 05:12
[ICEM] Export ICEM mesh to Gambit / Fluent romekr ANSYS Meshing & Geometry 1 November 26, 2011 13:11
how to pass bc from one zone to second zone Maria FLUENT 2 March 16, 2004 00:12


All times are GMT -4. The time now is 14:12.