EI、Scopus 收录

 引用本文: 段庆林 李锡夔. 不可压缩N-S方程的稳定分步算法及耦合离散[J]. 力学学报, 2007, 23(6): 749-759.
Qinglin Duan, Xikui Li. Pressure stabilized fractional step algorithm for incompressible N-S Equations and coupled finite element and meshfree discretization[J]. Chinese Journal of Theoretical and Applied Mechanics, 2007, 23(6): 749-759.
 Citation: Qinglin Duan, Xikui Li. Pressure stabilized fractional step algorithm for incompressible N-S Equations and coupled finite element and meshfree discretization[J]. Chinese Journal of Theoretical and Applied Mechanics, 2007, 23(6): 749-759.

Pressure stabilized fractional step algorithm for incompressible N-S Equations and coupled finite element and meshfree discretization

• 摘要: 在有限增量微积分(finite increment calculus,FIC)的理论框架下，通过引入一个附加变量，发展了压力稳定型分步算法，有效改善了经典分步算法的压力稳定性，同时还避免了标准FIC方法中存在的空间高阶导数的计算. 为保证数值方法同时具有较快的计算速度和较好的健壮性，发展了有限元与无网格的耦合空间离散方法. 该方案可在网格发生扭曲的区域采用无网格法空间离散以保证求解的精度和稳定性，而在网格质量较好的区域以及本质边界上保留使用有限元法空间离散以提高计算效率和便于施加本质边界条件. 方腔流考题的数值模拟结果突出地显示了所发展的压力稳定型分步算法比经典分步算法具有更好的压力稳定性，能够有效消除速度-压力插值空间违反LBB条件而导致的压力场的虚假数值振荡. 平面Poisseuille流动和一个典型型腔充填过程的数值模拟结果, 表明了发展的耦合离散方案相对于单一的有限元法和单一的无网格法在综合考虑计算效率和算法健壮性方面的突出优点.

Abstract: In virtue of an additional variable in the framework ofthe Finite Increment Calculus(FIC) theory, a pressure stabilized fractionalstep algorithm is developed in this paper with enhanced pressure stabilityin comparison with the classic one. In addition, the calculation of the highorder spatial derivatives which exists in the standard FIC procedure is alsoavoided. To ensure superior overall performance of the proposed numericalscheme in accuracy, efficiency and robustness, a coupled finite element andmeshfree method is developed for the spatial discretization andinterpolation approximation, in which the meshfree approximation is adoptedin the region where the mesh is distorted to preserve the accuracy androbustness of numerical solutions from the deterioration of the meshquality, while the finite element approximation is employed in the regionwhere the quality of the mesh is acceptable and on the boundaries whereessential boundary conditions of flow problems are imposed to ensure highcomputational efficiency and proper imposition of the essential boundaryconditions. Numerical results for the lid-driven cavity flow problemdemonstrate the better pressure stability of the proposed pressurestabilized fractional step algorithm than that of the classic one, and itscapability in removing the spurious oscillations in the resulting pressurefield induced by the incompatible interpolation approximations for thevelocity and pressure fields violating the LBB condition. The two exampleproblems, i.e. the plane Poisseuille flow and the injection molding problemsare illustrated to prominently demonstrate the superiority of the proposedcoupled finite element and meshfree method over the independent finiteelement and meshfree methods in the overall performance.

/

• 分享
• 用微信扫码二维码

分享至好友和朋友圈