基本方程组数值求解3.pptVIP

3．时间导数项的差分;时间导数项 在控制容积 上积分，得 （45） 这里假定 、 在整个控制容积中是均匀的，网格不随时间而变 在时间导数项的差分式子中，包含了两个时间层（在类时间坐标上是两个截面）的参数值 与此相对应，其它坐标方向的对流项、扩散项以及源项中的参数，在时间层上可以有不同的取法，因而可以得到显式、隐式和克兰克-尼科尔森（Crank-Nicholson）格式;1）显式差分格式 （46） 2）隐式差分格式 （47） 3）克兰克-尼科尔森格式 若对流项、扩散项和源项中的参数都取两个时间层上的算术平均值，可得差分方程 （48） 克兰克-尼科尔森格式的差分方程也需要与其它结点的差分方程联立求解。当时间步长较小时，其精度比隐式格式高。;4、源项的差分; 代入差分方程得 （51） 其中Sc和Sp可能也?? 的函数，因此式（50）也是拟线性关系。 第二种办法优于第一种办法 为了使拟线性代数方程组收敛，得到有意义的解，Sc和Sp必须满足 , （52）; 在满足条件（52）的基础上， 和 的确定有一定的任意性。一个方便的办法是用切线使其线性化，即令 （53） 这时 ， 按照切线法得到的 和 往往很复杂