- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
苎些竺墨呈竺竺鱼三垒!墨!!-计算物理
篓!墨墨:兰 !里竺塑!苎些竺墨呈竺竺鱼三垒!墨!!旦型!! :兰』鳖
[文章编号]1001—246x(2004)04—0316—05
DSMC方法的压力边界条件实现
王沫然, 王金库, 李志信
(清华大学工程力学系,北京100084)
[摘 要] 提出了一种实现DsMc压力边界条件的新方法,新方法不仅可以避免传统“通量法”造成的计算发
散问题,又较之“人口平均法”有更快的收敛速度.使用这种方法,对不同压差驱动下等壁温微通道内的气体
流动进行了模拟,并与传统的基于连续介质假设的滑移理论所得结果进行了比较.
[关键词] 压力边界条件;微流动;DsMc
[中图分类号] 041;TKl6;V211 [文献标识码] A
O 引言
近年来,微电子机械系统(MEMs)得到了特别的重视和发展,并带动了许多相关学科从传统的常规尺度
向微细观尺度发展¨’2J.MEMs发展对微尺度流动和换热机理研究的迫切需求推动了微尺度流动和微尺度换
热的发展,并使之成为了近期的研究热点¨’4-.对于大多数微流动系统,由于系统特征尺度(f)的减小,使得流
体(特别是气体)的平均自由程(A)与之相比不可忽略.根据努森数的定义‰=A/z,流动可以分为4个区:
当砌0.01时,连续介质假定不再成立,尽管滑移区的流动仍然可以通过壁面滑移修正使得连续介质运动
方程得到近似解,然而在过渡区和自由分子流区,基于连续介质的模拟方法会带来巨大的偏差.Birdbo提出
的直接模拟Monte
超高速稀薄气体的研究中得到了验证,近年来,一些学者也成功的将DsMc应用在MEMs的高速流动中¨’61.
在常规尺度的模拟中,DSMC主要用来求解一些高速或超高速的流动,边界条件通常采用无穷远处来流
的速度边界条件,由于在超音速流动中,下游的流动状况对上游没有影响,因此在出口边界上也可以采用均
匀速度的简单处理方法进行求解.这种边界的选取已经被绝大多数的DSMc研究者所采用.另一种成功的做
法是在出流边界上赋以“真空”边界条件,这种处理方法所得的结果也得到了实验的验证’8].然而,在MEMs
中所遇到的流动通常都是通道内亚音速的流动,而且实验所能测得的宏观量通常是压力和温度等,因此,采
了一种压力边界的处理方法,在入口边界上采用“通量法”来计算粒子的分布速度,这种方法的缺点是程序处
理上比较复杂,且容易发散;kou和Fan一11。发展了这种方法,采用在人口处统计平均的方法分布粒子,程序
处理比较简单,但收敛速度较慢.本文在此基础上,提出一种新的压力边界处理方法,新方法避免了发散的产
生,而且计算收敛性有所提高.本文还使用新方法对微通道内的压力驱动流进行了模拟,并在滑移区内与已
有的解析解进行了对比.
1 压力边界
在传统的DsMc模拟中,一般情况下已知条件为无穷远来流温度r。,来流数密度n。和来流速度y。,
其中在微尺度条件下,n。和y。很难通过实验给出,因此本文假设当前的边界条件为人口压力P。出口压力
P。和入口温度E。.
[收稿日期]2003一03—28;[修回日期】2003一lo—17
[基金项目]国家自然科学基金(批号2)及国家重点基础研究发展规划(批号:G1999033106)资助项目
[作者简介]王沫然(1977一),男,吉林省吉林市,博士,主要从事MEMs及微尺度流动换热模拟.
第4期 王沫然等:DsMc方法的压力边界条件实现 317
1.1入口压力边界
由于在入口处的温度和压力已知,因此只需要确定从人口界面进入计算区域的粒子速度.Ikegawa等凹1
和Nance等…3最早提出使用入口处通量来确定人口处分子的分布速度
㈨^=警, (1)
其中下标m表示单元格位置,A为单元格边界处面积,n为气体数密度,n为通过边界处在某一方向上的粒
子数通量,其具体可以表示为
~Z,
n:2生堕旦!=£!±当!!型型,
您可能关注的文档
最近下载
- 2024年初级会计师《经济法基础》三色速记手册.pdf VIP
- 2024年贵州省黔东南苗族侗族自治州凯里市鸭塘镇招聘社区工作者真题及参考答案详解.docx VIP
- 《小学语文课程与教学》第2章小学语文课程资源.ppt VIP
- 公众聚集场所防火措施之加强消防安全管理模板范本.docx VIP
- SPC控制图的系数和公式表.xls VIP
- 安徽省舒城中学2024-2025学年高一上学期入学检测英语试题(无答案).docx VIP
- 关于冯唐的经典语录冯唐经典语录100句.pdf VIP
- 学校、幼儿园校舍能力提升改造建设项目监理服务方案投标文件(技术方案).doc
- 三维柔性焊接工作台系统焊接变位机焊接自动化.PDF VIP
- 国际法院罗马规约.doc VIP
文档评论(0)