论文部分内容阅读
在处理突发水污染环境事件中,污染物在河流中的分布情况可以用对流扩散方程来描述。同样很多其他环境相关的问题也都可以转化为对流扩散方程的问题进行分析和解决。因此,对流扩散方程在环境监测以及对污染物的预测和处理领域有着十分重要的意义。但是,很多对流扩散方程问题难以找到解析解,需要对其进行数值求解,而对于突发性水污染事件而言,精确的通过数值计算得到污染物精确扩散位置以及浓度的同时,时效性也是不可或缺的。针对浓度对流扩散方程的数值求解问题,本文主要研究内容如下:文章的第一部分首先针对浓度对流扩散方程进行高精度离散,对内点构造两层八点隐格式,进而,构造与内点格式精度相匹配的边界层差分格式,对浓度对流扩散方程的时间和空间项分别进行相应阶数的泰勒展开,使用待定系数法求出差分格式的差分系数,得到浓度对流扩散方程内点以及边界的时间三阶,空间六阶精度隐式差分格式。进而,对一般情况下的一维高精度差分格式进行Von Neumann稳定性分析,随后对相应算例进行数值验证。最终证明了本文构造的所有格式均满足时间三阶,空间六阶的精度要求,且在一定条件下稳定,稳定性范围宽广,同时一定范围内可以高精度计算对流系数较大的对流占优扩散方程问题。文章的第二部分首先基于第一章构造的高精度差分格式对所得到的三对角方程组提出了一种新的并行计算方法。在p个计算机处理核心分组并行处理的基础上可以并行计算,使得整体并行计算的效率更高,数值算例表明:该并行方法简单易行,解决了隐格式的不易并行计算的问题,并且加速效果随着空间节点总数以及分组数的增加变得更加明显,加速比和方程分块数基本满足线性关系,在保持高精度求解的基础上实现了优良的并行效果。值得注意的是在求解过程中,组与组交叉的未知量可以形成块三对角方程,同样可以使用该方法进行并行计算,可以更大程度上的提高并行求解的效率。因此,本文提出的并行方法适用于二维乃至更高维度的对流扩散方程的并行计算,本文以二维对流扩散方程数值求解为例,给出了本方法和串行方法的计算时间对比分析。加速效果随着空间节点总数以及分组数的增加变得更加明显,且能够很好的保持求解过程的精度需求。文章第三部分对并行计算过程中使用的并行语句进行深入分析,揭示了其循环中的系统周转时间、循环控制和计算规模对于计算效果的影响,通过采用内存映射的方法,高效访问磁盘上由于太大而无法保留在内存中或需要花太长时间而无法加载的大数据集,解决了大型矩阵的数据通讯时间影响整提计算速度的问题。利用MEX混合编译和MATLAB的扩展特性,同时结合C语言进行编码,将计算中的大型循环计算使用C/C++和MATLAB混合编译来完成。高效提升求解大型三对角方程时的并行效果。文章第四部分给出了本文所构造的高精度差分格式在实际环境问题中的应用。分别以上游围油栏作为第一类固定边界,下游收油装置作为第一类移动边界,模拟对河道溢油事故的处理过程。通过采用本文构建的一维浓度对流扩散方程高精度差分格式,数值模拟了溢油发生时,围油栏和收油装置作处理装置时溢油浓度的变化。