VollerPrakash.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) 2017-2020 OpenCFD Ltd.
9-------------------------------------------------------------------------------
10License
11 This file is part of OpenFOAM.
12
13 OpenFOAM is free software: you can redistribute it and/or modify it
14 under the terms of the GNU General Public License as published by
15 the Free Software Foundation, either version 3 of the License, or
16 (at your option) any later version.
17
18 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 for more details.
22
23 You should have received a copy of the GNU General Public License
24 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26\*---------------------------------------------------------------------------*/
27
28#include "VollerPrakash.H"
29#include "phasePair.H"
31
32// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
33
34namespace Foam
35{
36namespace porousModels
37{
40 (
44 );
45}
46}
47
48
49// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
50
52(
53 const dictionary& dict,
54 const fvMesh& mesh
55)
56:
58 Cu_(dict.get<scalar>("Cu")),
59 solidPhase_(dict.get<word>("solidPhase"))
60{}
61
62
63// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
64
67{
68 const volScalarField& solidAlpha =
69 mesh_.lookupObject<volScalarField>(solidPhase_);
70
71 if (mesh_.time().outputTime() && debug)
72 {
74 (
75 "Su" + solidAlpha.name(),
76 Cu_*sqr(solidAlpha)/(pow3(1.0 - solidAlpha) + 1e-3)
77 );
78 Su.write();
79 }
80
81 return Cu_*sqr(solidAlpha)/(pow3(1.0 - solidAlpha) + 1e-3);
82}
83
84
85// ************************************************************************* //
Macros for easy insertion into run-time selection tables.
#define addToRunTimeSelectionTable(baseType, thisType, argNames)
Add to construction table with typeName as the key.
const word & name() const noexcept
Return the object name.
Definition: IOobjectI.H:65
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
Definition: dictionary.H:126
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:91
Porous model to apply a porous media on a solid phase.
Definition: VollerPrakash.H:97
virtual tmp< volScalarField > S() const
Aspect ratio.
Definition: VollerPrakash.C:66
A class for managing temporary objects.
Definition: tmp.H:65
A class for handling words, derived from Foam::string.
Definition: word.H:68
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
Definition: className.H:121
dynamicFvMesh & mesh
zeroField Su
Definition: alphaSuSp.H:1
Namespace for OpenFOAM.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
dimensionedScalar pow3(const dimensionedScalar &ds)
dictionary dict
volScalarField & e
Definition: createFields.H:11