Formatting
[foam-extend-3.2.git] / src / transportModels / incompressible / viscosityModels / BirdCarreau / BirdCarreau.C
blobb1544a0418072083bc24711dedfe251ee4e223dd
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | foam-extend: Open Source CFD
4    \\    /   O peration     | Version:     3.2
5     \\  /    A nd           | Web:         http://www.foam-extend.org
6      \\/     M anipulation  | For copyright notice see file Copyright
7 -------------------------------------------------------------------------------
8 License
9     This file is part of foam-extend.
11     foam-extend 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 3 of the License, or (at your
14     option) any later version.
16     foam-extend is distributed in the hope that it will be useful, but
17     WITHOUT ANY WARRANTY; without even the implied warranty of
18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
19     General Public License for more details.
21     You should have received a copy of the GNU General Public License
22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
24 \*---------------------------------------------------------------------------*/
26 #include "BirdCarreau.H"
27 #include "addToRunTimeSelectionTable.H"
28 #include "surfaceFields.H"
30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
32 namespace Foam
34 namespace viscosityModels
36     defineTypeNameAndDebug(BirdCarreau, 0);
37     addToRunTimeSelectionTable
38     (
39         viscosityModel,
40         BirdCarreau,
41         dictionary
42     );
47 // * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * * //
49 Foam::tmp<Foam::volScalarField>
50 Foam::viscosityModels::BirdCarreau::calcNu() const
52     return
53         nuInf_
54       + (nu0_ - nuInf_)
55        *pow(scalar(1) + sqr(k_*strainRate()), (n_ - 1.0)/2.0);
59 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
61 Foam::viscosityModels::BirdCarreau::BirdCarreau
63     const word& name,
64     const dictionary& viscosityProperties,
65     const volVectorField& U,
66     const surfaceScalarField& phi
69     viscosityModel(name, viscosityProperties, U, phi),
70     BirdCarreauCoeffs_(viscosityProperties.subDict(typeName + "Coeffs")),
71     nu0_(BirdCarreauCoeffs_.lookup("nu0")),
72     nuInf_(BirdCarreauCoeffs_.lookup("nuInf")),
73     k_(BirdCarreauCoeffs_.lookup("k")),
74     n_(BirdCarreauCoeffs_.lookup("n")),
75     nu_
76     (
77         IOobject
78         (
79             "nu",
80             U_.time().timeName(),
81             U_.db(),
82             IOobject::NO_READ,
83             IOobject::AUTO_WRITE
84         ),
85         calcNu()
86     )
90 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
92 bool Foam::viscosityModels::BirdCarreau::read
94     const dictionary& viscosityProperties
97     viscosityModel::read(viscosityProperties);
99     BirdCarreauCoeffs_ = viscosityProperties.subDict(typeName + "Coeffs");
101     BirdCarreauCoeffs_.lookup("nu0") >> nu0_;
102     BirdCarreauCoeffs_.lookup("nuInf") >> nuInf_;
103     BirdCarreauCoeffs_.lookup("k") >> k_;
104     BirdCarreauCoeffs_.lookup("n") >> n_;
106     return true;
110 // ************************************************************************* //