[Nek5000-users] Traveling gravity waves
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Fri Dec 7 10:29:49 CST 2012
I understand H.O.T. do not make big difference in the displacement amplitude, but how can you have dispersion without having at least two components such as cos(theta) and cos(2*theta)
you can estimate how much run time you need to see significant dispersion based on the phase speed of each component.
I also had a lot of problems in the past with pseudospectral codes where there are small inconsistencies/errors in the initial conditions.
Ammar
On Dec 6, 2012, at 4:05 PM, nek5000-users at lists.mcs.anl.gov wrote:
> Hi Ammar,
>
> as I am using a single deepwater wave with amplitude a = 0.01 and wavenumber k = 1 the steepness := k*a = 0.01 is quite low.
> Higher Order terms from Stokes Theory http://en.wikipedia.org/wiki/Stokes_wave can be neglected as they are small
> For example the third order solution for deep water http://en.wikipedia.org/wiki/Stokes_wave#Third-order_Stokes_wave_on_deep_water doesn't differ from the Airy Theory so the velocities are correct to order 3 in the pertubation series with respect to the steepness k*a
> Only the surface elevation differs in the order O( a²k +a³k) which should be neglectable as well.
>
> I want to simulate a simple wave and than goto steeper waves with higher nonlinearities (Stokes) in order to create a test case to see what resolution and parameters are needed in NEK5000
> and then step up to more complex sea states. I already simulated this wave with a pseudo-spectral method which stability strongly depends on accurate initial Conditions and could simulate many periods without problems.
>
> All the best, Arne
>
> PS: I'll let you know if 3rd Stokes made a difference
>
> Zitat von nek5000-users at lists.mcs.anl.gov:
>
>> I am not an expert on surface waves but just wondering if
>> you are initializing with one wave component such as this?
>>
>> My understanding is you should have multiple components to see something interesting. This can be simply achieved by modulating the above
>> form by a spatial localization function such as exp(-(x-x0)2) or add other components such as a1 cos(kx...)+ a2 cos(2kx...)+...etc
>>
>> Ammar
>>
>>
>>
>>
>>
>>
>> On Thu, Dec 6, 2012 at 6:38 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>> 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
>>
>>
>>
>>
>>
>> --
>> http://aelghany.googlepages.com/home
>>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20121207/729929f6/attachment.html>
More information about the Nek5000-users
mailing list