矩形空腔内Stokes流的状态空间有限元法  

The state space finite element method for Stokes flow in rectangular cavity

作　　者：孟俊苗 邓子辰[1,2] 王艳[1] 

机构地区：[1]西北工业大学工程力学系,西安710129 [2]大连理工大学工业装备结构分析国家重点实验室,大连116024

出　　处：《计算力学学报》2014年第2期205-211,共7页Chinese Journal of Computational Mechanics

基　　金：国家自然科学基金(10972182;11172239;10902089);国家基础研究计划973项目(2011CB610300);国家111引智计划(B07050);高校博士点基金(20106102110019);大连理工大学工业装备结构分析国家重点实验室开放基金(GZ0802)资助项目

摘　　要：基于Hellinger-Reissner二类变分原理,从平面Stokes流问题的平衡方程、连续性要求和边界条件出发,得到相应的Hamilton函数,建立Hamilton正则方程后,采用分离变量法对场变量进行离散求解:在x方向采用有限元插值,在y方向采用状态空间法给出控制坐标方向的解析解。计算过程中的指数矩阵均采用精细积分法求解,使得本文算法具有高效率、高精度、对步长不敏感的优点。通过对侧边自由液面边界条件的单板驱动矩形空腔Stokes流问题的求解,得到与文献相同的结果,从而验证了本文方法的有效性。本文旨在将弹性力学状态空间有限元法的思想引入到低雷诺数流体力学中,为Hamilton体系下研究复杂边界Stokes流问题提供新的途径。Based on the Hellinger-Reissner variational principle, the Hamilton canonical equation of the plane incompressible Stokes flow was derived from the equilibrium equations, continuity conditions and the force boundary conditions. By the separation of variables, the general finite element method was employed in x direction and was derived by the state space control method. In addition, the precise integra- tion method for the exponential matrix was employed in the calculation. The effectiveness of the state space finite element method is demonstrated by analyzing and comparing the simulation example in case of a single-lid driven cavity with free surface side wails. The study of this paper is to introduce the idea of the semi-analytical method into the low Reynolds number flow problems, and lay a foundation of further study on the Stokes flow with complex boundary in the Hamiltonian system.

关 键 词：Hellinger-Reissner变分原理 STOKES流 正则方程 状态空间法 精细积分 

分 类 号：O357.1[理学—流体力学]