【技术实现步骤摘要】
基于双线性变换的一维真空Crank-Nicolson完全匹配层 实现算法
本专利技术涉及数值仿真
，特别涉及一种基于双线性变换的一维真空 Crank-Nicolson完全匹配层实现算法。
技术介绍
近年来，时域有限差分方法（FDTD)作为一种计算电磁方法被广泛地应用于各种 时域的电磁仿真计算中，如天线，射频电路，光学器件和半导体等等。FDTD具有广泛的适用 性、适合并行计算、计算程序通用性等特点。 然而，随着科学研究的深入和各种越来越广泛应用的需求，其算法本身受 Courant FriedrichsLewy(CFL)数值稳定性条件的限制的缺陷越来越凸显出来。算法 本身所受数字稳定性条件限制：在计算过程中时间步长和空间步长必须满足CFL约束 条件，即
【技术保护点】
基于双线性变换的一维真空Crank‑Nicolson完全匹配层实现算法，其特征在于具体设置步骤： 步骤1：将频域中麦克斯韦旋度方程修正为带有拉伸坐标算子的麦克斯韦方程； 步骤2：将频域中修正后的一维麦克斯韦旋度方程在直角坐标系中表示； 步骤3：根据频域和z域的映射变换关系，将直角坐标系中的一维麦克斯韦方程变换到z域表示； 步骤4：根据双线性变换方法，将拉伸坐标变量由频域变换到z域，并代入到一维麦克斯韦方程的z域表达式，同时设置辅助变量； 步骤5：基于Crank‑Nicolson时域有限差分算法的时域展开形式，以及根据z域和时域的映射关系，将z域形式的直角坐标系中一维麦克斯韦旋度方程展开成时域有限差分的形式，同时也将z域形式的辅助变量变换为时域有限差分的形式； 步骤6：将时域有限差分形式的方程整理成求解的形式，结果产生一组电场和磁场耦合的方程，是一组隐式方程； 步骤7：将这组隐式方程进行去耦，即将磁场分量的方程代入到电场分量的方程中，同时将辅助变量的显式方程组代入； 步骤8：将代入磁场和辅助变量后的电场分量的迭代方程进行整理，整理后获得等式左边为三对角矩阵形式的系数的电场显式迭代方程； ...
【技术特征摘要】
1. 基于双线性变换的一维真空Crank-Nicolson完全匹配层实现算法，其特征在于具 体设置步骤： 步骤1 :将频域中麦克斯韦旋度方程修正为带有拉伸坐标算子的麦克斯韦方程； 步骤2 :将频域中修正后的一维麦克斯韦旋度方程在直角坐标系中表示； 步骤3 :根据频域和z域的映射变换关系，将直角坐标系中的一维麦克斯韦方程变换到 z域表示； 步骤4 :根据双线性变换方法，将拉伸坐标变量由频域变换到z域，并代入到一维麦克 斯韦方程的z域表达式，同时设置辅助变量； 步骤5 :基于Crank-Nicolson时域有限差分算法的时域展开形式，以及根据z域和时 域的映射关系，将z域形式的直角坐标系中一维麦克斯韦旋度方程展开成时域有限差分的 形式，同时也将z域形式的辅助变量变换为时域有限差分的形式； 步骤6 :将时域有限差分形式的方程整理成求解的形式，结果产生一组电场和磁场耦 合的方程，是一组隐式方程； 步骤7 :将这组隐式方程进行去耦，即将磁场分量的方程代入到电场分量的方程中，同 时将辅助变量的显式方程组代入； 步骤8:将代入磁场和辅助变量后的电场分量的迭代方程进行整理，整理后获得等式 左边为三对角矩阵形式的系数的电场显式迭代方程； 步骤9 :利用追赶法求解系数为三对角矩阵的电场迭代方程，得到电场分量的值； 步骤10 :将求解出的电场值代入到磁场的迭代方程中，求解出磁场分量的值； 步骤11 :将求解出的电场值和磁场值代入到辅助变量的迭代方程中，求解出辅助变量 的值；返回到步骤9,循环步骤9-11，从而在时间上迭代求解。2. 由权利1所述的基于双线性变换的一维真空Crank-Nicolson完全匹配层实现算 法，其特征在于：步骤4,利用双线性变换方法将拉伸坐标变量由频域变换到...
【专利技术属性】
技术研发人员：李建雄，于洋，史伟光，
申请(专利权)人：天津工业大学，
类型：发明
国别省市：天津;12
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。