第六章 变几何问题.docVIP

第六章 变几何问题：连续和移动边界 V.R.GUNDABALA, W.B.J ZIMMERMAN, A.F.ROUTH 1Department of Chemical and Process Engineering, University of Newcastle Street 2Department of Chemical Engineering, Cambridge University, Cambridge E-mail: w.zimmerman@shef.ac.uk 求解过程中，由于几何模型发生变化而导致区域网格改变时会出现几何结构的连续性问题。本章中，我们举两个例子——流道中由于不同尺寸孔板引起额外压力损失问题是稳态几何连续性的一个例子。从概念来看，这个问题与第5章贝纳尔问题中通过Rayleigh数进行的变量连续性没什么不同。第二个例子是液体中悬浮的乳胶颗粒的干膜。在这个问题中，考虑了两个移动前缘；一个通过坐标变换，另一个通过光滑界面模型。该技术对于之前模拟该问题采用的移动弱项方法有所改进。 1. 引言 1.1 几何连续性 我们已经有了一些参数连续性的例子——通过参数在一定范围内一系列的较小变化，以相邻参数值的解作为新参数求解初值。只要参数不通过歧点且参数步长足够小，即可保证从旧解到新解的平滑过渡。即使存在歧点，旧的分支仍然是一个可行解，参见我们在第5章中对贝纳尔对流问题的讨论。 几何连续性与参数连续性在一个重要方面存在性质上面的不同。在几何连续性中，区域几何结构的改变使得需要对网格重新划分。我们应该小心的区分几何连续性的改变。例如：在管流中，众所周知流动可由雷诺数表征： (1) 该无量纲数包含了流体密度ρ，入口速度U，直径D，和粘度μ，可用于描述流体的动力相似性。因此，管道直径改变引起充分发展流并不属于几何连续性的范围，而更接近于常规的参数改变问题。 与上一章类似，模拟的例子给出了稳态和瞬态的模型，本章我们将给出稳态几何连续性与瞬态几何连续性的例子。对于前者，不同的模型所用的几何结构不同，因此网格会略有不同。因此，一个参数到下一个参数的解是不兼容的（不同的网格自由度）。如果旧的解作为新几何结构的初始值，则需要将旧的解在保持一致性的情况下映射到新的区域。在这里的例子中，稳态模型PDE方程组是线性的，因此解可以直接由1FEM步获得。因此，将旧的解映射到新几何结构没有附加值。对于瞬态问题，随着移动的前缘，涉及到区域减小的问题。该区域在每一时间步后均发生改变，因此将上一时间步的解映射到新的区域是必要的。这导致在每一时间步后必须重新划分网格。对于自由边界问题这是关键的一步。例如膜流动和喷气流，边界的位置和速度场的解密切相关，边界应位于能够满足压力平衡的位置。 2维二维不可压、层流纳维-斯托克斯方程能以一系列标准方法求解（有限差分，有限元，谱元，格子-波尔兹曼以及多重网格技术），对于固体边界条件以及复杂几何形状已能够由商用标准模拟引擎进行计算。标准计算流体动力学包有两个标准引擎：(1) 能够满足复杂几何形状的网格生成器，(2) 能够求解包含压力项（如纳维-斯托克斯方程中对于连续性方程的拉格朗日乘子）通用传递方程组的偏微分方程求解器。这两步通常分别进行，先生成网格，后进行求解，FMELAB软件在这方面并无不同。 1.2 自由边界问题和任意拉格朗日-欧拉变换 计算流体动力学对于自由边界问题处理的并非很好。可以考虑耦合流动解和网格生成的迭代算法，但是自动化的标准计算包很难实现。Ruschak[1] 描述了使用网格自适应方法实现压力边界条件的标准方法。Goodwin和Schowalter[2]成功应用于对毛细管粘度喷射的处理之中对网格位置与流体方程和边界条件的同时求解，采用有限元方法对模型方程进行离散并采用牛顿迭代法求解。FIDAP，通常用于处理自由边界流体，使用迭代流体求解/椭圆网格重生成方法，而不是同时牛顿迭代。在缩小了的区域中的瞬态模型，我们采用坐标变换和平滑界面的方法以改变整个时间域上的几何域及两个移动边界条件。通常对于考虑多相流的模拟，由于平滑界面模型允许拓扑结构的改变，因此对于移动边界模型较为适用。本书第8章以水平集合方法求解液滴结合的问题。 原则上，COMSOL Multiphysics软件也能够采用ALE模式（任意拉格朗日-欧拉变换）求解后者问题。ALE模式可通过COMSOL Multiphysics中集成的方程来实现，随着移动网格位置剩余而增加。在所有的应用模式下，包含这些附加项的标准化方法均为“激活”网格。例如，允许给定自由边界是对上一版本FEMLAB的改进，但涉及到自由边界计算的问题较少，只有少数如考虑到表面张力引起的流动等问题。实际上，对FEMLAB软件包的这一改动并非主要针对这个问题。