马兰戈尼对流分析

分类: 分析 | 综合版 2026-04-06
CAE visualization for marangoni convection theory - technical simulation diagram
马兰戈尼对流分析

马兰戈尼对流的理论基础

概要

🧑🎓

老师！ 今天是马兰戈尼对流分析的话题，是吧？ 那是什么样的东西呢？


🎓

由表面张力的温度梯度驱动的对流。焊接熔融池的流动模式。影响晶体生长的质量。




支配方程式




$$ Ma = \frac{|\partial\sigma/\partial T|\Delta T L}{\mu\alpha} $$
$$ \tau = \frac{\partial\sigma}{\partial T}\frac{\partial T}{\partial s} $$




离散化方法

🧑🎓

这个方程在计算机上实际如何求解呢？


🎓

采用有限元法（FEM）进行空间离散化。组装单元刚性矩阵，构建整体刚性方程。


🎓

进行弱形式（变分形式）的转换，使用试函数和形状函数进行Galerkin法定式化。要素类型的选择（低阶要素 vs. 高阶要素、完全积分 vs. 降阶积分）直接关系到解的精度和计算成本的权衡。




矩阵求解算法

🧑🎓

矩阵求解算法，具体指什么呢？


🎓

直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）来求解线性方程组。对于大规模问题，预处理迭代法最为有效。



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

也就是说，在有限元法的地方偷工减料，最后会吃大亏，对吧。铭记于心！


商用工具中的实现

🧑🎓

那么，做马兰戈尼对流分析用什么软件呢？


工具名称开发商/现在主要文件格式
Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
COMSOL MultiphysicsCOMSOL AB.mph
Ansys Mechanical（原ANSYS Structural）Ansys Inc..cdb, .rst, .db, .ans, .mac

供应商系谱和产品整合的历程

🧑🎓

各家软件的发展史，是不是挺戏剧化的？



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


🧑🎓

听到这儿，我终于理解为什么开发历史很重要了！



COMSOL Multiphysics

🧑🎓

请告诉我关于"COMSOL Multiphysics"的内容！


🎓

1986年在瑞典成立。以MATLAB连接的FEMLAB开始，后改名为COMSOL。多场耦合分析见长。

现在的所属: COMSOL AB


🧑🎓

哇～，开发历史的故事，太有意思了！ 还能继续讲吗？


文件格式和互操作性

🧑🎓

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


格式扩展名类型概要
STEP.stp/.step中立CAD符合ISO 10303标准的3D CAD数据交换格式。形状+PMI支持。
IGES.igs/.iges中立CAD早期的CAD数据交换规范。曲面数据的兼容性存在问题。正在向STEP迁移。
CGNS.cgnsCFD数据CFD General Notation System。CFD结果的标准交换格式。
VTK.vtk/.vtu可视化Visualization Toolkit格式。用于ParaView等。
🎓

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


🧑🎓

原来…格式看起来很简单，其实深度可深了。


实务中的注意事项

🧑🎓

有没有教科书里没有的"实地知识"之类的？


🎓

网格收敛性的确认、边界条件的妥当性验证、材料参数的敏感性分析非常重要。


🎓
  • 网格依赖性验证: 至少用3个网格密度水平确认收敛性
  • 边界条件的妥当性: 设置有物理意义的约束条件
  • 结果的验证: 与理论解、实验数据、已知基准问题的比较



    • 🎓

    嗯，做得不错！ 实际动手是最好的学习。有不懂的地方随时问我。


    Coffee Break 闲话

    卡洛·马兰戈尼——1865年的发现与"光学错觉"论文

    系统研究马兰戈尼对流的首位科学家是意大利物理学家卡洛·马兰戈尼，他在博士论文（1865年）中奠定了基础。但最早研究"为什么向玻璃杯倒入葡萄酒时，液体会沿着内壁流下来"这个现象的是英国医生詹姆斯·汤姆逊（1855年），马兰戈尼将其系统化了。这一现象被称为"葡萄酒的眼泪（Tears of Wine）"。有趣的是，马兰戈尼在最初面对批评者"这种对流不是由表面张力引起，而是眼睛的光学错觉"的反驳时，一边辩论一边继续研究，这告诉我们科学发展从来不是一帆风顺的。

    马兰戈尼对流的数值计算方法

    数值方法的详细内容

    🧑🎓

    具体用什么算法来求解马兰戈尼对流分析呢？



    🧑🎓

    老师的说明清楚了！ 对马兰戈尼对流分析的迷茫消散了。


    离散化的表述



    🎓

    用形状函数 $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)| $$

    基础方程的离散形式


    🎓

    用数式表示就是这样。


    $$ Ma = \frac{|\partial\sigma/\partial T|\Delta T L}{\mu\alpha} $$
    $$ \tau = \frac{\partial\sigma}{\partial T}\frac{\partial T}{\partial s} $$

    🧑🎓

    嗯…只有式子的话还是有点看不明白…什么意思呢？


    🎓

    连续体的支配方程离散化后，得到以下代数方程组:



    $$ [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)$的阶数减少误差（光滑解的情况）


        🧑🎓

        原来…网格细化看起来简单，其实深度可深了。


        求解器设置的建议

        🧑🎓

        具体用什么算法来求解马兰戈尼对流分析呢？


        参数推荐值备注
        迭代法收敛判定$10^{-6}$残差范数标准
        预处理方法ILU(0) or AMG根据问题规模
        最大迭代次数1000不收敛时需重新设定
        内存模式In-core尽可能使用

        单体法

        将所有物理场作为一个连立方程组同时求解。对强耦合问题稳定，但实现复杂，内存消耗大。

        分割法（分离迭代法）

        各物理场独立求解，在界面处进行数据交换。实现容易，可利用现有求解器。适合弱耦合问题。

        界面数据转移

        最近邻法（最简单但精度低）、射影法（保守）、RBF插值（对非匹配网格有利）。平衡保守性和精度很关键。

        子迭代

        在各耦合步骤内进行充分反复，确保界面条件的一致性。残差标准基于各物理场的典型值进行缩放。

        Aitken缓和

        自动调整耦合迭代的缓和系数。防止过缓和导致的发散，加快收敛的自适应方法。

        稳定性条件

        注意附加质量效应（流固耦合中流体密度≈结构密度的情况）。不稳定时应用Robin型界面条件或IQN-ILS法。

        马兰戈尼对流的实务应用

        实践指南

        🧑🎓

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


        🎓

        说明马兰戈尼对流分析的实务性分析流程和注意事项。



        分析流程

        🧑🎓

        从第一步开始教我！ 应该从什么开始呢？


        🎓

        1. 前处理 (Pre-processing)

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

        🎓

        2. 求解 (Solving)

        • 求解器设置（求解方法、收敛标准、输出控制）
        • 作业提交和计算执行
        • 收敛性监测

        🎓

        3. 后处理 (Post-processing)

        • 结果可视化（位移、应力、其他物理量）
        • 结果验证和妥当性确认
        • 报告制作


        网格生成的最佳实践

        🧑🎓

        如何判断网格好坏呢？



        单元品质指标

        🧑🎓

        请告诉我关于"单元品质指标"的内容！


        指标理想值允许范围影响
        纵横比1.0< 5.0精度下降
        雅可比比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
            COMSOL MultiphysicsCOMSOL AB.mph
            Ansys Mechanical（原ANSYS Structural）Ansys Inc..cdb, .rst, .db, .ans, .mac

            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个水平以上的网格密度确认了收敛性吗
            • 验证了力的平衡（反力总和）吗
            • 确认了结果在物理上合理的范围内吗
            • 与已知理论解或基准问题进行了比较吗



              • 🎓

              嗯，做得不错！ 实际动手是最好的学习。有不懂的地方随时问我。


              Coffee Break 闲话

              硅单晶拉晶——半导体制造背后的马兰戈尼对流

              作为半导体晶圆原料的硅单晶通过Czochralski法（CZ法）制造。在1400℃以上加热的硅融液中，用种晶缓缓拉上来的工程中，融液内的马兰戈尼对流左右着晶体质量。不均匀的对流会导致氧和碳的不均匀分布（"应变纹"），成为半导体器件特性波动的原因。因此在直径300mm的大型炉中应用磁场（MCZ法）来抑制马兰戈尼对流——外部磁场通过洛伦兹力阻挡对流的原理。说"智能手机芯片的品质由马兰戈尼对流的控制决定"也不为过。

              马兰戈尼对流的软件比较

              商用工具比较

              🧑🎓

              有很多不同的软件吧？ 各自的特点告诉我！


              🎓

              详细说明支持马兰戈尼对流分析的主要商用CAE工具的功能比较和各产品的历史背景。


              🧑🎓

              老师的说明清楚了！ 对马兰戈尼对流分析的迷茫消散了。


              支持工具清单

              🧑🎓

              那么，做马兰戈尼对流分析用什么软件呢？


              工具名称开发商/现在主要文件格式
              Ansys FluentAnsys Inc..cas, .dat, .msh, .jou
              Simcenter STAR-CCM+Siemens Digital Industries Software.sim, .java, .csv
              COMSOL MultiphysicsCOMSOL AB.mph
              Ansys Mechanical（原ANSYS Structural）Ansys Inc..cdb, .rst, .db, .ans, .mac

              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


              🧑🎓

              听到这儿，我终于理解为什么开发历史很重要了！



              COMSOL Multiphysics

              🧑🎓

              请告诉我关于"COMSOL Multiphysics"的内容！


              🎓

              1986年在瑞典成立。以MATLAB连接的FEMLAB开始，后改名为COMSOL。多场耦合分析见长。

              现在的所属: COMSOL AB



              Ansys Mechanical（原ANSYS Structural）

              🧑🎓

              请告诉我关于"Ansys Mechanical"的内容！


              🎓

              1970年由Swanson Analysis Systems Inc. (SASI)开发。以APDL（Ansys Parametric Design Language）为基础。

              现在的所属: Ansys Inc.


              🧑🎓

              啊，原来如此！ 开发历史是这样的机制啊。


              功能比较矩阵

              🧑🎓

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


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

              转换时的风险

              🧑🎓

              转换时的风险，具体是怎么回事？


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

                • 🧑🎓

                啊，原来如此！ 不同工具间的转换是这样的机制啊。


                许可形式

                🧑🎓

                "许可形式"听说过，但可能还没真正理解…


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

                选择指南

                🧑🎓

                到底该选哪个，判断标准是什么？


                🎓

                在马兰戈尼对流分析工具的选择中，应考虑以下因素:


                🎓
                • 分析规模: 可扩展性到数万～数亿DOF
                • 物理模型: 所需本构模型·单元类型的对应状况
                • 工作流: 与CAD的连接、自动化的容易程度
                • 成本: 初期投资 + 年度维护 + 教育成本
                • 支持: 技术支持的质量和响应速度



                  • 🎓

                  嗯，做得不错！ 实际动手是最好的学习。有不懂的地方随时问我。


                  Coffee Break 闲话

                  能进行马兰戈尼分析的工具——选择其实很有限

                  "真正能分析"马兰戈尼对流的工具在通用CFD中其实很有限。ANSYS Fluent可以通过CSF模型为基础，用UDF添加表面张力温度系数来适应，但设置工作量很大。STAR-CCM+拥有"Marangoni stress"选项作为标准配置，在焊接·晶体生长领域的采用实绩众多。OpenFOAM在社区公开了interFoam拓展的马兰戈尼对应版（marangoniDyFoam），但稳定性验证需要自己进行。焊接专用的SIMUFACT Welding拥有整合热源·马兰戈尼·相变化的模型，在汽车零部件焊接设计中广泛使用。

                  马兰戈尼对流的前沿研究

                  前沿话题和研究动向

                  🧑🎓

                  马兰戈尼对流分析领域今后怎样进化呢？


                  🎓

                  一起看看马兰戈尼对流分析中最新的研究动向和先进方法。


                  🧑🎓

                  老师的说明清楚了！ 对马兰戈尼对流分析的迷茫消散了。


                  最新的数值方法

                  🧑🎓

                  接下来是最新数值方法的话题吧。内容是什么样的呢？



                  🧑🎓

                  嗯…只有式子的话还是有点看不明白…什么意思呢？


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


                    • 高性能计算（HPC）的对应


                    并行化方法概要适用求解器
                    MPI（区域分割）分布式内存型。大规模问题的标准全主要求解器
                    OpenMP共享内存型。节点内并行多数求解器
                    GPU（CUDA