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

3．时间导数项的差分 在时间坐标上，一个显著的特点是只有一阶导数项，没有高阶导数项，考虑定常二维边界层方程 其中在 方向也只有一阶导数，可称 为类时间坐标 此处的讨论同样适用于有类时间坐标的问题。 时间导数项 在控制容积 上积分，得 （45） 这里假定 、 在整个控制容积中是均匀的，网格不随时间而变 在时间导数项的差分式子中，包含了两个时间层（在类时间坐标上是两个截面）的参数值 与此相对应，其它坐标方向的对流项、扩散项以及源项中的参数，在时间层上可以有不同的取法，因而可以得到显式、隐式和克兰克-尼科尔森（Crank-Nicholson）格式 1）显式差分格式 （46） 2）隐式差分格式 （47） 3）克兰克-尼科尔森格式 若对流项、扩散项和源项中的参数都取两个时间层上的算术平均值，可得差分方程 （48） 克兰克-尼科尔森格式的差分方程也需要与其它结点的差分方程联立求解。当时间步长较小时，其精度比隐式格式高。 4、源项的差分 设源项在控制容积内均匀，并等于中心点的值，则有 （49） 微分方程离散化后得到的是一个拟线性代数方程，但源项往往是参数的非线性函数 为能够用解线性代数方程的方法求解差分方程，需要将源项线性化 源项线性化有两种办法: 一是用前一时