基于Fokker-Planck方程的等离子体模拟.doc
《基于Fokker-Planck方程的等离子体模拟.doc》由会员分享,可在线阅读,更多相关《基于Fokker-Planck方程的等离子体模拟.doc(14页珍藏版)》请在沃文网上搜索。
1、 基于Fokker-Planck方程的等离子体模拟ALI SHAJII,DANIEL SMITHMKS Instruments, 90 Industrial Way, Wilmington, MA, USA, 01887对于各个计算组织,等离子体的模拟一直是个极大的挑战,有很多不同近似程度的模拟计算方法。包括完整的动力学计算方法,流体近似方法和关于漂移扩散方程的方法。近几年来,有人用Fokker-Planck方程处理等离子体中的电子,同时把离子当作流体进行耦合计算,获得了很好的计算结果。本章我们将介绍基于通用Fokker-Planck方程的计算求解过程,并通过一个具体实例得到电容放电过程的电子
2、密度分布。希望通过该简单模型使读者对等离子放电建模过程有个初步的了解。1引言各种工业等离子体应用“过程”中都存在一个关键步骤23。历史上曾采用各种不同方法对等离子体进行简化建模,分别对应于不同层面问题所需准确性357。这些层面包括:l 完整的动力学模型(多组分Boltzmann方程)4;l 使用Monte-Carlo方法的颗粒模拟3;l Fokker-Planck近似17;l 多尺度流动模型(也被称作漂移扩散模型)3。出于种种原因,使得等离子体的建模和模拟非常困难。首先,最直接的使用多流体方程的模型不能反应相关的等离子体物理过程。其次,“水动力学”系数完全取决于研究的特定问题,不能作为纯气体或
3、液体的常数简单测量。最重要的一点是,完整的动力学模型包括Boltzmann方程,计算求解非常困难。对于完整动力学模型和流动模型之间的需求空白,通常采用Fokker-Planck (FP)近似或者Monte Carlo (MC)颗粒模拟。这两种方法可以在所需计算复杂度和捕获等离子体重要物理细节之间找到一个很好的平衡。本章的主要目的是展现用COMSOL Multiphysics求解FP方程的功能。为了对该问题给出一个整体认识,我们把侧重点集中在一个简单的例子上。特别是在第二节我们通过一个简单的例子对FP方程给出一个直观描述,将它用于布朗运动的颗粒模拟。本章最主要的贡献在于介绍了如何在外部电场情况下
4、对电子动力学过程进行建模。最后,在第四节中对如何使用COMSOL Multiphysics实现该模型给出了详细的讨论。2一维FP和Langevin方程对FP方程直接求导可以得到Langevin方程1。考虑浸没在流体中的“布朗”颗粒,如果颗粒足够小,会同时受到两种力。一是颗粒和流体介质间的粘性力,它会降低平均颗粒速度。二是颗粒与流体“分子”间随机碰撞的力。该布朗颗粒的运动控制方程如下1: (1)其中是阻尼系数,随机项表示颗粒与背景流体的连续碰撞。对于这个简单的例子,通常我们假设粘性力线性依赖于颗粒速度。同时,根据随机近似,Langevin力必须满足以下方程1:其中 1,表示整体平均,T是流体温度
5、,k是Boltzmann常数,m是布朗颗粒的质量。同时注意到力可以很容易的通过MATLAB函数randn(见第四节)实现。给定适当的颗粒数和初始条件后,使用MC方法可以很容易的求解方程(1)。为了对颗粒在时间t时的速度进行精确静态测量,颗粒数通常需要超过一百万。最简单的初始条件为。总之,Langevin方程描述了背景流体介质中布朗颗粒的运动。当然,如果没有随机力,颗粒的路径为。但是由于的存在,方程(1)对于很多颗粒算出的解通常是条分布曲线,分布的宽度由q决定1。FP方程为MC方法提供了另一种思路。包含随机力的Langevin方程可以等价与以下偏微分方程: (2)初始条件和边界条件为: (3)其
6、中w是颗粒在速度空间中的分布概率,也就是颗粒出现的概率,在无限大系统中,速度通常表示为wdv。初始条件(3)如图1所示。求解该偏微分方程等价于对无限个颗粒进行Monte-Carlo模拟,以获得t0时的概率分布函数。图2给出了两种方法的比较。MC方法针对初始条件为的情况,求解了20000个常微分方程。这还不是强非线性阻尼力和三维空间情况下的颗粒运动,即使如此,对于这个简单的例子,求解偏微分方程的时间和常微分方程一样。当颗粒痕迹的常微分方程是非线性时,或者需要从FEA模拟中查询外部力的值时,MC方法的计算时间就会变得非常大。图1 初始条件w图2 各时间步长下Fokker-Planck和Monte
7、Carlo方法的对比作为例子,我们考虑修改后的Langevin方程: (4)和相应的FP方程1: (5)其中,。这种情况下的分布函数如图3所示。即使计算一百万个颗粒,常微分方程的分布函数也存在很大的数值噪声,而FP结果就很平滑。这里出现的静态噪声是有限颗粒MC方法模拟的经典问题(也就是说数据的后处理是很精细、很关键的步骤)。图3 各时间步长下基于非线性Langevin方程的Fokker-Planck和Monte Carlo方法的对比下面,根据本节介绍的背景知识,我们来看一个基于FP方程的更为复杂的例子。3外加场中的电子动力学现在我们考虑电子流在电容两极板间外加电场作用下的动力学过程。一个极板保
8、持常电压,另一个极板接地。假设电子流从接地极板侧以一定速度引入,并在阳极侧积累。并且两极板间的整个区域充满了温度为T的静止流体。图4给出了该几何结构。图4 模拟域示意图模型主要假设为:l 忽略由于电子和背景介质碰撞所产生的离子;l 电子间的相互作用可以忽略;l 在一维情况下求解该问题。为了精确求解物理空间(y)和速度空间(vy)的电子分布,需要求解FP方程。在本例中我们还要计算电子密度、电子流和电容两极板间的静电场。在通常情况下,FP方程比静电场问题多需要一个空间维数,所以我们需要同时求解一个一维静电场和二维FP方程。电子分布的控制方程为: (6)对于该电场 (7)用标量形式电势和电子密度表示
9、为 (8)其中 (9)并且 (10)注意这里的ne和Fy都是y的函数。边界条件和初始条件为: (11) (12) (13) (14)其中m是电子质量,k是Boltzmann常数,T是背景流体温度,e是电子电荷,是自由空间介电常数,是阻尼系数。同时,无限大系统中,在速度vy和位置y处找到一个电子的概率为wdvdy。方程(6)和(8)是两个耦合的不同空间维数的偏微分方程。由于COMSOL Multiphysics允许不同空间维数的任意数量几何形状进行耦合,所以我们现在同时求解这两个方程。通过投影耦合变量方法将电子密度引入静电场问题中,通过导出耦合变量的方法将颗粒上的力引入FP问题中。概括起来讲,电
10、子流以速度v0进入“电容放电”区域,同时受到三种力:,和。这些力导致电子速度和密度分布沿y方向变化,所以与之前简单线性分布例子相比,当修改极板电势V时,不存在这种情况。图5给出了函数w的解。注意到初始时刻,当电子进入计算域时速度降低。这种非线性、非单调曲线是由如图6所示的阴极电子屏蔽所导致的。查看域中间位置附近的最大电子密度。图8表明电容的电子电势在最大ne处有最小值。图5 概率分布图6 电子密度图7 平均速度图8 平行极板间电子电势随位置的分布最后,给出电子流量:由于质量守恒,流量在y方向应该保持为常数。图7沿y方向的电子流量表明该过程中质量守恒。在下一节中我们将详细给出在COMSOL中如何
- 1.请仔细阅读文档,确保文档完整性,对于不预览、不比对内容而直接下载带来的问题本站不予受理。
- 2.下载的文档,不会出现我们的网址水印。
- 3、该文档所得收入(下载+内容+预览)归上传者、原创作者;如果您是本文档原作者,请点此认领!既往收益都归您。
下载文档到电脑,查找使用更方便
20 积分
下载 | 加入VIP,下载更划算! |
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 Fokker Planck 方程 等离子体 模拟