溶质运移理论水动力弥散方程数值解法.ppt

* 四、对流占主导地位时的数值方法 * 四、对流占主导地位时的数值方法 模型精确度 取决于网格距大小。可以在同一网格上计算速度场。如果采用较少的网格距，可以改善流场的描述。 坐标为 的特征点所在的单元，可用下标i,j确定 精确度受到每一网格内所选择的最初的特征点数目影响。 * 四、对流占主导地位时的数值方法 缺点： （1）实际程序冗长，需记录特征点的踪迹； （2）抽水点或出流边界的特征点需要消除，入渗边界须有特征点生成，不透水边界特征点须反射回来； （3）要通过消去特征点来避免滞留点附近特征点堆积； （4）特征点耗尽的单元，必须产生新的特征点。 * 四、对流占主导地位时的数值方法 动坐标系方法与网格变形方法 能消除过量现象，还保持较陡的浓度锋面。 对式 设坐标变换 得， 取空间步长 ，则 * 四、对流占主导地位时的数值方法 动坐标系方法与网格变形方法 在固定坐标系中取格距为 ，有 动坐标系中格点运动与固定坐标系格点运动重合。 设固定坐标系中格点xj在tn+1时刻与动坐标系中Xi重合，有 动坐标系下浓度结点的表达式 * 四、对流占主导地位时的数值方法 动坐标系方法与网格变形方法 式变成 对半无限长砂柱，当DL=0时 * 四、对流占主导地位时的数值方法 网格变形方法 设tn时刻有限元结点在空间中的分布是 对应浓度分布为 * 四、对流占主导地位时的数值方法 网格变形方法 主要步骤： 结点位置随时间变化，基函数及有限元方程的系数都依赖于时间，每推进一个Δ t都需要重新计算。二维问题还须进行畸形判断 * 四、对流占主导地位时的数值方法 结合动坐标的网格变形方法 一维对流-弥散在某个动点在t时刻位置x可写成 其中，x0是动点的初始位置，有 取时间的一阶偏导数 对流-弥散方程写成 若控制动点速度dx/dt 并让它接近流速u，则可转换成弥散为主的的方程式 * 四、对流占主导地位时的数值方法 若使用伽辽金法，基函数因依赖于结点位置而成为时间的函数。若用 表示随时间变化的基函数，则 对任意t，有 故 * 四、对流占主导地位时的数值方法 随机步行法 采用示踪剂描述污染物运移。只有被污染的特征点在流场中运动。每个特征点都给定一个不变的污染物质量，此质量的和等于排进含水层的总污染物质量。 取许多单个特征点轨迹（随机步行）的平均值，可得到一个有弥散的特征点的分布。要得到一个浓度分布，就要先叠加一个网格并算每个网格单元所含有的特征点数。 * 四、对流占主导地位时的数值方法 随机步行法-以一维运移为例 在x=0处，一个瞬时注入的、质量为Δ M的理想失踪及，在t时刻浓度分布 对一个固定时间t， 看做正态分布 * 四、对流占主导地位时的数值方法 随机步行法-以一维运移为例 在时间t=0和位置x=0处，放置具有污染物质量为ΔM/N的N个特征点，每一个特征点移动距离x，而在时刻t到达它们各自位置，即 Z是一个平均值为0，标准差为1的正态分布随机变量 所得的频率分布，通过一个标准化因子，有 对有限数量的颗粒，可近似去顶C(x,t)，将空间变量划分为长度Δx的若干区间，区间中点浓度 * 四、对流占主导地位时的数值方法 随机步行法-以一维运移为例 若孔隙平均流速是时间和位置的函数，那么时间为0开始的颗粒P在时间间隔Δt内的随机轨迹xp(t)为 * 四、对流占主导地位时的数值方法 随机步行法-二维流场 首先假定流动方向平行于x轴，得到 Z和Z’是正态分布的随机变量的两个值。对任意方向的流动，需考虑弥散张量特性。 * 四、对流占主导地位时的数值方法 随机步行法-二维流场 首先假定流动方向平行于x轴，得到 Z和Z’是正态分布的随机变量的两个值。对任意方向的流动，需考虑弥散张量特性。 对流运动可确定在水流流动方向及其正交方向的弥散运动 * 四、对流占主导地位时的数值方法 随机步行法 * 四、对流占主导地位时的数值方法 随机步行法 数值结果的质量主要取决于特征点数目N，要考虑网格的选择。大网格间距会产生极平均的结果，小网格间距会产生十分粗糙的分布。 在稳定流动和一个