微化工机械DSMC-SPH多尺度耦合算法研究

微化工机械DSMC-SPH多尺度耦合算法研究

论文摘要

目前,基于MEMS/NEMS(Micro/Nano Electro Mechanical Systems,简称MEMS/NEMS)技术的微器件在化工领域应用愈来愈广泛。随着对其功能和性能等要求的提高,对微化工器件的物理特性和工作机理的研究成为近年来各国学者关注的热点,是推动微化工机械发展的一个关键因素。特别是,对于微化工流体机械,如微反应器、微分离器和微泵等,由于其结构特征尺寸是从微观跨越到宏观,其内部流动呈现出多尺度效应,使得传统的单一尺度流体动力学分析方法已不再适用于这类对象。因此,针对微化工器件进行多尺度流动机理研究,并根据其结果对微化工器件的结构及参数进行优化,为新结构、新器件等开发提供基础理论依据,是当前微化工机械技术中亟待深入开展的重要研究课题。针对微化工器件内部流场的多尺度流动特征,基于微观到宏观的跨尺度耦合流体动力学分析,本文提出了一种基于直接蒙特卡洛方法(direct simulationMonte Carlo,简称DSMC)和光滑粒子动力学方法(smoothed particlehydrodynamics,简称SPH)耦合的多尺度方法——DSMC-SPH多尺度耦合算法,并应用该方法对微化工分离器件中的多尺度流动问题进行了数值模拟和方法论证。结果表明,该多尺度耦合算法除了拥有传统多尺度算法的“高效率”、“收敛性好”和“精度较高”等优点外,还可用于解决传统多尺度流动模拟中难以处理的“复杂几何边界”、“瞬态特性”、“热质传递”和“多组分”等技术瓶颈,更加适合于诠释微化工过程机械的流动特性和工作机理。本文的主要研究内容和创新成果在于:(1)为将微尺度的DSMC方法和宏观尺度的SPH方法耦合,分别改进了两种方法的边界条件处理。主要是针对DSMC方法,基于单元格压力分布,提出了新的一阶和二阶压力边界处理方法;通过新方法与解析解的结果比较,给出了方法的有效性论证。结果表明,新方法具有更好的适用性、收敛性和精度。其次,针对SPH方法,基于进出口设置边界粒子数据交换区,提出了适用于连续通道流动的边界处理方法,分别应用于Couette和Poiseuille流动的模拟计算,并与解析解进行对比,验证了新方法的有效性。经修正后,二者的基于压力耦合的边界条件处理方法,更能满足微化工流动的跨尺度耦合计算要求。(2)采用改进边界条件处理的DSMC方法,对微尺度流场中两种常见通道结构(微直通道和微射流通道)的流动特性和热质传递特性进行了模拟分析。对于微直通道流动,通过对压力场、速度场、温度场分布、质量流量和热通量等的研究,获得了变几何特征AR及壁面温度Tw等作用下的热质传递特性,总结出了微直通道流场中不同于宏观流动的热质传递规律;对于微射流通道流动,主要是研究了改变扩展区温度的情况下,其流动过程中的热质传递特性,给出了分别对壁面进行加热和冷却条件下,微射流场不同于宏观射流场的热质传递规律。研究方法和结果可以用于指导新型微射流结构反应器的设计与开发。(3)融合DSMC方法处理微流动和SPH方法处理瞬态连续流动的优点,原创性地提出了一种DSMC-SPH多尺度耦合算法。该耦合算法处理中,采用基于粒子模拟的DSMC方法求解微尺度流场的稀薄特性,同时采用基于粒子模拟的连续介质假设SPH方法求解宏观流场的连续特性,实现了由DSMC到SPH的跨尺度耦合计算;重点探讨了计算区域划分、算法界面处的耦合迭代方式以及界面信息交换过程等;构造了DSMC-SPH多尺度耦合算法的计算流程,研究了多尺度算法的收敛判断条件,编写完成了整个计算程序;结合文献中的标准算例及结果,对该方法进行了多尺度流场计算的有效性考核和方法论证,探讨了不同参数对算法收敛性及效率的影响等。结果表明,该算法及程序不仅拥有传统多尺度算法的“高效率”、“收敛性好”和“精度较高”等优点,而且更适用于解决“复杂几何边界”、“瞬态特性”和“热质传递”等实际化工流动问题。(4)采用上述DSMC-SPH多尺度耦合算法,对用于多组分气体分离的微筛以及用于燃料电池氢提纯的微分离膜结构等的多尺度流场进行模拟与分析。对于微筛结构,重点计算与分析了筛孔特征尺寸、筛孔壁面错位等流道几何结构参数变化对微筛内多尺度流场特性的影响;对于微分离膜结构,主要是计算与分析了多尺度流场中H2和CO不同的流动特性,包括各自的压力、流速以及多尺度流动效应。为了进一步探索更广义的多尺度通道流动特性,通过改变流场压力边界条件,探讨了不同压差驱动下,流场的多尺度流动规律;通过改变流场不同区域的特征尺寸,探讨了各自对多尺度通道流动的压力、速度以及质量流量等分布的影响。从而,比较全面地论证了该方法的有效性和准确性。

论文目录

  • 摘要
  • ABSTRACT
  • 目录
  • 1 绪论
  • 1.1 引言
  • 1.2 微化工技术的发展
  • 1.3 宏观-微观多尺度耦合方法国内外研究现状
  • 1.3.1 采用统一控制方程的方法
  • 1.3.2 宏观-微观分区模拟界面耦合方法
  • 1.4 本文主要研究内容
  • 2 DSMC和SPH多尺度耦合边界条件改进
  • 2.1 引言
  • 2.2 DSMC方法及其边界条件改进
  • 2.2.1 DSMC方法简介
  • 2.2.2 DSMC方法计算流程
  • 2.2.3 DSMC方法边界条件的处理方法改进
  • 2.3 SPH方法及其边界条件改进
  • 2.3.1 SPH方法的基本原理
  • 2.3.2 SPH方法的控制方程
  • 2.3.3 适用于宏观通道流动的SPH方法边界条件改进
  • 2.4 改进边界条件DSMC方法的验证
  • 2.4.1 基于滑移壁面条件的微通道流场解析解
  • 2.4.2 数值验证算例
  • 2.4.3 改进边界条件DSMC方法计算结果与解析解的比较
  • 2.4.4 改进边界条件DSMC方法计算结果与实验值的比较
  • 2.4.5 改进压力边界处理方法的收敛性讨论
  • 2.5 改进边界条件SPH方法的验证
  • 2.5.1 物理模型的描述
  • 2.5.2 SPH方法对应的随时间演化连续流场解析解
  • 2.5.3 数值算例
  • 2.5.4 改进边界条件SPH方法计算结果与解析解的比较
  • 2.6 本章小结
  • 3 改进边界条件DSMC方法应用于微流场热质传递特性研究
  • 3.1 引言
  • 3.2 微直通道流场的热质传递特性研究
  • 3.2.1 边界条件控制方程
  • 3.2.2 流场温度、质量流量及热通量的计算方法
  • 3.2.3 数值计算模型
  • 3.2.4 AR变化下的微通道热质传递特性
  • w变化条件下的微通道热质传递特性'>3.2.5 Tw变化条件下的微通道热质传递特性
  • 3.3 微射流结构流场的热质传递特性研究
  • 3.3.1 边界条件控制方程
  • 3.3.2 流场温度、质量流量及热通量计算
  • 3.3.3 数值计算模型
  • 3.3.4 应用于微射流结构流场的方法比较
  • 3.3.5 微射流结构流场的热质传递特性及规律
  • 3.4 本章小结
  • 4 DSMC-SPH多尺度耦合算法研究
  • 4.1 引言
  • 4.2 多尺度流动中不同流场区域内的数值方法
  • 4.3 计算区域划分
  • 4.4 界面交换处理
  • 4.4.1 SPH向DSMC(SPH→DSMC)的数据传递
  • 4.4.2 DSMC向SPH的(DSMC→SPH)数据传递
  • 4.5 收敛判据
  • 4.6 计算流程
  • 4.7 耦合算法有效性验证
  • 4.8 计算效率的影响因素
  • 4.8.1 耦合重叠区大小对计算效率的影响
  • 4.8.2 初始条件设置对计算效率的影响
  • 4.9 本章小结
  • 5 DSMC-SPH多尺度耦合算法的应用与讨论
  • 5.1 引言
  • 5.2 微过滤器件及其多尺度流动效应
  • 5.2.1 微筛器件
  • 5.2.2 微滤膜
  • 5.2.3 微过滤结构的多尺度流动效应
  • 5.3 可控微筛结构的多尺度流动分析
  • 5.3.1 筛孔特征尺寸改变的多尺度流动特性
  • 5.3.2 筛孔壁错向移位下的多尺度流动特征
  • 5.4 微滤膜分离过程的多尺度流动分析
  • 2和CO的多尺度流动特性'>5.4.1 H2和CO的多尺度流动特性
  • 5.4.2 不同压差驱动下的多尺度流场特性分析
  • 5.4.3 微滤膜流场结构对多尺度流动的影响分析
  • 5.5 本章小结
  • 6 总结与展望
  • 6.1 主要研究与结论
  • 6.2 主要研究创新点
  • 6.3 后续研究展望
  • 符号说明
  • 参考文献
  • 在读期间科研成果
  • 致谢
  • 相关论文文献

    标签:;  ;  ;  ;  ;  ;  ;  ;  

    微化工机械DSMC-SPH多尺度耦合算法研究
    下载Doc文档

    猜你喜欢