第五章地下水污染预测重点.pptVIP

主要内容 概述（了解） 近似解析法（了解） 数值法（熟悉） 数理统计法（熟悉） 灰色预测法（了解） 国外地下水模型软件 有限差分法的基本思路是按照时间步长和空间步长将时间和空间区域剖分成若干网格，用未知函数在网格节点上的值所构成的差分商近似代替所用偏微分方程中出现的各阶导数，从而把表示变量连续变化关系的偏微分方程离散为有限个代数方程，然后解此线性代数方程，以求出污染物浓度在各网格节点上不同时刻的解。 设f(x)为任一足够光滑的函数，把f(x)沿x的正向展开为泰勒级数： 上式称为f(x)的一阶向前差分 如果把f(x)对x的负向展开为泰勒级数： 上式称为f(x)的一阶向后差分 把上面两个差分方程相减，可得： 上式称为f(x)的一阶中心差分 同理可得f(x)的二阶导数的差分公式： 若把上述一元函数f(x)换为浓度函数，并按网格节点编号，则对x的一阶向前、向后和中心差分分别为： C(x,y,z,t)对x的二阶差分公式为： 设在无限含水层中，存在一维均匀流场，其渗透速度V=nu，流动方向为x轴正向。则一维弥散方程为： 对时间区域[0,T]和空间区域[0,L]都作等距剖分，设时间步长为△t，空间步长为△x，把第i个节点xi处在tn时刻的浓度记为Cin。 显示差分格式 将 及 差分公式中的浓度取为tn时刻的值，便可得到其显示差分格式： 显示差分格式的稳定条件 在应用显示差分格式时，剖分的步长必须满足： 隐示差分格式的稳定条件 在应用隐示差分格式时，剖分的步长必须满足： 隐示差分格式的矩阵形式 Crank-Nicolson差分格式 取显示差分格式和隐示差分格式的平均，便可得到C-N差分格式： 令A·C=δ 由此可得： b1=β1 b1w1=γ1 → w1=γ1/b1 b2+a2w1=β2 biwi=γi → wi=γi/bi (i=1……M-1) bi+aiwi-1=βi → bi=βi-aiwi-1 bM=βM’-aM’wM-1 由以上等式可知： b1y1=δ1’ → y1=δ1’ /b1 a2y1+b2y2=δ2 → y2=(δ2-a2y1)/ b2 aiyi-1+biyi=δi → yi=(δi-aiyi-1)/ bi aM’yM-1+bMyM=δM → yM=(δM-aM’yM-1)/ bM 至此可以求出所有的Y(y1, y2, …, yM)T值。 追：根据B·Y=δ，求出bi和wi，并求出yi； 赶：根据W·C=Y，求出Cin+1。 若多孔介质是各向同性的，流场为均匀的一维流场，渗透速度V=nu，流动方向取x轴的正向，与流速垂直为y方向，对流扩散方程可简化为： 将研究区域网格化 建立矩形网差分格式 隐式差分格式 于是，内节点的差分方程为 整理后可得 边界条件 初始条件 将网格取为正方形，即△x=△y=h，方程就可简化为： 交替方向差分格式（ADI） ADI法的实质是将二维隐式差分格式化为一维的隐式差分格式。 在第n层与第n+1层两层正中间引入一个过渡层n+1/2。 交替方向差分格式（ADI） 分两步走：在前半时间步长内，对x偏导用隐式表示，对y的偏导用显示表示；在后半时间步长内，对x偏导用显示表示，对y的偏导用隐式表示。 交替方向差分格式（ADI） 前半时间步长： 交替方向差分格式（ADI） 后半时间步长： 计算步骤： （1）根据初始条件或前一时段计算结果，可知tn时刻的浓度值，那么第一个方程中只有三个未知数（n+1/2时刻）。如果固定j值，对i的联立方程组就是三对角线的（追赶法），即沿x方向追赶，可求出过渡层（n+1/2）层上的全部浓度值。 计算步骤： （2）根据前一步求出的过渡层（n+1/2）层上的浓度值，再按同样的方法求出第二个方程中的（n+1）时刻的浓度值。即对于一个固定的i值，方程也是一个三对角方程组（追赶法），沿y方向追赶。 在有限差分法中，研究区域被一系列的网格节点所代替，在各个节点上以差分代替导数而得到求解相应偏微分方程的线性代数方程组。 有限单元法是把研究区域剖分为有限个子区域