1.一种等离子体中使用电流密度卷积完全匹配层的实现方法,其特征在于,包括以下步骤:
步骤1,输入模型文件;
输入的模型文件具体为:计算区域大小Nx×Ny,其中Nx为x方向的网格数,Ny为y方向的网格数;空间步长Δη,η=x或者η=y,x为横坐标,y为纵坐标;时间步长Δt;真空中的电导率σ,磁导率μ0,介电常数ε0;等离子体中的碰撞频率υ与等离子体中的电子密度ne;等离子体在计算区域中的位置;吸收边界层数NPML与相关参数κηmax,αηmax,σηmax;κηmax取整数,κηmax取值范围为[1,60];αηmax取值范围为[0,1);σηmax/σopt取值范围为(0,12];仿真计算时长Tf;加权拉盖尔多项式的阶数q,q≥0且为整数;时间尺度因子s,s取值范围为[109,1013];观测点;场源参数;
步骤2,初始化参数和设置参数;
初始化的参数具体包括:将整个计算区域的电磁场分量系数 整个计算区域的极化电流密度 整个计算区域的电磁场分量系数的和整个计算区域的极化电流密度的和 整个计算区域的辅助变量和 拉盖尔多项式 全部初始化为零;其中Fζ表示Ex,Ey,Hz,其中PML系数(C1η,C2η,C3,C4,C5,C6)初始化为C1η=1/(1+0.5ε0s),C2η=1,C3=ε0/μ0,C4=2/(ε0s),C5=0,C6=2;式中,e,m分别是电子的电量和质量;
设置的参数具体包括:
设置CFS-PML吸收边界的参数ση,κη,αη;具体为:ση=σηmax|η-η0|m/dm;
κη=1+(κηmax-1)|η-η0|m/dm;
αη=αηmax(d-|η-η0|)/d;
式中,η0为PML层与非PML截面位置,d是PML吸收边界的厚度,σηmax根据σopt来设置,σopt=(m+1)/150πΔη,m取值范围为[1,20],Δη取值范围为 λ为源的波长;
设置PML系数C1η,C2η和与等离子参数相关的系数C5,C6,具体为:C1η=1/(κηαη+ση+0.5κηε0s),C2η=(2αη/ε0s+1);
2
C5=2ene/(ms+2mυ),C6=s/(0.5s+υ);
步骤3,更新计算整个计算区域的y方向上电场分量系数具体为:
步骤3.1,给出电场分量系数 在计算区域的方程,如下所示:式中,i表示横坐标上的第i个计算网格,j表示纵坐标上的第j个计算网格;
步骤3.2,使用追赶法对整个计算区域的电场分量系数 进行求解;
步骤4,根据整个计算区域的y方向上电场分量系数 来更新计算整个计算区域的x方向上电场分量系数 步骤4具体为:
步骤4.1,给出电场分量系数 在计算区域的方程,如下所示:步骤4.2,使用追赶法求解系数为三对角的整个计算区域的电场分量系数步骤5,添加场源到磁场分量系数中,并更新计算整个计算区域的磁场分量系数;
其中,所添加场源的表达式为:
2 2
Imz(t)=sin(2πf(t-t0))×exp(-(t-t0) /τ);
式中,t0,τ为场源参数;
具体更新公式为:
步骤6,更新计算整个计算区域的极化电流密度 具体更新公式为:式中的 分别表示x方向上极化电流密度分量系数和y方向上极化电流密度分量系数;
步骤7,更新计算整个计算区域的电磁场分量系数的辅助变量,具体更新公式为:步骤8,更新计算观测点处的电磁场分量,具体按照以下公式更新计算:上式中U表示电磁场分量Ex,Ey,Hz,Uq表示q阶电磁场分量系数, 是q阶加权拉盖尔多项式, 是带有时间尺度因子s>0的扩展时间;
步骤9,将q+1赋值给q,并判断拉盖尔多项式的阶数q是否达到预设值,若未达到预设值,则返回步骤3;若达到预设值,则结束。