search

LHEEA / HOS-NWT

This project has migrated to GitLab: https://gitlab.com/lheea/HOS-NWT
https://gitlab.com/lheea/HOS-NWT
GNU General Public License v3.0
49 stars 25 forks source link

Usage of `dphis_add_rkramp` #19

Open rickyspaceguy opened 6 years ago

rickyspaceguy commented 6 years ago

What is the use of dphis_add_rkramp?

Line 151 resol_wmkr.f90 reads as

dphis_add_rkramp_loc(i3,i2) = - alphat*dpos1stdt(i3,i2) - alphatt*pos1st(i3,i2)

If this is velocity on the rhs then it should be- dphis_add_rkramp_loc(i3,i2) = - alpha*dpos1stdt(i3,i2) - alphat*pos1st(i3,i2)

Is it used anywhere?

gducrozet commented 6 years ago

This term is used for the time-marching of the wave maker boundary condition. Due to the presence of the time-ramp, there is a specific term to take into account, which is using this variable dphis_add_rkramp.

This is detailed in part II.2.4.2 of my PhD thesis (sorry, it is french). This term is used to compute the time integral in the Eq. (II.2.13).

As a consequence, this term is used in the subroutine solveHOS2, which is consequently at use in the time-integration runge_kutta.f90 and more specifically to compute the term a_add_mramp.

Hope this helps!

Guillaume.