[Nek5000-users] Traveling gravity waves

[nek5000-users at lists.mcs.anl.gov]

thanks Aleks. I tried that as well, but the velocity field still is
completely different to the initial Condition after one step (dt = 0.00025).

Arne

On 06.12.2012 13:45, nek5000-users at lists.mcs.anl.gov wrote:
> 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
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users