输运方程谱有限元耦合方法及在核工业中的应用.pdfVIP

输运方程谱有限元耦合方法 及其在核工业中的应用 梅立泉 西安交通大学理学院，710049 摘 要 本文研究求解Boltzmann方程的谱有限元耦合方法，给出了求解输运方程的谱标 准Galerkin有限元、谱流线扩散有限元、谱间断有限元耦合格式。研究这些格式的收敛 性，并将它们应用到了核测井领域。 考虑中子和1光子的输运过程所满足的Boltzmann输运方程 三皇!￡掣+五v圣(7，菇，E，t)+∑。(事，E)中(7，五，E，￡) 优 dE， ：，f…毪(彳，E，)，(矗7，E，^÷菇，E)心(7，五’，E一)d矗7 J4ⅡJ0 +寿s(7川E) (1) 源，巾(彳，五，厅)为粒子角通量。 将方程(1)中与方向有关的量：中子通量m．和反应指示函数，用球谐函数展开， ∞ n ／ (2) f(E’_÷厅，西’_5=j)=∑∑，儿(矽_F)％，。(目，妒)K。(p’，P，)’ n=Om=一n ∞ n 圣(寸，矗，F)=∑∑西。(7，E)z。(D，妒)， (3) n=0mtl J 然后将展开式代入方程(1)，并利用球谐函数的正交性可得方程的一种等价形式．取展 、 开式中有限项，即取n=0，1，2，…，N，则得到方程(1)的近似形式，一般称它为％近 似方程，再关于能量变量使用能量分群技术虬并在求解区域移上，我们赋以内流边 、 界条件得到，对每一个能群，‰近似方程有如下的向量形式 肿 || + 旧 哪=Q in口 Ⅳ (4) 刘 on 一 + 0口 Ⅲ 磋m q即 锄 唑曲” 等舭 10 x 其中G，G，巴，∑，M为2”+12N+I的矩阵，口，m为2”+1阶向量，而向量(n。，～，n：) 为边界彻的单位外法向．以下记_Ⅳ，_2N+I令{而：兀={71))为空间口上的有限元 剖分，定义有限元空间 ％=w∈H1∽)：妒IT∈Pr(T)VT∈^) ㈣、’ 【(L田6，置)十{(一D)皿“，互=(0，三)V兰∈(H)Ⅳ’ c圳惯：，三f?嚣筹