溶质运移理论水动力弥散方程的值解法-课件.pptVIP

地下水流数值模拟方法第六章水动力弥散方程的数值解法中国地质大学环境学院2019春

2一、有限差分法基本思想：将空间划分成若干网格，把时间分成若干小时段，每一个网格中心点（格点）处的未知变量（H或C）视为该网格平均值。利用差商代替微商。基本步骤：（1）剖分（2）建立差分方程组（3）求解

3一、有限差分法-导数的有限差分近似图中，去x轴上任意一点i，其坐标为在改点左右相聚为处分别取（i-1）和（i+1），其坐标分别为和以i为中心，泰勒展开C(x)

4一、有限差分法-导数的有限差分近似整理并略去余项（6-1）-（6-2），再除以略去余项（6-1）-（6-2），再除以略去余项：分别一阶导数的向前差分、向后差分及中心差分二阶中心差分

5一、有限差分法-一维水动力弥散的差分解法一维水动力弥散方程（6-7）（1）显格式式中仅有一个未知数，解得式（6-7）中的对流项取中心差分

6

可以证明，稳定性准则要求

即即。（1）（2）。由条件（2），格距要求很小；由（2）可知，鉴于较小，导致不能太大，将（2）代入（1）式中，得到

7若对流项改为向后差分解得