波的反射、折射与吸收
在等离子体物理理论与建模中，波与介质的相互作用是理解能量传输、加热以及诊断技术的核心。当电磁波或静电波在等离子体中传播时，由于等离子体自身具有色散特性且参数随空间分布不均匀，波在传播路径上不可避免地会发生反射、折射与吸收。掌握这三种物理过程的机制与建模方法，是进行准确等离子体模拟的前提。
波的反射发生在波遇到阻抗失配或到达截止区域时。在等离子体中，最典型的反射机制是由密度梯度引起的截止。
当波在等离子体中传播时，其色散关系通常可以表示为 $k^2 = \frac{\omega^2}{c^2} \varepsilon(\omega, x)$，其中 $\varepsilon$ 为等效介电常数。对于无碰撞冷等离子体，当波的频率 $\omega$ 等于局部等离子体频率 $\omega_{pe}$ 或上混合频率等特征频率时，等效介电常数趋近于零，波数 $k \to 0$。此时波的相速度趋于无穷大，群速度降为零，波无法继续向前传播，从而发生反射。这一空间位置被称为截止层。
在建模中，处理波的反射通常采用以下方法：
- WKB 近似：在缓变介质中，可以利用 WKB 近似追踪波包的轨迹。但在截止点附近 WKB 解失效，需要引入 Airy 函数进行连接处理。
- 全波模拟：对于波长与梯度尺度可比拟的情况，直接求解麦克斯韦方程组，数值格式本身能够自然地捕捉反射波与入射波的干涉图样。
波的折射与射线追踪
折射是指由于介质折射率的空间不均匀性导致波传播方向发生改变的现象。等离子体通常呈现边缘密度低、中心密度高的分布，这意味着边缘的折射率大，中心的折射率小，这种分布对波具有“散焦”作用。
在射频加热（如电子回旋共振加热 ECRH 或低杂波电流驱动 LHCD）的建模中，必须精确计算波的折射，以确保波束能够准确到达预定的共振层。此时，射线追踪是最常用的理论工具。
射线追踪基于几何光学近似，其核心方程组为：
- 哈密顿方程：$\frac{d\mathbf{r}}{dt} = \frac{\partial \omega}{\partial \mathbf{k}} = \mathbf{v}_g$
- 波矢演化方程：$\frac{d\mathbf{k}}{dt} = -\frac{\partial \omega}{\partial \mathbf{r}}$
通过在相空间 $(\mathbf{r}, \mathbf{k})$ 中积分上述方程，可以获得射线路径。在数值实现中，常采用四阶 Runge-Kutta 法进行求解。
示例：在托卡马克等离子体中，从低场侧发射的电子回旋波，由于极向磁场的存在和密度的梯度，射线在向内传播时会发生弯曲。如果不进行精确的射线追踪计算，波束可能会因为强烈的折射而偏离磁轴，导致加热效率大幅下降。
波的吸收与阻尼机制
波的吸收是等离子体加热和电流驱动的物理基础。波的能量通过不同的阻尼机制转化为等离子体粒子的动能。在理论建模中，吸收通常体现为波振幅随传播距离的指数衰减，即 $E \propto e^{-\gamma x}$，其中 $\gamma$ 为空间阻尼率。
主要的吸收机制包括：
- 碰撞吸收（欧姆吸收）：在低温或高密度等离子体边缘，电子与中性原子或离子的碰撞频率较高。波的电场驱动电子运动，电子通过碰撞将波的能量转化为热能。碰撞吸收率正比于电子-离子碰撞频率。
- 朗道阻尼：当波的相速度与带电粒子的热速度相近时，满足 $v_{ph} \approx v_{th}$ 的粒子会与波发生共振相互作用。速度略低于相速度的粒子被波加速吸收能量，速度略高的粒子减速释放能量。在麦克斯韦分布下，由于低速粒子数多于高速粒子数，净效应是波的能量被粒子吸收。这是离子回旋波（ICRH）和低杂波（LHW）加热的主要机制。
- 回旋阻尼：在磁化等离子体中，当波的频率接近离子的回旋频率（$\omega \approx \omega_{ci}$）或其谐波时，波的电场在粒子回旋轨道上的投影能够持续对粒子做正功，从而实现高效的共振吸收。这是离子回旋共振加热（ICRH）的核心机制。
在数值建模中，计算吸收率通常需要求解等离子体动理学方程（如 Vlasov 方程或 Fokker-Planck 方程），以获得包含速度分布函数的复数介电张量，进而通过其虚部提取阻尼系数。
综合建模与数值实现
在实际的等离子体物理建模中，反射、折射与吸收往往同时发生且相互耦合。例如，波在向高密度区传播时发生折射，到达截止层时发生反射，而在整个传播路径上都会经历碰撞或无碰撞阻尼。
构建完整的波传播模型需要遵循以下步骤：
- 确定色散关系：根据等离子体参数（密度、温度、磁场）选择冷等离子体或热等离子体模型，推导复数色散关系。
- 选择数值方法：
- 对于大尺度、长波情况，优先使用射线追踪法结合阻尼计算，计算效率高。
- 对于存在截止层、模式转换或强吸收区域的情况，必须采用全波 FDTD（时域有限差分）或 FEM（有限元）方法求解麦克斯韦方程组，并在方程中引入复数介电常数以同时涵盖反射与吸收。
- 能量守恒检验：在模拟结果中，入射波功率必须等于反射波功率、透射波功率与等离子体吸收功率之和，这是验证模型准确性的关键判据。
深刻理解波的反射、折射与吸收机制，并熟练运用相应的数学与数值工具，是开展等离子体波与介质相互作用研究的基石，也是优化聚变装置加热方案和工业等离子体源设计的关键所在。