其实也不难,你用grep大概能猜出个所以然来。
grep -rn "validate"
其他都是空函数只有一个函数有点意思
template<class BasicTurbulenceModel>
void Foam::ReynoldsStress<BasicTurbulencemModel>::validate()
{
correctNut();
}
然后我又grep 了correctNut()
常用的是kOmega 我看了下
template<class BasicTurbulenceModel>
void kOmega<BasicTurbulenceModel>::correctNut()
{
this->nut_ = k_/omega_;
correctBC;
BasicTurbulenceModel::correctNut():
}
我虽然没有一行行debug, 大概猜出来就是湍流初始化。最开始给定k omega 初始值计算出第一步的nut 初始值给k omega方程,即previous time step,欢迎指正。