chns方程推导简单的能量耗散律，分部积分向量形式，sav初简介

文章目录

chns方程
- 能量耗散律
- [🎯 推导目标](#🎯 推导目标)
- [📋 前提条件](#📋 前提条件)
- [🔍 逐步推导](#🔍 逐步推导)
- [🎉 最终结论](#🎉 最终结论)
- [💡 物理意义](#💡 物理意义)
分部积分的向量形式（并非格林公式）
- [你观察得非常仔细！这确实是一个容易混淆的地方。你之所以觉得它和格林公式不一样，是因为你看到的是**格林公式的向量形式**，而我们通常学的是**标量形式**。让我为你彻底拆解这个过程，你会发现它们本质上是相同的。](#你观察得非常仔细！这确实是一个容易混淆的地方。你之所以觉得它和格林公式不一样，是因为你看到的是格林公式的向量形式，而我们通常学的是标量形式。让我为你彻底拆解这个过程，你会发现它们本质上是相同的。)
- [🎯 核心问题：标量 vs 向量形式的格林公式](#🎯 核心问题：标量 vs 向量形式的格林公式)
- - **你熟悉的格林公式（标量形式）**
- [对于两个**标量函数** P P P 和 Q Q Q，二维格林公式是： ∬ D ( ∂ Q ∂ x − ∂ P ∂ y ) d x d y = ∮ ∂ D ( P d x + Q d y ) \iint_D \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dxdy = \oint_{\partial D} (Pdx + Qdy) ∬D(∂x∂Q−∂y∂P)dxdy=∮∂D(Pdx+Qdy) 或者更常用的**散度形式**： ∬ D ∇ ⋅ F d A = ∮ ∂ D F ⋅ n d s \iint_D \nabla \cdot \mathbf{F} dA = \oint_{\partial D} \mathbf{F} \cdot \mathbf{n} ds ∬D∇⋅FdA=∮∂DF⋅nds 其中 F = ( P , Q ) \mathbf{F} = (P, Q) F=(P,Q) 是一个**向量场**。](#对于两个标量函数 P P P 和 Q Q Q，二维格林公式是： ∬ D ( ∂ Q ∂ x − ∂ P ∂ y ) d x d y = ∮ ∂ D ( P d x + Q d y ) \iint_D \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dxdy = \oint_{\partial D} (Pdx + Qdy) ∬D(∂x∂Q−∂y∂P)dxdy=∮∂D(Pdx+Qdy) 或者更常用的散度形式： ∬ D ∇ ⋅ F d A = ∮ ∂ D F ⋅ n d s \iint_D \nabla \cdot \mathbf{F} dA = \oint_{\partial D} \mathbf{F} \cdot \mathbf{n} ds ∬D∇⋅FdA=∮∂DF⋅nds 其中 F = ( P , Q ) \mathbf{F} = (P, Q) F=(P,Q) 是一个向量场。)
- [🔍 你看到的分部积分：向量场点乘标量函数的梯度](#🔍 你看到的分部积分：向量场点乘标量函数的梯度)
- [📚 分步推导：从第一性原理出发](#📚 分步推导：从第一性原理出发)
- [将结果代回： ∫ Ω u ⋅ ∇ f d x = ∮ ∂ Ω f ( u ⋅ n ) d S − ∫ Ω f ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla f dx = \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=∮∂Ωf(u⋅n)dS−∫Ωf(∇⋅u)dx 在流体力学问题中，通常有**无穿透边界条件** u ⋅ n = 0 \mathbf{u} \cdot \mathbf{n} = 0 u⋅n=0（流体不流过边界）或者**固壁边界条件** u = 0 \mathbf{u} = \mathbf{0} u=0。在这两种情况下，边界项都为零： ∮ ∂ Ω f ( u ⋅ n ) d S = 0 \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS = 0 ∮∂Ωf(u⋅n)dS=0 因此，我们得到最终的分部积分公式： ∫ Ω u ⋅ ∇ f d x = − ∫ Ω f ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla f dx = - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=−∫Ωf(∇⋅u)dx](#将结果代回： ∫ Ω u ⋅ ∇ f d x = ∮ ∂ Ω f ( u ⋅ n ) d S − ∫ Ω f ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla f dx = \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=∮∂Ωf(u⋅n)dS−∫Ωf(∇⋅u)dx 在流体力学问题中，通常有无穿透边界条件 u ⋅ n = 0 \mathbf{u} \cdot \mathbf{n} = 0 u⋅n=0（流体不流过边界）或者固壁边界条件 u = 0 \mathbf{u} = \mathbf{0} u=0。在这两种情况下，边界项都为零： ∮ ∂ Ω f ( u ⋅ n ) d S = 0 \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS = 0 ∮∂Ωf(u⋅n)dS=0 因此，我们得到最终的分部积分公式： ∫ Ω u ⋅ ∇ f d x = − ∫ Ω f ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla f dx = - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=−∫Ωf(∇⋅u)dx)
- [🎨 应用到你的具体问题](#🎨 应用到你的具体问题)
- [在你的例子中， f = ∣ u ∣ 2 f = |\mathbf{u}|^2 f=∣u∣2，所以： ∫ Ω u ⋅ ∇ ( ∣ u ∣ 2 ) d x = − ∫ Ω ∣ u ∣ 2 ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla (|\mathbf{u}|^2) dx = - \int_{\Omega} |\mathbf{u}|^2 (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇(∣u∣2)dx=−∫Ω∣u∣2(∇⋅u)dx 两边乘以 − 1 2 -\frac{1}{2} −21，就得到了你看到的那个变换： − 1 2 ∫ Ω u ⋅ ∇ ( ∣ u ∣ 2 ) d x = 1 2 ∫ Ω ∣ u ∣ 2 ( ∇ ⋅ u ) d x -\frac{1}{2} \int_{\Omega} \mathbf{u} \cdot \nabla (|\mathbf{u}|^2) dx = \frac{1}{2} \int_{\Omega} |\mathbf{u}|^2 (\nabla \cdot \mathbf{u}) dx −21∫Ωu⋅∇(∣u∣2)dx=21∫Ω∣u∣2(∇⋅u)dx](#在你的例子中， f = ∣ u ∣ 2 f = |\mathbf{u}|^2 f=∣u∣2，所以： ∫ Ω u ⋅ ∇ ( ∣ u ∣ 2 ) d x = − ∫ Ω ∣ u ∣ 2 ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla (|\mathbf{u}|^2) dx = - \int_{\Omega} |\mathbf{u}|^2 (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇(∣u∣2)dx=−∫Ω∣u∣2(∇⋅u)dx 两边乘以 − 1 2 -\frac{1}{2} −21，就得到了你看到的那个变换： − 1 2 ∫ Ω u ⋅ ∇ ( ∣ u ∣ 2 ) d x = 1 2 ∫ Ω ∣ u ∣ 2 ( ∇ ⋅ u ) d x -\frac{1}{2} \int_{\Omega} \mathbf{u} \cdot \nabla (|\mathbf{u}|^2) dx = \frac{1}{2} \int_{\Omega} |\mathbf{u}|^2 (\nabla \cdot \mathbf{u}) dx −21∫Ωu⋅∇(∣u∣2)dx=21∫Ω∣u∣2(∇⋅u)dx)
- [💡 与一维分部积分的类比](#💡 与一维分部积分的类比)
- [🎯 总结](#🎯 总结)
sav的离散与求解

chns方程

能量耗散律

逐步推导这个能量耗散定律的证明过程。这是CHNS方程的一个重要性质，体现了系统的物理稳定性。

🎯 推导目标

证明： d E d t = − γ ∥ ∇ μ ∥ L 2 2 − ν ∥ ∇ u ∥ L 2 2 ≤ 0 \frac{dE}{dt} = -\gamma\|\nabla\mu\|{L^2}^2 - \nu\|\nabla u\|{L^2}^2 \leq 0 dtdE=−γ∥∇μ∥L22−ν∥∇u∥L22≤0

📋 前提条件

总能量定义 ： E ( ϕ , u ) = ∫ Ω [ 1 2 ( ∣ u ∣ 2 + λ ∣ ∇ ϕ ∣ 2 ) + λ F ( ϕ ) ] d x E(\phi,u) = \int_\Omega \left[\frac{1}{2}(|u|^2 + \lambda|\nabla\phi|^2) + \lambda F(\phi)\right]dx E(ϕ,u)=∫Ω[21(∣u∣2+λ∣∇ϕ∣2)+λF(ϕ)]dx
边界条件 ： ∂ ϕ ∂ n = ∂ μ ∂ n = 0 , u = 0 \frac{\partial\phi}{\partial n} = \frac{\partial\mu}{\partial n} = 0, u = 0 ∂n∂ϕ=∂n∂μ=0,u=0 (在边界上)

🔍 逐步推导

第一步：对能量函数求时间导数

d E d t = ∫ Ω [ u ⋅ ∂ u ∂ t + λ ∇ ϕ ⋅ ∂ ∂ t ( ∇ ϕ ) + λ F ′ ( ϕ ) ∂ ϕ ∂ t ] d x \frac{dE}{dt} = \int_\Omega \left[u \cdot \frac{\partial u}{\partial t} + \lambda\nabla\phi \cdot \frac{\partial}{\partial t}(\nabla\phi) + \lambda F'(\phi)\frac{\partial\phi}{\partial t}\right]dx dtdE=∫Ω[u⋅∂t∂u+λ∇ϕ⋅∂t∂(∇ϕ)+λF′(ϕ)∂t∂ϕ]dx
关键变换：

∂ ∂ t ∣ ∇ ϕ ∣ 2 = 2 ∇ ϕ ⋅ ∂ ∂ t ( ∇ ϕ ) \frac{\partial}{\partial t}|\nabla\phi|^2 = 2\nabla\phi \cdot \frac{\partial}{\partial t}(\nabla\phi) ∂t∂∣∇ϕ∣2=2∇ϕ⋅∂t∂(∇ϕ)
F ′ ( ϕ ) = g ( ϕ ) F'(\phi) = g(\phi) F′(ϕ)=g(ϕ) (根据方程1.1b)

第二步：利用控制方程替换时间导数

替换时间倒数！

从方程(1.1a)： ∂ ϕ ∂ t = γ Δ μ − u ⋅ ∇ ϕ \frac{\partial\phi}{\partial t} = \gamma\Delta\mu - u \cdot \nabla\phi ∂t∂ϕ=γΔμ−u⋅∇ϕ

从动量方程(1.1c)： ∂ u ∂ t = − u ⋅ ∇ u − ∇ p + ν Δ u + μ ∇ ϕ \frac{\partial u}{\partial t} = -u \cdot \nabla u - \nabla p + \nu\Delta u + \mu\nabla\phi ∂t∂u=−u⋅∇u−∇p+νΔu+μ∇ϕ
代入后得到 ：
d E d t = ∫ Ω [ u ⋅ ( − u ⋅ ∇ u − ∇ p + ν Δ u + μ ∇ ϕ ) + λ ∇ ϕ ⋅ ∇ ( γ Δ μ − u ⋅ ∇ ϕ ) + λ g ( ϕ ) ( γ Δ μ − u ⋅ ∇ ϕ ) ] d x \frac{dE}{dt} = \int_\Omega \left[u \cdot (-u \cdot \nabla u - \nabla p + \nu\Delta u + \mu\nabla\phi) + \lambda\nabla\phi \cdot \nabla(\gamma\Delta\mu - u \cdot \nabla\phi) + \lambda g(\phi)(\gamma\Delta\mu - u \cdot \nabla\phi)\right]dx dtdE=∫Ω[u⋅(−u⋅∇u−∇p+νΔu+μ∇ϕ)+λ∇ϕ⋅∇(γΔμ−u⋅∇ϕ)+λg(ϕ)(γΔμ−u⋅∇ϕ)]dx

第三步：分部积分处理关键项

处理流体动能项

∫ Ω u ⋅ ( − u ⋅ ∇ u ) d x = − 1 2 ∫ Ω u ⋅ ∇ ( ∣ u ∣ 2 ) d x = 1 2 ∫ Ω ∣ u ∣ 2 ( ∇ ⋅ u ) d x = 0 \int_\Omega u \cdot (-u \cdot \nabla u)dx = -\frac{1}{2}\int_\Omega u \cdot \nabla(|u|^2)dx = \frac{1}{2}\int_\Omega |u|^2(\nabla \cdot u)dx = 0 ∫Ωu⋅(−u⋅∇u)dx=−21∫Ωu⋅∇(∣u∣2)dx=21∫Ω∣u∣2(∇⋅u)dx=0
(利用 ∇ ⋅ u = 0 \nabla \cdot u = 0 ∇⋅u=0)
∫ Ω u ⋅ ( − ∇ p ) d x = − ∫ ∂ Ω p ( u ⋅ n ) d S + ∫ Ω p ( ∇ ⋅ u ) d x = 0 \int_\Omega u \cdot (-\nabla p)dx = -\int_{\partial\Omega} p(u \cdot n)dS + \int_\Omega p(\nabla \cdot u)dx = 0 ∫Ωu⋅(−∇p)dx=−∫∂Ωp(u⋅n)dS+∫Ωp(∇⋅u)dx=0
(利用边界条件 u = 0 u = 0 u=0和 ∇ ⋅ u = 0 \nabla \cdot u = 0 ∇⋅u=0)

处理粘性项

∫ Ω u ⋅ ( ν Δ u ) d x = ν ∫ ∂ Ω u ⋅ ∂ u ∂ n d S − ν ∫ Ω ∣ ∇ u ∣ 2 d x = − ν ∥ ∇ u ∥ L 2 2 \int_\Omega u \cdot (\nu\Delta u)dx = \nu\int_{\partial\Omega} u \cdot \frac{\partial u}{\partial n}dS - \nu\int_\Omega |\nabla u|^2dx = -\nu\|\nabla u\|_{L^2}^2 ∫Ωu⋅(νΔu)dx=ν∫∂Ωu⋅∂n∂udS−ν∫Ω∣∇u∣2dx=−ν∥∇u∥L22
(利用边界条件 u = 0 u = 0 u=0)

第四步：处理相场相关项

处理化学势耦合项

∫ Ω u ⋅ ( μ ∇ ϕ ) d x + ∫ Ω λ g ( ϕ ) ( − u ⋅ ∇ ϕ ) d x \int_\Omega u \cdot (\mu\nabla\phi)dx + \int_\Omega \lambda g(\phi)(-u \cdot \nabla\phi)dx ∫Ωu⋅(μ∇ϕ)dx+∫Ωλg(ϕ)(−u⋅∇ϕ)dx

利用方程(1.1b)： μ = − λ ( Δ ϕ − g ( ϕ ) ) \mu = -\lambda(\Delta\phi - g(\phi)) μ=−λ(Δϕ−g(ϕ))，代入后：
∫ Ω u ⋅ [ − λ ( Δ ϕ − g ( ϕ ) ) ∇ ϕ ] d x − λ ∫ Ω g ( ϕ ) ( u ⋅ ∇ ϕ ) d x \int_\Omega u \cdot [-\lambda(\Delta\phi - g(\phi))\nabla\phi]dx - \lambda\int_\Omega g(\phi)(u \cdot \nabla\phi)dx ∫Ωu⋅[−λ(Δϕ−g(ϕ))∇ϕ]dx−λ∫Ωg(ϕ)(u⋅∇ϕ)dx
关键观察：这两项相互抵消！

处理扩散项

∫ Ω λ ∇ ϕ ⋅ ∇ ( γ Δ μ ) d x + ∫ Ω λ g ( ϕ ) ( γ Δ μ ) d x \int_\Omega \lambda\nabla\phi \cdot \nabla(\gamma\Delta\mu)dx + \int_\Omega \lambda g(\phi)(\gamma\Delta\mu)dx ∫Ωλ∇ϕ⋅∇(γΔμ)dx+∫Ωλg(ϕ)(γΔμ)dx

对第一项分部积分：
λ γ ∫ ∂ Ω ∇ ϕ ⋅ ∂ μ ∂ n d S − λ γ ∫ Ω Δ ϕ ⋅ Δ μ d x \lambda\gamma\int_{\partial\Omega} \nabla\phi \cdot \frac{\partial\mu}{\partial n}dS - \lambda\gamma\int_\Omega \Delta\phi \cdot \Delta\mu dx λγ∫∂Ω∇ϕ⋅∂n∂μdS−λγ∫ΩΔϕ⋅Δμdx

利用边界条件 ∂ μ ∂ n = 0 \frac{\partial\mu}{\partial n} = 0 ∂n∂μ=0和方程(1.1b)：
= − λ γ ∫ Ω Δ ϕ ⋅ Δ μ d x = − γ ∫ Ω μ ⋅ Δ μ d x = -\lambda\gamma\int_\Omega \Delta\phi \cdot \Delta\mu dx = -\gamma\int_\Omega \mu \cdot \Delta\mu dx =−λγ∫ΩΔϕ⋅Δμdx=−γ∫Ωμ⋅Δμdx

再分部积分：
= − γ ∫ ∂ Ω μ ∂ μ ∂ n d S + γ ∫ Ω ∣ ∇ μ ∣ 2 d x = γ ∥ ∇ μ ∥ L 2 2 = -\gamma\int_{\partial\Omega} \mu\frac{\partial\mu}{\partial n}dS + \gamma\int_\Omega |\nabla\mu|^2dx = \gamma\|\nabla\mu\|_{L^2}^2 =−γ∫∂Ωμ∂n∂μdS+γ∫Ω∣∇μ∣2dx=γ∥∇μ∥L22

第五步：合并所有结果

将所有处理后的项合并：
d E d t = − ν ∥ ∇ u ∥ L 2 2 − γ ∥ ∇ μ ∥ L 2 2 \frac{dE}{dt} = -\nu\|\nabla u\|{L^2}^2 - \gamma\|\nabla\mu\|{L^2}^2 dtdE=−ν∥∇u∥L22−γ∥∇μ∥L22

🎉 最终结论

d E d t = − γ ∥ ∇ μ ∥ L 2 2 − ν ∥ ∇ u ∥ L 2 2 ≤ 0 \frac{dE}{dt} = -\gamma\|\nabla\mu\|{L^2}^2 - \nu\|\nabla u\|{L^2}^2 \leq 0 dtdE=−γ∥∇μ∥L22−ν∥∇u∥L22≤0

💡 物理意义

能量耗散：系统总能量随时间单调递减
稳定性保证：这个性质是数值格式设计的基础
物理一致性 ：反映了粘性耗散和相场扩散的不可逆性
这个推导展示了CHNS方程的优美数学结构，以及物理守恒律与数学分析的完美结合。

分部积分的向量形式（并非格林公式）

你观察得非常仔细！这确实是一个容易混淆的地方。你之所以觉得它和格林公式不一样，是因为你看到的是格林公式的向量形式，而我们通常学的是标量形式。

让我为你彻底拆解这个过程，你会发现它们本质上是相同的。

🎯 核心问题：标量 vs 向量形式的格林公式

你熟悉的格林公式（标量形式）

对于两个标量函数 P P P 和 Q Q Q，二维格林公式是：

∬ D ( ∂ Q ∂ x − ∂ P ∂ y ) d x d y = ∮ ∂ D ( P d x + Q d y ) \iint_D \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dxdy = \oint_{\partial D} (Pdx + Qdy) ∬D(∂x∂Q−∂y∂P)dxdy=∮∂D(Pdx+Qdy)
或者更常用的散度形式 ：
∬ D ∇ ⋅ F d A = ∮ ∂ D F ⋅ n d s \iint_D \nabla \cdot \mathbf{F} dA = \oint_{\partial D} \mathbf{F} \cdot \mathbf{n} ds ∬D∇⋅FdA=∮∂DF⋅nds
其中 F = ( P , Q ) \mathbf{F} = (P, Q) F=(P,Q) 是一个向量场。

🔍 你看到的分部积分：向量场点乘标量函数的梯度

我们要处理的积分是：
I = ∫ Ω u ⋅ ∇ f d x I = \int_{\Omega} \mathbf{u} \cdot \nabla f dx I=∫Ωu⋅∇fdx

其中：

u \mathbf{u} u 是向量场（如流体速度）
f f f 是标量函数 （在你的例子中， f = ∣ u ∣ 2 f = |\mathbf{u}|^2 f=∣u∣2）
这个形式不能直接套用标量格林公式，但我们可以用向量微积分的乘法法则 和散度定理来推导。

📚 分步推导：从第一性原理出发

第一步：使用乘法法则（关键技巧）

我们要构造一个散度项。利用标量与向量乘积的散度公式 ：
∇ ⋅ ( f u ) = ( ∇ f ) ⋅ u + f ( ∇ ⋅ u ) \nabla \cdot (f \mathbf{u}) = (\nabla f) \cdot \mathbf{u} + f (\nabla \cdot \mathbf{u}) ∇⋅(fu)=(∇f)⋅u+f(∇⋅u)

重新整理这个公式，得到我们需要的项：
u ⋅ ∇ f = ∇ ⋅ ( f u ) − f ( ∇ ⋅ u ) \mathbf{u} \cdot \nabla f = \nabla \cdot (f \mathbf{u}) - f (\nabla \cdot \mathbf{u}) u⋅∇f=∇⋅(fu)−f(∇⋅u)

第二步：对两边积分

∫ Ω u ⋅ ∇ f d x = ∫ Ω ∇ ⋅ ( f u ) d x − ∫ Ω f ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla f dx = \int_{\Omega} \nabla \cdot (f \mathbf{u}) dx - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=∫Ω∇⋅(fu)dx−∫Ωf(∇⋅u)dx

第三步：应用散度定理（高斯定理）

对右边的第一项应用散度定理：
∫ Ω ∇ ⋅ ( f u ) d x = ∮ ∂ Ω ( f u ) ⋅ n d S \int_{\Omega} \nabla \cdot (f \mathbf{u}) dx = \oint_{\partial \Omega} (f \mathbf{u}) \cdot \mathbf{n} dS ∫Ω∇⋅(fu)dx=∮∂Ω(fu)⋅ndS

其中 n \mathbf{n} n 是边界 ∂ Ω \partial \Omega ∂Ω 的单位外法向量。

第四步：合并并应用边界条件

将结果代回：

∫ Ω u ⋅ ∇ f d x = ∮ ∂ Ω f ( u ⋅ n ) d S − ∫ Ω f ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla f dx = \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=∮∂Ωf(u⋅n)dS−∫Ωf(∇⋅u)dx
在流体力学问题中，通常有无穿透边界条件 u ⋅ n = 0 \mathbf{u} \cdot \mathbf{n} = 0 u⋅n=0（流体不流过边界）或者固壁边界条件 u = 0 \mathbf{u} = \mathbf{0} u=0。在这两种情况下，边界项都为零：
∮ ∂ Ω f ( u ⋅ n ) d S = 0 \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS = 0 ∮∂Ωf(u⋅n)dS=0
因此，我们得到最终的分部积分公式：
∫ Ω u ⋅ ∇ f d x = − ∫ Ω f ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla f dx = - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=−∫Ωf(∇⋅u)dx

🎨 应用到你的具体问题

在你的例子中， f = ∣ u ∣ 2 f = |\mathbf{u}|^2 f=∣u∣2，所以：

∫ Ω u ⋅ ∇ ( ∣ u ∣ 2 ) d x = − ∫ Ω ∣ u ∣ 2 ( ∇ ⋅ u ) d x \int_{\Omega} \mathbf{u} \cdot \nabla (|\mathbf{u}|^2) dx = - \int_{\Omega} |\mathbf{u}|^2 (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇(∣u∣2)dx=−∫Ω∣u∣2(∇⋅u)dx
两边乘以 − 1 2 -\frac{1}{2} −21，就得到了你看到的那个变换：
− 1 2 ∫ Ω u ⋅ ∇ ( ∣ u ∣ 2 ) d x = 1 2 ∫ Ω ∣ u ∣ 2 ( ∇ ⋅ u ) d x -\frac{1}{2} \int_{\Omega} \mathbf{u} \cdot \nabla (|\mathbf{u}|^2) dx = \frac{1}{2} \int_{\Omega} |\mathbf{u}|^2 (\nabla \cdot \mathbf{u}) dx −21∫Ωu⋅∇(∣u∣2)dx=21∫Ω∣u∣2(∇⋅u)dx

💡 与一维分部积分的类比

这个公式其实是你熟悉的一维分部积分的自然推广。
一维情况 ：
∫ a b u ( x ) f ′ ( x ) d x = [ u ( x ) f ( x ) ] a b − ∫ a b u ′ ( x ) f ( x ) d x \int_a^b u(x) f'(x) dx = [u(x)f(x)]a^b - \int_a^b u'(x) f(x) dx ∫abu(x)f′(x)dx=[u(x)f(x)]ab−∫abu′(x)f(x)dx
向量情况 ：
∫ Ω u ⋅ ∇ f d x = ∮ ∂ Ω f ( u ⋅ n ) d S − ∫ Ω f ( ∇ ⋅ u ) d x \int{\Omega} \mathbf{u} \cdot \nabla f dx = \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS - \int_{\Omega} f (\nabla \cdot \mathbf{u}) dx ∫Ωu⋅∇fdx=∮∂Ωf(u⋅n)dS−∫Ωf(∇⋅u)dx

一维	向量
u ( x ) u(x) u(x)	u \mathbf{u} u (向量场)
f ′ ( x ) f'(x) f′(x)	∇ f \nabla f ∇f (梯度)
u ′ ( x ) u'(x) u′(x)	∇ ⋅ u \nabla \cdot \mathbf{u} ∇⋅u (散度)
[ u ( x ) f ( x ) ] a b [u(x)f(x)]_a^b [u(x)f(x)]ab	∮ ∂ Ω f ( u ⋅ n ) d S \oint_{\partial \Omega} f (\mathbf{u} \cdot \mathbf{n}) dS ∮∂Ωf(u⋅n)dS (边界项)

🎯 总结

你是对的 ：这个形式和你学的标量格林公式看起来不一样。
本质相同 ：它是格林公式（或散度定理）在向量场点乘标量梯度这一特定情况下的应用。
关键技巧 ：使用乘法法则 ∇ ⋅ ( f u ) = ( ∇ f ) ⋅ u + f ( ∇ ⋅ u ) \nabla \cdot (f \mathbf{u}) = (\nabla f) \cdot \mathbf{u} + f (\nabla \cdot \mathbf{u}) ∇⋅(fu)=(∇f)⋅u+f(∇⋅u) 来构造散度项。
边界条件 ：在流体问题中，边界项通常为零，使得公式变得简洁。
这个推导是流体力学和PDE分析中的基本技巧，掌握了它，你就能看懂很多能量估计和稳定性证明了。希望这个解释能帮你彻底理解！

chns方程 推导简单的能量耗散律，分部积分向量形式，sav初简介