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-2022 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 {
34 
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 // ************************************************************************* //
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
Generic GeometricBoundaryField class.
Generic GeometricField class.
const Boundary & boundaryField() const
Return const-reference to the boundary field.
Boundary & boundaryFieldRef()
Return a reference to the boundary field.
A HashTable with keys but without contents.
Definition: HashSet.H:62
An abstract base class with a fat-interface to all derived classes covering all possible ways in whic...
Definition: fvsPatchField.H:83
virtual bool coupled() const
Return true if this patch field is coupled.
virtual tmp< surfaceScalarField > interpolate(const volScalarField &vf) const
Return the face-interpolate of the given cell field.
Centred interpolation interpolation scheme class.
Definition: linear.H:53
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type and name.
const fvMesh & mesh() const
Return mesh reference.
A class for managing temporary objects.
Definition: tmp.H:55
T & ref() const
Return non-const reference or generate a fatal error.
Definition: tmpI.H:181
label patchi
Namespace for OpenFOAM.
dimensionedScalar sign(const dimensionedScalar &ds)
surfaceInterpolationScheme< scalar >::addMeshFluxConstructorToTable< interfaceCompressionNew > addinterfaceCompressionScalarMeshFluxConstructorToTable_
const wordHashSet compressionSchemes
defineTypeNameAndDebug(combustionModel, 0)