基于非对称NACA 2412翼型混合促进器强化真空膜蒸馏性能的研究

《Membranes》：Enhancing Vacuum Membrane Distillation Performance via Asymmetric NACA 2412 Airfoil-Based Mixing Promoters

　　

尽管用于低温海水淡化的真空膜蒸馏（VMD）是一项创新技术，但其工业应用仍受到浓度极化和温度极化的限制，而传统被动式间隔网在缓解这些问题的同时伴随着高水力能耗。本文旨在利用非对称NACA 2412弧形翼型作为被动混合器来削弱边界层限制。研究人员在向上和向下两种取向下评估了翼型攻角（AoA）的影响，气流模拟的雷诺数（Re）为1000和1500。结果表明，向上取向的翼型在引导气流方面非常有效，可在边界层内诱发流动分离和相干尾迹结构，同时保持较低的形状阻力及相应较小的水力学代价。在Re = 1500、AoA = 12°的向上设计中，渗透通量较零攻角基线（AoA = 0°）提高了14.6%，达到66.03 kg/m2h。这一增益的代价较为温和，通道压降仅增加约10 Pa（从13.33 Pa升至23.52 Pa）。这些结果表明，基于翼型的混合促进器在通量增强与水力学代价之间实现了有利平衡，使其成为VMD通道中一种有前景的被动设计策略。

论文解读：基于非对称NACA 2412翼型混合促进器强化真空膜蒸馏性能的数值研究

一、研究背景与问题

全球淡水资源短缺推动了对更高效、可持续的海水淡化与水净化技术的探索。真空膜蒸馏（VMD）作为膜蒸馏（MD）家族中极具前景的新技术，通过在渗透侧施加真空形成跨膜蒸气压差，可在远低于传统热法蒸馏的进料温度下运行，具有对非挥发性溶质近乎100%的截留率、较低显热损失以及可利用低品位工业余热或太阳能热能的优势。当前全球海水淡化产能约为每天9537万立方米，其中反渗透（RO）占68.7%，多级闪蒸（MSF）和多效蒸馏（MED）分别占17.6%和6.9%。与依赖电力驱动的RO不同，MD类技术可利用余热实现热电联产，在中东和北非（MENA）地区具有理想的应用条件。在DCMD、AGMD、SGMD和VMD四种构型中，VMD因渗透侧真空消除了冷凝侧阻力而具有最高的热力学驱动力。

然而，VMD的实际性能长期受制于两个相互关联的边界层现象：温度极化（TP）和浓度极化（CP）。水蒸气跨膜传输导致进料侧局部温度下降，非挥发性溶质在膜表面富集，两者均显著削弱了蒸气传输的驱动力，加剧膜污染并降低跨膜通量。传统被动混合方案多采用间隔网或静态障碍物，虽能缓解极化但代价是高水力学损失和泵送能耗。航空翼型（NACA）因其优化的升阻比和可控流动分离特性，有望在极小形状阻力代价下产生强烈局部涡量和混沌对流。然而，已有研究集中于对称翼型或平板插入件，尚无研究系统考察非对称弧形翼型在VMD通道中作为被动混合器的取向与攻角效应。为此，研究人员开展了本项二维计算流体力学（CFD）研究。

二、研究目的与主要技术方法

本研究旨在通过合理设计通道水力学条件，利用气动翼型轮廓实现VMD性能的显著提升。研究采用Ansys Fluent 2025 R2软件进行二维稳态k?ω SST湍流模型模拟。计算域为典型平板式VMD通道，通道长L = 60 mm、高H = 3 mm，内置弦长c = 5 mm的非对称NACA 2412翼型。研究在Re = 1000和1500两种湍流工况下，系统考察了翼型向上和向下两种取向在攻角（AoA）为5°、8°、12°时的性能，并以0°攻角为基线。进料以80 °C、海水盐度（35000 ppm）均匀速度入口进入，渗透侧真空压力为5 kPa。膜界面采用Knudsen-黏性扩散模型计算膜渗透系数，跨膜通量由Antoine方程修正盐度后的蒸气压差驱动。网格独立性研究表明，140000单元的M3网格为最优配置，y⁺值在所有工况下均小于1，通量数值不确定度约为±1.05%。模型经与Anqi等人（2016）的实验和CFD数据验证，在5个工况中3个偏差小于5.5%。

三、研究结果

3.1 网格独立性与模型验证

网格独立性研究显示，M3（140000单元）与M4（240000单元）网格间的通量差异小于0.5%，翼型表面y⁺≤0.205、膜表面y⁺≤0.119，充分解析了黏性底层。采用Celik等人（2008）的网格收敛指数（GCI）方法计算得到表观收敛阶p≈0.89，Richardson外推通量为56.99 kg/m²h，M3网格的离散不确定度为GCI≈1.05%。所有报告的攻角和取向驱动的通量变化均超出该不确定度带至少一个数量级，确认观察到的趋势是数值上充分分辨的而非离散化伪影。模型验证方面，与Anqi等人（2018）的空通道VMD数据对比，最小偏差2.20%（Re = 1000，T = 50 °C），最大偏差18.05%（Re = 200，T = 80 °C），5个工况中3个偏差在5.5%以内。

3.2 比较性能分析

在Re = 1500下，向上取向AoA = 12°时性能最优：渗透通量从基线的57.60 kg/m²h提升至66.03 kg/m²h（+14.6%），TPC从0.9319提升至0.9381，CPC仅从1.0017微变至1.0019，通道压降从13.33 Pa升至23.52 Pa。在Re = 1000下，向上取向12°时通量从51.74升至57.27 kg/m²h（+10.7%），TPC从0.9275升至0.9318。向下取向在Re = 1500、?12°时通量仅从57.60升至62.20 kg/m²h（+8.0%），且?5°和?8°时通量（56.80和56.45 kg/m²h）略低于基线；在Re = 1000下，?5°和?8°时通量（50.40和49.43 kg/m²h）同样低于基线（51.74 kg/m²h）。所有构型的绝对压降均保持在23.5 Pa以下，比传统间隔网和金属泡沫的数百帕至千帕级压降低2–3个数量级。研究还定义了优值数（Me）评估通量增益与压降代价的比例关系，所有构型Me < 1，反映了极低基线压降下比例式指标的固有局限，但绝对压降增量在VMD操作中微不足道。

3.3 水动力学机制

向上取向的翼型在吸力面与膜之间形成喷嘴效应，将加速流持续导向膜表面，产生连贯的分离尾迹结构，有效更新近壁流动并扰动热边界层。TPC从基线0.9319渐进提升至12°时的0.9381，对应14.6%的通量增益。向下取向因压力面朝向膜，流动被偏转远离膜面，近壁对流减弱。尽管在Re = 1500下向下取向的压降（14.46、16.86、23.38 Pa）与向上取向对应值相差约3%，但通量始终较低，且小攻角下低于基线。速度场分析显示，向上12°时膜面形成贯穿全通道的高速壁面射流（最高约0.21 m/s），而向下12°时高速带出现在底部非膜壁面，膜侧出现低速回流区。温度场云图进一步证实，向上取向使冷膜层显著变薄，12°时几乎在整个通道长度上维持膜面接近主流温度；向下取向则在翼型后方形成停滞冷流体聚集，冷膜层更厚且不均匀。

四、讨论与结论

研究表明，向上取向的NACA 2412翼型在Re = 1500、AoA = 12°时实现最大渗透通量66.03 kg/m²h，较基线提升14.6%，TPC从0.9319提升至0.9381，压降仅从13.33 Pa升至23.52 Pa。该最大压降比金属泡沫插入法在相近工况下报告的3332.5 Pa约低两个数量级。与Qamar等人（2023）采用NACA1410翼型间隔网在超滤中实现228%通量提升相比，本研究的14.6%增益虽较小，但源于完全不同的传输限制机制——VMD是热驱动过程，增益来自温度极化降低而非生物污垢控制。向下取向虽承受与向上取向几乎相同的压降代价，但通量增益较小，低攻角下甚至低于基线，证实了流动偏转方向而非翼型存在本身是决定性能的关键因素。CPC在所有构型中几乎不随攻角变化（1.0017–1.0023），表明浓度极化影响可忽略，通量增益完全由TPC改善驱动，因此该设计对高盐度进料预期仍有效，但翼型本身不能缓解浓度极化驱动的结垢，高结垢倾向进料仍需补充防垢措施。研究结论指出，向上取向的弧形NACA 2412翼型以极低压降代价实现了显著的边界层扰动，突破了以往被动混合技术在能量消耗与产水增益之间的权衡限制。未来工作需开展三维CFD模块级分析、实验验证（包括PIV流场测量、压降和通量直接测量）、扩展雷诺数范围、以及长期运行下污染、结垢和润湿行为的评估。