用一般形式的KFVS方法计算溃坝波问题.pdfVIP

·256· 数学、力学、物理学、高新技术研究进展——2008(12)卷 用一般形式的KFVS方法计算溃坝波问题} 陈秀荣 于加举 程 冰袁冬梅姜德民 青岛农业大学山东青岛，266109 摘要通过构造一个分布函数，利用无碰撞的Boltzmann方程，给出求解一般形式的KFVS 方法，计算了溃坝波问题，并与精确解进行了比较，结果表明该方法能很好地模拟 间断解问题。 关键词Boltzmann方程；流矢量分裂法；溃坝波 1．引言 vector 流矢量分裂(flux 和Warmingtll,VanLeeld2]等人提出的。他们将Euler方程中的流失量分成正、负两部分， 采用迎风、TVD等格式计算Euler方程中的对流项，通常的分裂方法是以方程的特征值 动。他假定气体分子的速度分布函数是局部平衡的Maxwell分布函数，然后分别计算出 fluxvector 沿正负两个方向运动分子的流失量。这就是分子运动论的流失量分裂(kinetic spliRing，KFVS)法。Reizl4]用算子分裂(operaterspliRing)法计算了分子运动的微观方程一 Leer ChapmanEnskog型分布函数，将KFVS方法用于计算Navier-Stokes方程。 o】等人用KFVS Nsdigais,9]提出了近平衡的离散速度方法，用于计算欧拉方程。汤华中【l 方法计算了一维激波管问题，徐昆IllJ等人用KFVS方法计算了非平底浅水方程。本文用 KFVS方法计算圣维南方程，并用KFVS方法模拟了溃坝波问题。 2．一维Boltzmann方程以及其与宏观变量间基本关系 一维无碰撞Boltzrnarm方程为 笪+c娑+丘箬：o ”7 8)1t 缸 mac ( fc-u)。 舯似^c)-赤e了是分子平懿分布醭㈡-河道过水断酾桃删子 速度；z—分子的位置；r是时间；Jj一水深；U一流速；C-作用在分子上的合外力；m． ’作者简介：陈秀荣(1978．)，女．山东菏泽人，硕士研究生．助教，主要研究方向：复杂系统建模 用一般形式的KFVS方法计算溃坝波问题 ·257· 分子质量，g．是重力加速度。 一维水流中微观变量厂与宏观变量即过水面积A，流速材，流量Q有如下基本关系： (2) 肛砒J．啦=Q；J_c2肚=Q2／A+圭驴矗2 其中甜=钐，B-河宽。圣维南方程和(1)的联系建立如下，在方程(1)两边同乘以(1，c)’， 对C在整个空间积分，由式(2)可得： 拍 塑：o 一·-‘· af 弧 塑。 a(Q2／A+g鼬2／2)