研究了应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)方法的直角网格技术.通过对边界三角形表面微元进行统计,建立直角网格缺陷单元与表面微元之间的映射关系,并进一步将缺陷单元区分为一级缺陷单元以及二级缺陷单元,同...研究了应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)方法的直角网格技术.通过对边界三角形表面微元进行统计,建立直角网格缺陷单元与表面微元之间的映射关系,并进一步将缺陷单元区分为一级缺陷单元以及二级缺陷单元,同时结合可变时间步长技术,从根本上解决了在DSMC方法中利用直角网格技术贴体光滑地匹配复杂几何型面的问题.给出了相应的自适应网格调整方法.基于以上策略,编制了通用DSMC计算程序RGD-Tool(Rarefied Gas Dynamic Tool).利用RGD-Tool程序对典型的圆球绕流以及无限薄斜挡板绕流两个算例进行了验证计算,结果表明所讨论直角网格技术是有效的.展开更多
NS-DSMC(Navier Stokes-Direct Simulation Monte Carlo)耦合方法是计算连续-稀薄跨流域流动的主要方法,应用过程中如何确定连续流域和稀薄流域的界面是此方法的关键问题之一,界面位置通常通过连续失效参数来判定.为合理选择连续失效参...NS-DSMC(Navier Stokes-Direct Simulation Monte Carlo)耦合方法是计算连续-稀薄跨流域流动的主要方法,应用过程中如何确定连续流域和稀薄流域的界面是此方法的关键问题之一,界面位置通常通过连续失效参数来判定.为合理选择连续失效参数,对目前广泛使用的两种连续失效参数KnQ和B参数进行了理论上的分析和比较,表明虽然他们建立的出发点不同,但在数学形式上具有一定相似性.通过圆柱绕流问题的数值试验进一步验证了两种连续失效参数具有一定对应关系.KnQ和B两种连续失效参数在NS-DSMC耦合方法中应用效果相差不大,但KnQ的阈值0.05适用范围较广,而B参数的阈值对不同流动问题会有变化.展开更多
研究了Voronoi网格技术并将其应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)计算.基于Dirichlet镶嵌与Voronoi图理论,Voronoi网格利用特征点表征网格单元,具有建立粒子与网格单元之间映射关系的独特算法,适合于DSMC方法...研究了Voronoi网格技术并将其应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)计算.基于Dirichlet镶嵌与Voronoi图理论,Voronoi网格利用特征点表征网格单元,具有建立粒子与网格单元之间映射关系的独特算法,适合于DSMC方法的统计特点.在剔除过于靠近边界的特征点以及必要情况下边界细化的基础上,通过区分由边界节点表征的非完整Voronoi网格单元以及由计算区域内镶嵌点表征的完整Voronoi网格单元,解决了Voronoi网格的二维边界匹配问题.Voronoi网格技术支持自适应DSMC计算.映射效率对比表明,Voronoi网格的DSMC计算效率高于三角形网格,低于多级直角网格.通过MEMS微喷管流动数值模拟,验证了Voronoi网格技术在DSMC方法中的有效性.展开更多
研究了适合在三维直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法计算中根据入口边界条件合理地布置进入计算区域内粒子参数的边界采样技术.提出了一种根据穿越入口边界、出口边界以及撞击壁面边界的粒子统计数据获取边界...研究了适合在三维直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法计算中根据入口边界条件合理地布置进入计算区域内粒子参数的边界采样技术.提出了一种根据穿越入口边界、出口边界以及撞击壁面边界的粒子统计数据获取边界位置处宏观流场参数的新型边界统计技术,弥补了DSMC计算中传统统计技术仅能够获取流场内部参数,不能够获取界面处流场参数的不足.这种界面统计技术具有普遍的适用性,不仅限于DSMC方法.以无穷远来流圆球绕流为例进行计算,验证了DSMC方法的边界采样与统计技术的有效性.展开更多
分别从分子运动论及连续流理论出发,对体积力驱动的微尺度平面泊肃叶(Poiseuille)流的横向分布特征进行了分析.分子水平模拟采用直接模拟蒙特卡罗(direct simulation Monte Carlo,DSMC)方法;连续流理论则主要考察了伯内特(Burnett)及超...分别从分子运动论及连续流理论出发,对体积力驱动的微尺度平面泊肃叶(Poiseuille)流的横向分布特征进行了分析.分子水平模拟采用直接模拟蒙特卡罗(direct simulation Monte Carlo,DSMC)方法;连续流理论则主要考察了伯内特(Burnett)及超伯内特(Super--Burnett)等高阶连续模型,在平行流假设下,获得一组高阶非线性常微分方程,补充完整的边界条件,并应用龙格--库塔(Runge--Kutta)方法求解.结果表明,即使对于过渡领域流动,高阶连续模型可以给出与DSMC结果完全相符的压力分布,而速度分布当努森(Knudsen)数约为0.2时即在壁面开始出现偏差;对于温度的横向分布,伯内特模型回复到纳维--斯托克斯(Navier--Stokes)水平,不能得到与DSMC一致的双峰结构,而超伯内特模型在滑移流动领域与DSMC定性相符,在过渡领域却仅能正确预测主流区温度分布,壁面附近差异明显;横向热流与纳维--斯托克斯模型预测接近,但机理上存在本质区别.本文结果提示选用连续模型时,不仅要根据流动参数来判断,还可以根据所关注的物理量来进行调整,适度扩大连续模型的适用范围.但即使采用高阶本构关系,连续模型仍然不能完全描述壁面附近区域的非平衡效应(如努森层效应),这是试图扩大连续模型适用范围时必然会遇到的困难.展开更多
文摘研究了应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)方法的直角网格技术.通过对边界三角形表面微元进行统计,建立直角网格缺陷单元与表面微元之间的映射关系,并进一步将缺陷单元区分为一级缺陷单元以及二级缺陷单元,同时结合可变时间步长技术,从根本上解决了在DSMC方法中利用直角网格技术贴体光滑地匹配复杂几何型面的问题.给出了相应的自适应网格调整方法.基于以上策略,编制了通用DSMC计算程序RGD-Tool(Rarefied Gas Dynamic Tool).利用RGD-Tool程序对典型的圆球绕流以及无限薄斜挡板绕流两个算例进行了验证计算,结果表明所讨论直角网格技术是有效的.
文摘NS-DSMC(Navier Stokes-Direct Simulation Monte Carlo)耦合方法是计算连续-稀薄跨流域流动的主要方法,应用过程中如何确定连续流域和稀薄流域的界面是此方法的关键问题之一,界面位置通常通过连续失效参数来判定.为合理选择连续失效参数,对目前广泛使用的两种连续失效参数KnQ和B参数进行了理论上的分析和比较,表明虽然他们建立的出发点不同,但在数学形式上具有一定相似性.通过圆柱绕流问题的数值试验进一步验证了两种连续失效参数具有一定对应关系.KnQ和B两种连续失效参数在NS-DSMC耦合方法中应用效果相差不大,但KnQ的阈值0.05适用范围较广,而B参数的阈值对不同流动问题会有变化.
文摘研究了Voronoi网格技术并将其应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)计算.基于Dirichlet镶嵌与Voronoi图理论,Voronoi网格利用特征点表征网格单元,具有建立粒子与网格单元之间映射关系的独特算法,适合于DSMC方法的统计特点.在剔除过于靠近边界的特征点以及必要情况下边界细化的基础上,通过区分由边界节点表征的非完整Voronoi网格单元以及由计算区域内镶嵌点表征的完整Voronoi网格单元,解决了Voronoi网格的二维边界匹配问题.Voronoi网格技术支持自适应DSMC计算.映射效率对比表明,Voronoi网格的DSMC计算效率高于三角形网格,低于多级直角网格.通过MEMS微喷管流动数值模拟,验证了Voronoi网格技术在DSMC方法中的有效性.
文摘研究了适合在三维直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法计算中根据入口边界条件合理地布置进入计算区域内粒子参数的边界采样技术.提出了一种根据穿越入口边界、出口边界以及撞击壁面边界的粒子统计数据获取边界位置处宏观流场参数的新型边界统计技术,弥补了DSMC计算中传统统计技术仅能够获取流场内部参数,不能够获取界面处流场参数的不足.这种界面统计技术具有普遍的适用性,不仅限于DSMC方法.以无穷远来流圆球绕流为例进行计算,验证了DSMC方法的边界采样与统计技术的有效性.
基金国家自然科学基金委员会与中国工程物理研究院联合基金资助项目(10976029)The project was supported by the National Natural Science Foundation of China-NSAF(10976029)
文摘分别从分子运动论及连续流理论出发,对体积力驱动的微尺度平面泊肃叶(Poiseuille)流的横向分布特征进行了分析.分子水平模拟采用直接模拟蒙特卡罗(direct simulation Monte Carlo,DSMC)方法;连续流理论则主要考察了伯内特(Burnett)及超伯内特(Super--Burnett)等高阶连续模型,在平行流假设下,获得一组高阶非线性常微分方程,补充完整的边界条件,并应用龙格--库塔(Runge--Kutta)方法求解.结果表明,即使对于过渡领域流动,高阶连续模型可以给出与DSMC结果完全相符的压力分布,而速度分布当努森(Knudsen)数约为0.2时即在壁面开始出现偏差;对于温度的横向分布,伯内特模型回复到纳维--斯托克斯(Navier--Stokes)水平,不能得到与DSMC一致的双峰结构,而超伯内特模型在滑移流动领域与DSMC定性相符,在过渡领域却仅能正确预测主流区温度分布,壁面附近差异明显;横向热流与纳维--斯托克斯模型预测接近,但机理上存在本质区别.本文结果提示选用连续模型时,不仅要根据流动参数来判断,还可以根据所关注的物理量来进行调整,适度扩大连续模型的适用范围.但即使采用高阶本构关系,连续模型仍然不能完全描述壁面附近区域的非平衡效应(如努森层效应),这是试图扩大连续模型适用范围时必然会遇到的困难.
