如何提高数值模拟计算性能（2）.PDFVIP

第二十一期：如何提高数值模拟计算性能（2 ） 如何提高数值模拟计算性能（2 ） 在第十一期讲义《如何提高数值模拟计算性能（1）》中讲解了如何读计算记录信 息，数值模拟计算的基本原理，如何确定最小时间步，以及如何检查时间步过小等内容。 本期内容将继续讲解控制方程、内迭代、外迭代、物质平衡误差、矩阵求解失败以及如 何读取牛顿迭代信息等。 在进行数值模拟运算的时候，有时会出现较差的计算性能，例如时间步过小、每 个时间步牛顿迭代过多、较多的截断（CUT ）、温度或压力值超过限值、物质平衡误差 很大、矩阵求解失败、警告（warning ）、错误（Error ）过多等问题，造成这些问题的 原因很多，我们要学会如何从众多的原因中找出最关键的原因，并及时调整，使模拟能 够继续重新顺利运行。 CMG软件的数值模拟计算是根据如下的流程进行： 1 第二十一期：如何提高数值模拟计算性能（2 ） 一、CMG中的控制方程 对于有某一个网格（单元体）的质量（和能量）守恒方程为： 组分i 的质量守恒方程为： 组分i 的能量守恒方程为： 总的未知数数目 = 总的方程数 未知数联立同步求解 需要用牛顿迭代对方程组进行线性化 2 第二十一期：如何提高数值模拟计算性能（2 ） 二、牛顿迭代（外迭代） 牛顿迭代，是将非线性偏微分方程组线性化的一种方法。牛顿迭代是牛顿在17世纪 提出的一种在实数域和复数域上近似求解方程的方法。解非线性方程f(x)=0 的牛顿法是 把非线性方程线性化的一种近似方法。把f(x)在x0点附近展开成泰勒级数 f(x) = f(x0)+(x －x0)f(x0)+(x －x0)^2*f(x0)/2! +… 取其线性部分，作为非线性方程f(x) = 0 的近似方程， 即泰勒展开的前两项，则有f(x0)+f(x0)(x －x0)-f(x)=0 设f(x0)≠0则其解为x1=x0 － f(x0)/f(x0) 这样，得到牛顿法的一个迭代序列：x(n+1)=x(n) －f(x(n))/f(x(n)) 。 通过牛顿迭代，将控制方程中的偏微分方程组进行线性化处理，最后形成线性的雅 克比（Jacobian ）矩阵，也叫方程组，即各个方程关于未知量均为一次的方程组。 3 第二十一期：如何提高数值模拟计算性能（2 ） 三、求解雅克比矩阵（外迭代） Δ Δ Δ Δ Δ Δ ⎛ f o 1 f o 1 f o 1 f o 1 f o 1 f o 1 ⎞ Δ Δ Δ Δ Δ0 Δ0 0 ⎛Δp ⎞ ⎛f ⎞ p s T p s T o1 o1 ⎜Δ o 1 Δ o 1 Δ 1 Δ o 2 Δ o 2 Δ 2 ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ f w 1 f w 1 f w 1 f w 1 f w 1 f w 1 ⎟ Δ