针对自由表面流动与弹性结构的流固耦合计算效率低、计算耗时长的问题,将流体体积法与基于结构-虚拟弹性体的快速动网格方法相结合,发展了一种适用于自由表面流动的高效流固耦合方法。使用流体体积(volume of fluid,VOF)法对流体自由表...针对自由表面流动与弹性结构的流固耦合计算效率低、计算耗时长的问题,将流体体积法与基于结构-虚拟弹性体的快速动网格方法相结合,发展了一种适用于自由表面流动的高效流固耦合方法。使用流体体积(volume of fluid,VOF)法对流体自由表面进行追踪;将流体域视为虚拟弹性体并构建结构-虚拟弹性体系统,以流固耦合界面的多相流体力为激励求解系统的动力学方程得到结构振动位移和流场网格变形;在每一个时间步内依次求解流体流动、结构变形和流场动网格,实现流固耦合计算。基于发展的方法计算了溃坝水流冲击下弹性挡板的流固耦合响应,得到了溃坝水流的自由液面和弹性挡板的运动行为,结果表明:自由液面演变和弹性挡板振动位移的计算结果与已有算法的结果吻合良好;在同等网格规模下,与已有算法相比本文方法可减少33.3%的计算时间;在水流冲击作用下,弹性挡板向冲击侧小幅弯曲。随后水流沿挡板左侧上升并形成射流,挡板向另一侧大幅弯曲。最后由于两侧流体的阻尼,挡板振幅逐渐衰减。展开更多
基于数值方法,采用流体体积函数模型(volume of fluid,VOF)对10-4g0和g0重力环境下水平方管内空气-水两相流和制冷剂R134a蒸汽-液体两相流进行数值模拟,分别得到泡状流、弹状流、搅混流和环状流4种典型流型,但两种混合物在流型上存在较...基于数值方法,采用流体体积函数模型(volume of fluid,VOF)对10-4g0和g0重力环境下水平方管内空气-水两相流和制冷剂R134a蒸汽-液体两相流进行数值模拟,分别得到泡状流、弹状流、搅混流和环状流4种典型流型,但两种混合物在流型上存在较大差异。通过对数值结果的统计分析,得到两种混合物在不同重力环境下的压降分布。结果显示,微重力下两种混合物的压降均大于常重力环境,且压降都随气、液速度的增大而增大;相同工况下,空气-水的压降大于R134a蒸汽-液体两相流的压降。将得到的压降数值结果与均相流模型、Friedel模型和Chisholm模型依次进行对比。重新根据分液相雷诺数(Reynolds)将流动分为层流区、过渡区和紊流区,并对Chisholm关系式进行了修正。结果显示,修正后的压降模型能较好地预测微重力环境下的气液两相流动压降。根据汽液两相流动特性,分析了发生以上现象的原因。展开更多
文摘针对自由表面流动与弹性结构的流固耦合计算效率低、计算耗时长的问题,将流体体积法与基于结构-虚拟弹性体的快速动网格方法相结合,发展了一种适用于自由表面流动的高效流固耦合方法。使用流体体积(volume of fluid,VOF)法对流体自由表面进行追踪;将流体域视为虚拟弹性体并构建结构-虚拟弹性体系统,以流固耦合界面的多相流体力为激励求解系统的动力学方程得到结构振动位移和流场网格变形;在每一个时间步内依次求解流体流动、结构变形和流场动网格,实现流固耦合计算。基于发展的方法计算了溃坝水流冲击下弹性挡板的流固耦合响应,得到了溃坝水流的自由液面和弹性挡板的运动行为,结果表明:自由液面演变和弹性挡板振动位移的计算结果与已有算法的结果吻合良好;在同等网格规模下,与已有算法相比本文方法可减少33.3%的计算时间;在水流冲击作用下,弹性挡板向冲击侧小幅弯曲。随后水流沿挡板左侧上升并形成射流,挡板向另一侧大幅弯曲。最后由于两侧流体的阻尼,挡板振幅逐渐衰减。
文摘基于数值方法,采用流体体积函数模型(volume of fluid,VOF)对10-4g0和g0重力环境下水平方管内空气-水两相流和制冷剂R134a蒸汽-液体两相流进行数值模拟,分别得到泡状流、弹状流、搅混流和环状流4种典型流型,但两种混合物在流型上存在较大差异。通过对数值结果的统计分析,得到两种混合物在不同重力环境下的压降分布。结果显示,微重力下两种混合物的压降均大于常重力环境,且压降都随气、液速度的增大而增大;相同工况下,空气-水的压降大于R134a蒸汽-液体两相流的压降。将得到的压降数值结果与均相流模型、Friedel模型和Chisholm模型依次进行对比。重新根据分液相雷诺数(Reynolds)将流动分为层流区、过渡区和紊流区,并对Chisholm关系式进行了修正。结果显示,修正后的压降模型能较好地预测微重力环境下的气液两相流动压降。根据汽液两相流动特性,分析了发生以上现象的原因。
