地下水污染问题的对流扩散方程在运动质点坐标系上的数值解法.pdfVIP

地下水污染问题的对流扩散方程 在运动质点坐标系上的数值解法 丁家平 (南京水利科学研究院，南京，210029) 摘要：本文给出一种运动质点坐标系上的欧拉一拉格朗日方法解关于地下水污染问题的对流扩散方程。模型网 格采用不规则四边形单元，引入运动质点局部坐标系处理对流扩散过程，该坐标系的两个主轴与水流流向及其 法向一致；对非网格结点上的变鼍采用次参数型单元二阶拉格朗El插值。该方法数值计算简捷，计算精度较高， 对复杂的模型轮廓具有良好的几何适应性，也较好地解决了对流扩散方程的数值弥散和振荡问题。该方法及所 编制的计算机程序经标准算例考核，并应用于某地下水污染实际问题。 关键词：对流扩散方程 欧拉一拉格朗日法 运动质点局部坐标系 次参数单元 对流扩散方程的大部分数值解法可以归结为欧拉法、拉格朗日法和欧拉也格朗El法。普通的欧拉法具有 固定的网格。但计算结果表明它有时不能适应陡峭的浓度峰面，会产生较明显的数值弥散和振荡。拉格朗日 法虽能适应陡峭的浓度峰面，但变动的网格容易产生变形和翻转等弊病。欧拉也格朗日法则在一定程度上兼 有上述两种方法的优点，又克服了它们的缺点。本文在已有的研究【1巧】的基础上加以改进，建立一种基于运 为：①在质点沿流线运动的过程中引入一个运动质点局部坐标系(SN)，它的两个主轴与水流流向及其法向是 一致的，计算格式简捷；②固定网格由不规则四边形构成，它比起正方形或矩形来，对轮廓有较好的几何适 应性；③运动质点浓度插值采用具有次参数单元概念[61的二阶拉格朗日多项式。 1 对流扩散方程的欧拉——拉格朗日解法 1．1通常的欧拉一拉格朗日法 不可压缩流体的二维对流扩散方程为 塞+圪塞+巧嘉=昙(巩塞+％争+瓦0【％瓦ac+％争 (·) 式中：C为浓度；％、巧为真实流速分量；Dxx，p∥、％、D为流体动力扩散系数张量的分量。 流场中任意运动质点其坐标随时间的变化曲线为 x=工(f) Y=y(f) (2) 根据流速的定义 ， 。 圪=五dx 0=石dx (3) 沿特征线，运动质点的浓度对时间的全导数(或随体导数)可以写作 ‘业：丝-II-丝垒+丝生 I‘．- 一==…■·一一 (4) c3t况cot Oy +dtOt， 将式(1)代人式(4)，得 面dc=昙(巩塞+％矽CqC+瓦a【％瓦ac+％瓦ac) (5) 这是对流扩散方程的一般的拉格朗日形式。解式(5)和特征线方程可以分别获得因扩散和对流引起的质点 作者简介：丁家平，1943年生，男，教授级高工，从事地下水水流和水质数模研究及坝工渗流试验、计算和安全控制分 析。E-mail：jpd@hotmail．com ·429· 浓度变化。 1．2运动质点坐标系(SN)上的欧拉啦格朗日法+ 引入运动质点局部坐标系(SN)，它是以九结点四边形单元ABCD的结点的上风点为坐标原点的。以九结 点四边形单元的中间点Q为例，将其上风点P作为坐标原点，以流线PQ及其法向作为二个主轴，构成运动 质点P的局部坐标系(SN)，见图l。根据欧拉也格朗日方法的基本原理，在以P点为原点的运动质点局部坐 标系(SN)上，可以写出类似式(5)的方程式