Stefan方程是描述挤压力、黏度、圆盘尺寸、间隙和接近速度之间关系的润滑理论关系式。 对于半径为 (R) 且相互接近的两个平行圆盘，在采用半间隙约定 (H(t)) 时，方程为

[ F(t)=\frac{3\pi \mu U(t)R^4}{8H(t)^3}, ]

其中，(F) 为法向力，(\mu) 为动力黏度，(U=-dH/dt) 为每个圆盘的接近速度，(H) 为圆盘之间的半间距。在教学型中试装置中，测量 (F)、(H) 和 (U)，即可让学生检验润滑理论，并根据力-位移数据估算黏度。

核心结论是间隙的强反立方关系：减小半间隙会使所需挤压力按 (H^{-3}) 增加，而增大圆盘半径会使力按 (R^4) 增加。当间隙狭窄、流动缓慢且黏性效应占主导时，该方程具有可靠性。

Stefan方程描述的内容

物理实验

实验使用两个彼此平行的圆盘，中间隔着一层薄薄的流体。圆盘相互靠近，迫使液体在逐渐变窄的间隙中沿径向向外流动。

该装置代表一种简化的挤压流动操作，与润滑、聚合物加工、涂布、成型和流变测量等领域相关。

关键变量

变量包括：

(R)：圆盘半径
(H(t))：圆盘之间的半间距
(U(t))：每个圆盘的接近速度
(F(t))：施加的法向力
(\mu)：动力黏度
(\rho)：流体密度
(r)：距圆盘中心的径向距离

流体层的总厚度为 (2H)。这一约定非常重要，因为其他参考资料通常将完整间隙定义为 (h=2H)，从而得到外观不同的数值系数。

基本假设

推导基于以下假设：

爬流，惯性可忽略
间隙狭窄，(H_0/R \ll 1)
流体为牛顿型且不可压缩
圆盘表面满足无滑移条件
圆盘平行且流动轴对称
每个圆盘上的接近速度均匀
边缘效应和流体泄漏可忽略

方程的推导方法

从径向挤压流动开始

随着圆盘相互靠近，流体无法穿过固体圆盘，因此只能沿径向向外流动。

由于间隙相对于圆盘半径很小，跨越间隙方向的速度梯度远大于径向速度梯度。这是润滑理论的核心近似。

简化动量方程

在爬流条件下，纳维-斯托克斯方程中的惯性项可以忽略。径向动量平衡方程简化为

[ 0=-\frac{\partial P}{\partial r} +\mu\frac{\partial^2 u_r}{\partial z^2}, ]

其中 (u_r) 为径向速度，(z) 为跨越流体层的坐标。

压力主要随半径变化，而径向速度则在狭窄间隙内发生显著变化。

应用圆盘边界条件

无滑移条件要求径向速度在两个圆盘表面均为零：

[ u_r=0 \quad \text{at} \quad z=-H \text{ and } z=H. ]

将简化后的动量方程沿间隙方向积分两次，可得到抛物线形径向速度分布。与该速度分布对应的体积流量必须等于相互靠近的圆盘所排开的流体量。

利用连续性方程

对于对称接近的两个圆盘，半间隙以速度 (U) 减小：

[ \frac{dH}{dt}=-U. ]

变化中的间隙以与局部圆盘面积成正比的速率排开流体。轴对称连续性方程随后将径向流动与间隙缩小速率联系起来。

将连续性方程与抛物线形速度分布相结合，可得到径向压力梯度

[ \frac{\partial P}{\partial r}

-\frac{3\mu U r}{2H^3}. ]

负号表示压力从中心向外缘逐渐降低。

积分得到压力场

取圆盘边缘压力为参考压力，即 (P(R)=0)，从 (r) 积分到 (R) 可得

[ P(r)=\frac{3\mu U}{4H^3}\left(R^2-r^2\right). ]

因此，压力在中心处最高，并向边缘呈二次曲线下降。

对压力积分求得力

总法向力等于压力在圆盘面积上的积分：

[ F=2\pi\int_0^R P(r),r,dr. ]

将压力分布代入可得

[ F

2\pi\int_0^R \frac{3\mu U}{4H^3} \left(R^2-r^2\right)r,dr, ]

因此

[ \boxed{ F=\frac{3\pi\mu U R^4}{8H^3} }. ]

这就是在上述半间隙和对称接近约定下的Stefan方程。

标度关系的重要性

间隙闭合时，力会快速增加

方程中包含以下因子

[ F\propto \frac{1}{H^3}. ]

如果 (U)、(R) 和 (\mu) 保持不变，将半间隙减半会使预测力增加八倍。

正是这种强依赖关系，使得间隙测量中的微小误差可能导致计算黏度或预测力出现较大误差。

圆盘半径的影响更为显著

力与半径的四次方成正比：

[ F\propto R^4. ]

较大的圆盘会产生明显更大的阻力，因为它排开的流体更多，且压力作用的径向距离更长。

黏度和接近速度呈线性影响

力与黏度和接近速度均成正比：

[ F\propto \mu U. ]

只要流体仍为牛顿型，且爬流假设仍然成立，黏度更高的液体或更快的压缩过程就需要更大的法向力。

在中试装置中应用方程

准备并检查装置

安装圆盘，确认其彼此平行，并测量圆盘半径。小心填充间隙，避免困入空气，因为空气的压缩性会使力的测量产生偏差。

应监测温度，因为黏度可能随温度显著变化。

控制运动

设定已知的接近速度曲线 (U(t))，或施加受控的位移历史 (H(t))。控制器应记录圆盘的实际位置，而不能仅依赖其指令位置。

半间隙根据测得的圆盘间距计算：

[ H(t)=\frac{\text{圆盘总间距}}{2}. ]

记录力和位移

在挤压步骤中，记录以下数据：

法向力 (F(t))
圆盘位置或间隙 (H(t))
接近速度 (U(t))
流体温度
时间

力-位移曲线尤其有用，因为它能够显示间隙变窄时阻力如何升高。

估算黏度

重新整理Stefan方程可得

[ \boxed{ \mu

\frac{8F H^3}{3\pi U R^4} }. ]

学生可以根据每个测量数据点计算瞬时表观黏度，也可以将模型拟合到完整的力-时间或力-间隙数据集。

对于恒定黏度，计算得到的数值在有效操作范围内应大致保持一致。

验证预测的标度关系

可以通过逐次改变一个参数来检验模型：

绘制 (F) 与 (U) 的关系图，以检验速度的线性依赖。
使用不同黏度的流体，绘制 (F) 与 (\mu) 的关系图。
绘制 (F) 与 (H^{-3}) 的关系图，以检验间隙依赖。
比较不同半径的圆盘，以检验 (R^4) 依赖关系。

双对数图可以帮助更容易地在实验中识别预测的幂律指数。

理解各种权衡因素

狭窄间隙近似可能失效

模型要求

[ \frac{H_0}{R}\ll 1. ]

如果相对于圆盘半径而言间隙不够小，径向效应和边缘效应就会变得重要，一维润滑近似的准确性也会降低。

高速测试中可能受到惯性影响

爬流条件要求

[ Re=\frac{H_0U_0}{\nu}\ll 1. ]

在接近速度较高或流体黏度较低时，惯性力可能变得不可忽略。施加的挤压力应满足以下所述的实际条件

[ F\ll \frac{\mu^2}{\rho} \left(\frac{R}{H_0}\right)^4, ]

这有助于在中试装置实验中减小惯性干扰。

边缘泄漏会降低测得的力

理想推导假设圆盘边界清晰明确，且圆盘边缘的压力等于参考压力。如果流体不均匀地逸出，压力场和测得的力都会低于预测值。

平行度不佳也可能造成同类误差，因为它会形成不均匀的间隙。

牛顿型行为可能不适用

Stefan方程采用恒定动力黏度。聚合物熔体、高浓度悬浮液、凝胶以及许多结构化流体可能表现出剪切变稀、剪切增稠、黏弹性或屈服应力特性。

对于这些材料，测得的力可能代表表观黏度或与速率相关的黏度，而不是单一的材料常数。

间隙测量尤其敏感

由于力随 (H^{-3}) 变化，微小的归零误差或圆盘间距不确定性都可能导致较大的黏度误差。

在采集数据前，应对装置的位移、力、圆盘平行度和初始间隙进行校准。

瞬态效应可能很重要

该方程描述的是准稳态假设下的黏性挤压响应。在运动循环的开始和结束阶段，执行器加速度、仪器响应时间、困入空气和流体惯性都可能造成瞬态偏差。

应在指令运动与测量运动一致的区域内解释数据。

如何将其应用于您的项目

最可靠的教学实验应结合受控运动、准确的间隙测量、温度控制以及在解释力数据前进行的无量纲检查。

如果您的主要目标是验证润滑理论：使用低黏度或性质明确的牛顿型流体、较小的间隙-半径比，并设置多种接近速度，以检验 (F\propto U) 和 (F\propto H^{-3})。
如果您的主要目标是测量黏度：重新整理方程计算 (\mu=8FH^3/(3\pi U R^4))，然后将结果与独立黏度计的测量值进行比较。
如果您的主要目标是研究尺度效应：系统地改变圆盘半径或初始间隙，并将测得的力与预测的 (R^4) 和 (H^{-3}) 依赖关系进行比较。
如果您的主要目标是聚合物或非牛顿加工：将结果视为表观黏度测量，并研究推算出的黏度如何随接近速度和间隙变化。
如果您的主要目标是获得具有说服力的数据：在拟合Stefan模型前，确认 (Re\ll1)、(H_0/R\ll1)、圆盘平行度、流体密闭性和温度稳定性。

只要几何结构和操作条件处于其假设适用范围内，Stefan方程就能在测得的挤压力与狭窄间隙流动中的黏性物理之间建立直接联系。

总结表：

主题	关键信息
方程	F = 3πμUR⁴/(8H³)
主要变量	F：力；μ：黏度；U：接近速度；R：圆盘半径；H：半间隙
核心标度关系	力随μ、U和R⁴增加，并随H⁻³减小时急剧上升
主要假设	爬流、狭窄间隙、牛顿型流体、无滑移、平行圆盘、泄漏可忽略
实验测量值	法向力、圆盘间距、接近速度、温度和时间
黏度计算	μ = 8FH³/(3πUR⁴)
验证测试	改变速度、黏度、间隙或圆盘半径，并比较测得的幂律趋势
常见误差来源	困入空气、平行度不佳、间隙测量误差、边缘泄漏、惯性和温度变化
教学应用	润滑理论、流变学、聚合物加工、涂布、成型和流体力学培训