新手关于su2的几个疑问

349516385

1、su2官网上说自己是：Unstructured meshes with median-dual control volumes (vertex-based) 这是啥意思？

2、su2空间积分时，里面有 Upwind_Residual 、 Viscous_Residual   、 Source_Residual  、  SetResidual_DualTime 等函数，这些函数都是计算残差？ 我看内部计算都是类似的，但是不知道差别

3、上面那些函数里对 LinSysRes  进行了操作，如：
    LinSysRes.SubtractBlock(iPoint, Res_Visc);
    LinSysRes.AddBlock(jPoint, Res_Visc);
  LinSysRes 是不是代表了系统的所有残差？LinSysSol 又是啥？

4、有句代码：  
    /*--- Compute and update residual ---*/
    numerics->ComputeResidual(Res_Visc, Jacobian_i, Jacobian_j, config);

   这个 Jacobian_i,和Jacobian_j 是啥？根据这个能算出残差？


5、最后是不是把 所有  Jacobian_i,和Jacobian_j  组成 Jacobian，Jacobian是隐式迭代用的？ 这个是怎么用的？

6、处理网格时貌似用了一个RCM算法对节点重编号，这是做什么？

也许有些问题很愚昧，但是确实不知道。。。谢啦

flowinflow

(1) vertex-based v.s. cell-based.
typical example is that fluent is cell-based and CFX is vertex-based.
for exact def. refer to blazek book.

(2) upwind is for convective flux, viscous is for viscous flux (physical diffusion)
source residual is for source term, such as transport eqn. e.g., SA turbulence

(3) res denotes residual, Sol denotes solution

(4) not sure

(5) this is the left hand side matrix for time stepping

(6) node reordering algorithm,  for improving the bandwitdh of the matrix during lin sol.

xici8080

1 关于median-dual control volumes ，请参考附件pdf 9-11页的内容
2-3 2楼朋友的已经说的很好了， 另外SetResidual_DualTime 是非定常计算时才出现的残值
4 是当前edge对应的i点和j点的Jacobian，SU2的隐式求解用的是Full Implicit，需要构造first order Jacobian，所以在计算残值的时候同时会计算对应的Jacobian
5 关于隐式求解和Jabobian构造可以参考Blazek书上时间推进那部分，再配合代码看会比较容易理解
6 RCM是非常经典的带宽缩减算法，一般是为了提高线性求解的效率