Frank.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | www.openfoam.com
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8  Copyright (C) 2014-2018 OpenFOAM Foundation
9  Copyright (C) 2020 OpenCFD Ltd.
10 -------------------------------------------------------------------------------
11 License
12  This file is part of OpenFOAM.
13 
14  OpenFOAM is free software: you can redistribute it and/or modify it
15  under the terms of the GNU General Public License as published by
16  the Free Software Foundation, either version 3 of the License, or
17  (at your option) any later version.
18 
19  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
20  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
21  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
22  for more details.
23 
24  You should have received a copy of the GNU General Public License
25  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
26 
27 \*---------------------------------------------------------------------------*/
28 
29 #include "Frank.H"
30 #include "phasePair.H"
32 
33 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
34 
35 namespace Foam
36 {
37 namespace wallLubricationModels
38 {
39  defineTypeNameAndDebug(Frank, 0);
41  (
42  wallLubricationModel,
43  Frank,
44  dictionary
45  );
46 }
47 }
48 
49 
50 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
51 
53 (
54  const dictionary& dict,
55  const phasePair& pair
56 )
57 :
59  Cwd_("Cwd", dimless, dict),
60  Cwc_("Cwc", dimless, dict),
61  p_(dict.get<scalar>("p"))
62 {}
63 
64 
65 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
66 
68 {}
69 
70 
71 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
72 
74 {
75  volVectorField Ur(pair_.Ur());
76 
77  const volVectorField& n(nWall());
78  const volScalarField& y(yWall());
79 
80  volScalarField Eo(pair_.Eo());
81  volScalarField yTilde(y/(Cwc_*pair_.dispersed().d()));
82 
83  return zeroGradWalls
84  (
85  (
86  pos0(Eo - 1)*neg(Eo - 5)*exp(-0.933*Eo + 0.179)
87  + pos0(Eo - 5)*neg(Eo - 33)*(0.00599*Eo - 0.0187)
88  + pos0(Eo - 33)*0.179
89  )
90  *max
91  (
93  (1 - yTilde)/(Cwd_*y*pow(yTilde, p_ - 1))
94  )
95  *pair_.continuous().rho()
96  *magSqr(Ur - (Ur & n)*n)
97  *n
98  );
99 }
100 
101 
102 // ************************************************************************* //
Foam::phasePair
Description for mass transfer between a pair of phases. The direction of the mass transfer is from th...
Definition: phasePair.H:53
Foam::dimLength
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
Definition: dimensionSets.H:52
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:61
Foam::wallLubricationModels::Frank::~Frank
virtual ~Frank()
Destructor.
Definition: Frank.C:67
Foam::pos0
dimensionedScalar pos0(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:188
Foam::exp
dimensionedScalar exp(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:261
Foam::magSqr
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
n
label n
Definition: TABSMDCalcMethod2.H:31
Foam::dimensionedScalar
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Definition: dimensionedScalarFwd.H:42
Foam::pow
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Definition: dimensionedScalar.C:75
Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:47
dict
dictionary dict
Definition: searchingEngine.H:14
Foam::dictionary
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:123
addToRunTimeSelectionTable.H
Macros for easy insertion into run-time selection tables.
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.C:33
Foam::wallLubricationModels::Frank::Frank
Frank(const dictionary &dict, const phasePair &pair)
Construct from components.
Definition: Frank.C:53
Foam::wallLubricationModels::defineTypeNameAndDebug
defineTypeNameAndDebug(Antal, 0)
Foam::wallLubricationModels::addToRunTimeSelectionTable
addToRunTimeSelectionTable(wallLubricationModel, Antal, dictionary)
Foam::GeometricField< vector, fvPatchField, volMesh >
Foam::neg
dimensionedScalar neg(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:199
Foam::wallLubricationModels::Frank::Fi
tmp< volVectorField > Fi() const
Return phase-intensive wall lubrication force.
Definition: Frank.C:73
Frank.H
Foam::dimless
const dimensionSet dimless
Dimensionless.
Definition: dimensionSets.C:189
Foam::wallLubricationModel
Definition: wallLubricationModel.H:56
y
scalar y
Definition: LISASMDCalcMethod1.H:14