$v_s = \dfrac{(\rho_p - \rho_f)\,d_p^2\,g\,C_c}{18\mu}$
Cunningham修正
$C_c = 1 + \dfrac{2\lambda}{d_p}\!\left(1.257 + 0.4\,e^{-1.1d_p/2\lambda}\right)$
扩散系数
$D_p = \dfrac{k_B T\,C_c}{3\pi\mu\,d_p}$
针对粒径0.01~100μm的粒子,实时计算并可视化沉降速度、扩散系数、斯托克斯数、Cunningham修正系数及肺部沉积率。
重力沉降速度(斯托克斯定律):这是计算粒子在静止空气中最终匀速下落速度的核心公式。当粒子密度大于空气时,它会沉降。
$$v_s = \dfrac{(\rho_p - \rho_f)\,d_p^2\,g\,C_c}{18\mu}$$其中,$v_s$是沉降速度,$\rho_p$是粒子密度,$\rho_f$是空气密度,$d_p$是粒子直径,$g$是重力加速度,$C_c$是Cunningham修正系数,$\mu$是空气的动力粘度。注意这里有$d_p^2$,所以粒径增大一倍,沉降速度会变为四倍!
Cunningham滑流修正系数:当粒子小到与空气分子的平均自由程(约65纳米)相当时,空气不再像连续流体,粒子会“打滑”,使得阻力变小,沉降变快。这个系数就是修正这个效应的。
$$C_c = 1 + \dfrac{2\lambda}{d_p}\!\left(1.257 + 0.4\,e^{-1.1d_p/2\lambda}\right)$$$\lambda$是空气分子的平均自由程。对于0.1微米(100纳米)左右的粒子,$C_c$可能大于2,意味着实际沉降速度是经典斯托克斯定律预测的两倍以上!
扩散系数:描述微小粒子因受空气分子随机碰撞(布朗运动)而扩散的剧烈程度。粒径越小,温度越高,扩散越强。
$$D_p = \dfrac{k_B T C_c}{3 \pi \mu d_p}$$$D_p$是扩散系数,$k_B$是玻尔兹曼常数,$T$是绝对温度。这个公式揭示了为什么捕捉纳米级粒子非常依赖扩散机制。
斯托克斯数:衡量粒子惯性大小的无量纲数,是预测粒子是否会偏离流线发生惯性碰撞的关键参数。
$$Stk = \frac{\rho_p d_p^2 U C_c}{18 \mu D}$$$U$是特征流速(如气流速度),$D$是特征尺寸(如纤维直径或气管直径)。$Stk$越大,粒子惯性越强,越容易撞击障碍物。
空气净化与口罩设计:高效滤网(HEPA)正是利用不同粒径粒子的不同动力学行为。对大颗粒(>1μm)主要靠筛分和惯性撞击拦截;对最难捕捉的0.1-0.3微米粒子(穿透窗口),则依赖布朗扩散使其吸附到超细纤维上。模拟器中的扩散系数和斯托克斯数就是设计的核心依据。
吸入式药物递送:治疗哮喘或肺部疾病的药物需要精准沉积在气管、支气管或肺泡。通过控制药物粒子的空气动力学直径(通常在1-5微米),利用其沉降和惯性撞击原理,可以使其在目标区域沉积。工程现场常见的是使用干粉吸入器。
工业粉尘控制与职业健康:在矿山、铸造车间,了解不同大小粉尘的沉降速度和扩散能力,是设计通风除尘系统的前提。例如,呼吸性粉尘(可进入肺泡,通常<10μm)需要高效的湿式除尘或袋式除尘器来捕获。
环境科学与流行病学:研究PM2.5、病毒气溶胶(如新冠病毒)在大气中的传输与沉降规律,评估其传播距离和风险。例如,大飞沫(>100μm)快速沉降,而核膜飞沫(<5μm)可悬浮很久并远距离扩散,这直接影响了社交距离和通风指南的制定。
首先,“粒子密度总可以默认为与水相同的1g/cm³”这种想法。这其实相当危险。例如,金属氧化物粉尘(密度约4g/cm³)与花粉(密度约0.5g/cm³)即使粒径同为10μm,其沉降速度也可相差数倍。由于它们在肺部的沉积行为也不同,因此第一步应是查明实际粒子的真实密度。
其次,容易认为“空气状态固定为标准条件(20℃, 1个大气压)”。在高海拔(气压降低)或高温环境中,空气的粘度μ和密度ρ_f会发生变化,从而影响计算结果。例如,在海拔2000米处空气密度约为标准条件的80%,浮力减小,沉降速度会略微加快。模拟实际环境时,切记调整温度与压力参数。
最后,“斯托克斯数Stk大于1就意味着很快沉降到地面”这种简单化解读。Stk是衡量粒子在气流“拐弯处”行为的指标。即使直线运动性强,若初始速度为零,沉降过程依然缓慢。例如,被通风口气流携带的大粒子(Stk>1)容易在管道弯头处撞击管壁,而同样的粒子在静止空气中则仅进行斯托克斯沉降。请养成习惯,将工具计算出的各参数与其所应用的流场类型结合起来思考。
某污染源排放PM2.5粒径1.2μm的气溶胶,环境温度20°C,大气压力101.3 kPa,气流速度0.5 m/s。输入参数后计算得:沉降速度vs=3.8×10⁻⁶ m/s,扩散系数D=1.95×10⁻⁵ cm²/s,斯托克斯数Stk=2.1×10⁻⁴(远小于1,粒子随气流运动),肺部细支气管沉积率约32%。若粒径增大至5μm,沉降速度跃升至7.2×10⁻⁵ m/s,肺部沉积率降至8%。