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., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00024 00025 Class 00026 Koch, Hill 00027 00028 Description 00029 Koch, Hill drag law 00030 based on Koch Hill 2001,"Inertial effects in suspensions and porous-media flows", Annual Review of fluid mechanics. 00031 00032 including interpolation of the velocity to the exact position 00033 00034 \*---------------------------------------------------------------------------*/ 00035 00036 #ifndef KochHillDrag2_H 00037 #define KochHillDrag2_H 00038 00039 #include "forceModel.H" 00040 #include "interpolationCellPoint.H" 00041 00042 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00043 00044 namespace Foam 00045 { 00046 00047 /*---------------------------------------------------------------------------*\ 00048 Class KochHillDrag2 Declaration 00049 \*---------------------------------------------------------------------------*/ 00050 00051 class KochHillDrag2 00052 : 00053 public forceModel 00054 { 00055 private: 00056 dictionary propsDict_; 00057 00058 bool verbous_; 00059 00060 word velFieldName_; 00061 00062 const volVectorField& U_; 00063 00064 word densityFieldName_; 00065 00066 const volScalarField& rho_; 00067 00068 word voidfractionFieldName_; 00069 00070 const volScalarField& voidfraction_; 00071 00072 bool interpolation_; // use interpolated field values 00073 00074 public: 00075 00076 //- Runtime type information 00077 TypeName("KochHillDrag2"); 00078 00079 00080 // Constructors 00081 00082 //- Construct from components 00083 KochHillDrag2 00084 ( 00085 const dictionary& dict, 00086 cfdemCloud& sm 00087 ); 00088 00089 // Destructor 00090 00091 ~KochHillDrag2(); 00092 00093 00094 // Member Functions 00095 void setForce 00096 ( 00097 double** const& mask, 00098 double**&, 00099 double**&, 00100 double**& 00101 ) const; 00102 }; 00103 00104 00105 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00106 00107 } // End namespace Foam 00108 00109 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00110 00111 #endif 00112 00113 // ************************************************************************* //