OpenFOAM logo
Open Source CFD Toolkit

CrossPowerLaw.H

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2005 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software; you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by the
00013     Free Software Foundation; either version 2 of the License, or (at your
00014     option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM; if not, write to the Free Software Foundation,
00023     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
00024 
00025 Class
00026     CrossPowerLaw
00027 
00028 Description
00029      Cross power law non-Newtonian viscosity model.
00030 
00031 SourceFiles
00032     CrossPowerLaw.C
00033 
00034 \*---------------------------------------------------------------------------*/
00035 
00036 #ifndef CrossPowerLaw_H
00037 #define CrossPowerLaw_H
00038 
00039 #include "viscosityModel.H"
00040 #include "dimensionedScalar.H"
00041 #include "volFields.H"
00042 
00043 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00044 
00045 namespace Foam
00046 {
00047 namespace viscosityModels
00048 {
00049 
00050 /*---------------------------------------------------------------------------*\
00051                            Class CrossPowerLaw Declaration
00052 \*---------------------------------------------------------------------------*/
00053 
00054 class CrossPowerLaw
00055 :
00056     public viscosityModel
00057 {
00058     // Private data
00059 
00060         dictionary CrossPowerLawCoeffs_;
00061 
00062         dimensionedScalar nu0_;
00063         dimensionedScalar nuInf_;
00064         dimensionedScalar m_;
00065         dimensionedScalar n_;
00066 
00067         volScalarField nu_;
00068 
00069     // Private Member Functions
00070 
00071         //- Calculate and return the laminar viscosity
00072         tmp<volScalarField> calcNu() const;
00073 
00074 
00075 public:
00076 
00077     //- Runtime type information
00078     TypeName("CrossPowerLaw");
00079 
00080 
00081     // Constructors
00082 
00083         //- Construct from components
00084         CrossPowerLaw
00085         (
00086             const word& name,
00087             const dictionary& viscosityProperties,
00088             const volVectorField& U,
00089             const surfaceScalarField& phi
00090         );
00091 
00092 
00093     // Destructor
00094 
00095         ~CrossPowerLaw()
00096         {}
00097 
00098 
00099     // Member Functions
00100 
00101         //- Return the laminar viscosity
00102         tmp<volScalarField> nu() const
00103         {
00104             return nu_;
00105         }
00106 
00107         //- Correct the laminar viscosity
00108         void correct()
00109         {
00110             nu_ = calcNu();
00111         }
00112 
00113         //- Read transportProperties dictionary
00114         bool read(const dictionary& viscosityProperties);
00115 };
00116 
00117 
00118 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00119 
00120 } // End namespace viscosityModels
00121 } // End namespace Foam
00122 
00123 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00124 
00125 #endif
00126 
00127 // ************************************************************************* //
For further information go to www.openfoam.org