用BGK格式计算不可压缩流场.pdfVIP

力学学.报 E ~吨。.6 ~-9.~1IbIIíl期 拥[.tIi]年 !Il月 A TA l\o!ECHANICA SIi\lC ∞l !EEZ歪倒gE但嘀岳Z缸~Z窗理旦 苏铭穗 • {清华大学飞程;力学系，北京!liJiD. 徐昆 M. Ghidaoui {曹海军司t量大字数学系，士本系，香德) 摘要简要介锦 BGK 方法的基本思想，数也方法及其对不可压缩流场rn值汁摩的推广.然后. 摊细介绍将r.rern方法应用于平丽方枢自匠的黯巢.将其结果与 æm 的生在巢进行了详细的比. 较，本方法淤，用来模拟后台阶流动，所得站到2与相应的实碰数据迹行比较.以上比校我明本方. 泌的有效性，同时也确定在计算不叮lli!iì流休时参她的正确选择 关键饲~泣，不可耐流动 E 号!言. 不可压结流场是当班体运功速度比声速小很多时的一种近似.在这样的近似条件下z 流体. 的密度被视为常数.看起来这种简化便问题变得简单了.事实上在求解不可压撞流体的 tðnQt;. ~方程时，由于密度是不变的，所以连蝶万程中不出现密度时时间的导数项，所以它不能. 与动量方程联功lI]解.但;沟满是连攒方惶，在:求解过程中必须保证速度的散度为零，这给数值. 计算带来了许多困难.为了克服这些困难，人们进行了不懈的努力.发展了诸如投影注 t(ïJmtmI lEmlIID动力修正法 (Abdallah， 1987 ~ )等以保证直度塌的散度始终为零.另外也有不少人重新. 将不可压缩流场视为低马赫茧的流动，豆豆在连续万程上附您一个人为的压力不定常项èJ.SJ;;; 1968(1)) ， 构造伪不应常法.类似的作法可见邓(tm!m等和王 (W缸ε、lPmII论文. 最近人们又发展了用格子气 [~mtll%tiil 1 ，1 方法用米计算不可压缩流场 [c.!m[fI等肉，我得. 了成功.但诸方法只能用于古牵连的不可压流场的计算，有一定的局限性 近年来，由圳等人 (Xu 等. 19!)3.1995， 1998[7~9J 发展了一种由 Bolt函副I芳草面衷面. 基于 rrr~ 模型的计算方法(简称 l~陆).其基本思想是从基于统计物耀中的批rr.m四n 方. 程出友，万程中的碰撞项来用[:rnl扎和~在iJtMj]年提出的~模型进行. 简化.以此为基础建立宏现量的辈分格式称为 BGK 锚式.可以证明.在f助理量连续变化的条. 件下，这样作等价于求解~mtm!j切怡5 万程. 对于物理量不连镜变化的条B下并不等价子求解­ ~方理.因此它可以求得比基于 rat!Ic.i!i..;，麟更好的结果.正是这个原因使得~方法得. 到越来越多的应用. 如它已被应用于水放运功. 磁流体问题， 多种抗体掺j昆，揣部等的数值模. 缸中去.这些应用表明该方法具有丰富的物理内涵.计算上更如灵活适应性更强，更稳键‘ • [IUI~:Iß;l收虫第…稿 1 心~收到修改稿 1J1 香熄科核大学弱了吃.;t~ i: (RGC97j98.HKUST6166j97P, RGC96/97.HKUSTrl 8/96盼及嗣依自然科学基金和国旅. 登计划资助项目 力. 学. 学. 报. ~ 瞌矗矗..监画画 本文的目的是将 ræm 方法应用于不可压缩流体流动的数值锁拟.过去读方法主要应用于. 计算可压缩有间航的