第六章变几何问题连续和移动边界.docVIP

第六章 变几何问题：和移动边界V.R.GUNDABALA, W.B.J ZIMMERMAN, A.F.ROUTH 1Department of Chemical and Process Engineering, University of Sheffield, Newcastle Street, Sheffield S1 3JD United Kingdom 2Department of Chemical Engineering, Cambridge University, Pembroke Street, Cambridge, CB2 3RA E-mail: w.zimmerman@shef.ac.uk 求解过程中，由于几何模型发生变化而导致区域网格改变时会出现几何结构的问题。本章中，我们举两个例子流道中由于不同尺寸孔板引起额外压力损失问题是稳态几何的一个例子。从概念来看，这个问题第5章贝纳尔问题中通过Rayleigh数进行。第二个例子是液体中悬浮的乳胶颗粒的干膜。在这个问题中，考虑了两个移动前缘；一个通过坐标变换，另一个通过光滑界面模型。该技术对于之前模拟该问题采用的移动弱项方法有所改进。1. 引言 1.1 几何连续性 我们已经有了一些参数连续性的例子通过参数在一定范围内一系列的较小变化，以相邻参数值的解作为新参数求解初值。只要参数不通过歧点且参数步长足够小，即可保证从旧解到新解的平滑过渡。即使存在歧点，旧的分支仍然是一个可行解，参见我们在第5章中对贝纳尔对流问题的讨论。几何连续性与参数连续性在一个重要方面存在性质上面的不同。在几何连续性中，区域几何结构的改变使得需要对网格重新划分。我们应该小心的区分几何连续性的改变。例如：在管流中，众所周知流动可由雷诺数表征： (1) 该无量纲数包含了流体密度ρ，入口速度U直径D和粘度μ可用于描述流体的动力相似性。因此，管道直径改变引起充分发展流并不属于几何连续性的范围，而更接近于常规的参数改变问题。 与上一章类似，模拟的例子给出了稳态和瞬态的模型，本章我们将给出稳态几何连续性与瞬态几何连续性的例子。对于前者，不同的模型所用的几何结构不同，因此网格会略有不同。因此，一个参数到下一个参数的解是不兼容的不同的网格自由度。如果旧的解作为新几何结构的初始值，则需要将旧的解在保持一致性的情况下映射到新的区域。在这里的例子中，稳态模型PDE方程组是线性的，因此解可以直接由1FEM步获得。因此，将旧的解映射到新几何结构没有附加值。对于瞬态问题，随着移动的前缘，涉及到区域减小的问题。该区域在每一时间步后均发生改变，因此将上一时间步的解映射到新的区域是必要的。这导致在每一时间步后必须重新划分网格。对于自由边界问题这是关键的一步。例如膜流动和喷气流，边界的位置和速度场的解密切相关，边界应位于能够满足压力平衡的位置。 二维不可压层流1.2 自由边界问题和任意拉格朗日-欧拉变换 计算流体动力学2. 稳态几何连续性：孔板引起的流道压降 ε时，保持恒定体积流率所需的额外压降。图中绘制了当板阻塞率为40%时所生成的典型网格。由于有必要分辨孔板周围的微小特征，孔板周围的网格必须保持足够致密，网格加密成为一个重要问题。我们将从几何参数改变来外推所需的附加压降，而不是直接求解高阻塞率问题。 图 1. 粘性流管道中孔板问题生成的网格。用于表示阻塞百分数的参数为ε＝0.4 尽管有可能考虑孔板周围任意雷诺数的计算，层流的主要影响与人为消去的雷诺数-斯托克斯方程相似。该现象的根由在于大多数的耗散发生在孔板缝隙之间，在该位置，流体加速，然而小的缝隙导致了较强的粘性摩擦，可近似根据斯托克斯方程模拟动量传递： (2) 其中，μ为流体粘度，ρ为E须满足： (3) 其中，Q为体积流率，W为功流。c为未知比例常数，该常数数值由Dugdale根据吸附能极值从理论上进行了计算，也能够通过测量Q和压降从实验上获得。在二维体系中，通过相似量纲分析可获得单位长度耗散损失和截面面积流率，从而可得： (4) Dugdale报道了通过糖浆实验确定常数c范围在3.17至3.30之间，他的理论计算结果为3.30。Bond[6]给出了孔板问题与管长为2ka的哈根-伯叔叶管流问题的相似性，其中a为孔板半径，压降等值为k＝0.631，表明c=3.21. 我们曾对管道中紧合粒子的拖曳较感兴趣。对于与(3)(4)相同的理论基础，管道中的颗粒由狭缝宽度决定。Zimmerman对于薄片[7]（宽边方向的运动）以及球体[8] 在圆柱管中的沉积做了相关研究，表明当颗粒半径较大时曳力增大较快（狭缝宽度为a）。对小颗粒半径(1-a)采用摄动法并将级数展开相加，有可能确定当颗粒接近管壁时问题的奇异性。通过对薄片宽