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 forceModel 00027 00028 Description 00029 00030 \*---------------------------------------------------------------------------*/ 00031 00032 #ifndef forceModel_H 00033 #define forceModel_H 00034 00035 #include "fvCFD.H" 00036 #include "cfdemCloud.H" 00037 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00038 00039 namespace Foam 00040 { 00041 00042 /*---------------------------------------------------------------------------*\ 00043 Class forceModel Declaration 00044 \*---------------------------------------------------------------------------*/ 00045 00046 class forceModel 00047 { 00048 00049 protected: 00050 00051 // Protected data 00052 const dictionary& dict_; 00053 00054 cfdemCloud& particleCloud_; 00055 00056 Switch treatExplicit_; // marker to treat force in implicit way (otherwise explicit) 00057 00058 Switch treatDEM_; // marker to use the force only on DEM side 00059 00060 mutable volVectorField impParticleForces_; // sum of implicit particle forces [N] 00061 00062 mutable volVectorField expParticleForces_; // sum of explicit particle forces [N] 00063 00064 bool coupleForce_; 00065 00066 const word modelType_; 00067 00068 public: 00069 00070 //- Runtime type information 00071 TypeName("forceModel"); 00072 00073 // Declare runtime constructor selection table 00074 00075 declareRunTimeSelectionTable 00076 ( 00077 autoPtr, 00078 forceModel, 00079 dictionary, 00080 ( 00081 const dictionary& dict, 00082 cfdemCloud& sm 00083 ), 00084 (dict,sm) 00085 ); 00086 00087 00088 // Constructors 00089 00090 //- Construct from components 00091 forceModel 00092 ( 00093 const dictionary& dict, 00094 cfdemCloud& sm 00095 ); 00096 00097 00098 // Destructor 00099 00100 virtual ~forceModel(); 00101 00102 00103 // Selector 00104 00105 static autoPtr<forceModel> New 00106 ( 00107 const dictionary& dict, 00108 cfdemCloud& sm, 00109 word forceType 00110 ); 00111 00112 00113 // Member Functions 00114 virtual void setForce(double** const& mask,double**&,double**&,double**&) const = 0; 00115 00116 //tmp<volScalarField> provideScalarField(); 00117 00118 virtual void manipulateScalarField(volScalarField&) const; 00119 00120 // Access 00121 word modelType(){ return modelType_; }; 00122 00123 inline volVectorField& impParticleForces() const { return impParticleForces_;}; 00124 00125 inline volVectorField& expParticleForces() const { return expParticleForces_;}; 00126 00127 inline const bool& coupleForce() const { return coupleForce_;}; 00128 }; 00129 00130 00131 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00132 00133 } // End namespace Foam 00134 00135 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00136 00137 #endif 00138 00139 // ************************************************************************* //