基于流体并行计算的二维柱体涡激振动数值模拟.pdfVIP

基于流体并行计算的二维柱体涡激振动数值模拟 郑德乾，顾明，张爱社 ( 同济大学 土木工程防灾国家重点实验室,上海 200092) 摘要：基于流体软件 Fluent 6 ，采用松耦合方法搭建了求解气动弹性问题的数值计 算平台，其中流体域可采用串行或并行计算，结构域采用 New-mark 法编程，流体 域和结构域界面的耦合采用动网格技术实现。然后，主要介绍了本文求解气动弹性 问题数值计算平台的计算流程和流体域采用并行计算时的特别处理方法。最后，采 用该计算平台进行了 Re=22000 时二维方柱横风向涡激振动的数值模拟，通过与已 有文献结果的比较验证了其有效性。 关键词：二维方柱；气动弹性；涡激振动；数值模拟 Numerical Simulation of Vortex-Induced Vibration for 2D Square Cylinder Based on Parallel Computation for Fluid Domain ZHENG Deqian, GU Ming, ZHANG Aishe (State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China) Abstract: Based on Fluent 6, numerical simulation code for solving aeroelastic problems is first presented by using loosely coupled method in this paper. The fluid domain can be solved by adopting serial or parallel computation method, and the structure domain is solved by adopting New-mark method, and finally the interface coupling between the fluid and structure domains is implemented through dynamic mesh technique. The calculation process and the special method to deal with the parallel computation in fluid domain are then introduced. Finally, numerical simulation of vortex-excited vibration of a 2D square cylinder at Reynolds number of 22000 is done. The computed results and the comparison with those in literatures could verify the validity of the present program. Keywords: 2D square cylinder; aeroelasticity; vortex-induced vibration; numerical simulation 1 引言 风作用在结构上时，就会引起结构的振动。当结构的振动足够大时，气动弹性作用可能 就不能忽略，此时就应该采用气动弹性力学模型和方法来分析所感兴趣的问题[1,2] 。气动弹 性问题是流固耦合(FSI，Fluid Structure Interaction) 问题的一个重要部分。流固耦合问题可以 出现在许多领域的工程问题中如航空工程、土木工程、机械工程等。数值求解流固耦合问题， 一般采用两种耦合解法：紧耦合法[3]和松耦合法[4,5] 。紧耦合方法是把流体和结构看作通过耦 合界面连接的单一连续介质，用统一的算子来描述控制方程。这种方法概念清晰，主要用于 流体－结构耦合场的分析。松耦合方法是在每一个时间步内分别对流体方程和结构方程依次 求解，然后通