看到官方的求解器，根本摸不清楚边界条件和有限体积是如何交互的。有人能解释一下这和我上述思路的差别和联系吗？
粗略看了下，这个求解器完全没用openfoam的思路来，采用的他自己的思路，就像你说的，他在开发这个求解器的时候，可能是有其他的代码（比如同组的Fortran）转过来的。并且完全是面向过程的思想。
如果你要用OpenFOAM求解的思路，还是看OpenFOAM的求解器比较好他这个就是把求解的过程自己实现了，然而没有使用OpenFOAM自带的方法。
他的main函数已经写的很清楚了：
int main(int argc, char *argv[]) { #include "setRootCase.H" #include "createTime.H" #include "createMesh.H" #include "setInputValues.H" #include "createFields.H" #include "readFluxScheme.H" /// Time step loop /// Posts the non-blocking send/recv of fields long int iter = 0; while( runTime.loop() ) { /// 构造通量 #include "constructFaceFlux.H" /// 有限体积离散 #include "sumFlux.H" /// 边界通量修正 #include "boundaryFlux.H" Info << "Iteration = " << ++iter << " "; /// 矩阵计算 #include "stateUpdateLTS.H" Info << " Max residue = " << rhoResidMax << endl; /// 输出结果 runTime.write(); } return 0; }