36 template<
class RdeltaTType,
class RhoType,
class SpType,
class SuType>
39 const RdeltaTType& rDeltaT,
61 mesh.
Vsc0()().field()*rho.oldTime().field()
62 *psi0*rDeltaT/mesh.
Vsc()().field()
65 )/(rho.field()*rDeltaT - Sp.field());
71 rho.oldTime().field()*psi0*rDeltaT
74 )/(rho.field()*rDeltaT - Sp.field());
81 template<
class RhoType>
93 template<
class RhoType,
class SpType,
class SuType>
105 if (fv::localEulerDdt::enabled(mesh))
107 const volScalarField& rDeltaT = fv::localEulerDdt::localRDeltaT(mesh);
118 template<
class RhoType,
class PsiMaxType,
class PsiMinType>
125 const PsiMaxType& psiMax,
126 const PsiMinType& psiMin
159 const PsiMaxType& psiMax,
160 const PsiMinType& psiMin
167 if (fv::localEulerDdt::enabled(mesh))
169 const volScalarField& rDeltaT = fv::localEulerDdt::localRDeltaT(mesh);
170 limit(rDeltaT, rho, psi, phi, phiPsi, Sp, Su, psiMax, psiMin,
false);
176 limit(rDeltaT, rho, psi, phi, phiPsi, Sp, Su, psiMax, psiMin,
false);
194 const RdeltaTType& rDeltaT,
201 const PsiMaxType& psiMax,
202 const PsiMinType& psiMin
212 const label nLimiterIter
217 const scalar smoothLimiter
222 const scalar extremaCoeff
235 const surfaceScalarField::Boundary& phiBDBf =
239 const surfaceScalarField::Boundary& phiCorrBf =
275 const label own = owner[facei];
276 const label nei = neighb[facei];
278 psiMaxn[own] =
max(psiMaxn[own], psiIf[nei]);
279 psiMinn[own] =
min(psiMinn[own], psiIf[nei]);
281 psiMaxn[nei] =
max(psiMaxn[nei], psiIf[own]);
282 psiMinn[nei] =
min(psiMinn[nei], psiIf[own]);
284 sumPhiBD[own] += phiBDIf[facei];
285 sumPhiBD[nei] -= phiBDIf[facei];
287 const scalar phiCorrf = phiCorrIf[facei];
291 sumPhip[own] += phiCorrf;
292 mSumPhim[nei] += phiCorrf;
296 mSumPhim[own] -= phiCorrf;
297 sumPhip[nei] -= phiCorrf;
315 const label pfCelli = pFaceCells[pFacei];
317 psiMaxn[pfCelli] =
max(psiMaxn[pfCelli], psiPNf[pFacei]);
318 psiMinn[pfCelli] =
min(psiMinn[pfCelli], psiPNf[pFacei]);
325 const label pfCelli = pFaceCells[pFacei];
327 psiMaxn[pfCelli] =
max(psiMaxn[pfCelli], psiPf[pFacei]);
328 psiMinn[pfCelli] =
min(psiMinn[pfCelli], psiPf[pFacei]);
335 const label pfCelli = pFaceCells[pFacei];
337 psiMaxn[pfCelli] =
max(psiMaxn[pfCelli], psiMax[pfCelli]);
338 psiMinn[pfCelli] =
min(psiMinn[pfCelli], psiMin[pfCelli]);
344 const label pfCelli = pFaceCells[pFacei];
346 sumPhiBD[pfCelli] += phiBDPf[pFacei];
348 const scalar phiCorrf = phiCorrPf[pFacei];
352 sumPhip[pfCelli] += phiCorrf;
356 mSumPhim[pfCelli] -= phiCorrf;
361 psiMaxn =
min(psiMaxn + extremaCoeff*(psiMax - psiMin), psiMax);
362 psiMinn =
max(psiMinn - extremaCoeff*(psiMax - psiMin), psiMin);
364 if (smoothLimiter > small)
367 min(smoothLimiter*psiIf + (1.0 - smoothLimiter)*psiMaxn, psiMax);
369 max(smoothLimiter*psiIf + (1.0 - smoothLimiter)*psiMinn, psiMin);
379 (rho.field()*rDeltaT - Sp.field())*psiMaxn
382 - (V0().field()*rDeltaT)*rho.oldTime().field()*psi0
389 - (rho.field()*rDeltaT - Sp.field())*psiMinn
391 + (V0().field()*rDeltaT)*rho.oldTime().field()*psi0
399 (rho.field()*rDeltaT - Sp.field())*psiMaxn
401 - (rho.oldTime().field()*rDeltaT)*psi0
409 - (rho.field()*rDeltaT - Sp.field())*psiMinn
410 + (rho.oldTime().field()*rDeltaT)*psi0
418 for (
int j=0; j<nLimiterIter; j++)
425 const label own = owner[facei];
426 const label nei = neighb[facei];
428 scalar lambdaPhiCorrf = lambdaIf[facei]*phiCorrIf[facei];
430 if (lambdaPhiCorrf > 0)
432 sumlPhip[own] += lambdaPhiCorrf;
433 mSumlPhim[nei] += lambdaPhiCorrf;
437 mSumlPhim[own] -= lambdaPhiCorrf;
438 sumlPhip[nei] -= lambdaPhiCorrf;
451 const label pfCelli = pFaceCells[pFacei];
452 const scalar lambdaPhiCorrf =
453 lambdaPf[pFacei]*phiCorrfPf[pFacei];
455 if (lambdaPhiCorrf > 0)
457 sumlPhip[pfCelli] += lambdaPhiCorrf;
461 mSumlPhim[pfCelli] -= lambdaPhiCorrf;
471 (sumlPhip[celli] + psiMaxn[celli])
472 /(mSumPhim[celli] + rootVSmall),
479 (mSumlPhim[celli] + psiMinn[celli])
480 /(sumPhip[celli] + rootVSmall),
490 if (phiCorrIf[facei] > 0)
492 lambdaIf[facei] =
min 495 min(lambdap[owner[facei]], lambdam[neighb[facei]])
500 lambdaIf[facei] =
min 503 min(lambdam[owner[facei]], lambdap[neighb[facei]])
525 const label pfCelli = pFaceCells[pFacei];
527 if (phiCorrfPf[pFacei] > 0)
530 min(lambdaPf[pFacei], lambdap[pfCelli]);
535 min(lambdaPf[pFacei], lambdam[pfCelli]);
557 const RdeltaTType& rDeltaT,
564 const PsiMaxType& psiMax,
565 const PsiMinType& psiMin,
566 const bool returnCorr
574 const surfaceScalarField::Boundary& phiPsiBf = phiPsi.
boundaryField();
582 phiBDPf = phiPsiBf[
patchi];
628 phiPsi = phiBD + lambda*phiCorr;
649 const PsiMaxType& psiMax,
650 const PsiMinType& psiMin,
656 if (fv::localEulerDdt::enabled(mesh))
658 const volScalarField& rDeltaT = fv::localEulerDdt::localRDeltaT(mesh);
659 limit(rDeltaT, rho, psi, phi, phiPsi, Sp, Su, psiMax, psiMin, rtnCorr);
664 limit(rDeltaT, rho, psi, phi, phiPsi, Sp, Su, psiMax, psiMin, rtnCorr);
669 template<
class SurfaceScalarFieldList>
682 const surfaceScalarField::Boundary& bfld =
683 phiPsiCorrs[0].boundaryField();
687 if (bfld[
patchi].coupled())
705 template<
class SurfaceScalarFieldList>
708 const SurfaceScalarFieldList& alphas,
709 SurfaceScalarFieldList& phiPsiCorrs,
725 limitSum(alphasInternal, phiPsiCorrsInternal, fixed);
728 const surfaceScalarField::Boundary& bfld =
729 phiPsiCorrs[0].boundaryField();
733 if (bfld[
patchi].coupled())
754 limitSum(alphasPatch, phiPsiCorrsPatch, fixed);
void limit(const RdeltaTType &rDeltaT, const RhoType &rho, const volScalarField &psi, const surfaceScalarField &phi, surfaceScalarField &phiPsi, const SpType &Sp, const SuType &Su, const PsiMaxType &psiMax, const PsiMinType &psiMin, const bool returnCorr)
Surface integrate surfaceField creating a volField. Surface sum a surfaceField creating a volField...
#define forAll(list, i)
Loop across all elements in list.
const GeometricField< Type, PatchField, GeoMesh > & oldTime() const
Return old time field.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
const word & name() const
Return name.
void surfaceIntegrate(Field< Type > &ivf, const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
bool moving() const
Is mesh moving.
virtual bool fixesValue() const
Return true if this patch field fixes a value.
A list of keyword definitions, which are a keyword followed by any number of values (e...
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
const Boundary & boundaryField() const
Return const-reference to the boundary field.
void size(const label)
Override size to be inconsistent with allocated storage.
Ostream & endl(Ostream &os)
Add newline and flush stream.
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
static word timeName(const scalar, const int precision=precision_)
Return time name of given scalar time.
const Time & time() const
Return the top-level database.
const dictionary & solverDict(const word &name) const
Return the solver controls dictionary for the given field.
const labelUList & neighbour() const
Internal face neighbour.
tmp< DimensionedField< scalar, volMesh > > Vsc() const
Return sub-cycle cell volumes.
void limitSum(UPtrList< scalarField > &phiPsiCorrs)
virtual bool coupled() const
Return true if this patch field is coupled.
A class representing the concept of a field of 0 used to avoid unnecessary manipulations for objects ...
void limiter(scalarField &allLambda, const RdeltaTType &rDeltaT, const RhoType &rho, const volScalarField &psi, const surfaceScalarField &phiBD, const surfaceScalarField &phiCorr, const SpType &Sp, const SuType &Su, const PsiMaxType &psiMax, const PsiMinType &psiMin)
virtual tmp< Field< Type > > patchNeighbourField() const
Return patchField on the opposite patch of a coupled patch.
scalar deltaTValue() const
Return time step value.
virtual bool coupled() const
Return true if this patch field is coupled.
A templated 1D list of pointers to objects of type <T>, where the size of the array is known and used...
A 1D vector of objects of type <T>, where the size of the vector is known and can be used for subscri...
tmp< DimensionedField< scalar, volMesh > > Vsc0() const
Return sub-cycl old-time cell volumes.
const Mesh & mesh() const
Return mesh.
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
const labelUList & owner() const
Internal face owner.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
T lookupOrDefault(const word &, const T &, bool recursive=false, bool patternMatch=true) const
Find and return a T,.
Boundary & boundaryFieldRef()
Return a reference to the boundary field.
Upwind differencing scheme class.
Mesh data needed to do the Finite Volume discretisation.
dimensionedScalar lambda(laminarTransport.lookup("lambda"))
void correctBoundaryConditions()
Correct boundary field.
MULES: Multidimensional universal limiter for explicit solution.
const volScalarField & psi
A class for managing temporary objects.
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
Specialization of GeometricField which holds slices of given complete fields in a form that they act ...
void explicitSolve(const RdeltaTType &rDeltaT, const RhoType &rho, volScalarField &psi, const surfaceScalarField &phiPsi, const SpType &Sp, const SuType &Su)
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
An abstract base class with a fat-interface to all derived classes covering all possible ways in whic...
const fvBoundaryMesh & boundary() const
Return reference to boundary mesh.