摘要：针对流固耦合问题,发展了基于浸入边界-多松弛时间格子玻尔兹曼通量求解法（immersed boundary method multi-relaxation-time lattice Boltzmann flux solver,IB-MRT-LBFS）的弱耦合算法.依据多尺度Chapman-Enskog展开,建立不可压宏观方程状态变量和通量与格子玻尔兹曼方程中粒子密度分布函数之间的关系;采用强制浸入边界法处理流固界面使固壁表面满足无滑移边界条件,根据修正的速度求解动量方程力源项;结构运动方程采用四阶龙格-库塔法求解.格子模型与浸入边界法的引入使流固耦合计算可以在笛卡尔网格下进行,无需生成贴体网格及运用动网格技术,简化了计算过程.数值模拟了单圆柱横向涡激振动、单圆柱及串列双圆柱双自由度涡激振动问题.结果表明,IB-MRT-LBFS能够准确预测圆柱涡激振动的锁定区间、振动响应、受力情况以及捕捉尾流场结构形态,验证了该算法在求解流固耦合问题的有效性和可行性.

关键词：流固耦合 浸入边界法 格子玻尔兹曼通量求解法 涡激振动 

单位：哈尔滨工业大学能源科学与工程学院; 哈尔滨150001