带电粒子输运Fokker-Planck方程的有限元数值解.pdfVIP

带电粒子输运Fokker—Planck方程的 有限元数值解 王尚武银燕王同权 (国防科技大学应用物理系，长沙，410973) 的数值解，并由此分别计算了a粒子对背景等离子体中的离子、电子的能量沉积率随时空的变化。 一、物理方程 程： (1) 其中：^为产生带电粒子的独立源项，满足I正口2dvJ2rcdlt=n肼T(∞)Dr． A(。)=矗Sn一≯习m．，(即) B‘?)=童”西去卜暴讹)+务。] c㈤=‘砉”西嚣吲即) 旷(糕)”， 昂=堡筝坠 如(x)=习2。xe一，2咖一争一， 方程(1)之定解条件： o 『，(r，口，乒，￡)It=。2 j，(r=0，口，芦，￡)=y(r=0，口，一产，t)(u0) 1，(r=R，口,／1，￡)I删=0 (2) l，(r，口，P，￡)I。、=，(r，_。p，￡)I。％=0 =、F—P方程有限元数值解法 用“(≈2 分戒G群，其中第g群的速度范围为[％一1，2，％+1，2]，用算符： L：上fl^+锄卜n口2dv q。～ 。、．1n 口5△“’△口；，、△f^=fl+t一“t△u；={(_玉1在一口；一l，2) 一343 向的第”z个区间共同组成(r，卢)平面上的～个长方形单元e． 将单元e内任一点(r，∥)处的分布函数，(r，产)用四个节点上的函数值用双线性插值表示： ，(r．F)=髓T，。((r，卢)∈e) 其中研7=五i圭甄[Cldltc2dl，c2d2,Cld2]，于=[^，尼·f3，^]7 而 △F。=‰～Pm～1．△h=“一rH～I， ‘1=Fm一∥，C2=／a～户㈣1dt=rn—r。d2=，一rn-l 从而得到单元e内的离散方程，再用Galerkln方法在单元e内的离散方程两边乘以列矩阵勋，并作 rt rR r2出，得到有限元方程： 积分l 2nd／al 2 g／g-1 (g=1，，…．G)(3) 三。j。群_‘+1／2+豆筑‘“在+e巩{}1在]=手[a％’‘+E：】 ‘L 其中， 一=』2嘶『r2d廊r A；=箍(铀一(1-d-)一％h‘ 一 磁= 一岢·{z嘶弘励毒[c1--p2，警] +学[{z呻』r2瓣誓+{z呻』r2旃掣警] 毽=一遗(站-础+％)旧 戢=丽1，研Atk洲。)啦。{2nd／a!rzd前 根据能量边界条件，(3)可以写成以下矩阵形式 玩q ∥“41 日 ～ A§BlQ 弱 Iml (4) A台 B毛 魄“1刊 军曩E毛