二维河道移动边界的流矢量处理法.pdf

141 数学技术应用科学 Science Mathemali璐Tedmch硒舢耻ied 二维河道移动边界的流矢量处理法· 程冰m，周晓阳2 (1．莱阳农学院，山东青岛266109；2．华中科技大学数学系，湖北武汉430074) 摘要：本文从Bolt硼锄方程出发采用二维有限体积ⅪⅣS方法，给出了模拟河道洪水演进时于湿交界处移动边 界的一种新方法，并成功地模拟了大坝下游为干河道的二维溃坝模型． 关键词：Boltzmann方程；浅水方程；有限体积KFVS方法；移动边界 中图分类号：rvl31．4文献标识码：A文章编号：7—5416—2336—9(2006)一0141—04 明渠非恒定流是水力学中最常见的工程实际问 1 Boltzm 方程及Boltmnann变量与宏 题之一，它一般指流速和水深随时间变化的具有自由 观变量之间的关系 表面的水流．水流不恒定现象的产生机制很复杂，其 中之一就是流场固壁的运动变化使流场条件改变，因 二维无碰撞Bohzmann方程可写为如下形式： 此流体与河岸交界处的移动边界问题成为水力计算 菩+气箬+勺蒡+鲁老+鲁苦=。㈣ 中常见的边界条件．目前处理这一问题最常用的方法 是假设干地区域被一层水深很小的薄水层所覆盖再 其中，为粒子速度的Maxwell分布函数，Cx，c，为粒子 沿茗，Y方向的速度z‘，秽为水流在x,y方向的速度，见， 利用变动的网格进行计算【1剑，这种处理方法简单实 用，但是它要求计算域!必须覆盖到洪水可能到达的所 E为x,y方向合外力，m为粒子质量 将式(1)乘以币=(1，q，cT)7然后在粒子速度空 有地方，这在计算中造成很多额外计算量，而且，对于 间进行积分，可得二维明渠浅水方程： 某些问题，如溃坝波问题，这种处理方法可能会导致 完全不同的物理现象． 掣+丛+照：S (2)、7 函 Ox 03 近年来，一些学者利用Bhatrlagar—Gross—Krook 其中 U=(^，hu，hv)7 于1954年提出的BGK波尔兹曼方程为许多流体力学 S=(O，醇(‰一强)，gh(哂一％))1 问题建立了数值模型旧4】，取得了很好的效果．这种 模型采用自底向上的建模的思想，从分子运动论出 F=o＼＼c,fdc，dc,，f＼船ldc，， 发，对微观变量即分子速度概率分布函数求矩得到微 观变量与宏观变量之间的关系．xu等人最先利用 f fc，讹咄)7=(地，^扩+譬，huv)7 BGK形式的Bohzmann方程建立了可压缩大气流动 数值模型1；邓家泉第一个应用BGK—Boltzmaan方 G=t＼＼cfdc,dc,，＼＼c,cfdc,dcr， 程建立了适应水流、水质的数值模型怕1；我们曾利用 借助微观分子速度的符号构造宏观通矢量分裂的格 f fc弘，嵋)7=(hv，huv，^矿+学)7 式一有限体积KFVS格式处理了一维干坝的移动 边界问题；本文拟将此方法推广到二维情形． 上式中g为重力加速度，