00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043 #ifndef LESmodel_H
00044 #define LESmodel_H
00045
00046 #include "LESdelta.H"
00047 #include "fvm.H"
00048 #include "fvc.H"
00049 #include "fvMatrices.H"
00050 #include "incompressible/transportModel/transportModel.H"
00051 #include "wallFvPatch.H"
00052 #include "bound.H"
00053 #include "autoPtr.H"
00054 #include "runTimeSelectionTables.H"
00055
00056
00057
00058 namespace Foam
00059 {
00060
00061
00062
00063
00064
00065 class LESmodel
00066 :
00067 public IOdictionary
00068 {
00069
00070 protected:
00071
00072
00073
00074 const Time& runTime_;
00075 const fvMesh& mesh_;
00076
00077
00078 private:
00079
00080
00081
00082 const volVectorField& U_;
00083 const surfaceScalarField& phi_;
00084
00085 transportModel& transport_;
00086
00087 dictionary LESmodelProperties_;
00088
00089 dimensionedScalar k0_;
00090
00091 autoPtr<LESdelta> delta_;
00092
00093
00094
00095
00096
00097 LESmodel(const LESmodel&);
00098 LESmodel& operator=(const LESmodel&);
00099
00100
00101 public:
00102
00103
00104 TypeName("LESmodel");
00105
00106
00107
00108
00109 declareRunTimeSelectionTable
00110 (
00111 autoPtr,
00112 LESmodel,
00113 dictionary,
00114 (
00115 const volVectorField& U,
00116 const surfaceScalarField& phi,
00117 transportModel& transport
00118 ),
00119 (U, phi, transport)
00120 );
00121
00122
00123
00124
00125
00126 LESmodel
00127 (
00128 const word& type,
00129 const volVectorField& U,
00130 const surfaceScalarField& phi,
00131 transportModel& transport
00132 );
00133
00134
00135
00136
00137
00138 static autoPtr<LESmodel> New
00139 (
00140 const volVectorField& U,
00141 const surfaceScalarField& phi,
00142 transportModel& transport
00143 );
00144
00145
00146
00147
00148 virtual ~LESmodel()
00149 {}
00150
00151
00152
00153
00154
00155
00156
00157 inline const volVectorField& U() const
00158 {
00159 return U_;
00160 }
00161
00162
00163 inline const surfaceScalarField& phi() const
00164 {
00165 return phi_;
00166 }
00167
00168
00169 inline transportModel& transport() const
00170 {
00171 return transport_;
00172 }
00173
00174
00175
00176 inline const dictionary& LESmodelProperties()
00177 {
00178 return LESmodelProperties_;
00179 }
00180
00181
00182 inline const volScalarField& delta() const
00183 {
00184 return delta_();
00185 }
00186
00187
00188 inline const volScalarField& nu() const
00189 {
00190 return transport_.nu();
00191 }
00192
00193
00194 const dimensionedScalar& k0() const
00195 {
00196 return k0_;
00197 }
00198
00199
00200 dimensionedScalar& k0()
00201 {
00202 return k0_;
00203 }
00204
00205
00206
00207 virtual tmp<volTensorField> B() const = 0;
00208
00209
00210 virtual tmp<volScalarField> k() const = 0;
00211
00212
00213 virtual tmp<volScalarField> epsilon() const = 0;
00214
00215
00216 virtual tmp<volScalarField> nuSgs() const = 0;
00217
00218
00219 virtual tmp<volScalarField> nuEff() const
00220 {
00221 return tmp<volScalarField>
00222 (
00223 new volScalarField("nuEff", nuSgs() + nu())
00224 );
00225 }
00226
00227
00228
00229 virtual tmp<fvVectorMatrix> divB(volVectorField& U) const = 0;
00230
00231
00232 virtual void correct(const tmp<volTensorField>& gradU);
00233
00234
00235
00236
00237 void correct();
00238
00239
00240 virtual bool read() = 0;
00241 };
00242
00243
00244
00245
00246 }
00247
00248
00249
00250 #endif
00251
00252