Euler方法中多物质界面计算的模糊处理研究.pdfVIP

第八届中国青年运筹信息管理学者大套论文集 桂林，2006年8月18-22日，第780-788页 Euler方法中多物质界面计算的模糊处理4 陈龙伟 云南财经大学计算机科学系，昆明650221 摘要针对三维欧拉(Eulcr)方法中多物质界面处理的难题，基于流体网格(VOF)提出模糊 界面处理方法，给…压力计算、输运体积计算、输运的具体算法，并对三维三种介质的爆炸场 进行¨r数值模拟，说明模糊界面方法简单、适用。 关键词欧拉方法，混合网格，模糊方法，界面处理，数值模拟 1引言 糊理论的研究和实际应用蓬勃发展，取得了许多成果。近年来，模糊数学在数值计算方面 也有广泛的应用。 在七十年代中期以前，计算流体力学一般采用的是有限差分方法，之后有限元方法被 引入。按照所选用的空间坐标分类，无论有限元方法还是有限差分方法都可以采用 法的代表。 LagTange坐标随物质一起运动和变形，处理自由面和物质界面是直观而自然的，同 时由于一个网格或单元始终对应～块物质团，因此能较准确地描述不同部分丰才料的不同应 力历程，允许对不同部分的材料采用不同的本构关系。其缺点是对物质大变形会发生畸变、 扭曲从而导致诗算不稳定。另外在处理层裂、挤凿等材料破坏问题时，需要引入滑移面等 特别解决办法。然而，在Lagrange法中速度离散后的值一般取在网格角点上的值，这意 味着假定网格角点处的值连续，引入断裂和滑移面后相邻网格要分离或相对运动，这样网 格角点处的速度就不能再连续了。对这一问题的处理有许多办法．如主从界面法 (master—slave 运产生间断面那样自然“1。 Euler坐标是固定的空间坐标系，计算网格固定于几何空间，网格是不变形的，物质 通过网格边界流进流出，即使发生大畸变，也不会引起如Lagrange计算中的误差灾变和 不能使计算进行下去的困难。材料变形也不直接影响时间步长的计算。因此，Euler方法 能处理大变形问题。但是，采用空间坐标系难于跟踪物质的变形，也不能识别材料界面的 位形。在计算过程中有混合格问题，不能精确处理界面，使得Euler方法的发展受到限制。 混合格问题是Euler方法面临的主要问题之一”1。 ’云南财经大学科研基金资助(2004299201) Euler方法中多物质界面计算的模糊处理 Euler方法一般采用分裂算法。在体力不起作用的问题中，格元物理量的局部变化取 决于通过格元界面的对流、压力和应力偏量作用等三部分效应。分裂算法是对格元物理量 的先后按照三部分效应逐步计算相应的改变量。由于计算对流与应力作用引起的变形是分 r3】 开进行的，而且此种变形计算并没有跟踪物质，因此必然会带来误差。1982年McGlaun 开始在CSQ程序中，在计算对流的一步中放弃使用l阶精度格式，改用二阶精度格式的van Leer的守恒格式”3，对提高计算精度较为有益。近几年来，几乎所有的Euler形程序均在 这方面作了改进。在计算多种材料占据的混合网格或者材料末占满的自由面格中材料的对 流时，Euler型程序有一共同的缺点，就是引进了数值上的扩散，特别是自由面附近的材 和Elrick51利用本格元和相邻格元中各种材料所占的体积分数值，比较方便而又精细地确 定格元中材料界面的大致位置(用斜线段表示)，从而改进了通过格元侧面的输运量的计 算，大大减少了数值扩散。 这两种方法各有优缺点，人们想到把两种方法的优点结合起来在Euler法中吸收 表的ALE法(hrbitraryLagrangianEulerrian)为典型。当今的趋势是研制兼有两者优 点的程序，如CEL方法等。目前，运用Lagrange方法的流体动力学程序引入了许多新的 算法，吸收了Euler方法的一些思想，得到了很大发展，所以目前比较流行，相比之下， 运用Euler方法的程序的发展较为落后。因为Euler法在处理两种以上的混合网格时输运 了一些算法，在一定程度上解决了网格畸变等问题，但与Euler法相比，不可能象Euler 法那样自然的反映大变形。另外，由于一些高精度的界面跟踪方法和输运量计算方法的引 入，Euler法的应用已经有了很大的发展“1。 Euler方法中的多物质界面描述，以前人们所采用的方法，大多数属于确定性数学方 法或模型，如Level