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-2019 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 {
35  defineTypeNameAndDebug(motionDirectionalDiffusivity, 0);
36 
38  (
39  motionDiffusivity,
40  motionDirectionalDiffusivity,
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().timeName(),
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",
95  (n & cmptMultiply(fvc::interpolate(D), n))
96  );
97  }
98  else
99  {
101  (
103  (
104  IOobject
105  (
106  "faceDiffusivity",
107  mesh().time().timeName(),
108  mesh(),
111  ),
112  mesh(),
114  )
115  );
116  }
117 }
118 
119 
120 // ************************************************************************* //
virtual tmp< surfaceScalarField > operator()() const
Return diffusivity field.
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Definition: Istream.H:57
const Type & lookupObject(const word &name) const
Lookup and return the object of the given Type.
Macros for easy insertion into run-time selection tables.
static tmp< GeometricField< scalar, fvsPatchField, surfaceMesh > > New(const word &name, const Mesh &, const dimensionSet &, const word &patchFieldType=fvsPatchField< scalar >::calculatedType())
Return a temporary field constructed from name, mesh, dimensionSet.
dynamicFvMesh & mesh
word timeName
Definition: getTimeIndex.H:3
dimensioned< Type > cmptMultiply(const dimensioned< Type > &, const dimensioned< Type > &)
addToRunTimeSelectionTable(ensightPart, ensightPartCells, istream)
Abstract base class for cell-centre mesh motion diffusivity.
defineTypeNameAndDebug(combustionModel, 0)
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
motionDirectionalDiffusivity(const fvMesh &mesh, Istream &mdData)
Construct for the given fvMesh and data Istream.
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
Definition: dimensionSets.H:47
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
dimensioned< scalar > mag(const dimensioned< Type > &)
label n
mesh Sf()
A class for managing temporary objects.
Definition: PtrList.H:53
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:92
Namespace for OpenFOAM.
const dimensionSet dimVelocity