不同输送量(动量、能量和质量)的守恒定律可以从一些简单的原理推导出来。举个例子,假设一个在系统中守恒的量 Φ,其通量由矢量 j = (jx, jy, jz) 给出。通过尺寸为 Δx、Δy 和 Δz 的微小体积单元,可以得到 φ 的平衡,其中的通量 j 是单位面积、单位时间的量。产生项或消耗项 Rs 是单位体积、单位时间的量。

根据上图,可以得到以下平衡方程:

$$ \left. \Delta x\Delta y\Delta z\frac{\varphi _{t+\Delta t}-\varphi _t}{\Delta t}=-\left( \Delta y\Delta z\left( j_{x,x+\Delta x}-j_{x,x} \right) +\Delta x\Delta z\left( j_{y\cdot y+\Delta y}-j_{y\cdot y} \right) \right) +\Delta x\Delta y\left( j_{z,z+\Delta z}-j_{z,z} \right) \right) +\,\,\Delta x\Delta y\Delta zR_s $$

其中,jx,x 表示 x 位置 x 方向上的通量,jx,x+Δx 是 x+Δx 位置 x 方向上 的通量矢量。下标 t 和 t+Δt 表示相应时间的状态。该方程表明,如果没有产生或消耗(Rs = 0),则流入或流出体积单元的净通量必须通过该物理量在一段时间内的积累或消耗进行平衡。此外,如果在一段时间内没有积累或消耗,那么进入体积单元的通量必须完全平衡流出该体积单元的通量,使进出该体积单元的总通量为零。

将以上方程除以体积 ΔxΔyΔz,并使 Δx、Δy 和 Δz 都接近零,可以得到 φ 的以下方程:

$\frac{\varphi_{t+\Delta t}-\varphi_{t}}{\Delta t}=-\underbrace{\left(\frac{\partial j_{x}}{\partial x}+\frac{\partial j_{y}}{\partial y}+\frac{\partial j_{z}}{\partial z}\right)}_{\nabla \cdot \mathrm{j}}+R_{s}$

使 Δt 接近 0,可得:

$\frac{\partial \varphi}{\partial t}+\nabla \cdot \mathbf{j}-R_{s}=0$

为了遵循守恒定律,所研究系统中的连续介质(如流体或固体)中的每一微小体积单元都必须满足这一平衡方程。在传递现象的建模与仿真中,这个偏微分方程——“偏”是因为它表示为每次相对于一个自变量的变化而变化(x、y、z 和 t 为自变量)——能够描述动量、能量和质量的守恒定律。

在动量守恒中,守恒量是一个矢量,而通量项则以张量形式表示,其中包含应力张量。将动量守恒、动量通量的本构方程以及不可压缩牛顿流体的质量守恒相结合,可以得到纳维-斯托克斯方程。这些方程是流体流动(CFD)建模的基础,它们的解描述了流动流体的速度和压力场。如果守恒量为能量,则系统中的传热方程也可通过上述守恒方程导出。

最后,我们来看看质量传递。假设我们要研究同时存在传递和反应过程的流体的成分。为此,我们可以定义和求解流体中每种物质的质量守恒方程。每种物质 i 的浓度 ci 均保持守恒,其通量表示为 Ni。使用上面的守恒方程可以得到以下适用于每种物质的方程:

$\frac{\partial c_{i}}{\partial t}+\nabla \cdot \mathbf{N}_{i}-R_{i}=0$

如果我们进一步假设通量由扩散给出,并且菲克定律可定义该通量(例如溶剂中的溶质通量)的本构关系,则可以得到扩散-反应方程,常用于模拟可忽略流体流动的反应系统:

$\frac{\partial c_{i}}{\partial t}+\nabla \cdot\left(-D_{i} \nabla c_{i}\right)-R_{i}=0$

在这一方程中,Di 表示溶液中物质 i 的扩散系数。

如果发生对流,也就是整个溶液的净输送,则可以得到输运方程,常用于存在流体流动的反应系统:

$\frac{\partial c_{i}}{\partial t}+\nabla \cdot \underbrace{\left(-D_{i} \nabla c_{i}+c_{i} \mathbf{u}\right)}_{\mathbf{N}_{i}}-R_{i}=0$

在这个方程中,u 表示速度矢量。如果溶液受外加电场 E 作用,并且存在离子,则可以得到 Nernst-Planck 方程,可用于电化学系统:

$\frac{\partial c_{i}}{\partial t}+\nabla \cdot \underbrace{\left(-D_{i} \nabla c_{i}+c_{i} \mathbf{u}+z_{i} u_{i} c_{i} \mathbf{E}\right)}_{\mathbf{N}_{i}}-R_{i}=0$

在这个方程中,zi 表示物质 i 的化合价,ui 表示物质 i 的迁移率,通过 Nernst-Einstein 关系直接与扩散率相关。通量矢量中的第三项称为迁移项。

综上所述,定义模型方程所依据的原理非常简单,问题的关键是定义守恒定律以及用于产生通量的关系。为此,我们需要在不同的条件下反复求解给定系统的方程,然后研究计算结果,从而了解系统中的传递现象。

最后修改:2023 年 08 月 20 日
如果觉得我的文章对你有用,请随意赞赏