CFD计算过程发散诸多原因剖析.pptVIP

CFD计算过程发散诸多原因分析 　　cfd 计算是将描述物理问题的偏微分方程转化为代数方程组求解，从而得到离散空间上指定点上的值，而其他位置的值通过插值来完成。这本质是将非线性方程线性化并求解的过程，我是这样理解的，对于物理问题真解随时间的变化是一条曲线，而求解过程中不断的求解曲线的斜率，并实现时间的上的步进(可以理解为1阶导数)。如果在指定的时间步长内斜率变化很小，这种步进是可以满足方程要求的，并可以得到将来任意时刻的近似解。然而，如果变化比较大时，实际的变化过程将不能通过线性化方法来处理，这时候再使用这种方法就会出现问题。随着时间的推移，数值解越来越偏离真解，最终引起发散。对于我们求解的流动问题，有下面几种可能会引起计算过程的发散 　　1)时间步长过大 　　这点很容易理解，时间步长过长可能会造成时间段内，斜率变化较大，从而会造成数值解偏离真解。一般而言，对于显式方法要求在指定的时间步长内流体流动不能超过一个网格(库朗数限制)。 　　对于NS方程SIMPLE系列算法而言，需要通过压力来显式修正表面流率和速度，如果修正次数过少，通常会造成修正后的速度和真实速度偏离较大，下一时间步求解速度的时候会得到一个不真实的表面流率，从而影响下一个时刻的速度求解，进而影响压力，长期计算会影响造成误差积累设置发散。因此，压力修正次数要有一定的限制。一般而言对于2d问题，2~3次能够满足要求，而对