[发明专利]一种等离子体中使用电流密度卷积完全匹配层的实现方法

摘要

本发明公开了一种等离子体中使用电流密度卷积完全匹配层的实现方法，包括：输入模型文件；初始化参数、设置PML系数和吸收边界参数；分别更新计算整个计算区域y方向上和x方向上电场分量系数电场分量系数添加场源到磁场分量系数，并更新计算整个计算区域的磁场分量系数；更新计算整个计算区域的极化电流密度更新计算整个计算区域的电磁场分量系数的辅助变量；更新计算观测点处电磁场分量；将q+1赋值给q，并判断拉盖尔多项式的阶数q是否达到预设值，若未达到预设值，返回步骤3；若达到预设值，则结束。本发明的一种等离子体中使用电流密度卷积完全匹配层的实现方法，计算速度快，内存消耗小，且对于低频和凋落波具有很好的吸收效果。

主权项

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_ζ表示E_x,E_y,H_z，其中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；式中，η₀为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)；C5＝2e2ne/(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，添加场源到磁场分量系数中，并更新计算整个计算区域的磁场分量系数；其中，所添加场源的表达式为：Imz(t)＝sin(2πf(t‑t0))×exp(‑(t‑t0)2/τ2)；式中，t0，τ为场源参数；具体更新公式为：步骤6，更新计算整个计算区域的极化电流密度具体更新公式为：式中的分别表示x方向上极化电流密度分量系数和y方向上极化电流密度分量系数；步骤7，更新计算整个计算区域的电磁场分量系数的辅助变量，具体更新公式为：步骤8，更新计算观测点处的电磁场分量，具体按照以下公式更新计算：上式中U表示电磁场分量E_x,E_y,H_z，U^q表示q阶电磁场分量系数，是q阶加权拉盖尔多项式，是带有时间尺度因子s>0的扩展时间；步骤9，将q+1赋值给q，并判断拉盖尔多项式的阶数q是否达到预设值，若未达到预设值，则返回步骤3；若达到预设值，则结束。