68 const scalar kd =
k()*depth(), ka =
k()*amplitude(t);
70 const scalar
T = deep() ? 1 :
tanh(kd);
72 const scalar B22 = (3/
sqr(T) - 1)/T/4;
81 + (1/
k())*
sqr(ka)*B22*
cos(2*angle(t, u, x));
92 const scalar kd =
k()*depth(), ka =
k()*amplitude(t);
94 const scalar A22ByA11 = deep() ? 0 : 0.375/
pow3(
sinh(kd));
98 const scalar A11 = 1/
sinh(kd);
99 Info<<
"A22 = " << A22ByA11*A11 <<
endl;
104 + celerity()*
sqr(ka)*A22ByA11*vi(2, t, u, xz);
dimensionedScalar tanh(const dimensionedScalar &ds)
A list of keyword definitions, which are a keyword followed by any number of values (e...
defineTypeNameAndDebug(Airy, 0)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Stokes2(const objectRegistry &db, const dictionary &dict)
Construct from a database and a dictionary.
Ostream & endl(Ostream &os)
Add newline and flush stream.
label k
Boltzmann constant.
Macros for easy insertion into run-time selection tables.
virtual ~Stokes2()
Destructor.
dimensionedScalar cos(const dimensionedScalar &ds)
virtual tmp< vector2DField > velocity(const scalar t, const scalar u, const vector2DField &xz) const
Get the wave velocity at a given time, mean velocity and local.
virtual tmp< scalarField > elevation(const scalar t, const scalar u, const scalarField &x) const
Get the wave elevation at a given time, mean velocity and local.
addToRunTimeSelectionTable(waveModel, Airy, objectRegistry)
virtual tmp< vector2DField > velocity(const scalar t, const scalar u, const vector2DField &xz) const
Get the wave velocity at a given time, mean velocity and local.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensionedScalar sinh(const dimensionedScalar &ds)
A class for managing temporary objects.
Registry of regIOobjects.
virtual tmp< scalarField > elevation(const scalar t, const scalar u, const scalarField &x) const
Get the wave elevation at a given time, mean velocity and local.