读书报告--DSMC方法中分子碰撞的计算

1、. 1 / 8 dsmc 方法中分子碰撞的计算摘要dsmc 方法是依赖物理的概率模拟方法。从宏观参量到细观速度分布函数的水平上，该方法均能得到实验的支持。经过30 多年的开展 ,dsmc方法在航天、真空系统等离子体材料加工、微电子刻蚀、微机电系统、化工等21 世纪高技术开展的前沿领域获得了广泛的应用。dsmc 方法通常用有限个模拟分子代替大量真实气体分子，通过随机抽样并跟踪模拟分子轨迹实现求解真实气体流动问题。其主要特点是模拟分子的迁移运动与碰撞解耦。本文主要从分子碰撞模型、 分子碰撞的计算两方面详细阐述先关研究进展。1. 引言dsmc 方法通常称试验粒子为模拟分子或模拟粒子，它的根本做法可表

2、述为：用有限个模拟分子代替大量的真实气体分子，通过随机抽样模拟分子并跟踪模拟分子的运动轨迹来到达求解真实气体流动问题的目的。dsmc 方法的主要特点在于将模拟分子的迁移运动与碰撞作用解耦处理。在每个时间步长， 首先认为每个模拟分子作匀速直线运动得到新的位置坐标，如果模拟分子与边界发生相互作用那么进展相应处理；然后计算模拟分子之间的碰撞。dsmc 方法采用几率论而不是决定论的方法计算模拟分子间的碰撞，因而能够大大减少计算时间。 整个计算是一个时间过程， 从非定常向定常开展。 最后对各模拟分子的物理量统计平均得到所需的宏观流动参量。2. 分子碰撞模型的研究进展2.1 分子碰撞模型的研究dsmc 方

3、法中，最简单的分子模型是硬球模型，它由bird 于 1976 年提出，随后，还提出了逆幂律分子模型以与maxwell 分子模型。由于真实气体流动条件下，分子的碰撞截面是随着相对速度而改变的，这样才能保证 与温度 t 的关系与实际气体中一致。 而上述任何一种分子模型都不同时具有有限的碰撞截面又能实现粘性与温度的依赖关系， 为此，bird 提出了变径硬球模型 variable hard sphere，简称 vhs ，在此模型假设下，分子具有和硬球一样的均匀散射几率，但它的碰撞截面像逆幂律分子一样是分子间相对速度的函数，这样导出的粘性系数随温度的一定幂次变化， 满足真实气体所表现出来的规律。但由于

4、vhs模型没有考虑气体分子的散射非对称性， 这导致在考虑多组元混合气体中扩散起重要作. 2 / 8 用的情况下， vhs模型会给出与实际偏离较大的结果。koura 和 matsumoto提出了变径软球模型 variable soft sphe简称 vss ，克制了 vhs模型的这一缺陷。虽然 vhs与 vss模型用逆幂率作用势取代了硬球分子， 但它用单纯的排斥力来描述分子间的作用。 对于再入流场， 温度变化很大， 粘性对温度的依赖靠一个方幂已不能完全描述。 分子间的相互作用， 除了近距离的排斥力外在较大分子间距处吸引力也显现出来。为了在dsmc 模拟中能够同时再现分子间的排斥与吸引的作用势，

5、hassan和 hash提出了概念化硬球模型ghs ，是对 vhs模型和 vss模型概念的一种推广， 包含了同时具有吸引和排斥的分子间作用力。它的散射率与硬球模型一样，但其总碰撞截面与碰撞中的相对动能反映出了吸引排斥作用势。结 合 vss 模型 的 散 射 率 与 ghs的碰 撞面 ，文献 建议 了概 念 化 软 球 模型generalized soft sphere，简称 gss 。用 gss模型计算出的粘性系数、自扩散系数和扩散系数在整个温度变化围给出与实验数据相符很好的结果。而vhs ，vss与 ghs模型在低温和过高温度下那么与实验数据有一定的差距，gss模型在低温和极性分子的情况下是

6、特别有意义的。文献在前人根底上，提出了gss-3分子碰撞模型，能够给出正确的气体输运性质。2.2 碰撞抽样几率函数的研究在气体分子碰撞对的抽样方面， bird 首先提出了“时间计数器方法 tc ，该方法具有较高的效率而且同时保证了在较大网格分子数下正确的碰撞频率。但在极度非平衡的情况下， 如强激波的前缘问题， 由于碰撞概率与网格推进的时间关系，1)2/(rrcncnt，偶然中选但具有很小几率的碰撞，即较小的rrc，会使网格推进的时间ct过大，导致计算碰撞对的减少，从而使碰撞频率发生畸变。koura 给出了“零碰撞方法，可以克制tc方法的这一缺陷且不消耗较多的时间，但是不利于计算机的向量化计算。

7、bird 建议的非时间计数器法ntc 在保持了 tc法效率的同时， 克制了 tc法中的困难， 并且可以在程序碰撞过程开场之前事先算出网格单元的碰撞数，而且使得dsmc 方法能够实现向量化计算，所以这种方法在今天被广泛的应用。2.3 能激发与化学反响模拟的研究dsmc 模拟中，一般采用 larsen-bergnakke l-b模型模拟分子在碰撞中能量的交换与分配， 保证能的激发和松弛速率与实验所给出的结果一致。但是由于振动自由度为温度的函数， 在采用模型分配碰撞中的各种能量时存在奇异性，因此文献开展了一种推广的积累分布取舍联合抽样法，给出了有单奇异和双奇异性分布的抽样方法，克制了l-b 模型的困

8、难。在 dsmc 方法的模拟中，化学反响通常是与分子自由度松弛过程耦合在一起的。bird认为分子与分子碰撞时，以一定的几率发生化学反响， 并提出了处理化学反响的位阻因子概念，其值代表了. 3 / 8 弹性碰撞导致化学反响的几率。通常引入反响截面与碰撞截面的比值tr表示弹性碰撞导致化学反响发生的几率。文献建议了一个依赖空间取向的化学反响模型，它从一个双原子分子和另一个粒子发生碰撞导致其化学键断裂而产生离解或交换反响的微观判据出发，得到了化学反响速率常数的表达。3.dsmc 中碰撞的数值计算方法3.1 仿真分子碰撞对抽样方法1bird 的 time-counter 方法由分子气体动力学可知同组元气

9、体的分子平均碰撞频率为gnt式中符号“表示平均值。上式说明对于两个特定的气体分子，它们可能发生碰撞的几率与分子对碰撞截面t和相对运动速率g的乘积成正比，亦即( )coltpgg因而网格单元模拟分子对的碰撞几率函数可表示成max( )()tcoltgpgg于是对应于硬球分子模型， 抽样几率函数的表达形式为max( )/colpgg g，对应于负幂率分子模型，表达形式为1 4/1 4/max( )colgpgg。另一方面，还可得到在时间步长mt各网格单元中模拟分子间应发生的碰撞次数tn 的表达式为12tmtmnn ng t式中mn是网格单元中模拟分子的总数目，n是网格气体分子的数密度， 因子1/

10、2是由于每次碰撞都有两个分子参加而得到的权因子。对于硬球分子模型，为常数，tn 的计算式为12tmtmnn ng t. 4 / 8 对于负幂律模型，tn 的计算式为1 4/12tmtmnn ngt如果按照上述表达式中的任何一式计算tn ，都势必需要计算一个以气体分子相对运动速度g为自变量的函数的平均值。 由于计算这些平均值所需的计算机时是与单元所具有的分子数的平方成正比，因而需要消耗大量机时。 为了解决这一困难， bird 提出了 time-counter scheme，即所谓的“超前时间累加方法。这个方法是在每个单元设置一个计时器，当模拟分子在网格单元发生了一次碰撞，就在计时器上累加一个 “

11、超前时间cit 。对于不同的分子作用势模型， “超前时间cit的计算式也各不一样，对于硬球分子和负幂律分子，cit的计算式分别为2cimttn ng1 4/2cimttn ng一旦计时器的累加时间大于mt ，亦即cimitt，就停止该网格单元模拟分子的碰撞计算。2no time counter (ntc) 抽样方法在 dsmc 中，tc方法具有较高的效率同时保证了在较大网格分子数下正确的碰撞频率， 但在极度非平衡的情况下， 如强激波的前缘问题， 由于碰撞概率与网格推进的时间关系，1(/ 2)crrtnnc，偶然中选但具有很小几率的碰撞，即较小的rrc，会使网格推进的时间ct过大，导致计算碰撞对

12、的减少，从而使碰撞频率发生畸变。 bird 建议的 ntc法在保持了 tc法效率的同时，又克制了tc法中的困难， 并且可以在程序碰撞过程开场之前事先算出网格单元的碰撞数，这里介绍一下 ntc法的思想。考虑体积为cv 的网格单元，每一个模拟分子所代表的真实分子数为nf ，在t时间组成碰撞对的两个分子发生碰撞的几率p等于碰撞截面r以相对速度rc扫过的体积与网格体积之比，即. 5 / 8 /nrrcpfct v在上式中，rc 随碰撞对的选择变化，r是rc 的函数，其它量不依赖于碰撞对的选择。由于网格单元的真实分子数是cnv，平均模拟分子数/cnnnvf，因此在dsmc 中，通过式 (2.9) 需要计

13、算网格可能碰撞对的数目为(1)/2n n。由于碰撞几率p在一般情况下是非常小的量， 碰撞数(1)/ 2n n几乎正比于网格中分子数n的平方，因此如果按p的计算式逐个考察所有(1)/ 2n n个碰撞对，那么计算效率必定非常低。 为了提高效率 bird 建议仅需考虑所有碰撞对(1)/ 2n n的一局部并同时将p按一样的分数加以放大，如果该分数能够使最大的概率为1那么计算效率为最高，因此理想的分数即为最大碰撞概率maxmax()/nrrcpfct v通过maxp与2/ 2n的乘积可以得到所需考察的碰撞对数。然而在多数情况下网格单元数n是波动量，由于平方的平均不同于平均的平方，因此2n 应当被瞬时量n

14、与时间平均或系统平均n 的乘积代替。由以上讨论可得ntc方法中需要考察的分子碰撞对数与碰撞几率的表达式maxmax()22nrrcfctnnnnpvmax()rrrrcpc该方法到多组元混合气体的推广是直接简单的，在混合气体中， 网格单元组元p与组元q碰撞需要考察的数量与碰撞几率分别为maxmax()22pqpqnrrpqcn nn nfctpvmax()rrrrpqcpc3.2 碰撞计算由于在碰撞过程中不存在外力， 碰撞的双分子体系遵循质量、 动量和能量守. 6 / 8 恒定律。记ab为分子a和分子的b简并质量，g为碰撞分子的相对运动速度，g为碰撞分子的质心速度，那么abababm mmma

15、bgaabbabmmgmm根据质量和动量守恒定律有gg其中上标 “表示碰撞后的物理量。 可见气体分子间无论发生何种类型的碰撞，分子体系的质心速度都将保持不变。无反响碰撞分子的能量守恒方程为intintintintrelabrelab其中212relabg为碰撞分子的相对平动动能。从上式可见，如果碰撞分子之间不发生化学反响，分子能量仅在相对平动能与分子能之间传递。记123(,)nn nn为碰撞后分子相对运动速度g 的方向单位矢量， 由于相对运动速率g 可由能量守恒方程确定得到， 于是只需确定n。考虑到分子碰撞的随机性之后， bird 根据所提出的“可变硬球模型构造了确定n的方法为11cosn21

16、2sincosn312sinsinn其中，1和2为 euler 角，由随机抽样确定。其抽样方法为12cos11r222 r式中，1r、2r为随机数。由此可见bird 散射模型是一种球对称模型，认为相对运动速度方向是以球对称形式散射的。. 7 / 8 由能量守恒方程可确定得到碰撞后分子相对运动速度的模值，然后抽样得到相对运动速度的方向矢量，于是可得碰撞后单个分子的运动速度为abaaggmabbbggm3.3 分子碰撞中能量的交换在 dsmc 中，处理分子碰撞中能量交换的常用模型是由larsen 和 bergnakke引入的唯象论模型。其中心思想是假设碰撞中的相对平动能和能遵守能量守恒，碰撞后的能

17、按照动能和能组合的平衡分布取值，而能量松弛过程的速率通过调节弹性碰撞和非弹性碰撞的比例加以确定，使其满足实验的结果。 这里给出混合气体中不同组分的分子间larsen-bergnakke 模型在 vhs 或 vss 模型下的实现过程。在混合气体中组分1 与组分 2 碰撞对的总能为,1,2iiie组分 1 和组分 2 的能分布函数满足121,12,1,11,22,2,2()exp()()exp()iiiiiifektfekt组分 1 的能,1i到,1,1iid，组分 2 的能,2i到,2,2iid的比分为121122i,1,1,1,2()exp()iiiiieeddkt先固定,1i，此时,2iid

18、ed，对,1i从 0 积分到ie 得到总能从ie 到ide 的比分正比于1212()exp()iiiiief e deedekt如果对于不同组分可以定义出一个粘性系数的温度幂指数1,2，那么碰撞中的平动能t的分布函数可以写为1,232()exp()tttfkt. 8 / 8 由上述方程可以得到在碰撞中的总能量ctiee 中平动能t的分布函数正比于121,23122()exp()ctcteekt，由于t由ce所决定，因此为常数，所以指数项是常数，碰撞中的平动能分布正比于121,23122()tcte在碰撞中总能量,1,2ctitiiee保持不变。在非弹性碰撞中，碰撞后的平动能量*t和能从下式用取舍法进展取样121,23*122()()(1)tttcccfeee碰撞后的能*ictee在两个分子间的分布按照下式用取舍法进展取样12*11,1,1,122*()()(