对于通道内两层不相溶液体的稳态、充分发展流动，控制方程可简化为每个相各一个动量方程，并通过壁面条件和液-液界面条件相互耦合。 对于不可压缩牛顿流体，质量守恒要求速度场无散，而动量守恒则由纳维-斯托克斯方程描述。在常用的平行板近似下，每个相中的轴向速度沿通道横向呈二次分布，积分常数由壁面处的无滑移条件以及界面处的速度连续和剪切应力连续条件确定。

速度分布由每个相中的压力梯度驱动动量平衡决定。黏度差异通过各相特定的动量方程影响速度分布，而界面条件则确保液体之间的运动和应力传递符合物理规律。

必须求解的内容

各相中的连续性

对于密度恒定的不可压缩相，连续性方程为

[ \nabla \cdot \mathbf{v}^{(i)} = 0, \qquad i=1,2 ]

其中，(\mathbf{v}^{(i)}) 是相 (i) 的速度场。

对于稳态、充分发展、单向通道流动，

[ \mathbf{v}^{(i)} = \left(v_x^{(i)}(y),0,0\right) ]

由于轴向速度仅沿横向变化，连续性方程会自动得到满足。

动量守恒

每个相的一般不可压缩牛顿流体动量方程为

[ \rho_i \frac{D\mathbf{v}^{(i)}}{Dt}

\rho_i\mathbf{g} -\nabla P^{(i)} +\mu_i\nabla^2\mathbf{v}^{(i)} ]

其中，(\rho_i)、(\mu_i) 和 (P^{(i)}) 分别为相 (i) 的密度、动力黏度和压力。

该方程包括瞬态加速度和对流加速度、重力、压力作用力以及黏性动量扩散。

充分发展的轴向流动

当流动为稳态、充分发展且主要由轴向压力梯度驱动时，轴向动量方程变为

[ 0

-\frac{dP^{(i)}}{dx} + \mu_i\frac{d^2v_x^{(i)}}{dy^2} ]

或者等价地，对于每个液层，

[ \frac{d^2v_x^{(i)}}{dy^2}

\frac{1}{\mu_i}\frac{dP^{(i)}}{dx} ]

所得速度分布为

[ v_x^{(i)}(y)

\frac{1}{\mu_i}\frac{dP^{(i)}}{dx}\frac{y^2}{2} + a_i y+b_i ]

其中，(a_i) 和 (b_i) 是由边界条件和界面条件确定的积分常数。

固体壁面处的边界条件

无滑移条件

在静止固体壁面处，流体切向速度等于壁面速度。对于参考构型中的两层几何结构，

[ v_x^{(1)}(H_1)=0 ]

以及

[ v_x^{(2)}(-H_2)=0 ]

其中，相 1 占据界面与上壁面之间的区域，相 2 占据下壁面与界面之间的区域。

这些条件提供了确定积分常数所需四个方程中的两个。

不可穿透条件

在不可渗透的静止壁面处，法向速度为零：

[ \mathbf{n}\cdot\mathbf{v}^{(i)}=0 ]

该条件防止任一液体穿过通道壁面。在理想化的单向流动解中，这已经通过不存在壁面法向速度分量体现出来。

移动壁面或多孔壁面

如果壁面发生移动，切向条件变为

[ \mathbf{v}t=\mathbf{V}{w,t} ]

而不是零。多孔壁面或反应壁面需要单独的法向通量条件，因为流体或物种可能穿过边界。

液-液界面处的条件

界面运动学条件

对于位于 (y=0) 的静止界面，法向速度必须与界面位置相容。对于不可渗透且不发生反应的界面，

[ \mathbf{n}\cdot\mathbf{v}^{(1)}

\mathbf{n}\cdot\mathbf{v}^{(2)} ]

对于平行且充分发展的构型，两个相的法向速度均为零。

如果界面发生移动或变形，其速度必须满足相应的运动学条件。此时，界面位置也会成为待求解的未知量。

切向速度连续性

对于两个普通的、不相溶的牛顿液体，在不存在界面滑移的情况下，

[ v_x^{(1)}(0)=v_x^{(2)}(0) ]

该条件确保液层在相互接触的边界处具有相同的切向速度。

它提供了确定积分常数所需的第三个关系式。

切向应力连续性

当界面不存在由表面张力梯度产生的切向应力时，黏性剪切应力连续：

[ \mu_1\frac{dv_x^{(1)}}{dy}(0)

\mu_2\frac{dv_x^{(2)}}{dy}(0) ]

当两相黏度不同时，这一点十分重要。速度梯度通常会在界面两侧发生变化，从而保持剪切应力平衡。

它为两条二次速度分布中的四个积分常数提供了第四个关系式。

法向应力和表面张力

完整的两相模型还需要在界面处满足法向应力条件。其一般形式为

[ \mathbf{n}\cdot \left(\boldsymbol{\tau}^{(1)}-\boldsymbol{\tau}^{(2)}\right) \cdot\mathbf{n}

\gamma\kappa ]

其中采用适当的符号约定，(\boldsymbol{\tau}) 为总应力张量，(\gamma) 为界面张力，(\kappa) 为界面曲率。

对于表面张力恒定的平坦静止界面，曲率为零，法向应力条件可简化为法向应力相等。当界面发生变形、通道不是平面结构或毛细力控制相分布时，该条件尤为重要。

马兰戈尼效应

如果界面张力随温度或组成变化，切向应力就不能简单地视为连续。其跃变由表面张力梯度平衡：

[ \mathbf{t}\cdot \left(\boldsymbol{\tau}^{(1)}-\boldsymbol{\tau}^{(2)}\right) \cdot\mathbf{n}

\nabla_s\gamma\cdot\mathbf{t} ]

其中，(\nabla_s) 表示表面梯度。

对于等温且界面处组成均匀的液体，通常忽略该项。

通道入口和出口处的条件

指定入口流动

中试装置模型通常会指定入口速度分布、各相流量或总容积流量。充分发展的入口流动可以表示为

[ v_x^{(i)}(y)=v_{\text{in}}^{(i)}(y) ]

也可以通过各相流量表示：

[ Q_i=\int_{A_i}v_x^{(i)},dA ]

入口设定必须与相分布和界面位置保持一致。

指定压力梯度

对于较长的充分发展通道，压力梯度可能是控制输入：

[ \frac{dP^{(i)}}{dx} ]

在许多并流通道构型中，两相共用相同的轴向压力梯度，但由于黏度和液层厚度不同，最终形成的速度分布也不同。

出口条件

在出口处，模型可以指定压力、充分发展流动条件或适当的对流流出条件。具体选择取决于中试装置是作为压力驱动系统、流量控制系统，还是与下游设备相连的系统进行建模。

出口条件必须避免对正在发展中的界面或速度场施加人为约束。

黏度差异为何重要

速度分布曲率不同

每个速度分布的曲率与下式成正比：

[ \frac{1}{\mu_i}\frac{dP^{(i)}}{dx} ]

在相同压力梯度下，黏度较低的相会比黏度较高的相产生更大的速度曲率。

应力平衡控制梯度变化

尽管切向速度在界面处连续，但速度梯度不一定连续。应力条件要求

[ \mu_1\left.\frac{dv_x^{(1)}}{dy}\right|_{0}

\mu_2\left.\frac{dv_x^{(2)}}{dy}\right|_{0} ]

因此，黏度比会改变动量在两层液体之间的分配方式。

对单元操作的影响

所得速度场会影响停留时间分布、界面更新、压降、相含率、传质面积和分离性能。

这些影响对于液-液萃取模块、层流反应器、膜接触设备以及其他通道式中试设备都十分重要。

理解各种权衡

充分发展假设

只有在入口效应变得可以忽略之后，二次分布解才有效。在入口、分配器、弯头、收缩段和相分流设备附近，轴向发展和二次流动可能要求使用完整的多维方程。

平界面假设

将界面固定在 (y=0) 排除了界面变形和毛细演化。当流量较高、通道尺寸较小、界面张力较低或两相性质产生显著压力差时，这种假设可能不再适用。

层流假设

简化后的常微分方程描述的是黏性充分发展流动，最直接适用于层流条件。如果通道进入过渡流或湍流状态，可能需要额外的封闭关系和湍流模型。

雷诺数应基于适当的水力特征尺寸，例如对于相关通道几何结构，可采用水力平均直径

[ d_m=\frac{4A}{p} ]

计算。

界面无滑移假设

速度连续性假设界面滑移可以忽略。在某些特殊系统中，表面活性剂、污染物、界面流变性或微尺度效应可能改变界面条件，因此需要使用界面本构模型。

物性恒定假设

参考方程假设黏度和密度恒定。强烈的温度梯度或组成梯度可能使 (\mu_i) 和 (\rho_i) 发生变化，从而将动量方程与能量平衡和物种平衡耦合起来。

根据目标做出正确选择

在描述两相通道式中试装置时，请采用以下模型选择原则：

如果主要关注充分发展压力驱动流动：求解各相特定的简化动量方程，在壁面处施加无滑移条件，并在界面处施加切向速度连续和切向应力连续条件。
如果主要关注界面位置或相含率：增加界面运动学条件和法向应力平衡，并纳入表面张力和曲率。
如果主要关注入口、分配器或转弯区域的行为：保留完整的瞬态或稳态多维纳维-斯托克斯方程，并施加符合物理规律的入口、出口、壁面和界面条件。
如果主要关注传质或反应：将流动方程与物种守恒方程耦合，并在入口、壁面和液-液界面处施加浓度或通量条件。
如果主要关注放大和压降预测：纳入适当的水力直径、雷诺数流动状态、壁面粗糙度、各相流量以及物性变化。

可靠的两相中试装置模型不仅由微分方程定义，还取决于能够准确表示实际设备运行方式的边界条件和界面条件。

总结表：

模型要素	控制要求	目的
连续性	`∇ · v⁽ⁱ⁾ = 0`	确保每个相中的不可压缩质量守恒
动量	`ρᵢ Dv⁽ⁱ⁾/Dt = ρᵢg − ∇P⁽ⁱ⁾ + μᵢ∇²v⁽ⁱ⁾`	描述压力、重力、惯性和黏性效应
充分发展流动	`0 = −dP⁽ⁱ⁾/dx + μᵢ d²vₓ⁽ⁱ⁾/dy²`	在每个液层中产生二次速度分布
固体壁面	无滑移和不可穿透条件	定义流体在静止、移动或多孔壁面处的运动
液-液界面	法向速度相容性和切向速度连续性	维持物理上相容的相边界
剪切应力	`μ₁ dvₓ⁽¹⁾/dy = μ₂ dvₓ⁽²⁾/dy`	平衡跨界面的动量传递
法向应力	应力跃变由 `γκ` 平衡	在界面变形时考虑曲率和表面张力
马兰戈尼效应	切向应力跃变由 `∇ₛγ` 平衡	表示由温度或组成引起的表面张力梯度
入口和出口	指定流量、速度、压力或充分发展流动条件	使模型与实际中试装置运行相匹配
扩展模型	能量、物种、湍流和变物性方程	支持传质、反应、放大以及非理想运行条件的模拟

使用 LABPARK 设计可靠的两相通道式中试装置。我们的教育与职业单元操作中试装置支持化学工程、生物过程与生物技术以及环境与水处理应用。我们帮助大学、研究机构和企业在真实运行条件下评估压降、相含率、停留时间、传质、分离性能和放大效果。联系 LABPARK，讨论您的目标流体、流动状态、通道构型以及培训或研究需求。我们的团队可以帮助您选择实用的中试装置系统，提高实验重复性，支持实践教学，并为工业开发生成有价值的工艺数据。