isoAdvector是一种新的几何VOF(Volume of Fluid)方法,虽然克服了传统几何VOF方法难以适用于三维空间任意多面体网格的缺点,但不能直接用来模拟涉及动网格技术的液舱晃荡.为此,引入了运动通量修正,并提出了面-界面交线运动修正,使得修...isoAdvector是一种新的几何VOF(Volume of Fluid)方法,虽然克服了传统几何VOF方法难以适用于三维空间任意多面体网格的缺点,但不能直接用来模拟涉及动网格技术的液舱晃荡.为此,引入了运动通量修正,并提出了面-界面交线运动修正,使得修正后的isoAdvector方法可以应用到液舱晃荡的模拟中.基于不同的VOF方法对非共振、共振受迫晃荡和单次冲击波面进行数值模拟,并将模拟结果与试验结果以及解析解进行了比较.结果表明:相对于代数VOF方法,采用修正后的isoAdvector方法获得的自由液面位置和整体水动力载荷精度更高;捕捉的波面没有褶皱,能够较好地模拟波面的翻卷和破碎.此外,提出了界面厚度的估计方法,分析了自由液面波高精度提高的原因.展开更多
针对自由表面流动与弹性结构的流固耦合计算效率低、计算耗时长的问题,将流体体积法与基于结构-虚拟弹性体的快速动网格方法相结合,发展了一种适用于自由表面流动的高效流固耦合方法。使用流体体积(volume of fluid,VOF)法对流体自由表...针对自由表面流动与弹性结构的流固耦合计算效率低、计算耗时长的问题,将流体体积法与基于结构-虚拟弹性体的快速动网格方法相结合,发展了一种适用于自由表面流动的高效流固耦合方法。使用流体体积(volume of fluid,VOF)法对流体自由表面进行追踪;将流体域视为虚拟弹性体并构建结构-虚拟弹性体系统,以流固耦合界面的多相流体力为激励求解系统的动力学方程得到结构振动位移和流场网格变形;在每一个时间步内依次求解流体流动、结构变形和流场动网格,实现流固耦合计算。基于发展的方法计算了溃坝水流冲击下弹性挡板的流固耦合响应,得到了溃坝水流的自由液面和弹性挡板的运动行为,结果表明:自由液面演变和弹性挡板振动位移的计算结果与已有算法的结果吻合良好;在同等网格规模下,与已有算法相比本文方法可减少33.3%的计算时间;在水流冲击作用下,弹性挡板向冲击侧小幅弯曲。随后水流沿挡板左侧上升并形成射流,挡板向另一侧大幅弯曲。最后由于两侧流体的阻尼,挡板振幅逐渐衰减。展开更多
基于大涡模拟(large eddy simulation,LES)和流体体积(volume of fluid,VOF)法对同轴双通道气-液射流进行了3D数值模拟,用LES法求解气液相的流场,VOF法捕捉运动的气液界面.数值模拟结果表明,随着韦伯数(We)的增大,液体破裂出现了轴对称...基于大涡模拟(large eddy simulation,LES)和流体体积(volume of fluid,VOF)法对同轴双通道气-液射流进行了3D数值模拟,用LES法求解气液相的流场,VOF法捕捉运动的气液界面.数值模拟结果表明,随着韦伯数(We)的增大,液体破裂出现了轴对称模式和非轴对称模式,以破裂长度、未扰动长度及破裂形态为表征液体射流的特征量与实验值进行了对比,得到了较好的一致性.展开更多
文摘isoAdvector是一种新的几何VOF(Volume of Fluid)方法,虽然克服了传统几何VOF方法难以适用于三维空间任意多面体网格的缺点,但不能直接用来模拟涉及动网格技术的液舱晃荡.为此,引入了运动通量修正,并提出了面-界面交线运动修正,使得修正后的isoAdvector方法可以应用到液舱晃荡的模拟中.基于不同的VOF方法对非共振、共振受迫晃荡和单次冲击波面进行数值模拟,并将模拟结果与试验结果以及解析解进行了比较.结果表明:相对于代数VOF方法,采用修正后的isoAdvector方法获得的自由液面位置和整体水动力载荷精度更高;捕捉的波面没有褶皱,能够较好地模拟波面的翻卷和破碎.此外,提出了界面厚度的估计方法,分析了自由液面波高精度提高的原因.
文摘针对自由表面流动与弹性结构的流固耦合计算效率低、计算耗时长的问题,将流体体积法与基于结构-虚拟弹性体的快速动网格方法相结合,发展了一种适用于自由表面流动的高效流固耦合方法。使用流体体积(volume of fluid,VOF)法对流体自由表面进行追踪;将流体域视为虚拟弹性体并构建结构-虚拟弹性体系统,以流固耦合界面的多相流体力为激励求解系统的动力学方程得到结构振动位移和流场网格变形;在每一个时间步内依次求解流体流动、结构变形和流场动网格,实现流固耦合计算。基于发展的方法计算了溃坝水流冲击下弹性挡板的流固耦合响应,得到了溃坝水流的自由液面和弹性挡板的运动行为,结果表明:自由液面演变和弹性挡板振动位移的计算结果与已有算法的结果吻合良好;在同等网格规模下,与已有算法相比本文方法可减少33.3%的计算时间;在水流冲击作用下,弹性挡板向冲击侧小幅弯曲。随后水流沿挡板左侧上升并形成射流,挡板向另一侧大幅弯曲。最后由于两侧流体的阻尼,挡板振幅逐渐衰减。
文摘基于大涡模拟(large eddy simulation,LES)和流体体积(volume of fluid,VOF)法对同轴双通道气-液射流进行了3D数值模拟,用LES法求解气液相的流场,VOF法捕捉运动的气液界面.数值模拟结果表明,随着韦伯数(We)的增大,液体破裂出现了轴对称模式和非轴对称模式,以破裂长度、未扰动长度及破裂形态为表征液体射流的特征量与实验值进行了对比,得到了较好的一致性.
