interfaceCompression.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2020-2021 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "interfaceCompression.H"
27 #include "linear.H"
28 
29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
30 
31 namespace Foam
32 {
33  defineTypeNameAndDebug(interfaceCompressionNew, 0);
34 
35  surfaceInterpolationScheme<scalar>::
36  addMeshFluxConstructorToTable<interfaceCompressionNew>
38 
40  {
41  "interfaceCompression",
42  "noInterfaceCompression",
43  "PLIC",
44  "PLICU",
45  "MPLIC",
46  "MPLICU"
47  };
48 }
49 
50 
53 (
54  const volScalarField& vf
55 ) const
56 {
57  const surfaceScalarField& nHatf =
59  (
60  "nHatf"
61  );
62 
63  const surfaceScalarField vff
64  (
66  );
67 
69  (
70  cAlpha_*sign(phi_)*vff*(1 - vff)*nHatf/mesh().magSf()
71  );
72 
74 
75  // Do not compress interface at non-coupled boundary faces
76  // (inlets, outlets etc.)
78  {
79  fvsPatchScalarField& vfcp = vfcBf[patchi];
80 
81  if (!vfcp.coupled())
82  {
83  vfcp == 0;
84  }
85  }
86 
87  tmp<surfaceScalarField> tvff(tScheme_().interpolate(vf) + vfc);
88  tvff.ref().maxMin(0, 1);
89 
90  return tvff;
91 }
92 
93 
94 // ************************************************************************* //
dimensionedScalar sign(const dimensionedScalar &ds)
A HashTable with keys but without contents.
Definition: HashSet.H:59
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
virtual tmp< surfaceScalarField > interpolate(const volScalarField &vf) const
Return the face-interpolate of the given cell field.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
Centred interpolation interpolation scheme class.
Definition: linear.H:50
bool interpolate(const vector &p1, const vector &p2, const vector &o, vector &n, scalar l)
Definition: curveTools.C:75
T & ref() const
Return non-const reference or generate a fatal error.
Definition: tmpI.H:181
const wordHashSet compressionSchemes
void maxMin(const dimensioned< Type > &minDt, const dimensioned< Type > &maxDt)
fvMesh & mesh
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
virtual bool coupled() const
Return true if this patch field is coupled.
defineTypeNameAndDebug(combustionModel, 0)
label patchi
Boundary & boundaryFieldRef()
Return a reference to the boundary field.
A class for managing temporary objects.
Definition: PtrList.H:53
surfaceInterpolationScheme< scalar >::addMeshFluxConstructorToTable< interfaceCompressionNew > addinterfaceCompressionScalarMeshFluxConstructorToTable_
An abstract base class with a fat-interface to all derived classes covering all possible ways in whic...
Definition: fvsPatchField.H:65
Namespace for OpenFOAM.