[zhxter]

Hi, foamers.
I am now using the waves2foam and have some confusion about the parameter Tsoft.

For the stokesfirst wave theory (see stokesfirst.C)


// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //


scalar stokesFirst::factor(const scalar& time) const
{
scalar factor(1.0);
if (Tsoft_ > 0.0)
{
factor = Foam::sin(2*PI_/(4.0*Tsoft_)*Foam::min(Tsoft_, time));
}

return factor;
}


scalar stokesFirst::eta
(
const point& x,
const scalar& time
) const
{
scalar eta = H_/2.0*Foam::cos(omega_*time - (k_ & x) + phi_)
*factor(time) + seaLevel_;
return eta;
}

Why the defintion of wave elevation eta is multiplied by the factor(time) ?? I can understand it when the value of factor(time) equals to 1. But if Tsoft>0, then the value of factor(time) is no longer equal to 1. In such cases, the value of eta is not correct??

[ngj]

Good evening,

It has been some time since my last reply, so I will try to wrap everything up:

@Xiantoo: Try plotting the function factor(const scalar& time) and you will see that it goes from zero to one over the interval Tsoft.

Kind regards,

Niels