66 scalar phaseTot = Kx*
x + Ky*
y - omega*t +
phase;
68 return H*0.5*
cos(phaseTot);
83 for (
int i=1; i<=100; i++)
106 scalar
k =
sqrt(Kx*Kx + Ky*Ky);
107 scalar phaseTot = Kx*
x + Ky*
y - omega*t +
phase;
111 scalar v = u*
sin(waveAngle_);
112 u *=
cos(waveAngle_);
127 const scalar waveKx = waveK*
cos(waveAngle_);
128 const scalar waveKy = waveK*
sin(waveAngle_);
145 level[paddlei] = waterDepthRef_ + tCoeff*eta;
159 const scalar waveKx = waveK*
cos(waveAngle_);
160 const scalar waveKy = waveK*
sin(waveAngle_);
170 setPaddlePropeties(level, facei, fraction, z);
174 const label paddlei = faceToPaddle_[facei];
190 U_[facei] = fraction*
Uf*tCoeff;
219 if (regularWaveModel::readDict(overrideDict))
221 waveLength_ = waveLength(waterDepthRef_, wavePeriod_);
234 os <<
" Wave type: " << waveType() <<
nl;
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
InfoProxy< ensightCells > info() const
Return info proxy.
Mesh data needed to do the Finite Volume discretisation.
Single incompressible phase derived from the phase-fraction. Used as part of the multiPhaseMixture fo...
A patch is a list of labels that address the faces in the global face list.
Base class for waveModels.
virtual void setVelocity(const scalar t, const scalar tCoeff, const scalarField &level)
Calculate the wave model velocity.
virtual void setLevel(const scalar t, const scalar tCoeff, scalarField &level) const
Set the water level.
virtual scalar waveLength(const scalar h, const scalar T) const
Return the wavelength.
virtual vector UfBase(const scalar H, const scalar h, const scalar Kx, const scalar x, const scalar Ky, const scalar y, const scalar omega, const scalar t, const scalar phase, const scalar z) const
Wave velocity.
virtual bool readDict(const dictionary &overrideDict)
Read from dictionary.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
volScalarField H(IOobject("H", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimLength, Zero))
autoPtr< surfaceVectorField > Uf
OBJstream os(runTime.globalPath()/outputName)
constexpr scalar pi(M_PI)
constexpr scalar twoPi(2 *M_PI)
Different types of constants.
dimensionedScalar cosh(const dimensionedScalar &ds)
dimensionedScalar sin(const dimensionedScalar &ds)
dimensionedScalar tanh(const dimensionedScalar &ds)
dimensionedScalar sinh(const dimensionedScalar &ds)
void readFields(const typename GeoFieldType::Mesh &mesh, const IOobjectList &objects, const wordHashSet &selectedFields, LIFOStack< regIOobject * > &storedObjects)
Read the selected GeometricFields of the templated type.
dimensionedScalar sqrt(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedScalar cos(const dimensionedScalar &ds)
constexpr char nl
The newline '\n' character (0x0a)
#define forAll(list, i)
Loop across all elements in list.
const vector L(dict.get< vector >("L"))