粘性散逸
粘性散逸的理论基础
什么是粘性散逸
老师，粘性散逸就是"粘性产生热"吗，能更准确地解释一下这个现象吗？
流体的运动能量被粘性力不可逆地转化为内部能量（热）的现象。在分子水平上，流体层之间的剪切导致分子摩擦，运动能量散逸。
日常生活中的例子包括：
- 轮胎橡胶在行驶过程中受热（粘弹性体变形散逸）
- 航天飞机大气层再入时的高温（冲击波内的粘性散逸）
- 高分子注射成型时浇口部分局部加热
- 水库泄洪时水温略微升高
水库放水时温度也会升高吗！
理论上是的。从高度 $h = 100\,\text{m}$ 坠落的水的位置能量完全转化为热时：
这与Joule在19世纪的实验测量值基本相符。
粘性散逸函数的推导
能量方程中的粘性散逸项写成：
对于不可压缩的牛顿流体，这是：
张量记法为：
重要的是 $\Phi \geq 0$ 始终成立。这是热力学第二定律的要求，粘性散逸总是将运动能量转化为热，这是一个单向过程。
也就是说这是一个"不可逆过程"。
Brinkman数——粘性散逸的重要性指标
判断粘性散逸是否可以忽略的无量纲数是Brinkman数。
其中 $\mu$ 是粘度，$U$ 是代表速度，$k$ 是热导率，$\Delta T$ 是代表温度差。
| Br 的值 | 解释 | 例子 |
|---|---|---|
| $\text{Br} \ll 1$ | 粘性散逸可忽略 | 普通水的管内流 |
| $\text{Br} \sim O(1)$ | 应考虑粘性散逸 | 高分子加工、润滑油膜 |
| $\text{Br} \gg 1$ | 粘性散逸占主导 | 超高速流、航天器再入 |
只有当Br数很大时才需要在能量方程中加入散逸项吧。
对。不必要地加入散逸项会增加计算成本，也会恶化收敛性，所以用Br数进行事先判断很重要。
粘性散逸的发现史——从Joule加热到流体摩擦热的演进
将流体的粘性散逸（流动能量转化为热的现象）用流体力学方程表示的是乔治·加布里埃尔·斯托克斯（G.G. Stokes）。在1845年的论文"On the Theories of the Internal Friction of Fluids in Motion"中，他数学上描述了粘性应力张量与能量散逸的关系，明确了N-S方程粘性项的物理意义。在此之前，Joule在1843年展示了固体电阻中的电热转换（焦耳加热），而流体的摩擦热可以看作是其流体力学版本。高速飞行体（弹道导弹、航天器再入）的空气动力加热就源于这种粘性散逸，150年前斯托克斯的理论直接应用于阿波罗太空舱的热防护罩设计。
粘性散逸的数值计算方法
粘性散逸在能量方程中的应用
在CFD中计算粘性散逸时，要修改哪个方程？
将粘性散逸函数 $\Phi$ 作为源项加入能量方程。
在有限体积法的离散化中，$\Phi$ 作为体积积分来评估：
在每个单元的中心计算速度梯度，计算剪切速率，然后用 $\Phi = \mu |\dot{\gamma}|^2$ 来评估能量方程中的散逸源项。
动量方程不需要改变吗？
基本上是的。但是，如果粘性散逸导致的温度升高影响粘度（温度相关粘度），动量方程就会通过粘度反馈。在这种情况下，需要耦合求解能量方程和动量方程。
湍流场中的粘性散逸
在湍流场中，粘性散逸的处理与层流不同。在RANS（Reynolds平均）框架中：
$\overline{\Phi}_{\text{mean}}$ 是平均速度场引起的散逸，$\varepsilon$ 是湍流能量的散逸率（即 $k$-$\varepsilon$ 模型中的 $\varepsilon$）。
实际上，在大多数湍流问题中，$\varepsilon$ 作为能量方程的源项被自动考虑。湍流运动能量方程的散逸项正是这个：
$P_k$ 是湍流能量的生成项，$\varepsilon$ 是散逸率。最终 $\varepsilon$ 转化为热。
原来 $k$-$\varepsilon$ 模型中的 $\varepsilon$ 就是粘性散逸！
完全正确。湍流全部运动能量的级联的最后阶段就是粘性散逸（$\varepsilon$）。根据Kolmogorov的标度关系，散逸发生的最小尺度是：
数值精度注意事项
在粘性散逸数值计算中要注意的事项总结如下。
| 项目 | 注意点 | 对策 |
|---|---|---|
| 速度梯度精度 | $\Phi$ 与速度梯度的平方成正比，梯度精度直接影响 | 使用二阶或更高阶方案 |
| 网格分辨率 | 壁面附近的剪切最大 | 充分细化壁面网格 |
| 与数值散逸的区别 | 迎风差分的数值粘性也产生非物理散逸 | 用高阶精度方案降低数值散逸 |
| 非牛顿流体 | $\Phi = \eta(\dot{\gamma})|\dot{\gamma}|^2$ 取决于模型 | 确认实现的散逸模型与粘度模型一致 |
不能混淆数值散逸和物理粘性散逸。
特别是在LES（大涡模拟）中，SGS（小网格尺度）模型引入了额外的散逸。物理散逸与模型散逸的分离对精度至关重要。
粘性散逸的数值处理——高速流体的能量守恒精度保证
在高马赫数流和高粘性流体的分析中，准确求解能量方程中的粘性散逸项Φ=τ:∇u对保证能量守恒很关键。这一项与速度的平方成正比，在流速较高的区域迅速增大。数值上的问题是：在网格较粗的区域，速度梯度被低估，散逸也会被低估，导致能量守恒崩溃。实务上应先用"Brinkman数Br=μU²/(kΔT)"预评价散逸的重要性，当Br>0.1时务必启用散逸。在包含散逸项的计算中，温度残差和速度残差相互影响，因此需要设置更严格的收敛标准（残差≤1e-6）。
粘性散逸的实务应用
判断是否需要考虑粘性散逸
在实务中，如何判断是否应启用粘性散逸？
先进行Brinkman数的初步估算。除此之外，以下情况应考虑粘性散逸。
| 应用领域 | 典型Br | 散逸的影响 |
|---|---|---|
| 高分子注射成型 | $1 \sim 100$ | 浇口温度上升50-100K以上 |
| 轴承·润滑膜 | $0.1 \sim 10$ | 油膜温度控制性能 |
| 挤出成型（模具内） | $1 \sim 50$ | 螺杆旋转产生热 |
| 微通道流 | $0.01 \sim 1$ | 通道径小时Br增大 |
| 超声速/极超声速流 | $\gg 1$ | 冲击波内绝热温度上升 |
| 一般水管道 | $\ll 0.01$ | 可忽略 |
高分子加工中基本是必需的。
正是。在注射成型中，通过浇口时剪切速率达到 $10^4 \sim 10^5\,\text{s}^{-1}$，粘度为 $10^2 \sim 10^3\,\text{Pa}\cdot\text{s}$，所以粘性散逸的温度上升可达几十度。忽视这个温度上升会导致充填图案预测完全不对。
Ansys Fluent
在Fluent中启用粘性散逸的步骤：
1. Models > Energy 打开
2. Energy Dialog > Options 中勾选"Viscous Dissipation"
3. 这样就会自动在能量方程中加入 $\Phi$ 源项
用TUI命令：/define/models/energy yes yes yes（最后的yes是viscous dissipation）
OpenFOAM
OpenFOAM中的处理因使用的求解器不同而异。
- buoyantPimpleFoam：能量方程包含粘性散逸
- chtMultiRegionFoam：流体区域的能量方程支持
- 自定义求解器可用
fvm::Sp(mu*magSqr(symm(fvc::grad(U))), he)这样的项添加
STAR-CCM+
在STAR-CCM+中：
1. Physics > Energy Model 启用
2. Physics > Viscous Dissipation Model 启用
3. 自动集成到能量方程中
验证问题：Couette流中的粘性散逸
用有解析解的问题来验证。在上壁以速度 $U$ 移动的Couette流中，定常温度分布是：
上下壁都固定为温度 $T_0$ 时，最高温度在中央（$y = H/2$）：
与解析解对比，就能验证软件的粘性散逸实现是否正确。
完全正确。在2D平行板中进行足够精细的网格计算（$H$ 方向至少50个单元），验证与解析解的一致性后再进行实际问题。对水（$\mu = 10^{-3}\,\text{Pa}\cdot\text{s}$，$k = 0.6\,\text{W/(m}\cdot\text{K)}$）在 $U = 10\,\text{m/s}$，$H = 1\,\text{mm}$ 时，$\Delta T \approx 0.02\,\text{K}$ 非常小。对高分子熔体（$\mu = 1000\,\text{Pa}\cdot\text{s}$，$k = 0.2$），$\Delta T \approx 62,500\,\text{K}$ 在计算上看起来很大，但实际上粘度会因剪切加热而下降，存在反馈效应。
粘性散逸使咖啡升温——搅拌的意外物理
用勺子搅拌咖啡时，温度会略微升高。这就是粘性散逸的真实体现。速度梯度越大，散逸越大，转化为热的能量越多。实际上，用1000rpm的工业搅拌器搅拌高粘性膏体（η≈100 Pa·s）时，一分钟内温度可升高数摄氏度。食品和制药行业的工艺设计中，"散逸热导致的品质劣化"是一个严重问题，搅拌叶片的最优化设计本质上就是对粘性散逸分布的控制。看似无关的话题，实际上这个项扮演着关键角色。
粘性散逸的软件对比
支持粘性散逸的工具对比
哪些软件对粘性散逸计算支持最好？
粘性散逸本身几乎所有CFD求解器都支持，但在非牛顿流体耦合和温度相关粘度处理上有差别。
| 工具 | 粘性散逸 | 非牛顿耦合 | 温度相关粘度 | 特别说明 |
|---|---|---|---|---|
| Ansys Fluent | 单次复选启用 | 完全支持 | Arrhenius、WLF等 | 可用UDF完全自定义 |
| STAR-CCM+ | 模型选择启用 | 完全支持 | 各种内置模型 | Field Function灵活定义 |
| OpenFOAM | 求解器依赖 | 支持（需设置） | UDF等价物实现 | 源代码修改完全自由 |
| COMSOL | Weak形式添加 | 完全支持 | 任意公式输入 | 多物理耦合容易 |
| Ansys Polyflow | 标准内置 | 包括粘弹性 | 各种模型 | 高分子加工最充实 |
| Moldflow | 自动考虑 | 注射成型专用 | 金具DB链接 | 注射成型中默认ON |
高分子加工领域的选择
在粘性散逸特别重要的高分子加工领域，专用软件的优势明显。
Ansys Polyflow 能准确计算粘弹性流体的模具膨胀（die swell）和粘性散逸导致的温度分布。支持自由表面追踪、吹塑成型、热成型等。
Moldflow 专门用于注射成型，粘性散逸默认包含在能量方程中。标准使用Cross-WLF模型，材料数据库丰富：
Moldflow的材料数据库到底包含多少种材料？
超过10,000种塑料材料被预设。每种材料的Cross-WLF参数、PVT特性、热物性都包括在内，选择材料等级就能进行包含粘性散逸的精确仿真。
通用CFD的选择
对于润滑问题或一般热流体问题中的粘性散逸，Fluent或STAR-CCM+ 足以。选择准则如下。
| 用途 | 推荐工具 | 原因 |
|---|---|---|
| 注射成型 | Moldflow、Moldex3D | 材料DB、模具冷却耦合 |
| 挤出、吹塑 | Polyflow | 粘弹性、自由表面 |
| 润滑、轴承 | Fluent、STAR-CCM+ | 薄膜流处理 |
| 微通道 | COMSOL、Fluent | 多物理耦合 |
| 极超声速 | Fluent（Density-based） | 压缩性流 + 散逸 |
用途不同选择就不同，很清楚。
粘性散逸模型的工具实现差异——Fluent和CFX的默认设置陷阱
ANSYS公司的Fluent和CFX在粘性散逸处理上有差异。Fluent的"Energy Equation"面板中"Viscous Heating"选项在不可压流体中默认OFF，导致高速、高粘性流中容易漏掉散逸而犯错。而CFX在"Heat Transfer"设置中明确区分"Thermal Energy"（无散逸）和"Total Energy"（含散逸），选择错误少。OpenFOAM的buoyantSimpleFoam默认包含粘性散逸，但很多用户通过viscous标志关闭它，需确认设置意图。各工具的最可靠方法都是计算简单Couette流（解析解已知），验证能量守恒后再进行实际问题。
粘性散逸的前沿研究
粘弹性流体中的散逸
粘弹性流体中粘性散逸会如何变化？
在粘弹性流体中，能量散逸分为"粘性部分"和"弹性部分"。弹性能暂时储存，之后才散逸或作为功回收。
Oldroyd-B模型中的散逸函数是：
Peters & Baaijens (1997) 提出了粘弹性流体散逸的准确表述，方法是分离弹性储存和不可逆散逸。在注射成型浇口等处，弹性效应对温度场的影响可达几十K。
弹性能的储存与散逸的区分对粘弹性流体很关键。
极超音速流中的粘性散逸
在Ma > 5的极超音速流中，通过冲击波时温度上升达数千K。
正常冲击波的温度比是：
Ma = 10的情况（空气），$T_2/T_1 \approx 20$。入口 $T_1 = 250\,\text{K}$ 则冲击波后约5000 K，此时开始发生化学反应（分子解离、电离）。
在这个量级，普通Navier-Stokes方程框架不够，还需要：
- 高温气体效应：$c_p(T)$、$\mu(T)$、$k(T)$ 的温度相关性
- 化学反应：$\text{N}_2 \rightleftharpoons 2\text{N}$、$\text{O}_2 \rightleftharpoons 2\text{O}$ 的解离反应
- 辐射传热：高温气体的辐射
需要耦合求解所有这些。NASA的US3D、Dplr，以及带有化学动力学的Ansys Fluent（Density-based solver）被使用。
大气层再入本来就是粘性散逸的极端情形。
微观尺度下的粘性散逸
在微通道（$D_h < 100\,\mu\text{m}$）中，宏观尺度上可忽略的粘性散逸变得重要。
Brinkman数 $\text{Br} \propto U^2 / \Delta T$，但微通道中流速大、外部加热/冷却导致的 $\Delta T$ 小，所以Br变得不能忽略。
实验上也发现，微通道内摩擦系数与宏观尺度理论值（$f = 64/\text{Re}$）的偏差，一个原因就是粘性散逸导致的粘度变化。
研究进展
粘性散逸研究的最前沿是什么？
古典概念还有这么多研究空间。
超音速发动机内的粘性散逸——发热量超乎想象
超燃冲压发动机（极超音速燃烧冲压喷气引擎）中，进气流速达到5～10马赫。此时粘性散逸的加热极其巨大，燃烧室进口处空气温度已跳升至1000K以上。若不考虑散逸项 Φ=μ(∂u_i/∂x_j+∂u_j/∂x_i)∂u_i/∂x_j 进行分析，入口温度会被低估100～200K，导致燃料点火时序设计完全错误。尖端的超音速流体分析中，粘性散逸不是"细微修正"而是主要物理。
粘性散逸故障排除
粘性散逸计算的典型问题
启用粘性散逸后计算出问题的例子给我讲讲。
总结粘性散逸的典型问题模式。
1. 温度非物理性上升并发散
现象：局部温度变为数万K并发散
原因：粘性散逸 → 温度上升 → 粘度下降（温度相关粘度情况）→ 剪切速率增大 → 进一步散逸加大，形成正反馈循环
对策：
- 设置粘度的下限值（物理合理范围内）
- 减小时间步（非定常计算）
- 降低松弛系数（能量方程：0.7-0.8，粘度：0.5-0.7）
- 先以等温（粘度固定）收敛流场，再启用能量方程
2. 壁面附近出现温度尖峰
现象：壁面第一层单元温度异常高
原因：网格太粗，壁面的剪切速率被过高估计。因为 $\Phi \propto \mu (\partial u/\partial y)^2$，壁面第一层的速度梯度高估会导致散逸剧增。
对策：
- 充分细化壁面附近网格
- 棱柱层的增长率控制在1.2以下
- 确认 $y^+$ 值，与壁函数模型一致
网格质量直接影响散逸精度。
3. 散逸量大于理论值
现象：Couette流或Poiseuille流中，散逸导致的温度上升是解析解的好几倍
原因：数值粘性（人工散逸）叠加到物理散逸上。特别是一阶迎风差分（First Order Upwind）最明显。
对策：
- 使用二阶或更高阶的离散格式（Second Order Upwind、QUICK、Central Differencing）
- 充分细分网格，用Richardson外推验证网格收敛
- OpenFOAM中在
fvSchemes中使用Gauss linear或Gauss limitedLinear 1
4. 启用粘性散逸后收敛恶化
现象：没有散逸时300次迭代收敛，启用后5000次还不收敛
原因：高Br数时，能量方程和动量方程耦合强，分离法的交替求解收敛难
对策：
- 使用耦合求解器（Coupled Solver）
- 能量方程的松弛系数分步提升（初始0.5，稳定后0.8-0.9）
- 启用伪时间步（Pseudo Transient）
Br数越大耦合越强，收敛越难。
检查清单
粘性散逸计算的品质保证检查项：
- 事先估算Brinkman数，把握散逸的影响程度
- 用有解析解的问题（Couette流、Poiseuille流）验证
- 验证网格收敛，散逸量对网格无依赖
- 数值格式是二阶或以上吗
- 使用温度相关粘度时，粘度上下限是否合理设置
- 绝热壁 vs. 等温壁的边界条件与问题相符吗
从Br数估算开始，用验证问题确认，然后进实际问题。
完全正确。粘性散逸是物理上正确的现象，但数值人工散逸的混淆需要系统的验证。
CFD中温度无限上升与粘性散逸的二重计账及求解器收敛问题
压缩流CFD中出现"温度随迭代无限上升不收敛"的症状，可能是粘性散逸项的二重计账。Fluent中同时启用"Viscous Heating"和"Compressibility Effect"时，在某些求解器设置下散逸贡献会被重复累加，这是已知的实现问题。低速不可压流体中启用Viscous Heating会导致理论上应忽略的微小散逸累积，能量守恒崩溃。排查方法：①关闭能量方程，检查速度场是否收敛；②计算Brinkman数评价散逸重要性；③查阅手册"Energy Source Terms"实现规范。
相关主题
细节
错误