📚 课程目录第1章 COMSOL Multiphysics 入门

1. COMSOL Multiphysics 入门

COMSOL Multiphysics 是一款以偏微分方程（PDE）为底层求解核心的多物理场仿真平台。与许多“一个模块对应一类物理场”的仿真软件不同，COMSOL 的组织逻辑是：几何与网格是公共基础，物理场接口负责把物理定律翻译成方程，多物理场耦合节点负责把不同方程之间的相互作用显式连接起来，求解器再统一离散求解。理解这条主线，是入门阶段最重要的一件事。本章围绕界面布局、模型开发器、内置物理场接口与耦合概念展开，帮助读者建立“从模型树到物理方程”的整体认知。

一、界面布局：桌面端与模型开发器

启动 COMSOL 后，默认进入“新建”窗口，可选择模型向导（Model Wizard）或空模型。进入建模界面后，主窗口通常由以下几个区域构成：

模型开发器（Model Builder）：位于左侧，是模型的核心导航树，按“全局定义 → 组件 → 几何 → 材料 → 物理场 → 网格 → 研究 → 结果”的顺序组织。
设置窗口（Settings）：位于中右部，显示当前选中节点的具体参数，是输入物理量、边界条件、材料属性的主要位置。
图形窗口（Graphics）：显示几何、网格、物理场边界选择与结果云图，支持旋转、缩放、剖切与选择。
信息与进度窗口：显示求解日志、警告、错误与计算进度。
工具栏与快捷按钮：提供网格剖分、求解、绘图等常用操作。

模型开发器是 COMSOL 的“骨架”。它并非简单的文件夹，而是一棵有依赖顺序的树：几何必须先于网格，网格必须先于物理场离散，物理场必须先于研究步骤。若顺序颠倒，节点会显示为灰色或报错。因此，初学者应养成“自上而下、逐级展开”的操作习惯。

二、内置物理场接口：从物理定律到方程

COMSOL 内置了大量物理场接口，常见类别包括：

结构力学：固体力学、壳、梁、膜、压电、流固耦合等。
流体流动：层流、湍流、多相流、多孔介质流、管流等。
传热：固体传热、流体传热、辐射、相变传热等。
电磁：静电、静磁、电磁波、射频、等离子体等。
化学与传质：稀物质传递、浓物质传递、化学反应工程等。
声学：压力声学、热黏性声学、声固耦合等。

每个物理场接口本质上是一组预置的 PDE 模板。例如“固体传热”接口默认求解能量守恒方程，用户只需指定材料导热系数、密度、比热容，以及边界上的温度、热通量或对流换热系数。接口中的“域”“边界”“边”“点”节点，对应不同维度的方程项与边界条件。理解这一点，就能明白为什么 COMSOL 的物理场设置总是围绕“选择几何实体 + 指定物理量”展开。

三、多物理场耦合：显式连接不同物理场

多物理场耦合是 COMSOL 的核心优势。耦合并非自动发生，而是通过“多物理场”节点显式建立。常见耦合方式包括：

单向耦合：一个物理场的输出作为另一个物理场的输入，例如热膨胀中温度场影响结构变形，但变形不显著影响温度。
双向耦合：两个物理场相互影响，例如流固耦合中流体压力使固体变形，固体变形又改变流场。
共享变量耦合：不同物理场共享同一温度、压力或电势变量，例如焦耳热中电流产生热量，温度又改变电导率。

在模型开发器中，耦合通常表现为“多物理场”分支下的节点，如“热膨胀”“流固耦合”“焦耳热”“热电效应”等。每个耦合节点会要求用户指定参与耦合的物理场接口和几何实体。若缺少耦合节点，即使两个物理场同时存在，也可能各自独立求解，得不到真正的耦合结果。

四、一个最小示例：热膨胀耦合

为了把上述概念落到操作层面，下面给出一个可在 COMSOL 中复现的最小示例。对象是一根长度 0.2 m、截面 0.01 m × 0.01 m 的铝制杆件，左端固定，初始温度 20 ℃，右端施加 100 ℃ 温度边界。目标是观察热膨胀引起的轴向位移。

模型向导选择：
  空间维度：三维
  物理场：固体力学 + 固体传热
  研究：稳态

全局定义 → 参数：
  L  = 0.2   [m]     杆长
  T0 = 20    [degC]  初始温度
  T1 = 100   [degC]  右端温度

几何：长方体，尺寸 0.2 × 0.01 × 0.01 m

材料：铝
  导热系数 k   = 237   [W/(m·K)]
  密度 rho     = 2700  [kg/m^3]
  比热容 Cp    = 900   [J/(kg·K)]
  杨氏模量 E   = 70e9  [Pa]
  泊松比 nu    = 0.33
  热膨胀系数 α = 23e-6 [1/K]

固体传热：
  温度边界：左端 T = T0，右端 T = T1

固体力学：
  固定约束：左端面
  其余面自由

多物理场 → 热膨胀：
  温度场选择：固体传热（T）
  结构场选择：固体力学（u）

网格：自由四面体，单元尺寸“常规”

研究：稳态求解

求解后，在结果中绘制“固体力学 → 位移场 → 总位移”云图，可观察到右端位移最大。若以线弹性、小变形假设估算，轴向热应变约为 α·ΔT = 23e-6 × 80 = 1.84e-3，对应自由伸长量约 0.2 × 1.84e-3 = 3.68e-4 m，即约 0.368 mm。由于左端固定，右端位移应接近该量级。这里给出的是基于材料参数的估算值，实际仿真结果会因网格、边界约束和材料模型略有差异，应以软件计算结果为准。

判断耦合是否生效的实用方法：在“多物理场”节点下查看是否已添加“热膨胀”；在“固体力学”设置中检查“热膨胀”子节点是否引用了正确的温度变量；求解后对比“仅传热”与“传热+热膨胀”两种研究的位移结果，若后者位移明显非零，说明耦合已参与求解。

五、入门阶段最容易混淆的三个问题

第一，物理场接口与多物理场耦合节点的区别。物理场接口定义“有哪些方程”，耦合节点定义“方程之间如何交换变量”。只添加两个物理场而不添加耦合节点，通常得到的是两个独立问题，而非耦合问题。

第二，材料属性与物理场边界条件的分工。材料属性描述“介质本身的性质”，如导热系数、弹性模量；边界条件描述“外部环境或约束”，如固定约束、温度、热通量。二者不能互相替代。

第三，网格与求解精度的关系。网格越细，结果通常越接近收敛解，但计算量也越大。入门阶段建议先用“常规”网格跑通流程，再针对应力集中、边界层等区域做局部细化，而不是一开始就全局加密。

常见报错提示：若模型开发器中“物理场”节点显示为灰色，通常是因为几何尚未创建或组件维度不一致；若求解时报“未定义材料属性”，应检查材料节点是否已覆盖所有域；若耦合结果为零，优先检查耦合节点是否选择了正确的物理场接口和几何实体。

掌握以上内容后，读者应能独立完成“新建模型 → 选择物理场 → 建立几何 → 指定材料 → 设置边界条件 → 添加耦合 → 剖分网格 → 求解 → 查看结果”的完整链路。这条链路是后续所有章节的共同基础，也是判断一个 COMSOL 模型是否结构清晰、可复现、可扩展的关键依据。