DOI: 10.3390/fluids10010004 | Fluids 2025, 10, 41 仿真背景与价值圆柱绕流是流体力学中最为经典又最具挑战的基准问题之一。在亚临界雷诺数Re130,000条件下流场同时包含层流分离、剪切层Kelvin-Helmholtz不稳定性、转捩与湍流涡街等多个复杂机制对数值方法的精度和分辨率提出了极高要求。本文基于发表于《Fluids 2025, 10, 4》的研究同时采用商用CFD平台Ansys Fluent OpenFOAM以有限体积法FVM大涡模拟LES对亚临界圆柱绕流进行系统仿真。研究涵盖了非结构六面体网格HM2430万单元与曲线正交O型网格OM1360万单元两种拓扑以及k方程亚格子模型TKE与其动态版本dTKE并在单精度SP、双精度DP和混合精度SPDP三种算术模式下进行了全面对比。特别引入了Lyapunov稳定性理论分析SP与DP求解的动态系统稳定性差异这在CFD工程实践中尚属前沿探索对计算资源优化和GPU加速具有重要工程意义。图1 计算网格非结构六面体HM左与曲线正交O型OM右2 仿真建模方法与关键技术计算网格图1采用两种网格拓扑——非结构六面体HM30D×20D×10D圆柱周向分辨率Nθ184、展向Nz640、首层y≤4×10⁻⁴和曲线正交O型OM444×444×64、首层y1.25×10⁻⁴均为中粗网格约10-25M单元。数学模型基于Favre滤波可压缩Navier-Stokes方程SGS亚格子封闭采用k方程涡粘模型TKE及其动态版本dTKE时间推进采用隐式Euler法BDF-2空间离散采用有界中心差分CDS-2二阶迎风SOU。求解器方面Ansys Fluent使用代数多重网格AMGGauss-Seidel迭代OpenFOAM采用对称Gauss-Seidel几何多重网格GAMG。总运行矩阵含7个算例覆盖SP/DP/SPDP三种精度和TKE/dTKE两种SGS模型。入流Ma0.05、Re130,000出口波传输边界展向周期性壁面无滑移等温。统计收敛标准为τconv4时间平均取Nvs25个涡脱落周期。图2 升力系数与阻力系数时间历程OF-dTKE vs AF64-dTKE图3 Q准则涡量可视化Q10⁴OF-dTKE vs AF64-dTKE图4 圆柱表面压力系数a与摩擦系数b分布3 仿真结果对比分析积分参数对比表2图2、4、5各算例的阻力系数⟨Cd⟩在0.94–1.34范围与实验值1.2Cantwell Coles, Achenbach, Schewe相比OF-TKE/dTKE略高≈1.34AF-dTKE略低≈0.94但整体LES均值1.19±20%与实验一致。升力脉动⟨C′l⟩0.27–0.53Strouhal数St0.18–0.21与实验St0.18–0.20吻合良好。OF-TKE/dTKE的升力PDF呈双峰分布AF-dTKE呈高斯正态分布与Schewe实验一致。 压力与摩擦系数图4OF-TKE/dTKE与Cantwell Coles实验及Cao Tamura数值结果在压力系数分布上吻合良好AF-dTKE更吻合Achenbach实验数据。分离角预测θs83–85°略大于实验77–78°但与其他LES结果81–94°一致。 中心线速度及湍流统计图5OF与AF在远尾迹区x/D2均满意再现实验数据近尾迹区AF的再循环区长度Lr0.68D略长于OF的0.55–0.58D和实验0.4–0.5D。垂向速度脉动方面AF-dTKE略低15%OF-dTKE在近尾迹区偏差约30%。 能谱分析图6一维能量谱在全波数范围清晰呈现−5/3惯性子区幂律与Ong Wallace和Parnaudeau实验数据良好吻合f/fvs≈1之后各曲线开始缓慢衰减说明数值耗散可控。 单双精度对比图7Lyapunov分析SP与DP的积分参数和局部流场特性无显著差异Lyapunov散度曲线呈现典型线性增长区→饱和区的双区特征。SP动态系统Lyapunov稳定预测时间tp≈1.5τconv饱和后两系统达到各自统计稳态。这意味着在ANSYS Fluent框架下单精度LES完全可用于工程级圆柱绕流仿真在保持物理精度的前提下节省约10-20%的CPU内存和计算时间为GPU加速通常优化于SP提供了数值依据。图5 中心线平均流向速度a及流向/垂向RMS脉动分布b,c图6 近尾迹垂向速度一维能量谱含−5/3幂律及剪切层KH不稳定频率图7 Lyapunov散度曲线SPAF32与DPAF64速度场对比4 总结与服务本文通过Ansys Fluent OpenFOAM双平台的LES仿真系统评估了Re130,000圆柱绕流中网格拓扑、SGS模型及浮点精度对结果的影响。研究表明1OF和AF在积分参数和流动物理上基本一致差异主要源于数值扩散水平2单精度算术在Fluent中可实现与双精度同等的LES精度Lyapunov稳定3混合精度SPDP是兼顾精度与效率的折中方案。我方团队具备Fluent与OpenFOAM双平台LES/DES/WMLES的完整仿真能力可为钝体绕流、气动优化、风工程等复杂湍流问题提供高精度CFD解决方案。 LES圆柱绕流Fluent-OpenFOAM单双精度对比
DOI: 10.3390/fluids10010004 | Fluids 2025, 10, 41 仿真背景与价值圆柱绕流是流体力学中最为经典又最具挑战的基准问题之一。在亚临界雷诺数Re130,000条件下流场同时包含层流分离、剪切层Kelvin-Helmholtz不稳定性、转捩与湍流涡街等多个复杂机制对数值方法的精度和分辨率提出了极高要求。本文基于发表于《Fluids 2025, 10, 4》的研究同时采用商用CFD平台Ansys Fluent OpenFOAM以有限体积法FVM大涡模拟LES对亚临界圆柱绕流进行系统仿真。研究涵盖了非结构六面体网格HM2430万单元与曲线正交O型网格OM1360万单元两种拓扑以及k方程亚格子模型TKE与其动态版本dTKE并在单精度SP、双精度DP和混合精度SPDP三种算术模式下进行了全面对比。特别引入了Lyapunov稳定性理论分析SP与DP求解的动态系统稳定性差异这在CFD工程实践中尚属前沿探索对计算资源优化和GPU加速具有重要工程意义。图1 计算网格非结构六面体HM左与曲线正交O型OM右2 仿真建模方法与关键技术计算网格图1采用两种网格拓扑——非结构六面体HM30D×20D×10D圆柱周向分辨率Nθ184、展向Nz640、首层y≤4×10⁻⁴和曲线正交O型OM444×444×64、首层y1.25×10⁻⁴均为中粗网格约10-25M单元。数学模型基于Favre滤波可压缩Navier-Stokes方程SGS亚格子封闭采用k方程涡粘模型TKE及其动态版本dTKE时间推进采用隐式Euler法BDF-2空间离散采用有界中心差分CDS-2二阶迎风SOU。求解器方面Ansys Fluent使用代数多重网格AMGGauss-Seidel迭代OpenFOAM采用对称Gauss-Seidel几何多重网格GAMG。总运行矩阵含7个算例覆盖SP/DP/SPDP三种精度和TKE/dTKE两种SGS模型。入流Ma0.05、Re130,000出口波传输边界展向周期性壁面无滑移等温。统计收敛标准为τconv4时间平均取Nvs25个涡脱落周期。图2 升力系数与阻力系数时间历程OF-dTKE vs AF64-dTKE图3 Q准则涡量可视化Q10⁴OF-dTKE vs AF64-dTKE图4 圆柱表面压力系数a与摩擦系数b分布3 仿真结果对比分析积分参数对比表2图2、4、5各算例的阻力系数⟨Cd⟩在0.94–1.34范围与实验值1.2Cantwell Coles, Achenbach, Schewe相比OF-TKE/dTKE略高≈1.34AF-dTKE略低≈0.94但整体LES均值1.19±20%与实验一致。升力脉动⟨C′l⟩0.27–0.53Strouhal数St0.18–0.21与实验St0.18–0.20吻合良好。OF-TKE/dTKE的升力PDF呈双峰分布AF-dTKE呈高斯正态分布与Schewe实验一致。 压力与摩擦系数图4OF-TKE/dTKE与Cantwell Coles实验及Cao Tamura数值结果在压力系数分布上吻合良好AF-dTKE更吻合Achenbach实验数据。分离角预测θs83–85°略大于实验77–78°但与其他LES结果81–94°一致。 中心线速度及湍流统计图5OF与AF在远尾迹区x/D2均满意再现实验数据近尾迹区AF的再循环区长度Lr0.68D略长于OF的0.55–0.58D和实验0.4–0.5D。垂向速度脉动方面AF-dTKE略低15%OF-dTKE在近尾迹区偏差约30%。 能谱分析图6一维能量谱在全波数范围清晰呈现−5/3惯性子区幂律与Ong Wallace和Parnaudeau实验数据良好吻合f/fvs≈1之后各曲线开始缓慢衰减说明数值耗散可控。 单双精度对比图7Lyapunov分析SP与DP的积分参数和局部流场特性无显著差异Lyapunov散度曲线呈现典型线性增长区→饱和区的双区特征。SP动态系统Lyapunov稳定预测时间tp≈1.5τconv饱和后两系统达到各自统计稳态。这意味着在ANSYS Fluent框架下单精度LES完全可用于工程级圆柱绕流仿真在保持物理精度的前提下节省约10-20%的CPU内存和计算时间为GPU加速通常优化于SP提供了数值依据。图5 中心线平均流向速度a及流向/垂向RMS脉动分布b,c图6 近尾迹垂向速度一维能量谱含−5/3幂律及剪切层KH不稳定频率图7 Lyapunov散度曲线SPAF32与DPAF64速度场对比4 总结与服务本文通过Ansys Fluent OpenFOAM双平台的LES仿真系统评估了Re130,000圆柱绕流中网格拓扑、SGS模型及浮点精度对结果的影响。研究表明1OF和AF在积分参数和流动物理上基本一致差异主要源于数值扩散水平2单精度算术在Fluent中可实现与双精度同等的LES精度Lyapunov稳定3混合精度SPDP是兼顾精度与效率的折中方案。我方团队具备Fluent与OpenFOAM双平台LES/DES/WMLES的完整仿真能力可为钝体绕流、气动优化、风工程等复杂湍流问题提供高精度CFD解决方案。