流固耦合理论研究述评2—流固耦合理论基础及关键问题.pdfVIP

第七届全国现代结构工程学术研讨会 流固耦合理论研究述评(2) ’一流固耦合理论基础及关键问题 袁行飞1 苏波2聂国隽3j钱若军2 (1浙江大学土木工程学院，杭州310027)(2同济大学士术工程学院．上海200092)(3同济大学航空航天与力学学院，上海200092X 摘要：流固耦台(FsI)是目前很多领域研究的热点和也是难点之一．目前，国内研究埘此相对滞后或比较片面。为在倒一 软件平台上实现多种耦合算法，本文对流固耦合理论进行了研究。内容包括FsI研究进展：流体力学基车理论和流体计算力 学力学(cFD)有限元法；耦台界面条件、边界追踪、网格更新、以及大型非线性方程的求解等诸多问题；特别是对目前广为 关注和讨论的强耦合和弱耦合算法作了理论上的探讨和区分．本文为流圊耦合理论数值模拟提供理论和方法的指导。 关键词：流固耦合．强耦合，弱耦合，边界条件，札E法 三、流体力学基本理论 (一)纳维一斯托克斯O卜s)方程组 流体运动虽然复杂，但必须满足三大守恒定律，即质量守恒律、动量平衡律、能量守恒律。采用欧拉 描述方法，对流体微元体蚴dz应用这三大定律，可得最根本的流体运动控制方程，1即N·s方程组。其表 达形式可以是守恒型或非守恒型，限于篇幅，这里只写出N．s方程组守恒形式‘11： 挈Ⅳ(pV)圳质量蛔 (3la) 等+V．(pw—t)=^动量守．回 (31b) (3．1c) 曼譬+v．(pvE—T．v+q)：fa．v+矿(能量守恒) 出 式中，f为时间，p为流体密度，v为速度矢量，f。为流体域的体积力矢量，f为应力张量，E 为总能量，矿为热量生成速率。 导系数，g为重量加速度。 对式(3．1)中的各方程，可以用～个通式表示，并可写出求解该方程的源程序，对于不同的妒值，只 要重复调用该程序，就可以求解不同类型的流体流动及传热问题。 n保守形式： 710 工业建筑2007增刊 第七届全国现代结构工程学术研讨会 巫笋砌(删珈(珊删)+s (32) 式cfl矿为通用变量，可以表示速度u，温度丁等求解变量；r为广义扩散系数：S为广义源项。 2)非保守形式： 若将式(2．2)的瞬态项和对流项中的物理量从微分符号中移出，便可得到通用控制方程的非守恒形式： P(警+“警+V等十w警)=折(rg蒯矿)+sP【盖州盖¨茜+w舌J-批(rg删纠c，圆(3’3’ 需要注意的是N—s方程本身并不闭合，对之求解需要外加本构方程和状态方程。 (二)湍流模型 流体运动状态可分为层流和湍流。湍流(tllrbulence)是一种高度复杂的非稳态、带旋转的不规则流动。 在湍流中流体的各种物理参数，如速度、压力、温度等都随时间与空间发生随机变化。湍流的计算流体力 学模型大致可分为四类： 1、直接数值模拟DNS(DirectNumcIicalsi删lation) 从理论上讲，N．s方程可以描述所有的运动。如直接根据N_s方程来计算湍流．湍流模拟方法称为直接 数值模拟。但目前计算机的存储能力和运算速度，使得DNs模型难以直接应用于实际工程。 Navier．Stokcs 2、雷诺平均N—S方程队NS(Reynoldsavemged Equad蚰) 此方法考虑物理量的瞬时值庐，由时均值≯及脉动值∥组成： 争=季+∥ 时流体的应力与应变率的本构方程，引入湍流粘性系数“，将湍湍流脉动附