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

1、地下水流数值模拟方法地下水流数值模拟方法第六章 水动力弥散方程的数值解法2一、有限差分法一、有限差分法3一、有限差分法一、有限差分法- -导数的有限差分近似导数的有限差分近似xxi ixx1-x1 - i ix1x1ii4一、有限差分法一、有限差分法- -导数的有限差分近似导数的有限差分近似x2分别一阶导数的向前差分、向后差分及中心差分分别一阶导数的向前差分、向后差分及中心差分2x二阶中心差分二阶中心差分5一、有限差分法一、有限差分法- -一维水动力弥散的差一维水动力弥散的差分解法分解法（6-76-7）1niC11122221()2()()2nnnnLLLiiiiDtu tDtDtu tCCC

2、Cxxxxx（1显格式式中仅有一个未知数，解得 1111122()2nnnnnnniiiiiiiLCCCCCCCDutxx式6-7中的对流项取中心差分 6可以证明，稳定性准则要求可以证明，稳定性准则要求221()LDtx2()2LxtD 22()Lu tDtxx()22LLLLDa uDxau 即 即 。（1）（2）。2Lxa xt22()222LLLLLDa uxDatDuu由条件由条件2），格距要求很小；由，格距要求很小；由2可知，鉴于可知，鉴于较小，导致较小，导致不能太大，将不能太大，将2代入代入1式中，得到式中，得到 7111122()nnnnnnniiiiiiiLCCCCCCCDut

3、xx11122212()()()nnnnLLLiiiiDtu tu tDtDtCCCCxxxxx2(1)22(1)LLLxxxDtDDuxuuxx 若对流项改为向后差分若对流项改为向后差分 解得 稳定性要求稳定性要求不难看出，稳定性限制比对流项取中心差分有所改善。不难看出，稳定性限制比对流项取中心差分有所改善。8（2 2隐式格式隐式格式111111111122()2nnnnnnniiiiiiiLCCCCCCCDutxx 1111122221()2()()2nnnnLLLiiiiDtu tDtDtu tCCCCxxxxx整理后得到整理后得到 隐格式是无条件稳定的。9（2 2Grank-Nicol

4、sonGrank-Nicolson格式格式 整理后得到整理后得到 取隐式和显示的平均，即取隐式和显示的平均，即Grank-Nicolson格式格式xCCxCCCDtCninniLn2u-221C, 1in12, 1i, in1, i1ni,，（），xCCxCCCDninniiL2u-21, 1i1n121, 1i1,1 -n11, 121,21, 122-t-122t-niLniLniLCxtuxDCxtDCxtuxDniLniLniLCxtuxDCxtDCxtuxD, 1 -2,2, 122-t-122t10一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（

5、4-564-56）（1显格式式中仅有一个未知数式4-56中的对流项取中心差分 2, 1i, in ,1, i1nj,i,2CxCCCDtCnjnjjiLnj，xCCyCCCDnjjinjnjjiT2u-2, 1, in ,12, 1, i, in , 1，11一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法化简后，有njiTLnjiLnjiCyDxDCxDC,22, 121,t2t21x2tu-tnjiTnjiTnjiLCyDCyDCxD, 1 -,2, 1,2, 1 -2ttx2tu-t涉及以i，j为中心的5个网格点在tn时刻的已知浓度12一、有限差分法一、

6、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（2隐格式式4-56中右端的对流项取中心差分，右端个C的时阶均取n+1水平 21, 1i, 1, in , 11, i1nj,i,2CxCCCDtCnjnjjiLnj，xCCyCCCDnjjinjnjjiT2u-21, i1n ,121, 1, i1, in , 1，13一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（2隐格式整理后收敛且无条件稳定涉及以i，j为中心的5个网格点在tn+1时刻的未知浓度14一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（3Gra

7、nk-Nicolson格式将隐式格式的两式相加除以215一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（3Grank-Nicolson格式整理得16一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（4交替方向隐式法ADI法）分两次对三对角矩阵求逆，将一个t分成两个t/2计算 第一个半时间步，对x方向的偏导数采用隐式差分，对y方向的偏导数采用显示差分。2, 1, i, in , 1221, 1i21, i21n ,1 -22yCCCDxCCCDnjnjjiTnjnjjiL，2tu-j , i21nj21j , 1i21-njn

8、iniCCxCC，17一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（4交替方向隐式法ADI法）整理得 第二个半时间步，对y方向的偏导数采用隐式差分，对x方向的偏导数采用显示差分。21, 1, i1, i1n , 1221, 1i21, i21n ,122yCCCDxCCCDnjnjjiTnjnjjiL，18一、有限差分法一、有限差分法- -二维水动力弥散的差二维水动力弥散的差分解法分解法（4交替方向隐式法ADI法）整理得收敛且无条件稳定2tu-j , i21nj21j , 1i21-njniniCCxCC，19一、有限差分法一、有限差分法- -求解差分方程

9、的计算求解差分方程的计算机程序举例机程序举例算例算例 见教材见教材P81-85P81-8520二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法（6-6-2626）21二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法C（6-6-2929）22二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法23二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法（6-6-3232）代入代入6-296-29）24二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法25二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法基

10、函数NLx，y形状如同一顶高等于1、有4条直线斜边和4条下凹型曲边的尖顶斗笠26二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法27二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法28二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法29二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法（6-6-3434）30二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法nnmmjjiiyCyCyCyC,x,x,x,x、 nnmmjjiiyyyy,x,x,x,x、31二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金

11、有限单元法txC, y,eyNyNyNyNeeee, x, x, x, xnmji、（6-6-3636）32二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法LDxy,4neL33二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法34二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法35二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法36二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法37二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法1t 1tC lC B 2tt C Cttt与38二、有限

12、单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法39二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法算例：单元弥散矩阵与整体弥散矩阵的关系算例：单元弥散矩阵与整体弥散矩阵的关系设研究区划分为设研究区划分为8 8个单元，个单元，1616个结点，四周为隔水边个结点，四周为隔水边界编号如下：界编号如下： 40二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法算例：单元弥散矩阵与整体弥散矩阵的关系算例：单元弥散矩阵与整体弥散矩阵的关系计算每个单元弥散矩阵，按双下标已知，放置单元弥计算每个单元弥散矩阵，按双下标已知，放置单元弥散矩阵并叠加，得到总体弥散矩阵。散

13、矩阵并叠加，得到总体弥散矩阵。41二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法42二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法43二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法积分求解可用高积分求解可用高斯求积方法，详斯求积方法，详见见P97-105P97-10544二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法45二、有限单元法二、有限单元法- -伽辽金有限单元法伽辽金有限单元法算例：二维水动力弥散的伽辽金法计算机程序算例：二维水动力弥散的伽辽金法计算机程序 数学模型写成数学模型写成46三、过程与数值弥散三、

14、过程与数值弥散意义。意义。47三、过程与数值弥散三、过程与数值弥散48三、过程与数值弥散三、过程与数值弥散49三、过程与数值弥散三、过程与数值弥散相当于弥散系数的弥散项相当于弥散系数的弥散项50三、过程与数值弥散三、过程与数值弥散此条件下的数值弥散系数此条件下的数值弥散系数LeeLLDxuPPDxuDD即22n51三、过程与数值弥散三、过程与数值弥散eP2eP2ePePePeP控制控制xx的大小从而减少数值弥散的大小从而减少数值弥散52三、过程与数值弥散三、过程与数值弥散53三、过程与数值弥散三、过程与数值弥散结果也造成数值弥散结果也造成数值弥散54四、对流占主导地位时的数值方法四、对流占主导

15、地位时的数值方法对流项的有限差分化过程中乘上一个适当的对流项的有限差分化过程中乘上一个适当的权因子，则波动和过量得到改善或消除权因子，则波动和过量得到改善或消除55四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法2/1i2/1ix,x56四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法57四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法58四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法59四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法2/1ix2C1iiC2/1ix60四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法6

16、1四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法项用沿特征线的示踪质点的运动描述项用沿特征线的示踪质点的运动描述o0t0eH或0div tyt ,x62四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法63四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法64四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法 ttppyx、65四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法特征点。特征点。66四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法o取空间步长取空间步长，那么，那么oo67四、对流占主导地位时的数值方法四、对流占主导

17、地位时的数值方法o动坐标系下浓度结点的表达式动坐标系下浓度结点的表达式68四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法69四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法nnXCnX70四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法o71四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法o若控制动点速度若控制动点速度dx/dt dx/dt 并让它接近流速并让它接近流速u u，则可转换成弥散为主的的方程式，则可转换成弥散为主的的方程式72四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法t ,xNi txNtCtxCNii,i173四

18、、对流占主导地位时的数值方法四、对流占主导地位时的数值方法o74四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法t4exp4,20LLDutxtDCtxC75四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法o 对有限数量的颗粒，可近似去顶对有限数量的颗粒，可近似去顶C(x,t)，将，将空间变量划分为长度空间变量划分为长度x的若干区间，区间中的若干区间，区间中点浓度点浓度N76四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法77四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法78四、对流占主导地位时的数值方法四、对流占主导地位时的数值方法对流运动可确定在水流流动方对流运动可确定在水流