使用有限元方法对各向异性等离子体鞘层中TE波散射进行数值模拟：验证与内部场分析

《Optik》：Numerical Modeling of TE-Wave Scattering from an Anisotropic Plasma Sheath Using the Finite Element Method: Validation and Internal Field Analysis

电磁波与等离子体相互作用机理及数值建模研究。本研究基于有限元素法构建高精度数值模型，系统分析了轴对称导电圆柱表面冷等离子体鞘层（各向异性、径向不均匀）对H波（TE波）的散射特性，重点揭示近场轴向磁场分布规律。通过网格收敛性和解析基准验证（NRMSE≤1.08%），建立了包含等离子介电张量正确定义与多物理场耦合的可靠仿真框架，发现散射共振源于混合模态的局域激发，为可调谐隐身器件设计提供了理论支撑。

　　

阿尼斯·索尔塔尼（Anis Soltani）

伊朗阿拉克大学（Arak University, Arak, Iran）

摘要

在这项研究中，我们提出了一个精确的数值模型，用于分析经典的H波（TE波）散射问题。该问题涉及一个完全导电的圆柱体，周围环绕着一种冷的非各向同性等离子体鞘层，该等离子体鞘层具有径向不均匀性。我们采用有限元方法（FEM）对这一现象进行了分析。尽管Rusch和Yeh在1967年首次研究了这种结构的远场散射特性，但控制共振现象的内部物理机制仍大部分未被探索。在这里，我们在COMSOL Multiphysics中开发了一个高保真的色散模型，通过严格的网格研究和对等离子体介电张量的正确处理确保了模型的收敛性。通过与分析基准的定量验证，我们发现该模型的性能非常出色，对于薄等离子体情况，归一化均方根误差（NRMSE）低至1.08%。这项工作的主要创新在于对等离子体鞘层内部轴向磁场分布（$H_z$）的分析和可视化。这种近场分析表明，尖锐的散射共振是由在非各向同性层内强烈限制的混合鞘模式激发引起的。通过将密度剖面梯度与驻波模式相关联，本研究加深了对波-等离子体相互作用的物理理解，并为设计可调的隐身装置提供了坚实的框架。

引言

电磁波与等离子体介质的相互作用是电磁学中的一个基本课题，对航空航天技术和国防应用具有关键意义。它在卫星通信[2]、确保航天器再入时的信号稳定性[3][4]以及优化雷达吸收涂层[5][6]方面尤为重要。在高速飞行条件下，车辆周围会形成一层电离等离子体。由于电子密度的空间变化以及外部磁场（例如地球地磁场）的存在，这层等离子体鞘层会变成一个非均匀且各向异性的介质[7][8]，这极大地影响了入射波的雷达截面（RCS）和散射特性[9][10]。

Rusch和Yeh（1967年）[1]在这个领域建立了开创性的基准研究，他们为一个涂有非各向同性等离子体鞘层的无限长导电圆柱体推导出了H波（TE波）散射的解析解。他们的发现为后续的许多研究[12][13]奠定了基础。后来的发展将这一框架扩展到了各向同性等离子体涂层[11]、斜入射配置[16][17]，以及磁活性鞘层中的模态或衍射方法[18][19][20]。最近，数值技术的进步——如FDTD模拟碰撞磁化等离子体[4]和与基于等离子体的天线相关的建模[14][15][21]——加深了经典散射理论与应用波工程之间的联系。

与等离子体电磁学并行发展的是，现代光子学也独立强调了各向异性和可调环境在先进波操控中的关键作用。值得注意的例子包括利用双折射特性的超高重复率微光纤谐振器[22]、基于γ-MnO₂材料的双芯光纤器件[23]、宽带光调制器[24]，以及实现太赫兹重复率的同步多波长孤子光纤激光器[25]。尽管这些光子技术在不同的频率域内工作，但它们都与等离子体鞘层有共同的物理基础：即依赖于材料各向异性和模态工程来控制波的传播。

受到这些发展趋势的启发，本文采用了一个严谨且可重复的计算框架重新审视了经典的等离子体散射基准问题。主要目标是首先验证Rusch和Yeh[1]报告的关键散射特性，其次提供对等离子体鞘层内部轴向磁场分布（$H_z$）的新分析——这是原始研究中缺失的部分。这种内部视角阐明了共振状态下场集中和驻波形成的机制，明确了它们与在等离子体[16][18][19][20]以及复杂各向异性介质[22][23][24][25]中观察到的混合模式特征之间的联系。

我们对鞘层厚度、磁偏压和电子密度剖面（抛物线型与准平面型）进行了系统的参数探索，以量化它们对前向/后向散射平衡和共振锐度的影晌。重点在于透明地报告控制方程、旋光张量的处理方法以及网格统计信息，以便能够进行准确的复制[4][14][21]。这项工作的创新之处在于将经典基准问题扩展到现代FEM框架中，提供了定量验证，并引入了将宏观散射响应与混合鞘模式激发联系起来的内部场分析。

章节摘录

理论与方程

在这项研究中，我们分析了一个由半径为、长度无限的完全导电圆柱体组成的经典配置，该圆柱体被一个冷且磁化的等离子体鞘层包围，该等离子体鞘层在径向上不均匀分布，并向外延伸至半径为$a$（见图1）。

该结构被一个在+x方向上传播的时谐平面H波（TE波）照射，其磁场矢量与圆柱体轴线平行（$H_{\mathrm{inc}}$）。假设时间依赖性为

数值实现与验证

为了确保有限元模型的准确性和可重复性，我们使用COMSOL Multiphysics进行了专门的数值验证。这包括网格收敛性研究和与Rusch和Yeh[1]建立的解析基准的定量比较。

虽然本研究的结果主要使用归一化无量纲参数来表示以确保可扩展性（例如$k_{0}a$$k_{p}a$$k_{g}a$），但模拟是使用具体的物理值进行的

结果与讨论

为了准确分析散射问题，我们基于有限元方法（FEM）开发了一个二维计算模型，并使用COMSOL Multiphysics进行了模拟。在该模型中，物理配置——包括内部导电圆柱体、非均匀且各向异性的等离子体鞘层以及由完美匹配层（PML）界定的周围自由空间区域——被仔细模拟。入射激励被定义为平面波

结论

本研究使用有限元方法（FEM）对一个涂有各向异性、径向不均匀等离子体鞘层的导电圆柱体的电磁散射进行了严格的数值分析。与以往主要关注远场雷达截面（RCS）模式的解析研究不同，本研究将研究范围扩展到了近场领域，揭示了控制波-等离子体相互作用的物理机制。色散FEM模型的准确性得到了验证

CRediT作者贡献声明

阿尼斯·索尔塔尼（Anis Soltani）：撰写——审稿与编辑、原始草稿撰写、可视化、验证、软件开发、方法论研究、数据分析、概念化。

利益冲突声明

作者声明他们没有已知的可能会影响本文所述工作的财务利益或个人关系。