Forward compatibility: flex
[foam-extend-3.2.git] / src / dynamicMesh / meshMotion / fvMotionSolver / motionDiffusivity / inverseDistance / inverseDistanceDiffusivity.C
blob55061b04607440e544606cda9f710781917ef20b
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | foam-extend: Open Source CFD
4    \\    /   O peration     | Version:     3.2
5     \\  /    A nd           | Web:         http://www.foam-extend.org
6      \\/     M anipulation  | For copyright notice see file Copyright
7 -------------------------------------------------------------------------------
8 License
9     This file is part of foam-extend.
11     foam-extend is free software: you can redistribute it and/or modify it
12     under the terms of the GNU General Public License as published by the
13     Free Software Foundation, either version 3 of the License, or (at your
14     option) any later version.
16     foam-extend is distributed in the hope that it will be useful, but
17     WITHOUT ANY WARRANTY; without even the implied warranty of
18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
19     General Public License for more details.
21     You should have received a copy of the GNU General Public License
22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
24 \*---------------------------------------------------------------------------*/
26 #include "inverseDistanceDiffusivity.H"
27 #include "addToRunTimeSelectionTable.H"
28 #include "patchWave.H"
29 #include "HashSet.H"
30 #include "surfaceInterpolate.H"
31 #include "zeroGradientFvPatchFields.H"
33 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
35 namespace Foam
37     defineTypeNameAndDebug(inverseDistanceDiffusivity, 0);
39     addToRunTimeSelectionTable
40     (
41         motionDiffusivity,
42         inverseDistanceDiffusivity,
43         Istream
44     );
48 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
50 Foam::inverseDistanceDiffusivity::inverseDistanceDiffusivity
52     const fvMotionSolver& mSolver,
53     Istream& mdData
56     uniformDiffusivity(mSolver, mdData),
57     patchNames_(mdData)
59     correct();
63 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
65 Foam::inverseDistanceDiffusivity::~inverseDistanceDiffusivity()
69 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
71 Foam::tmp<Foam::scalarField> Foam::inverseDistanceDiffusivity::y() const
73     const polyMesh& mesh = mSolver().mesh();
75     labelHashSet patchSet(mesh.boundaryMesh().patchSet(patchNames_));
77     if (patchSet.size())
78     {
79         return tmp<scalarField>
80         (
81             new scalarField(patchWave(mesh, patchSet, false).distance())
82         );
83     }
84     else
85     {
86         return tmp<scalarField>(new scalarField(mesh.nCells(), 1.0));
87     }
91 void Foam::inverseDistanceDiffusivity::correct()
93     const fvMesh& mesh = mSolver().mesh();
95     volScalarField y_
96     (
97         IOobject
98         (
99             "y",
100             mesh.time().timeName(),
101             mesh
102         ),
103         mesh,
104         dimless,
105         zeroGradientFvPatchScalarField::typeName
106     );
107     y_.internalField() = y();
108     y_.correctBoundaryConditions();
110     faceDiffusivity_ = 1.0/fvc::interpolate(y_);
114 // ************************************************************************* //