清华大学计算固体力学第七次讲义ALE公式.pptVIP

作为一个使用修正的弹性方程的ALE网格更新的例子，考虑一个位于一个矩形流体域中的圆柱绕流的有限元网格。在圆柱沿着y方向移动了一个位移0.25w后的网格见图b。矩形域的边界保持固定。由图可见，圆柱附近的网格精度在更新后的网格中得以保持，并且无明显的单元扭曲发生。 网格更新的例子 6 网格更新算法 7 Petrov-Galerkin方法 伽辽金(Galerkin)方法，是利用满足位移和应力边界条件 的函数寻求积分方程的解答。 关键问题是增加粘性，消除不稳定项，保证数值稳定性。 建立Petrov-Galerkin方法的迎风流线(Streamline Upwind Petrov-Galerkin－SUPG)公式。对流－扩散方程是一个有用的方法，它对应于动量方程的线性化。对于离散的稳态对流—扩散方程，将得到闭合解答。将证明当网格参数(已知的Peclet数)超过临界值时，这个解答在空间是振荡的。通过建立P-G方法以消除这些振荡，即纠正不稳定性。 在一维中，离散方程类似于标准迎风方程。然而在多维中，它们提供了沿着流线引导迎风项的一致理论框架。 迎风格式的基本思想：当用差分方程求解偏微分方程时，利用特征线方向一侧的单边差商来代替空间偏导数。目的是保证数值稳定性。如一阶线性常系数双曲型方程： 7 Petrov-Galerkin方法 在特征线方向一侧的单