42
double REDUCE::GRAD(double gradSI)
43
{ return (eta*gradSI);}
44
// reduced body force
45
double REDUCE::F(double fSI)
46
{ return ((lamda*eta/gammaIA)*fSI);}
47
48
double REDUCE::KW(double kwSI)
49
50
51
52
53
54
55
56
57
58
59
{ return ((nu/(lamda*lamda))*kwSI);}
60
// reduced velocity
61
double REDUCE::V(double vSI)
62
{ return ((eta*nu/(lamda*gammaIA))*vSI);}
63
// reduced divergence
64
double REDUCE::DIV(double divSI)
65
{ return (eta*divSI);}
66
// reduced time
67
double REDUCE::T(double tSI)
68
{ return ((lamda*gammaIA/(eta*eta*nu))*tSI);}
69
// reduced volumetric heat capacity
70
double REDUCE::C(double cSI)
71
{ return ((lamda*theta0/gammaIA)*cSI);}
72
73
double REDUCE::KH(double khSI)
74
{ return ((nu*theta0/(gammaIA*gammaIA))*khSI);}
75
// reduced heat flux
76
double REDUCE::QH(double qhSI)
77
{ return (eta*nu/(gammaIA*gammaIA)*qhSI);}
78
79
//
80
//
81
//
82
//
SOIL FUNCTIONS
83
//
84
85
Tsoil::Tsoil(){
86
InitWATER(0.42,0.02,1.0e-8);
87
InitBC(0.36,2.6,11.196);
88
InitTherm(0.52,2.32,3.42);
89
InitStress();
90
}
91
92
void Tsoil::InitWATER(double wsat, double wd, double ksat){
93
// saturated volumetric water content
94
//
WSAT = 0.42;
95
WSAT = wsat;
19