[Nek5000-users] Traveling gravity waves

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Dec 6 06:45:30 CST 2012


Hi Arne,

Just one comment that ffx is a force per unit volume so you do not need to multiple it by density, i.e. param(1)

Aleks

----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Thursday, December 6, 2012 6:38:42 AM
Subject: [Nek5000-users] Traveling gravity waves

Hi Neks,

I am trying to set up a traveling gravity wave and started from the
std.wv example.

- In the useric I prescribed the the Initial velocity by using linear
wave theory (Airy)
 
http://en.wikipedia.org/wiki/Airy_wave_theory#Solution_for_a_progressive_monochromatic_wave
  wave amplitude
- the initial mesh was generate with genbox with periodic sides, bottom
wall and a free surface
- the surface tension was excluded by setting sigma & We = 0
- reverted density & viscosity to dimensional  variables (SI Units) in
the .rea to get correct velocities
- set the gravitational force in userf ffy = -9.81*param(1)

But instead of a traveling wave the system reverts to a sort of standing
wave immediately

Do I have set the initial velocity of the moving mesh as well? And if so
are wx,wy the correct variables to do so?
Did I apply the dimensions correctly?

I attached the .rea & .usr, any help with my Problem would be greatly
appreciated.

Thanks
Arne

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users


More information about the Nek5000-users mailing list