motionDirectionalDiffusivity.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) 2011-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 
27 #include "surfaceInterpolate.H"
30 
31 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
32 
33 namespace Foam
34 {
36 
38  (
41  Istream
42  );
43 }
44 
45 
46 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
47 
49 (
50  const fvMesh& mesh,
51  Istream& mdData
52 )
53 :
54  motionDiffusivity(mesh),
55  diffusivityVector_(mdData)
56 {}
57 
58 
59 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
60 
62 {}
63 
64 
65 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
66 
69 {
70  if (mesh().foundObject<volVectorField>("cellMotionU"))
71  {
72  const volVectorField& cellMotionU =
73  mesh().lookupObject<volVectorField>("cellMotionU");
74 
76  (
77  IOobject
78  (
79  "D",
80  mesh().time().name(),
81  mesh()
82  ),
83  diffusivityVector_.y()*vector::one
84  + (diffusivityVector_.x() - diffusivityVector_.y())*cellMotionU
85  /(mag(cellMotionU) + dimensionedScalar(dimVelocity, small)),
86  zeroGradientFvPatchVectorField::typeName
87  );
88  D.correctBoundaryConditions();
89 
90  const surfaceVectorField n(mesh().Sf()/mesh().magSf());
91 
93  (
94  "faceDiffusivity",
96  );
97  }
98  else
99  {
101  (
103  (
104  IOobject
105  (
106  "faceDiffusivity",
107  mesh().time().name(),
108  mesh(),
111  ),
112  mesh(),
114  )
115  );
116  }
117 }
118 
119 
120 // ************************************************************************* //
static const Foam::dimensionedScalar D("D", Foam::dimTemperature, 257.14)
label n
Macros for easy insertion into run-time selection tables.
Generic GeometricField class.
static tmp< GeometricField< Type, PatchField, GeoMesh > > New(const word &name, const Internal &, const PtrList< PatchField< Type >> &, const HashPtrTable< Source > &=HashPtrTable< Source >())
Return a temporary field constructed from name,.
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:99
An Istream is an abstract base class for all input systems (streams, files, token lists etc)....
Definition: Istream.H:60
static const Form one
Definition: VectorSpace.H:119
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:99
Abstract base class for cell-centre mesh motion diffusivity.
MotionDirectional finite volume mesh motion diffusivity.
motionDirectionalDiffusivity(const fvMesh &mesh, Istream &mdData)
Construct for the given fvMesh and data Istream.
virtual tmp< surfaceScalarField > operator()() const
Return diffusivity field.
A class for managing temporary objects.
Definition: tmp.H:55
static tmp< SurfaceField< Type > > interpolate(const VolField< Type > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
Namespace for OpenFOAM.
addToRunTimeSelectionTable(polyPatch, mergedCyclicPolyPatch, word)
word name(const bool)
Return a word representation of a bool.
Definition: boolIO.C:39
const dimensionSet dimless
dimensioned< scalar > mag(const dimensioned< Type > &)
defineTypeNameAndDebug(combustionModel, 0)
dimensioned< Type > cmptMultiply(const dimensioned< Type > &, const dimensioned< Type > &)
const dimensionSet dimVelocity
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.