63 RootDeltaT(
sqrt(DeltaT)),
66 Alpha(UOdict.
lookup<scalar>(
"UOalpha")),
67 Sigma(UOdict.
lookup<scalar>(
"UOsigma")),
68 Kupper(UOdict.
lookup<scalar>(
"UOKupper")),
69 Klower(UOdict.
lookup<scalar>(
"UOKlower")),
74 scalar sqrKupper =
sqr(Kupper);
75 scalar sqrKlower =
sqr(Klower) + small;
80 if ((sqrK =
magSqr(K[i])) < sqrKupper && sqrK > sqrKlower)
82 UOfield[i] = Scale*Sigma*WeinerProcess();
104 scalar sqrKupper =
sqr(Kupper);
105 scalar sqrKlower =
sqr(Klower) + small;
110 if ((sqrK =
magSqr(K[i])) < sqrKupper && sqrK > sqrKlower)
114 (1.0 - Alpha*DeltaT)*UOfield[i]
115 + Scale*Sigma*WeinerProcess();
119 Info<<
" Number of forced K = " << count <<
nl;
#define forAll(list, i)
Loop across all elements in list.
FvWallInfoData< WallInfo, label > label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
scalar scalarNormal()
Advance the state and return a scalar sample from a normal.
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensionedSymmTensor sqr(const dimensionedVector &dv)
void size(const label)
Override size to be inconsistent with allocated storage.
const complexVectorField & newField()
dimensionedScalar sqrt(const dimensionedScalar &ds)
label count(const ListType &l, typename ListType::const_reference x)
Count the number of occurrences of a value in a list.
CGAL::Exact_predicates_exact_constructions_kernel K
Vector< complex > complexVector
complexVector obtained from generic Vector.
Calculate the wavenumber vector field corresponding to the space vector field of a finite volume mesh...
dimensioned< scalar > magSqr(const dimensioned< Type > &)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
static const direction dim
Dimensionality of space.
UOprocess(const Kmesh &kmesh, const scalar deltaT, const dictionary &)
Construct from wavenumber mesh and timestep.
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.