基于高效流固耦合方法的自由表面流动问题分析.pdfVIP

2025年第25卷第4期科学技术与工程ISSN1671—1815

2025,25(4):01325-07ScienceTechnologyandEngineeringCN11—4688/T

DOI:10.12404/j.issn.1671-1815.2401914

引用格式:温喆然,徐自力,薛杰,等.基于高效流固耦合方法的自由表面流动问题分析[J].科学技术与工程,2025,25(4):1325-1331.

WenZheran,XuZili,XueJie,etal.Analysisoffreesurfaceflowproblemsbasedonanefficientfluid-structurecouplingmethod[J].Science

TechnologyandEngineering,2025,25(4):1325-1331.

力学

基于高效流固耦合方法的自由表面流动问题分析

11∗222

温喆然,徐自力,薛杰,王旭阳,王珺

(1.西安交通大学航天航空学院复杂服役环境重大装备结构强度与寿命全国重点实验室,西安710049;

2.西安航天动力研究所航天液体动力全国重点实验室,西安710100)

摘要针对自由表面流动与弹性结构的流固耦合计算效率低、计算耗时长的问题,将流体体积法与基于结构-虚拟弹性体的

快速动网格方法相结合,发展了一种适用于自由表面流动的高效流固耦合方法。使用流体体积(volumeoffluid,VOF)法对流

体自由表面进行追踪;将流体域视为虚拟弹性体并构建结构-虚拟弹性体系统,以流固耦合界面的多相流体力为激励求解系统

的动力学方程得到结构振动位移和流场网格变形;在每一个时间步内依次求解流体流动、结构变形和流场动网格,实现流固

耦合计算。基于发展的方法计算了溃坝水流冲击下弹性挡板的流固耦合响应,得到了溃坝水流的自由液面和弹性挡板的运

动行为,结果表明:自由液面演变和弹性挡板振动位移的计算结果与已有算法的结果吻合良好;在同等网格规模下,与已有算

法相比本文方法可减少33.3%的计算时间;在水流冲击作用下,弹性挡板向冲击侧小幅弯曲。随后水流沿挡板左侧上升并形

成射流,挡板向另一侧大幅弯曲。最后由于两侧流体的阻尼,挡板振幅逐渐衰减。

关键词流固耦合;自由表面流动;流体体积(VOF)法;快速动网格方法;溃坝水流

中图法分类号O353.4;文献标志码A

AnalysisofFreeSurfaceFlowProblemsBasedonanEfficient

Fluid-structureCouplingMethod

11∗222

WENZhe-ran,XUZi-li,XUEJie,WANGXu-yang,WANGJun

(1.StateKeyLaboratoryforStrengthandVibrationofMechanicalStructures,SchoolofAerospaceEngineering,

XianJiaotongUniversity,Xian710049,China;

2.NationalKeyLaboratoryofAerospaceLiquidPropulsion,XianAerospacePropulsionInstitute,Xian710100,C