可压缩_不可压缩流体交界面高精度数值方法的研究.pdf

《可压缩_不可压缩流体交界面高精度数值方法的研究.pdf》由会员分享，可在线阅读，更多相关《可压缩_不可压缩流体交界面高精度数值方法的研究.pdf（147页珍藏版）》请在三一文库上搜索。

1、中国科学技术大学 博士学位论文 可压缩/不可压缩流体交界面高精度数值方法的研究 姓名：马东军 申请学位级别：博士 专业：流体力学 指导教师：尹协远;孙德军 2002.6.1 摘要 本文主要研究流体交界面高精度数值模拟方法。重点发展了V o l u m e O f - F l u i d 方 法和L e v e l S e t 方法在可压缩流动和不可压缩流动界面问题中的应用，如R i c h t m y e r - M e s h k o v 不稳定性以及表面张力诱导的液滴动力学等此外，本文还考虑了高效并行 计算方法在界面数值模拟中的应用，并讨论了浅水自由面问题的数值模拟方法 首先本文发展了一套

2、简单实用的混合型多介质P P M 方法，结合S a u r e l S h y u e 的 计算模型和经典的P P M 方法数值模拟高密度比的多介质可压缩流动问题同时本 文讨论了将L e v e l S e t 界面跟踪G h o s t F l u i d 方法同高精度的流场计算格式( T V D - A C W E N O A C ) 相结合应用于一般状态方程的可压缩多介质流动模拟算例表明两 种方法可以有效应用于可压缩多介质交界面的数值模拟，相对于多介质P P M 方法， L e v e l S e t 界面跟踪方法能更好地无数值耗散地确定交界面的位置，适合于高分辨率交 界面计算，但是后者

3、的健壮性不如前者，对于状态方程差异较大且高密度比的情况还 不能够很好地处理本文还探讨了计算一般状态方程多介质可压缩流动问题的数值方 法，本文应用L e v e l S e t 技术捕捉界面位置，在界面附近采用守恒数值离散，并用双波 近似求解一般状态方程R i e m a n n 问题一维算例表明，该方法能准确地模拟交界面的 位置，可处理任意状态方程的多介质问题 一本文基于M P I 消息传递并行环境发展了多介质P P M 方法和L e v e l S e t 界面追踪 G 脯s t ，F l u i d 方法的并行算法，以进行大规模界面不稳定性的数值模拟由于两种方法 均为显式格式，本文利用粗粒

4、度数据并行，可达到很高的并行效率本文还将两种并 行计算方法应用于R i c h t m y e r M e s h k o v 不稳定性现象的数值研究，讨论了二维气液交 界面的两种典型情况：重流体或轻流体中激波加速受正弦扰动的交界面，并清晰地得 到了明显的钉状结构和泡状结构，以及后期发展的蘑菇状交界面三维结果同二维情 况演化趋势极为相似，区别在于定量的交界面扰动强度和扰动增长率的不同此外。 本文还考虑了周向和轴向内聚爆炸波诱导的柱形R _ M 不稳定性 本文还介绍了两种高精度的不可压缩流动界面模拟方法：V O F 方法和L e v e l S e t 方法，以及它们的相关改进，并将它们归结到一

5、个统一的特征相函数框架下；通过不 可压动量方程粘性项的半隐处理，提高了计算方法的稳定性；同时借用了多重网格方 法，可以准确高效地处理变密度问题的压力P o i s s o n 方程，因而可以模拟高密度比和 高粘度比的实际问题数值算例表明两种方法都能够正确地模拟交界面的演化以及拓 扑结构变化；L e v e l S e t 方法在处理表面张力问题的时候较V O F 方法精度更高一些。 也更加准确一些；但是V O F 方法的体积守恒保持得要比L e v e l S e t 方法好得多 最后，本文将计算气体动力学中的混合有限体积有限元方法推广到一般浅水流 动的数值模拟，通过基于水面高度的浅水方程来解

6、决由于底坡和摩阻源项带来的通量 梯度和源项不平衡问题，使用简单的中心差源项离散得到满足精确守恒C 性质的流场 解同时本文将非结构网格自适应方法应用到浅水方程的模拟中来，并仔细讨论了一 种基于截断误差的误差指示器数值算例表明，本非结构网格自适应混合有限体积 有限元方法是一种高效算法，能够很好地捕捉流场的细节，适用于复杂几何区域的一 般浅水自由面流动y A b s t r a c t T h i st h e s i s p r e s e n t s t h en u m e r i c a li n v e s t i g a t i o no ft h eh i g hr e s o l u

7、 t i o ni n t e r f 靠 c i a ls i m u l a t i o nm e t h o d s e s p e c i a l l yV o l u m e O f - F l u i dm e t h o d sa n dL e v e l S e tm e t h o d s ， a n dt h e i ra p p l i c a t i o n si nc o m p r e s s i b l ea n di n c o m p r e s s i b l ef l o w s s u c ha sR i c h t m y e r - M e s h

8、k o v ( R - M ) i n s t a b i l i t y ，d r o pd y n a m i c si n d u c e db ys u r f a c et e n s i o n ，a n ds oo n F u r - t h e r m o r e ，w ec o n s i d e rt h eh i g hp e r f o r m a n c ep a r a l l e lc o m p u t i n gi nt h es i m u l a t i o no f i n t e r f a c i a lf l o w sa n dt h en u

9、m e r i c a lm e t h o df o rs h a l l o ww a t e rf r e e - s u r f a c ep r o b l e m s F i r s t l y ，as i m p l ef l u i d m i x t u r et y p eP P Ma l g o r i t h mf o rm u l t i c o m p o n e n tf l o w ， w h i c hc o m b i n e st h ef l u i d m i x t u r em o d e lo fS a u r e l S h y u ea n

10、dt h ec l a s s i c a lP a r a b o l i c P i e c e w i s e M e t h o d ( P P M ) ，i sd e v e l o p e dt os i m u l a t et h ec o m p r e s s i b l eh i g hd e n s i t y r a t i o m u l t i f l u i df l o w s T h eG h o s t F l u i dm e t h o dw i t h1 e v e l s e ti n t e r f a c et r a c k i n ga n

11、 dh i g h r e s o l u t i o nn u m e r i c a ls c h e m e s u c ha sT V D A Ca n dW E N O A Cs c h e m e i sa l s ou t i l i z e d t os i m u l a t et h em u l t i f l u i df l o w sw i t hv a r i o u se q u a t i o no fs t a t e ( E O S ) N u m e r i c a lr e - s u i t ss h o wt h a tb o t hm e t h

12、 o d sc a nr e s o l v et h em u l t i f l u i di n t e r f a c i a lf l o w sc o r r e c t l y I n c o n t r a s tw i t hm u l t i c o m p o n e n tP P Mm e t h o d ，t h eL e v e l S e tm e t h o di sm o r ea c c u r a t e a tl o e a l i z a t i o no ft h ei n t e r f a c ea n dd on o ts m e a ro u

13、 tt h ed e n s i t yp r o f i l e b u tt h eP P M m e t h o di s s i m p l et op r o g r a ma n dm o r er o b u s tf o rh i g hd e n s i t y ，r a t i of l o wa n dl a r g e d i f f e r e n c eo ft h eE O Sa tb o t hs i d e so ft h ei n t e r f a c eI na d d i t i o n ，w ed i s c u s sah i g h r e s

14、o l u t i o ni n t e r f a c et r a c k i n gm e t h o dt os i m u l a t em u l t i f l u i df l o w sw i t hg e n e r a lE O S I t c a p t u r e st h ei n t e r f a c eu s i n gL e v e l S e tt e c h n i q u ew i t hc o n s e r v a t i v en u m e r i c a ld i s c r e t i z a - t i o na ti n t e r f

15、 a c e T h er e s u l t ss h o w st h a tt h ec o m p u t a t i o ni sh i g ho r d e ra c c u r a t ea n d n on u m e r i c a lo s c i l l a t i o na n ds m e a r i n ga tm u l t i f l u i di n t e r f a c e T h e p a r a l l e la l g o r i t h m sf o rm u l t i c o m p o n e n tP P Mm e t h o da n

16、dL e v e l S e tG h o s t F l u i dm e t h o d ，b a s e do nt h eM P Ip a r a l l e le n v i r o n m e n t ，a r ed e v e l o p e dt oi m p l e m e n t t h el a r g es c a l ep a r a l l e ls i m u l a t i o n T h eh i g hp a r a l l e le f f i c i e n c yi sg a i n e db yt h i c k g r a i n t y p ed

17、 a t ap a r a l l e lm o d e T h et w op a r a l l e lm e t h o da r ea p p l i e df o rn u m e r i c a ls t u d y o fR M I a n dt h eo b v i o u ss p i k ea n db u b b l es t r u c t u r ea n dt h ef l n a lm u s h r o o ms h a p e c a nb er e s o l v e d T h eq u a l i t a t i v eb e h a v i o ro

18、f2 DR M Ii st h es a m ea st h a to f3 D b u t t h eq u a n t i t a t i v eb e h a v i o ri sd i f i e r e n t F u r t h e r m o r e w ed i s c u s st h eR M Id r i v e nb yt h e i m p l o d i n gc y l i n d r i c a ls h o c kw i t ha z i m u t h a la n da x i a ld i s t u r b e di n t e r f a c e T

19、 h eh i g hr e s o l u t i o ni n t e r f a c i a ln u m e r i c a lm e t h o d sf o ri n c o m p r e s s i b l ef l o w s V O F m e t h o da n dL e v e l S e tm e t h o d ，a n ds o m ei m p r o v e m e n to ft h e s em e t h o d sa r ea l s od i s c u s s e di nt h i sd i s s e r t a t i o n T h eg

20、 o v e r n i n ge q u a t i o n sw i t hu n i f o r i l lp h a s ef u n c t i o na r e p r e s e n t e de s p e c i a l l yf o rt h ec a p i l l a r yf o r c e s T h em o r es t a b l es c h e m ec a nb eo b t a i n e d t h r o u g h a s e m i i m p l i c i tS t o k e ss o l v e r ，a n dam u l t i g

21、 r i da l g o r i t h m f o rP o i s s o ne q u a t i o ni n v a r i a b l ed e n s i t yf l o w si su t i l i z e df o rh i g h v i s c o s i t yr a t i oa n d h i g h d e n s i t yr a t i op r a c t i c a l p r o b l e m s N u m e r i c a lr e s u l t ss h o wt h a tb o t hm e t h o d s c a nd e 甜

22、w i t ht o p o l o g yc h a n g e s L e v e l S e tm e t h o dh a st h eh i g h e rr e s o l u t i o nf o rs u r f a c et e n s i o nc o m p u t a t i o nt h a nV O F m e t h o d b u tr i g o r o u sm a s sc o n s e r v a t i o ni se l u s i v ei nt h eL e v e l S e tm e t h o dw h i l ea 1 1 t h ea

23、 l g o r i t h m si nV O Fm e t h o df u l l yc o n s e r v et h et o t a lm a s s A tl a s t Av e r t e x c e n t r e df i n i t e - v o l u m e f i n i t e e l e m e n tm e t h o di sd e v e l o p e df o rs o l v - i n g 2 Ds h a l l o ww a t e r e q u a t i o n sw i t h s o u r c et e r m sw r i

24、t t e ni nas u r f a c ee l e v a t i o n s p l i t t i n g f o r m T h em e t h o di si m p l e m e n t e do nt h ea d a p t i v eu n s t r u c t u r e dg r i d sb a s e do na ni m p r o v e de r r o ri n d i c a t o r R e s u l t si n d i c a t et h em e t h o di s a c c u r a t e e m c i e n ta n

25、 ds u i t a b l e f o rs t m l l o ww a t e rf r e e s u r f a c ep r o b l e m sw i t hc o m p l e xt o p o g r a p h y I I 第1 章绪论 1 1研究交界面及自由面问题的意义 自然界中存在着各种各样的交界面和自由面问题，比如水中的气泡在浮力的作用下上升直到 和水面作用然后破碎和消失；在沸腾流动中众多气泡不断浮起并发生融合并附带有复杂的热传递 过程；空气中液滴下落到水面引起的小液滴飞溅和水面荡漾的现象；还有超新星爆炸中的交界面 不稳定性现象，等等同时在各种工业实际应用中也存

26、在着广泛的运动界面，例如水下气体高速射 流以及水下爆炸过程，晶体生长中枝晶的形成、演化和发展，惯性约束核聚变中的R a y l e i g h - T a y l o r 不稳定性和R i c h t m y e r - M e s h k o v 不稳定性现象。还有电雾化、金属的铸造和熔焊等众多的应用在 本文中主要考虑的是比较基础的互不掺混可能含表面张力作用的多流体交界面和自由面问题 对于交界面和自由面问题的研究及其数值模拟有着重要意义：( 1 ) 交界面问题涉及到高精度 的数值算法的发展早期的L a g r a n g e 方法精度都比较低，而且难于处理交界面的复杂演化，特别 是当交界面发

27、生融合或破裂等拓扑结构剧烈变化的情况现代的方法一方面可以处理复杂交界面 演化，同时也更注重高精度地捕捉交界面的位置交界面上一般存在密度间断等数学上比较奇异 的同题，如何准确地模拟交界面间断一直是近几十年来研究的重点课题之一( 2 ) 交界面同题也 有着广泛的实际应用背景比如R a y l e i g h - T a y l o r 不稳定性和R i c h t m y e r M e s h k o v 不稳定性对于湍 流机理的研究有着重要的作用综上所述，交界面和自由面问题的研究及其数值模拟有着多方面 的重要学术意义和广泛的实际应用价值，其研究内容相当的丰富博大，其中还要很多问题尚待研 究并且

28、极具挑战性 本章将在5 1 2 节中首先概要介绍交界面数值模拟方法的研究历史和发展现状，然后在1 3 节 和1 4 节分别介绍可压缩流动和不可压缩流动中的交界面问题的数值研究，以及其中重要的模拟 方法、实际应用和尚需解决的问题关于大规模并行计算及其在交界面问题中的应用将在1 5 节 中加以介绍，最后在1 6 节详细列出了本文的主要研究方向和工作内容 1 2 交界面数值模拟方法概述 近几十年来，交界面的数值模拟由于其深刻的现实意义和广阔的应用前景而被广泛研究【8 9 ， 1 6 4 主要的交界面数值模拟方法可以分为两大类 8 9 ，5 3 ，1 7 2 ，1 6 6 】：L a g r a n

29、g e 界面跟踪方法 ( F r o n t T r a c k i n g ) 和E u l e r 界面捕捉方法( F r o n t C a p t u r i n g ) 前者能准确跟踪交界面的演化。间断 处理无数值耗散，但难于处理界面拓扑结构的变化( 如界面融合和破碎) ，且拓展到高维情况时算 法非常复杂后者计算简单，格式统一，易于处理界面拓扑改变，但不能精确定位多介质交界面， 不适合j 二对界面位置要求很高的数值模拟 交界面数值模拟方法并不局限于上述的两大类方法，可以根据所处理问题的不同而采用不同 的方法比较常用的还有格子气( L a t t i c e G a s ) 与B o

30、l t z m a n n 格子气( B o l t z m a n n L a t t i c e - G a s ) 方法 ( B e n z i1 9 9 2 1 2 1 ，R o t h m a n Z a l e s k i1 0 9 4 1 1 5 9 1 ，R o t h m a n Z a l e s k i1 9 9 7 1 1 6 0 ) ，相场( P h a s e - F i e l d ) 方 法( A n t a n o v s k i i1 9 9 5 1 1 1 】，J a c q m i n1 9 9 6 1 9 2 ，J a c q m i n1 9 9 9

31、 1 9 3 ) 等 2 第1 章绪论 1 2 1L a g r a n g e 界面跟踪( F r o n t - T r a c k i n g ) 法 L a g r a n g e 界面跟踪方法通过参数表示的闭合曲线来表征界面位置，界面是固定在这组闭合曲 线坐标系下的，界面随着表征曲线的运动而运动，所以叫做L a g r s z t g e 方法一般地，这组闭合曲 线是一些标记点，由这些标记点的运动来跟踪交界面位置的变化，所以叫做F r o n t - T r a c k i n g 方法 边界积分方法( B o u n d a r yI n t e g r a lM e t h o

32、d ) J 2 0 5 ，1 4 6 ，8 6 】也可以归入到L a g r a n g e 界面跟踪方法中来 【3 2 】L a g r a n g e 方法可以给出活动交界面的高阶准确近似，并可以避免数值扩散，而一般数值扩 散会使本来锐利的界面变得很平滑而模糊了界面的准确位置 格子类方法是一类重要的L a g r a n g e 界面跟踪方法( 图1 1 ) 其中P a r t i c l e - I n - C e 1 ( P I C ) 方法 ( H a r l o w1 9 6 4 1 7 7 ) 是最基本的格子类方法，P I C 方法把连续的流场看作为有限个质量集中的流体 质点系，

33、它们是分布在E u l e r 网格内的、具有L a g r a n g e 特征的离散点通过对这些离散点的计算、 追踪和调整，实现流场的数值模拟由G e n t r y , M a r t i n 和D a l y 提出的F L u i d I n C e l I ( F L I C ) 方法( 1 9 6 6 1 5 9 1 ) 是P I C 方法的一种重要改进，避免了质点的引入，同时精细地考虑了网格界面的输运计 算M a r k e r A n d - C e l l ( M A C ) 方法( H a r l o w W e l c h1 9 6 5 1 7 8 1 ) 采用无质量的标

34、记点( M a r k e r ) 来代 替P I C 方法中的质点系，并以之确定自由面的形状这种方法是许多现代界面追踪方法的基础 图1 1 ：L a g r g m g i e a x 界面跟踪方法( 图中用标记点显式跟踪交界面的演化) L a g r a n g i a n 界面跟踪( F r o n t - T r a c k i n g ) 方法的优点很明显， 出于是显式地跟踪界面的位置，所以界面跟踪方法模拟交界面的位置较其它方法更为精确 界面跟踪方法在交界面上使用精确的间断跳跃条件，不引入数值扩散，因而交界面没有弥散， 这对于某些宏观不混的多介质可压缩流动的模拟非常重要 由于使用了大

35、量的标记点，可以有子网格尺度的分辨率如图1 1 所示不稳定性的涡卷起一 般的界面方法在一个网格内是难以得到比交界面更精细的结构的，而界面跟踪方法可以 由于界面跟踪方法能更准确地得到交界面的位置，因而也能够更精确地计算交界面的平均曲 率和表面张力 但是这种复杂的计算方法缺点也是相当多的： 界面跟踪方法的算法非常复杂，必须考虑到各种情况，比如标记点的重新分布等 显式跟踪复杂的交界面是很困难的。特别是如果出现奇异性。界面处理的算法将会相当复杂， 甚至不能够处理界面奇异性问题而出现错误的结果 有时需要对表征交界面的M a r k e r 标记点进行专门的数值离散，因为交界面的处理是独立于流 体计算网格

36、的处理的，比如可压缩流界面处理的F r o n t - T r a c k i n g 方法就需要在交界面上求解 复杂的多维R i e m a n n 问题 难于处理拓扑结构的变化例如当两个气泡互相融合的时候。就必须考虑交界面标记点的重 1 2 交界面数值模拟方法概述 3 新分布；当一个气泡分裂的时候就需要考虑加入新的标记点；甚至当两个气泡非常接近时， 还必须考虑标记点是标记的哪一个气泡 推广到高维情况非常复杂，因为不仅是空间维数的增加而使工作量大大增加，而且以上的各 项缺点会在高维下更加明显和错综复杂 详细的有关F r o n t - T r a c k i n g 方法的实现和发展可以参考

37、G l i m m 等【3 4 ，6 6 ，6 3 ，65 ，n y 黯v a - s o n 等【1 9 3 ，1 9 6 ，1 9 5 ，和S h y y 等【1 7 1 ，1 7 2 的文章另一类使用L a g r a n g e 观点的粒子方法是 S P H ( S m o o t h e d - P a r t i d e - H y d r o d y n a m i c 8 ) 方法1 6 0 ，1 3 0 ，它也可以归结为无网格( M e s h l e s s ) 方法 17 】这种方法实现简单，适合复杂几何区域的计算，原本是从天体物理中发展起来并用于可压缩 流动模拟，近来也

38、越来越广泛地应用于不可压缩流动界面问题中【1 3 1 ，1 3 3 ，1 3 2 ；但它的精度尚不 够高，而且还有很多具体问题待解决 1 2 2E u l e r 界面捕捉( F r o n t C a p t u r i n g ) 不同于L a g r a n g e 界面跟踪方法显式地跟踪界面位置，E u l e r 界面捕捉方法是将交界面隐式地 包含在一系列特征函数中，通过这些函数的输运来表征界砸的运动，再从这些函数中获得交界面 的位置及其运动 1 V o l u m e - O f - F l u i d ( V O F ) 方法 V o l u m e O f - F l u i

39、d ( V O F ) 方法是最重要的运动界面追踪数值模拟方法之一，他的主要特点是将 运动界面在空间网格内定义成一种流体体积函数，并构造这种流体体积函数的发展方程，通过与 流体运动耦合的流体体积输运，精细地确定该运动界面的位置、形状和变化。从而达到界面追踪 的目的 在二维情况，交界面可以看成是连续的、分段光滑的直线段构成的如果已知界面网格以及相 邻网格的目标流体体积分数，交界面的位置可以通过重构来获得，就是在每个界面网格内求出一 个交界面线段的近似最简单的V O F 方法是简单直线界面计算( N o h W o o d w a r d ，1 9 7 6 1 1 3 9 ) 和 S O L A

40、V O F 算法( H i r t N i c h o l l s ，1 9 8 1 1 8 2 ) ，它们用平行。】：网格的直线来近似界面( S i m p l eL i n e I n t e r f a c eC a l c u l a t i o n ，S L I C ) ，因而在界面重构上仅一阶精度，对界面的处理还相当粗糙( 见示意 图1 2 a ) L a f a u r i ee ta 1 ( 1 9 9 4 ) 1 0 5 对这种方法做了有限的改进更精确的界面重构技术是采用 分段线性的直线和几何构造方法来重构界面，一般可以归为P L I C 方法( P i e c e w i

41、s eL i n e a rI n t e d a c e C o n s t r u c t i o n ) ( 图1 2 b ) 其中的主要界面重构算法有；A s h g r i z P o o 的F L A I R 方法( 1 9 9 1 【1 2 ) ， P a r k e r Y o u n g s 的方法( 1 9 9 2 【1 4 2 ) ，R i d e r K o t h e 的方法( 1 9 9 5 【1 5 5 ) ，P i l l i o dJ r P u c k e t t 的方法( 1 9 9 7 9 T ) 。以及G u e y f f i e r 的方法( 1

42、9 9 9 7 4 ) 等 图12 ：V O F 两种界面重构方法( a ) 一阶S L I C ( b ) 二阶P L I C 4 第1 章绪论 V O F 方法得到广泛应用，主要有以下几个方法的优点， 能够处理界面拓扑结构的复杂变化，比如交界面的融合与分裂因为V O F 方法是隐式处理交 界面。所以不必像F r o n b - D a c k i n g 方法那样复杂的算法 V O F 方法能够简单自然地将二维算法拓展到三维去。没有本质的变化 V O F 方法能够保持体积分数输运过程的守恒性，因此整个计算的流体体积( 或面积) 是守恒 的 由于体积分数的更新只依赖于相邻计算网格的体积分数值

43、，因而V O F 方法能够很容易地实 现并行算法，特别是区域分I 睥算法 但是V O F 方法对于计算交界面的平均曲率不够精确，因而当应用到依赖于曲率的问题中时会 存在较大的误差或者需要进行繁琐的细节处理 2 L e v e l - S e t 方法 L e v e l - S e t 方法首先是由O s h e r S e t h i a n ( 1 9 8 8 ) 1 4 1 提出并应用于几何界面追踪同题这种方 法的主要思想是将n 一1 维的曲线或曲面用一个n 维光滑函数的零等值面，它的一个重要特点是 可以将间断问题转化为光滑问题来加以解决一个典型应用是将运动界面( 一般是一个密度、压力

44、以及状态方程的间断面) 定义为一个光滑的l e v e l s e t 函数的零等值面，毋( x t ) = 0 然后交界面的 运动就是该零等值面的运动，零等值面的位置也就是交界面的位置( 如图1 3 所示) o o b 1 图l3 ：L e v e l - S e tM e t h o d a ) 实线部分为定义的l e v e l s e t 函数。纵轴为时间。阴影部分为初始时刻两个分离的 交界丽；b ) 随着时问演化。界面拓扑结构变化融合为一个交界面如阴影部分所示 L e v e l S e t 方法有着许多优点t L e v e l S e t 方法实现简单，同F r o n t -

45、T r a c k i n g 方法比较，L e v e l S e t 方法不需要引入专门的数 据结构来显式描述交界面的位置及其运动，而是隐含在l e v e l - s e t 函数的时间演化中了，这同 V O F 方法隐式处理界面的思想相同 L e v e l S e t 方法可以方便地处理各种复杂拓扑结构的变化，所以这些交界面的融合和分裂等奇 1 3 可压缩流动界面问题的数值研究 5 异现象都可以由l e v e l s e t 函数自然地解决，避免了F r o n t T r a c k i n g 等方法复杂的处理方法 L e v e l S e t 方法适合于交界面曲率依赖问题的

46、解决根据光滑的L e v e l - S e t 函数可以很方便地计 算交界面的平均曲率以及外法线方向。可以比较准确地计算表面张力等 另一方面L e v e l - S e t 方法也有其不足和尚需改进的地方； L e v e l - S e t 方法对于界面的处理没有F r o n t - T r a c k i n g 方法的精度高。例如对于激波同交界面的 相互作用、考虑表面张力的两相流动等 为了处理界面的间断以及计算平均曲率。需要求L e v e l - S e t 函数具有距离函数的性质，即需要 对L e v e l S e t 函数进行重新初始化( R e i n i t i a l

47、 i z a t i o n ) ，这会带来额外的计算量 实际上，L e v e l S e t 方法本身并不能在输运过程中保持流体体积守恒，将在一定程度上使计算 结果失真 N a r r o wB a n dL e v e l S e tm e t h o d 【5 ，l o o 可以将计算量由O ( N 2 ) 减少到O ( k N ) ，这里假设 为沿一个方向的网格点数，且为N a r r o wB a n d 带宽内的网格点数由于该方法可以大大降低计 算量而越来越得到广泛的研究和应用 1 3 可压缩流动界面问题的数值研究 复杂的多介质可压缩流动问题的数值模拟是近十几年来计算流体力学的重

48、要研究方向之一， 特别是考虑互不掺混的多流体界面问题的高精度计算方法类似于交界面模拟方法，多介质可压缩 流动的计算方法 4 】也可以根据其使用的界面处理方法分为两大类：F r o n t - T r a c k i n g 方法和F r o n t C a p t u r i n g 方法前者能准确跟踪交界面的演化，间断处理无数值耗散，但需要显式地在交界面上 处理界面跳跃条件，难于处理界面的拓扑结构的变化( 如界面融合和破碎) ，G l i m m 等在这方面做 出了系统的工作 3 4 ，6 3 ，6 6 ，65 】，主要是进行可压缩流动如R a y l e i g h T a y l o r 和R i c h t m y e r - M e s h k o v 不稳定性的数值模拟，但是其算法非常复杂后者计算简单，格式统一，易于处理界面拓扑变化， 但不能精确定位多介质交界面，一般地不适合于对界面位置要求很高的数值模拟界面捕捉方法 由于简单容易实现而得到了广泛的应用，早期的界面捕捉方法是借用多组分流动模拟中基于质量 分数的方法来描述不同流体，但是这样会导致在多流体交界面上的速度和压力振荡【2 1 为防止 在界面附近的压力振荡而提出了众多的数值方法，例如基于7 模型和L e v e l _ S e