00001
00002
00003
00004
00005
00006
00007 #ifndef AMROC_PROBLEM_H
00008 #define AMROC_PROBLEM_H
00009
00010 #define DIM 3
00011
00012 #include "LBMProblem.h"
00013 #include "LBMD3Q19Thermal.h"
00014
00015 typedef double DataType;
00016 typedef LBMD3Q19Thermal<DataType> LBMType;
00017
00018 #define OWN_LBMSCHEME
00019
00020
00021 #include "LBMStdProblem.h"
00022
00023
00024
00025
00026
00027 class LBMSpecific : public LBMType {
00028 typedef LBMType base;
00029 public:
00030 LBMSpecific() : base(), Tmp(1.), deltaTp(1.), gp(1.), betap(1.), Diffp(1.),GasRho(1.), Gasnu(1.), Gasp0(1.), Qvol(0.) {}
00031
00032 virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {
00033 base::register_at(Ctrl,prefix);
00034 RegisterAt(base::LocCtrl,"Tm",Tmp);
00035 RegisterAt(base::LocCtrl,"deltaT",deltaTp);
00036 RegisterAt(base::LocCtrl,"g",gp);
00037 RegisterAt(base::LocCtrl,"beta",betap);
00038 RegisterAt(base::LocCtrl,"Diff",Diffp);
00039 RegisterAt(base::LocCtrl,"Gas_rho",GasRho);
00040 RegisterAt(base::LocCtrl,"Gas_nu",Gasnu);
00041 RegisterAt(base::LocCtrl,"Gas_p0",Gasp0);
00042 RegisterAt(base::LocCtrl,"Qvol", Qvol);
00043 }
00044
00045
00046 virtual void SetupData(GridHierarchy* gh, const int& ghosts) {
00047 base::SetupData(gh,ghosts);
00048
00049 BBox wb = base::GH().wholebbox();
00050 DCoords lc = base::GH().worldCoords(wb.lower(),wb.stepsize());
00051 DCoords uc = base::GH().worldCoords(wb.upper()+wb.stepsize(),wb.stepsize());
00052
00053
00054 DataType Hp=uc(1)-lc(1);
00055
00056 DataType Speed = base::SpeedUp();
00057
00058 DataType viscp = Gasnu;
00059 DataType Pr = viscp/Diffp;
00060
00061 DataType tScale = base::TimeScale();
00062 DataType lScale = base::LengthScale();
00063 DataType THp = Tmp+deltaTp/DataType(2.);
00064 DataType TCp = Tmp-deltaTp/DataType(2.);
00065
00066
00067 DataType GasCsp = base::LatticeSpeedOfSound();
00068
00069 base::SetGas(Gasp0,GasRho,viscp,GasCsp);
00070 base::SetThermalGas(TCp,THp,Diffp,gp,betap);
00071
00072
00073 std::cout << "D3Q19Thermal: Gas_rho=" << base::GasDensity() << " Gas_Cs=" << base::GasSpeedofSound()
00074 << " Gas_nu=" << base::GasViscosity() << " Gas_nut=" << base::GasViscosityT() << std::endl;
00075 std::cout << "\nControl Parameters : " << std::endl;
00076 std::cout << "Pr = " << Pr << std::endl;
00077 std::cout << "\nTemperature : " << std::endl;
00078 std::cout << "TH = " << THp << " TC = " << TCp << std::endl;
00079 std::cout << "\nFluid values : " << std::endl;
00080 std::cout << "viscp = " << viscp << " Diffp = " << Diffp << " gp = " << gp << " betap = " << betap << std::endl;
00081 std::cout << "\nGeometry : " << std::endl;
00082 std::cout << "Channel Width = " << Hp << std::endl;
00083 std::cout << "\nLattice Units : " << std::endl;
00084 std::cout << "(Amroc) Omega = " << base::Omega(base::TimeScale()) << " (Amroc) OmegaT = " << base::OmegaT(base::TimeScale()) << std::endl;
00085 std::cout << "(Amroc) Gas_nu = " << base::LatticeViscosity(base::Omega(base::TimeScale())) << " (Amroc) Gas_nuT = " << base::LatticeViscosityT(base::OmegaT(base::TimeScale())) << std::endl;
00086
00087
00088
00089 std::cout << "\nScaling Factors: " << std::endl;
00090 std::cout << "LengthScale = " << base::LengthScale() << std::endl;
00091 std::cout << "TimeScale = " << base::TimeScale() << std::endl;
00092 std::cout << "VelocityScale = " << base::VelocityScale() << std::endl;
00093 std::cout << "SpeedUp = " << base::SpeedUp() << std::endl;
00094
00095 }
00096
00097 inline virtual void Collision(MicroType &f, const DataType dt) const {
00098
00099 MacroType q = MacroVariables(f);
00100 MicroType feq = Equilibrium(q);
00101
00102 DataType force = DataType(0.5)*dt/S0*gbeta*(q(4)-Tref);
00103 DataType omega = Omega(dt), omegaT = OmegaT(dt);
00104
00105 BBox wb = GH().wholebbox();
00106 DCoords lc = GH().worldCoords(wb.lower(),wb.stepsize());
00107 DCoords uc = GH().worldCoords(wb.upper()+wb.stepsize(),wb.stepsize());
00108
00109
00110 DataType width=uc(0)-lc(0);
00111 DataType bodyforce = 8.*LatticeViscosity(Omega(dt))*Qvol/width/width;
00112 bodyforce = bodyforce/10.;
00113
00114
00115 for (int i=0; i<=25; i++){
00116 if(i<19){
00117 f(i)=(DataType(1.)-omega)*f(i) + omega*feq(i);
00118 if(i==3) f(3) += (force+bodyforce);
00119 if(i==4) f(4) -= (force+bodyforce);
00120 if(i==7) f(7) += bodyforce;
00121 if(i==8) f(8) -= bodyforce;
00122 if(i==9) f(9) -= bodyforce;
00123 if(i==10) f(10) += bodyforce;
00124 if(i==11) f(11) += bodyforce;
00125 if(i==12) f(12) -= bodyforce;
00126 if(i==13) f(13) += bodyforce;
00127 if(i==14) f(14) -= bodyforce;
00128 }
00129 else f(i)=(DataType(1.)-omegaT)*f(i) + omegaT*feq(i);
00130 }
00131
00132 }
00133
00134 protected:
00135 DataType Tmp, deltaTp, gp, betap, Diffp;
00136 DataType GasRho, Gasnu, Gasp0, Qvol;
00137 };
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147
00148
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237
00238
00239
00240
00241
00242
00243
00244
00245
00246
00247
00248
00249
00250
00251
00252
00253
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264
00265
00266
00267
00268
00269
00270
00271
00272
00273
00274
00275
00276
00277
00278
00279
00280
00281
00282
00283
00284
00285
00286
00287
00288
00289
00290
00291
00292
00293
00294
00295
00296
00297 #endif