反应扩散系统熵与熵衰减计算的正确性验证请求

阿华AIGC实验室

2026-4-22

我想要解决下面这个习题：

我们考虑两个分子$U_1$可逆反应生成一个分子$U_2$的例子：

$$
\begin{pmatrix}
2 & 0\
0 & 1
\end{pmatrix}
\begin{pmatrix}
U_1 \
U_2
\end{pmatrix}
\longrightarrow
\begin{pmatrix}
0 & 1\
2 & 0
\end{pmatrix}
\begin{pmatrix}
U_1 \
U_2
\end{pmatrix}
$$

如果物质$U_i$（i=1,2）在有界区域$\Omega$内以恒定扩散系数$d_i$扩散，反应速率用质量作用定律建模且与温度无关，我们得到如下系统：

$$
\begin{cases}
\partial _t u_1 - d_1 \Delta u_1=-2(u_1 ^2-u_2)\
\partial _t u_2 - d_2 \Delta u_2=u_1 ^2-u_2
\end{cases}
$$

同时满足齐次Neumann边界条件$n(x)\nabla _{x} u_i(t,x)=0$在$\partial\Omega$上，其中$n(x)$是$\Omega$的单位外法向量。

需要计算该系统的熵以及对应的熵衰减。

我对这个主题不太熟悉，不确定自己的方法是否正确，我的计算步骤如下：

我定义熵为：
$$
\epsilon (t)=\int_{\Omega} (-u_1 \ln(u_1)-u_2 \ln(2)) dx
$$

为了分析熵衰减，我计算了$\epsilon$的时间导数：
$$
\frac{d\epsilon}{dt}=-\int_{\Omega}(\partial _t u_1 \ln(u_1)+\partial _t u_2 \ln(u_2))dx
$$

然后代入系统方程：
$$
\frac{d\epsilon}{dt}=\int_{\Omega}\left(d_1\Delta u_1 \ln(u_1)+d_2\Delta u_2 \ln(u_2)+2(u_1 ^2-u_2)\ln(u_1)+(u_1 ^2-u_2)\ln(u_2)\right)dx
$$

接下来我利用Neumann边界条件和分部积分：
$$
\frac{d\epsilon}{dt}=\int_{\Omega}\left(-d_1\nabla (u_1 \nabla \ln(u_1))-d_2 \nabla (u_2 \nabla \ln(u_2))+2(u_1 ^2-u_2)\ln(u_1)+(u_1 ^2-u_2)\ln(u_2)\right)dx
$$

之后应用散度定理处理前两项：
$$
\frac{d\epsilon}{dt}=-\int_{\Omega}(2(u_1 ^2-u_2)\ln(u_1)+(u_1 ^2-u_2)\ln(u_2))dx-\int_{\partial\Omega} (d_1 u_1\nabla (\ln(u_1))\cdot n+d_2 u_2 \nabla \ln(u_2)\cdot n )ds
$$

最后得到：
$$
\frac{d\epsilon}{dt}=-\int_{\Omega}[2(u_1 ^2-u_2)\ln(u_1)+(u_1 ^2-u_2)\ln(u_2)]dx
$$

请问我的解法正确吗？有没有哪里出错了？谢谢！

分析与修正

你的整体思路方向是对的，但在几个关键步骤上存在错误，下面逐一说明：

熵定义的笔误
你初始定义熵时写了$-u_2 \ln(2)$，但后续求导却用了$\ln(u_2)$，这显然是笔误。正确的熵函数应该符合热力学中常用的形式，对于浓度场的熵应该定义为：
$$
\epsilon(t) = \int_{\Omega} \left(-u_1 \ln u_1 - u_2 \ln u_2\right) dx
$$
这样才能和后续的求导步骤自洽。
分部积分的错误推导
处理$\int \Delta u_i \ln u_i dx$时，你错误地写成了$-\nabla(u_i \nabla \ln u_i)$，正确的分部积分过程应该是：
$$
\int_{\Omega} \Delta u_i \ln u_i dx = \int_{\partial\Omega} \nabla u_i \cdot n \ln u_i ds - \int_{\Omega} \nabla u_i \cdot \nabla \ln u_i dx
$$
结合齐次Neumann边界条件$\nabla u_i \cdot n = 0$，边界项直接消失。而$\nabla \ln u_i = \frac{\nabla u_i}{u_i}$，因此$\nabla u_i \cdot \nabla \ln u_i = \frac{|\nabla u_i|^2}{u_i}$，这是一个非负项（因为浓度$u_i>0$，梯度平方非负）。
漏掉扩散项的耗散贡献
你最终的结果完全丢掉了扩散项的积分，这是严重的疏漏。扩散过程是耗散过程，必然会贡献熵衰减的一部分。正确的熵时间导数应该包含扩散和反应两部分的耗散项：
$$
\frac{d\epsilon}{dt}=-\int_{\Omega}\left(d_1 \frac{|\nabla u_1|^2}{u_1} + d_2 \frac{|\nabla u_2|^2}{u_2}\right)dx - \int_{\Omega}(u_1^2 - u_2)\ln\left(u_1^2 u_2\right)dx
$$
反应项的整理与熵衰减验证
反应项可以整理为$(u_1^2 - u_2)\ln(u_1^2 u_2)$，结合平衡态条件$u_1^2 = u_2$时该项为0；当系统偏离平衡时，该项的符号始终保证整个积分非正，和扩散项一起满足熵衰减（$\frac{d\epsilon}{dt} \leq 0$）的热力学要求。