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 #ifndef CoEulerDdtScheme_H
00043 #define CoEulerDdtScheme_H
00044
00045 #include "ddtScheme.H"
00046
00047
00048
00049 namespace Foam
00050 {
00051
00052
00053
00054 namespace fv
00055 {
00056
00057
00058
00059
00060
00061 template<class Type>
00062 class CoEulerDdtScheme
00063 :
00064 public fv::ddtScheme<Type>
00065 {
00066
00067
00068 scalar maxCo_;
00069
00070
00071
00072
00073
00074 CoEulerDdtScheme(const CoEulerDdtScheme&);
00075
00076
00077 void operator=(const CoEulerDdtScheme&);
00078
00079
00080 tmp<volScalarField> CorDeltaT() const;
00081
00082
00083 tmp<surfaceScalarField> CofrDeltaT() const;
00084
00085
00086 public:
00087
00088
00089 TypeName("CoEuler");
00090
00091
00092
00093
00094
00095 CoEulerDdtScheme(const fvMesh& mesh, Istream& is)
00096 :
00097 ddtScheme<Type>(mesh, is),
00098 maxCo_(readScalar(is))
00099 {}
00100
00101
00102
00103
00104
00105 const fvMesh& mesh() const
00106 {
00107 return fv::ddtScheme<Type>::mesh();
00108 }
00109
00110 tmp<GeometricField<Type, fvPatchField, volMesh> > fvcDdt
00111 (
00112 const dimensioned<Type>&
00113 );
00114
00115 tmp<GeometricField<Type, fvPatchField, volMesh> > fvcDdt
00116 (
00117 const GeometricField<Type, fvPatchField, volMesh>&
00118 );
00119
00120 tmp<GeometricField<Type, fvPatchField, volMesh> > fvcDdt
00121 (
00122 const dimensionedScalar&,
00123 const GeometricField<Type, fvPatchField, volMesh>&
00124 );
00125
00126 tmp<GeometricField<Type, fvPatchField, volMesh> > fvcDdt
00127 (
00128 const volScalarField&,
00129 const GeometricField<Type, fvPatchField, volMesh>&
00130 );
00131
00132 tmp<fvMatrix<Type> > fvmDdt
00133 (
00134 GeometricField<Type, fvPatchField, volMesh>&
00135 );
00136
00137 tmp<fvMatrix<Type> > fvmDdt
00138 (
00139 const dimensionedScalar&,
00140 GeometricField<Type, fvPatchField, volMesh>&
00141 );
00142
00143 tmp<fvMatrix<Type> > fvmDdt
00144 (
00145 const volScalarField&,
00146 GeometricField<Type, fvPatchField, volMesh>&
00147 );
00148
00149 typedef typename ddtScheme<Type>::fluxFieldType fluxFieldType;
00150
00151 tmp<fluxFieldType> fvcDdtPhiCorr
00152 (
00153 const volScalarField& rA,
00154 const GeometricField<Type, fvPatchField, volMesh>& U,
00155 const fluxFieldType& phi
00156 );
00157
00158 tmp<fluxFieldType> fvcDdtPhiCorr
00159 (
00160 const volScalarField& rA,
00161 const volScalarField& rho,
00162 const GeometricField<Type, fvPatchField, volMesh>& U,
00163 const fluxFieldType& phi
00164 );
00165
00166 tmp<surfaceScalarField> meshPhi
00167 (
00168 const GeometricField<Type, fvPatchField, volMesh>&
00169 );
00170 };
00171
00172
00173 template<>
00174 tmp<surfaceScalarField> CoEulerDdtScheme<scalar>::fvcDdtPhiCorr
00175 (
00176 const volScalarField& rA,
00177 const volScalarField& U,
00178 const surfaceScalarField& phi
00179 );
00180
00181
00182 template<>
00183 tmp<surfaceScalarField> CoEulerDdtScheme<scalar>::fvcDdtPhiCorr
00184 (
00185 const volScalarField& rA,
00186 const volScalarField& rho,
00187 const volScalarField& U,
00188 const surfaceScalarField& phi
00189 );
00190
00191
00192
00193
00194 }
00195
00196
00197
00198 }
00199
00200
00201
00202 #ifdef NoRepository
00203 # include "CoEulerDdtScheme.C"
00204 #endif
00205
00206
00207
00208 #endif
00209
00210