壁面模型LES（WMLES）

分类: 流体分析（CFD） | 整合版 2026-04-06
CAE visualization for wmles theory - technical simulation diagram
壁面模型LES（WMLES）

壁面模型LES（WMLES）的理论基础

概要

🧑🎓

老师！今天就是壁面模型LES（WMLES）的内容，对吧？具体是什么呢？


🎓

通过对壁面近傍进行建模来削减LES的计算成本。适用于高Reynolds数流动。




支配方程式




$$ \tau_w=f(u_{LES}, y_{match}) $$
$$ y_{match}/\delta \sim 0.1 $$



🧑🎓

原来如此……壁面模型看似简单，但实际上内涵很深啊。


离散化方法

🧑🎓

这些方程具体怎样用计算机求解呢？


🎓

我们使用有限单元法（FEM）进行空间离散化。组装单元刚度矩阵，构建全局刚度方程。


🎓

我们将微分方程转换为弱形式（变分形式），使用Galerkin法进行定式化，采用试函数和形状函数。选择单元类型（低阶单元 vs. 高次单元、完全积分 vs. 低阶积分）与解精度和计算成本的折衷息息相关。




矩阵求解算法

🧑🎓

矩阵求解算法具体是怎样的呢？


🎓

直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）来求解联立方程。对于大规模问题，前处理迭代法很有效。



求解法分类内存使用量适用规模
LU分解直接法O(n²)小～中规模
Cholesky分解直接法（对称正定）O(n²)小～中规模
PCG法迭代法O(n)大规模
GMRES法迭代法O(n·m)大规模·非对称
AMG前处理前处理O(n)超大规模
🧑🎓

也就是说，在有限单元法这个环节偷工减料的话，后来就要吃苦头。我铭记于心了！


商用工具中的实现

🧑🎓

那做壁面模型LES（WMLES）需要什么样的软件呢？


工具名开发者/现属主要文件格式
Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
OpenFOAM开源软件（OpenCFD/ESI、OpenFOAM Foundation）字典文件（blockMeshDict等）, .foam
Ansys CFXAnsys Inc..cfx, .def, .res, .ccl

供应商的系统和产品整合的经过

🧑🎓

各个软件的发展历程，是不是有很多故事呢？



Ansys Fluent

🧑🎓

下面说说Ansys Fluent的内容吧。是什么呢？


🎓

由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。

现在的所属: Ansys Inc.



Simcenter STAR-CCM+

🧑🎓

下面说说Simcenter STAR的内容吧。是什么呢？


🎓

由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌。多面体网格是其特点。

现在的所属: Siemens Digital Industries Software


🧑🎓

听到这里，我终于明白为什么供应商系统那么重要了！



OpenFOAM

🧑🎓

OpenFOAM具体是怎样的呢？


🎓

源自Imperial College London的开源CFD。OpenCFD Ltd（ESI Group旗下）和The OpenFOAM Foundation并行开发。

现在的所属: 开源软件（OpenCFD/ESI、OpenFOAM Foundation）


🧑🎓

哦，供应商系统的话题超有趣啊！请继续讲述。


文件格式和互操作性

🧑🎓

在不同软件间传递数据时有什么要注意的吗？


格式扩展名种类概要
CGNS.cgnsCFD数据CFD General Notation System。CFD结果的标准交换格式。
VTK.vtk/.vtu可视化Visualization Toolkit格式。用于ParaView等软件。
🎓

在不同求解器间转换模型时，要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高阶单元或特殊单元（内聚单元、用户定义单元等）在求解器间通常无法直接转换。


🧑🎓

原来如此……文件格式看似简单，但实际上内涵很深啊。


实务中的注意事项

🧑🎓

教科书里没有的"现场智慧"之类的东西有吗？


🎓

网格收敛性的确认、边界条件的合理性验证、材料参数的灵敏度分析非常重要。


🎓
  • 网格依赖性验证：至少用3个网格密度等级确认收敛性
  • 边界条件的合理性：设置物理意义明确的约束条件
  • 结果验证：与理论解、实验数据、已知基准问题的对比


    • 🧑🎓

    哇，壁面模型LES（WMLES）好深啊……但是听了老师的讲解，总算理清头绪了！


    🎓

    很好！走在正确的路上。实际上手操作是最好的学习方式，有不懂的随时问我啊。


    Coffee Break 趣谈

    WMLES的"$Re^{2.6}$之壁"——为什么壁面分辨LES行不通

    壁面分辨LES（Wall-Resolved LES）的网格点数为$N \sim Re^{13/7} \approx Re^{1.86}$数量级，但考虑到时间步长，实际计算成本约为$Re^{2.6}$数量级。$Re = 10^7$的客运飞机翼，即使用当今最强的超级计算机也需数十年来完成计算。WMLES通过对壁面近傍的细小尺度进行建模，大幅降低这个指数——"不是物理上的妥协，而是工程上的理性"——这就是WMLES存在的理由。

    壁面模型LES（WMLES）的数值计算方法

    数值方法的详细

    🧑🎓

    具体怎样用什么算法来求解壁面模型LES（WMLES）呢？



    🧑🎓

    也就是说，在壁面模型的地方偷工减料的话，后来就要吃苦头。我铭记于心了！


    离散化的定式化



    🎓

    用形状函数$N_i$来近似未知量：



    $$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$




    🎓

    用数式表达如下。


    $$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$

    基础方程式的离散形


    🎓

    用数式表达如下。


    $$ \tau_w=f(u_{LES}, y_{match}) $$
    $$ y_{match}/\delta \sim 0.1 $$

    🧑🎓

    嗯……光看式子的话不太能理解……这表示什么呢？


    🎓

    连续体的支配方程进行离散化后，可得以下代数方程组：



    $$ [K]\{u\} = \{F\} $$


    🎓

    其中$[K]$是全局刚度矩阵（或等价的系统矩阵），$\{u\}$是未知节点变量矢量，$\{F\}$是外力矢量。


    🧑🎓

    啊，这样啊！连续体的支配方程是这样一种机制呀。


    单元技术

    🧑🎓

    听说过"单元技术"这个词，但可能理解得不够准确……


    单元类型阶次节点数(3D)精度计算成本
    四面体1次线性4低（剪切锁定）低
    四面体2次二次10高中
    六面体1次线性8中中
    六面体2次二次20非常高高
    棱柱线性/二次6/15中～高中

    积分方案

    🧑🎓

    积分方案具体是怎样的呢？


    🎓
    • 完全积分：对所有项进行精确积分。有刚度高估的倾向（锁定）
    • 低阶积分：削减积分点数。提高计算效率，但有沙漏模式出现的风险
    • 选择性低阶积分 (B-bar法)：将体积项和偏差项分开积分。避免锁定

      • 🧑🎓

      听到这里，我终于明白为什么单元类型那么重要了！


      收敛性和稳定性

      🧑🎓

      不收敛的话首先要检查什么？


      🎓
      • h-细分：细分网格（减小单元尺寸h）来提高精度
      • p-细分：提高单元的多项式阶次来提高精度
      • hp-细分：同时优化h和p

        • 🎓

        收敛速度：二次单元的误差按$O(h^2)$数量级减小（光滑解的情况）


        🧑🎓

        原来如此……细分网格看似简单，但实际上内涵很深啊。


        求解器设置的推荐事项

        🧑🎓

        具体怎样用什么算法来求解壁面模型LES（WMLES）呢？


        参数推荐值说明
        迭代法的收敛判定$10^{-6}$残差范数标准
        前处理方法ILU(0) or AMG取决于问题规模
        最大迭代次数1000不收敛时需重新调整设置
        内存模式In-core尽可能

        迎风差分（Upwind）

        一次迎风：数值扩散大但稳定。二次迎风：精度提高但有振荡风险。高Reynolds数流动中是必要的。

        中心差分（Central Differencing）

        二次精度，但当Pe数 > 2时会出现数值振荡。低Reynolds数扩散主导流动中适用。

        TVD方案（MUSCL、QUICK等）

        通过限流器函数抑制数值振荡，同时维持高精度。对于冲击波和陡峭梯度的捕捉有效。

        有限体积法 vs 有限单元法

        FVM：自然满足守恒律。是CFD的主流。FEM：在复杂形状和多物理耦合中有优势。SPH等无网格法也在发展中。

        CFL条件（Courant数）

        显式法：CFL ≤ 1是稳定条件。隐式法：CFL > 1也稳定，但影响精度和迭代次数。LES：推荐CFL ≈ 1。物理意义：1个时间步长内信息不会传播超过1个单元。

        残差监测

        连续方程·动量·能量各残差下降3～4个数量级即可判为收敛。质量守恒的残差尤其重要。

        松弛系数

        压力：0.2～0.3、速度：0.5～0.7是一般的初值。发散时降低松弛系数。收敛后可提高来加速。

        非定常计算的内部迭代

        每个时间步内进行迭代至定常解。内部迭代数：5～20次是目安。若残差在时间步间波动，需重新检查时间步长。

        壁面模型LES（WMLES）的实务适用

        实践指南

        🧑🎓

        老师，关于"实践指南"，请教我！


        🎓

        介绍壁面模型LES（WMLES）的实务分析流程和注意事项。


        🧑🎓

        也就是说，在壁面模型的地方偷工减料的话，后来就要吃苦头。我铭记于心了！


        分析流程

        🧑🎓

        从第一步开始教我！从哪里开始好呢？


        🎓

        1. 前处理 (Pre-processing)

        • 导入CAD数据并简化形状
        • 定义材料特性
        • 网格生成（单元类型·尺寸的确定）
        • 设置边界条件和荷载条件

        🎓

        2. 求解 (Solving)

        • 求解器设置（求解法、收敛标准、输出控制）
        • 投入作业并执行计算
        • 收敛监测

        🎓

        3. 后处理 (Post-processing)

        • 结果可视化（位移、应力、其他物理量）
        • 结果验证与合理性确认
        • 报告编制


        网格生成的最佳实践

        🧑🎓

        网格的好坏怎样判断呢？



        单元品质指标

        🧑🎓

        关于"单元品质指标"，请教我！


        指标理想值允许范围影响
        纵横比1.0< 5.0精度下降
        Jacobian比1.0> 0.3单元退化
        翘曲0°< 15°精度下降
        斜率0°< 45°收敛性恶化
        锥形比0< 0.5精度下降

        网格密度的确定

        🧑🎓

        网格密度的确定具体是怎样的呢？


        🎓
        • 应力集中部：至少配置3层以上的单元
        • 应力梯度大的区域：单元尺寸减小至周围的1/3～1/5
        • 荷载施加点近傍：局部细分
        • 远方区域：粗网格以保证计算效率


          • 边界条件的设置指南

          🧑🎓

          听说边界条件这里搞错了，后面全都完蛋……


          🎓
          • 注意过约束：刚体运动的约束仅6个自由度
          • 对称条件的活用：削减计算规模
          • 荷载的等价分配：集中荷载 vs. 分布荷载的选择

            • 🧑🎓

            啊，这样啊！注意过约束就是这样一种机制呀。


            按商用工具的实现步骤

            🧑🎓

            各种各样的软件呀。分别有什么特点？


            工具名开发者/现属主要文件格式
            Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
            Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
            OpenFOAM开源软件（OpenCFD/ESI、OpenFOAM Foundation）字典文件（blockMeshDict等）, .foam
            Ansys CFXAnsys Inc..cfx, .def, .res, .ccl

            Ansys Fluent

            🧑🎓

            下面说说Ansys Fluent的内容吧。是什么呢？


            🎓

            由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。

            现在的所属: Ansys Inc.



            Simcenter STAR-CCM+

            🧑🎓

            下面说说Simcenter STAR的内容吧。是什么呢？


            🎓

            由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌。多面体网格是其特点。

            现在的所属: Siemens Digital Industries Software


            🧑🎓

            老师的讲解很清楚！关于工具名的疑问总算散开了。


            常见失败和对策

            🧑🎓

            初学者容易犯什么错误？想事先知道！


            症状原因对策
            计算不收敛网格品质不良、不恰当的边界条件改善网格、重新检查拘束条件
            应力异常大应力特异点、网格依赖回避特异点、局部网格细分
            位移非现实材料常数错误、单位制不统一确认输入数据
            计算时间过长不必要的细分、低效的求解法网格优化、并行计算

            质量保证检查清单

            🧑🎓

            教科书里没有的"现场智慧"之类的东西有吗？


            🎓
            • 用3个网格密度等级以上确认了网格收敛性吗
            • 验证了力的平衡（反力合计）吗
            • 确认了结果在物理合理范围内吗
            • 与已知理论解或基准问题做过比较吗


              • 🧑🎓

              哇，壁面模型LES（WMLES）好深啊……但是听了老师的讲解，总算理清头绪了！


              🎓

              很好！走在正确的路上。实际上手操作是最好的学习方式，有不懂的随时问我啊。


              Coffee Break 趣谈

              风力发电涡轮的WMLES——向实机规模的挑战

              在风力发电涡轮叶片（翼弦长数米，$Re \sim 10^7$）的CFD中，WMLES正受到关注。传统的RANS在评估叶片剥离·失速时存在局限，年平均发电量预测误差可超5%。使用WMLES后，剥离始发角和再附着的预测精度改善，翼型优化的可信度上升。计算成本是RANS的100～1000倍，但在数十亿元级的风力发电项目中，仅将设计失误减少1%就足以抵消成本——这就是现实中WMLES采用不断推进的原因。

              壁面模型LES（WMLES）的软件比较

              商用工具比较

              🧑🎓

              各种各样的软件呀。分别有什么特点？


              🎓

              介绍支持壁面模型LES（WMLES）的主要商用CAE工具的功能比较及各产品的历史背景。


              🧑🎓

              也就是说，在壁面模型的地方偷工减料的话，后来就要吃苦头。我铭记于心了！


              支持工具列表

              🧑🎓

              那做壁面模型LES（WMLES）需要什么样的软件呢？


              工具名开发者/现属主要文件格式
              Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
              Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
              OpenFOAM开源软件（OpenCFD/ESI、OpenFOAM Foundation）字典文件（blockMeshDict等）, .foam
              Ansys CFXAnsys Inc..cfx, .def, .res, .ccl

              Ansys Fluent

              🧑🎓

              下面说说Ansys Fluent的内容吧。是什么呢？


              🎓

              由Fluent Inc.开发。2006年被Ansys收购。基于非结构网格的通用CFD求解器。

              现在的所属: Ansys Inc.



              Simcenter STAR-CCM+

              🧑🎓

              下面说说Simcenter STAR的内容吧。是什么呢？


              🎓

              由CD-adapco开发。2016年被Siemens收购，整合到Simcenter品牌。多面体网格是其特点。

              现在的所属: Siemens Digital Industries Software


              🧑🎓

              听到这里，我终于明白为什么供应商系统那么重要了！



              OpenFOAM

              🧑🎓

              OpenFOAM具体是怎样的呢？


              🎓

              源自Imperial College London的开源CFD。OpenCFD Ltd（ESI Group旗下）和The OpenFOAM Foundation并行开发。

              现在的所属: 开源软件（OpenCFD/ESI、OpenFOAM Foundation）



              Ansys CFX

              🧑🎓

              关于"Ansys CFX"，请教我！


              🎓

              由AEA Technology (UK) 开发的CFX。2003年被Ansys收购。结合型求解器是其特点。

              现在的所属: Ansys Inc.


              🧑🎓

              啊，这样啊！供应商系统就是这样一种机制呀。


              功能比较矩阵

              🧑🎓

              预算和时间都有限，性价比最强的是哪个？


              功能FluentStar-CCM+OpenFOAMCFX
              基本功能○○○○
              高级功能○○○△
              自动化/脚本○○○○
              并行计算○○○○
              GPU支持△△△○

              转换时的风险

              🧑🎓

              转换时的风险具体是怎样的呢？


              🎓
              • 单元类型的不兼容：求解器固有单元无法用中立格式表达
              • 材料模型的差异：同名但内部实现有所不同的情况存在
              • 边界条件的重新定义：多数情况下需要手动重新设置
              • 结果数据的比较：输出变量的定义（节点值 vs. 单元值、积分点值）存在差异

                • 🧑🎓

                啊，这样啊！不同工具间的模型转换就是这样一种机制呀。


                许可证形式

                🧑🎓

                听说过"许可证形式"这个词，但可能理解得不够准确……


                工具许可证特点
                商用FEA节点锁定/浮动价格高但有官方支持
                OpenFOAMGPL无偿但支持需付费
                COMSOL节点锁定/浮动按模块购买
                Code_AsterGPLEDF开发的开源求解器

                选择指南

                🧑🎓

                总之该选哪个，能教我判断标准吗？


                🎓

                在壁面模型LES（WMLES）工具选择中应考虑以下方面：


                🎓
                • 分析规模：数万～数亿自由度的可扩展性
                • 物理模型：必要的本构关系·单元类型的支持情况
                • 工作流：与CAD的整合、自动化的易用性
                • 成本：初期投资 + 年度维护 + 教育成本
                • 支持：技术支持的质量与响应速度


                  • 🧑🎓

                  哇，壁面模型LES（WMLES）好深啊……但是听了老师的讲解，总算理清头绪了！


                  🎓

                  很好！走在正确的路上。实际上手操作是最好的学习方式，有不懂的随时问我啊。


                  Coffee Break 趣谈

                  壁面模型LES工具——Ansys Fluent WMLES vs OpenFOAM dynamicKEqn

                  壁面模型LES（WMLES）的商用工具主要有Ansys Fluent（WMLES模式、σ-模型＋美国壁面模型）和Charles（Cadence）、OpenFOAM（dynamicKEqn SGS模型）等。Fluent自2019年以来强化了WMLES专用求解器，y+<1不必要的壁面模型设置和HPC并行效率得到改善。OpenFOAM的dynamicKEqn作为动态Smagorinsky模型实现了dynamic Lagrangian法，具有灵活的定制性。Charles是专门针对X-LES（壁面模型＋LES）的商用工具，在航空航天领域评价很高。HPC计算成本（WMLES: RANS×100～1000）是商业应用的主要障壁。

                  壁面模型LES（WMLES）的前沿研究

                  前沿话题和研究趋势

                  🧑🎓

                  壁面模型LES（WMLES）领域今后会如何发展呢？


                  🎓

                  介绍壁面模型LES（WMLES）中的最新研究动向和先进方法。


                  🧑🎓

                  也就是说，在壁面模型的地方偷工减料的话，后来就要吃苦头。我铭记于心了！


                  最新的数值方法

                  🧑🎓

                  下面说说最新的数值方法吧。是什么呢？



                  🧑🎓

                  光看式子的话不太能理解……这表示什么呢？


                  🎓
                  • 等参数分析 (IGA)：直接使用NURBS基函数，实现CAD-CAE间的无缝协作
                  • 粒子法 (SPH, MPM)：无网格方法追踪大变形·破坏
                  • 相场法 (Phase-Field)：界面的隐式表示实现复杂界面追踪
                  • 机器学习支持：替代模型、物理信息神经网络 (PINN)


                    • 高性能计算 (HPC) 的应对


                    并行化方法概要支持求解器