基于OpenFOAM的浅水流动和物质输移扩散模拟.pdf

中国环境与生态水力学2008 基于OpenFoAM的浅水流动和 物质输移扩散模拟 王智勇 朱德军 陈永灿 (清华大学水利水电工程系，北京 100084) 摘 要 在openFOAM基础上进行二次开发，建立了适用于平面二维浅水流动 和物质输移扩散过程模拟的数值模型，利用类似PIs0算法的“速度一水位”预 测校正法处理“速度一水位”耦合问题，采用水深平均的矗一￡紊流模型封闭紊 动扩散项。将模型应用于侧向排放以及长江三峡水库涪陵磷肥厂污水排放问题进 行模型验证，计算结果与实测数据吻合良好。研究表明，OpenFOAM开源、易 于快速开发、用户群广的特点使其在水动力与水环境数值模拟研究中具有良好的 应用前景。 关键词 CFD 浅水流动 面向对象 0penFoAM “速度一水位”预测校正 1 引言 传统的CFD程序大多采用面向过程的程序设计方法，在这种开发方式下常常需要 面临底层代码实现细节，数据的传递较为复杂，代码重用性和维护性较差[1]。近年来， 面向对象的程序设计方法因其代码的高度重用性和易维护性，以及适合解决大规模复杂 Field And 问题的优势而在CFD领域逐渐得到了应用。0penFOAM(OpenOperation 言编写，提供了大量类库和一系列CFD求解器。作为开源软件，其支持二次开发以及 面向对象的程序设计方案使用户可以根据自身需要很方便地开发针对某特定问题的求解 器，因此目前已拥有较多的用户群．广泛应用于各种流体计算以及偏微分方程求解等科 学领域‘引。 自然界中的水流如宽阔河流、湖泊、海湾、河口等多为浅水流动，其水平方向上的几 何、运动尺度要远大于垂直方向上的尺度。在一维模型不能提供足够信息和三维模型计算 量过大的情况下，沿水深积分得到的平面二维模型成为模拟浅水流动的首选。0pen一 基金项目：国家重点基础研究发展计划(2006CB403304) 作者简介：王智勇(1983一)，男，河北承德人，博士研究生，主要研究方向：湿地水环境，E—mail： s．thu．edu．cn wang—zy@mail · ·39 FOAM并未提供平面二维浅水流动求解器．但其采用任意形状非结构网格的特点以及针 对流体力学控制方程通用形式的面向对象求解方案非常适用于各种流体力学数值模型的开 发。本文基于openFOAM平台开发了二维浅水流动模型和物质输移扩散模型，并通过算 例对模型进行验证。 2浅水流动和物质输移扩散模型的开发 2．1 openFOAM中的PDE求解体系 流体力学问题可以用偏微分方程(PDE)来描述，通常具有下面的形式： 旦—}譬+甲．(PuQ)一可．(矿甲o)一s，o+s，(1) 式中：u为流速；p为密度；Q为流体的某种属性，可以是标量(如浓度c)也可以是向 量(如流速u)。 积法进行求解。其基本思路是：①利用有限体积离散，计算 域被划分为互不重叠但充满整个计算域n的离散单元V．，如 图l所示；②在每个单元中对控制方程进行积分，并利用高 斯定理将散度项转化为单元界面上的通量积分t则问题最终 转化为如何利用网格中心点的变量插值求得界面上的通量； ③对瞬时流动，以上转换过程将式(”离散为关于时间￡的 圈1 0penFOAM非结构 常微分方程，需要再利用有限差分方法进行时间离散；④将 网格示意图Ⅲ