1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2008 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
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 the
13 Free Software Foundation; either version 2 of the License, or (at your
14 option) any later version.
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
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM; if not, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
25 \*---------------------------------------------------------------------------*/
27 #include "JohnsonJacksonFrictionalStress.H"
28 #include "addToRunTimeSelectionTable.H"
30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
34 defineTypeNameAndDebug(JohnsonJacksonFrictionalStress, 0);
36 addToRunTimeSelectionTable
38 frictionalStressModel,
39 JohnsonJacksonFrictionalStress,
45 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
47 Foam::JohnsonJacksonFrictionalStress::JohnsonJacksonFrictionalStress
49 const dictionary& dict
52 frictionalStressModel(dict)
56 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
58 Foam::JohnsonJacksonFrictionalStress::~JohnsonJacksonFrictionalStress()
62 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
64 Foam::tmp<Foam::volScalarField> Foam::JohnsonJacksonFrictionalStress::
67 const volScalarField& alpha,
68 const dimensionedScalar& alphaMinFriction,
69 const dimensionedScalar& alphaMax,
70 const dimensionedScalar& Fr,
71 const dimensionedScalar& eta,
72 const dimensionedScalar& p
77 Fr*pow(max(alpha - alphaMinFriction, scalar(0)), eta)
78 /pow(max(alphaMax - alpha, scalar(5.0e-2)), p);
82 Foam::tmp<Foam::volScalarField> Foam::JohnsonJacksonFrictionalStress::
83 frictionalPressurePrime
85 const volScalarField& alpha,
86 const dimensionedScalar& alphaMinFriction,
87 const dimensionedScalar& alphaMax,
88 const dimensionedScalar& Fr,
89 const dimensionedScalar& eta,
90 const dimensionedScalar& p
95 eta*pow(max(alpha - alphaMinFriction, scalar(0)), eta - 1.0)
96 *(alphaMax-alpha) + p*pow(max(alpha - alphaMinFriction, scalar(0)), eta)
97 )/pow(max(alphaMax - alpha, scalar(5.0e-2)), p + 1.0);
101 Foam::tmp<Foam::volScalarField> Foam::JohnsonJacksonFrictionalStress::muf
103 const volScalarField& alpha,
104 const dimensionedScalar& alphaMax,
105 const volScalarField& pf,
106 const volTensorField& D,
107 const dimensionedScalar& phi
110 return dimensionedScalar("0.5", dimTime, 0.5)*pf*sin(phi);
114 // ************************************************************************* //