质量守恒定律的数学表达

质量守恒定律是固体力学与连续介质力学的基石之一,它表明在任意闭合系统中,质量既不会凭空产生,也不会凭空消失。该定律在数学上的严谨表达为质量守恒方程,它既可以写成积分形式,也可以写成微分形式。本文围绕质量守恒定律的数学表达展开,系统阐述其推导、适用范围以及典型应用示例,帮助读者在实际工程分析中准确使用。

质量守恒定律的基本形式

在任意固定控制体 (V)(体积随时间不变)内,质量守恒可以表述为:

[
\frac{d}{dt}\int_{V}\rho , dV = -\int_{\partial V}\rho \mathbf{v}\cdot\mathbf{n}, dS
]

  • (\rho(\mathbf{x},t)):材料密度(kg·m(^{-3}))
  • (\mathbf{v}(\mathbf{x},t)):质点速度向量(m·s(^{-1}))
  • (\mathbf{n}):控制面外法向量
  • (\partial V):控制体的边界面

左侧是控制体内质量随时间的变化率,右侧是通过控制面流出的质量通量。该式是质量守恒的积分形式。

连续介质中的质量守恒

在固体力学中,材料通常被视为连续介质。对任意微小体元 (dV) 进行体积分守恒的局部化,可得到微分形式:

[
\frac{\partial \rho}{\partial t} + \nabla!\cdot(\rho \mathbf{v}) = 0
]

该方程的含义是:在任意时刻、任意空间点,密度的局部变化率与质量通量的散度相互抵消。

常见特例

  • 不可压材料(如多数金属)
    (\rho = \rho_0 = \text{常数}),则方程化简为
    [
    \nabla!\cdot \mathbf{v}=0
    ]
    表示速度场是无散的(体积保持)。

  • 可压材料(如橡胶、土体)
    必须保留 (\partial \rho/\partial t) 项,才能描述密度随应力或温度的变化。

积分形式与微分形式的转换

利用高斯散度定理((\displaystyle \int_{V}\nabla!\cdot\mathbf{A},dV = \int_{\partial V}\mathbf{A}\cdot\mathbf{n},dS)),可将积分形式转化为微分形式:

[
\begin{aligned}
\frac{d}{dt}\int_{V}\rho , dV
&= \int_{V}\frac{\partial \rho}{\partial t}, dV \
&= -\int_{\partial V}\rho \mathbf{v}\cdot\mathbf{n}, dS \
&= -\int_{V}\nabla!\cdot(\rho \mathbf{v}), dV
\end{aligned}
]

因为控制体 (V) 任意,内部积分的被积函数必须相等,遂得微分形式。

示例:一维稳态流动中的质量守恒

考虑一根截面积恒定的管道,流体沿 (x) 方向稳态流动((\partial/\partial t = 0)),密度随 (x) 变化。质量守恒微分方程简化为:

[
\frac{d}{dx}(\rho u)=0
]

其中 (u(x)) 为流速。积分得到:

[
\rho(x) , u(x) = \text{常数}= \dot{m}
]

(\dot{m}) 为质量流率(kg·s(^{-1})),在整个管道内保持不变。若已知入口密度 (\rho_1) 与入口速度 (u_1),则任意截面上的速度可通过

[
u(x)=\frac{\rho_1 u_1}{\rho(x)}
]

直接求得。该结果在管道设计、流体输送等工程中经常使用。

常见误区与注意事项

  • 把积分形式误当作局部方程
    积分形式描述的是整体质量守恒,只有在控制体无限小(微分化)时才等价于微分形式。

  • 忽略密度随温度/应力的变化
    对可压材料或高温环境,必须保留 (\partial \rho/\partial t) 项,否则会低估体积变化。

  • 误用不可压假设
    在高速冲击或大变形问题中,材料的体积变化显著,直接使用 (\nabla!\cdot\mathbf{v}=0) 会导致错误的应力预测。

小结

质量守恒定律的数学表达分为积分形式和微分形式,两者通过高斯散度定理相互转化。对不可压材料可简化为无散度条件,对可压材料则需完整保留密度变化项。掌握这些表达式并正确选择适用范围,是进行固体力学分析、流体力学计算以及多物理场耦合时的基本前提。通过本文的推导与示例,读者应能够在实际工程问题中自如地建立质量守恒方程,并避免常见的使用误区。