基于Ritz与Kim方法的发育等离子体湍流相干模间非线性能量转移：IMPED模拟与实验验证

《Journal of Plasma Physics》：Nonlinear energy transfer analysis between coherent modes in developing plasma turbulence

　　

等离子体中波动物理参数的各个谱分量之间的能量转移由非线性相互作用引起，但传统线性谱方法通常无法捕捉这些效应。在反镜等离子体实验装置（Inverse Mirror Plasma Experimental Device, IMPED）中测量的等离子体密度涨落具有多种不稳定性之间非线性模相互作用的特征，例如瑞利–泰勒（Rayleigh–Taylor, RT）模和漂移波（drift-wave, DW）模。本文详细研究了这些模之间由非线性波相互作用（通过二次耦合过程）导致的能量转移。研究人员探索了单场模型中已有的计算方法，如Ritz方法和Kim方法，以理解相干模（coherent modes）之间非线性相互作用的动力学。两种方法均被应用，并在模拟以及来自IMPED的等离子体相干模实验数据中得到验证。研究人员发现，这些方法的有效性和适用性取决于数据的统计性质，特别是高阶矩如峰度（kurtosis），并取决于空间平稳性。使用这些方法进行的能量转移分析主要揭示了能量从RT模以及DW与RT的混合模向相对低频DW模的转移，证明了该方法能够量化相干模之间的谱能量输运。

研究背景。等离子体湍流（plasma turbulence）通常从小振幅涨落在线性不稳定性驱动下增长开始，随后非线性模耦合使谱展宽，形成向耗散尺度级联，并最终在增长、非线性转移和耗散之间趋于统计稳态。传统线性谱方法基于二阶统计，如自功率谱和交叉功率谱，假设模之间统计独立，不能捕捉非线性动力学产生的相位相干相互作用，因此对谱分量间能量交换的认识有限。三阶统计的双谱（bispectrum）和双相干（bicoherence）可识别相位相关波三元组，但主要量化耦合存在性，不直接给出能量流方向。Ritz & Powers等提出的功率传递函数和能量传递函数可量化净非线性能量交换，但Ritz方法采用Millionshchikov近似（Millionshchikov approximation），将四阶矩用二阶矩平方代替，对非高斯和高峰度数据产生误差，并在完全发展湍流中产生非物理阻尼。Kim等提出的方法保留四阶矩，将数据分解为理想和非理想分量，并施加空间平稳闭合，可用于高斯和非高斯数据，但要求输入输出谱趋势相似。物理上，三波耦合主要来自漂移波（drift-wave, DW）的E×B二次对流非线性，RT模则与有效重力、磁曲率或离心力作用于不利密度梯度有关。IMPED密度涨落已有RT和DW等模非线性相互作用特征。因此，研究人员在模拟和IMPED实验数据中系统比较Ritz方法和Kim方法，重点考察峰度等统计性质和空间平稳性对方法适用性的影响，并量化相干模间能量转移。

研究人员开展了单场模型的非线性能量转移分析，并用非线性黑箱模型进行模拟基准测试；实验数据来自IMPED装置，在350 G和4×10^-4 mbar条件下，对两个径向位置的密度涨落进行分析。结论是Ritz方法仅对近高斯、低峰度数据可靠，Kim方法对非高斯数据更稳健；实验显示能量主要从RT模以及DW与RT混合模向相对低频DW模转移。该工作的意义在于明确了两类能量传递函数方法的适用条件，并证明可用其量化发育等离子体湍流中相干模间的谱能量输运。

关键技术方法。研究人员在单场非线性波耦合框架下采用功率谱、交叉功率谱、自/交叉双谱和双相干表征二次相位耦合；以Ritz方法和Kim方法估计线性传递函数L_f和二次传递函数Q_f^f₁,f₂，进而计算功率传递函数T_f(f₁, f₂)与能量传递函数W_NL^f，并用W_mis和Δ_im评估能量守恒与空间平稳性。模拟中构建非线性黑箱，将零均值白高斯噪声串联通过5个或10个黑箱，取迭代5/10的输入输出对，使用M=2500、N=64，并用M>3500检查收敛。实验数据来自IMPED装置，使用三探针Langmuir探针测量密度涨落?和浮动电位涨落φ?_f，采样1 MHz、时长10 s、约10⁷个样本；两探针极向间距Δx=6.4 mm，经交叉相关函数（cross-correlation function, CCF）确定输入输出方向，在r=2.24 cm和r=5.76 cm分析，取N=2048、M=4882，频率分辨率约244 Hz。

研究结果。4.1 流动传播方向估计：通过CCF分析，在r=5.76 cm处看到正时间延迟峰，确认波从X传播到Y，从而确定输入和输出信号。4.2 实验数据选择：依据偏度、峰度和空间平稳性选择两个径向位置；r=2.24 cm处分布近高斯，偏度0.042、超额峰度0.093；r=5.76 cm处偏离高斯，偏度0.205、超额峰度0.961，并呈多峰分布。4.3 分析和结果如下。

4.3.1 径向位置r=5.76 cm的非线性能量转移分析：自功率谱和自双相干显示三个主模，2.5 kHz为DW模，8.3 kHz为DW与RT混合模，10.8 kHz为RT模；8.3 kHz与2.5 kHz耦合到10.8 kHz，满足f₁+f₂=f，并满足k_θ,1+k_θ,2=k_θ；21 kHz为10.8 kHz的二次谐波，与8–12.5 kHz频带强耦合。

4.3.1.1 使用Kim方法分析：功率传递函数T_f(f₁, f₂)和能量传递函数W_NL^f显示，10.8 kHz RT模与8.3 kHz混合模通过(10.8 kHz)+(-8.3 kHz)→2.5 kHz将能量传给2.5 kHz DW模；2.5 kHz处为正峰，8.3 kHz和10.8 kHz处为负峰。10.8 kHz还向21 kHz、10.2 kHz和31.8 kHz转移能量；10.8 kHz与10.2 kHz耦合到21 kHz；21 kHz向8.5–10.6 kHz带转移能量。W_mis=0.456。结论是Kim方法在该位置可靠，RT模和混合模向低频DW模及其他模转移能量。

4.3.1.2 使用Ritz方法分析：同一数据给出与Kim方法不同甚至相反的趋势，如2.5 kHz+8.3 kHz→10.8 kHz；12.7 kHz→21 kHz+(-8.3 kHz)；10.8 kHz→31.8 kHz+(-21 kHz)。W_mis=-0.199，结果更嘈杂。由于该位置超额峰度约0.96，超过模拟给出的约0.40阈值，Ritz方法失效，结果不可物理采信。

4.3.2 径向位置r=2.24 cm的非线性能量转移分析：自双相干显示2.5 kHz为DW模，8.8 kHz为混合模，11.3 kHz为RT模；8.8 kHz与2.5 kHz耦合到11.3 kHz，另有较弱的6.3 kHz与2.5 kHz耦合到8.8 kHz；极向波数匹配不严格，与较弱耦合及波数估计不确定性有关。

4.3.2.1 使用Kim方法分析：11.3 kHz RT模通过11.3 kHz→8.8 kHz+2.5 kHz向8.8 kHz混合模和2.5 kHz DW模转移能量；6.3 kHz+2.5 kHz→8.8 kHz使8.8 kHz增益；8.8 kHz+(-6.3 kHz)→2.5 kHz使2.5 kHz增益。W_NL^f在11.3 kHz为负峰，在8.8 kHz和2.5 kHz为正峰，净能量从11.3 kHz RT模流向8.8 kHz混合模和2.5 kHz DW模。W_mis=-0.18。

4.3.2.2 使用Ritz方法分析：Ritz方法结果与Kim方法定性一致，主导耦合相同，W_NL^f趋势相同，同样显示11.3 kHz RT模向8.8 kHz混合模和2.5 kHz DW模转移能量。W_mis=-0.06。该位置超额峰度约0.09，低于约0.40阈值，因此Ritz方法可靠，并与Kim方法一致。

4.3.3 增长率的讨论：在r=2.24 cm，原文报告约10.8 kHz模具有小正增长率，因而不稳定；约2.5 kHz和8.3 kHz模强阻尼。10.8 kHz模线性增长获得能量，并通过非线性相互作用把能量再分配给其他模或背景平衡；较低频模通过非线性相互作用接收能量。高频带正增长率表明这些模不稳定，其能量输运到低频带或背景平衡。该位置Δ_im=0.06，虽小但表明系统仍处于发育湍流状态，Kim方法仍有效。

讨论部分总结：研究人员将Ritz方法和Kim方法先在发育等离子体湍流模型的解析传递函数上做基准测试，再用于IMPED实验密度涨落。Ritz方法对近高斯、低偏度和低峰度涨落表现较好；峰度较大时，因四阶矩近似误差，传递函数估计恶化。Kim方法不采用Millionshchikov近似，保留四阶矩，因而对高峰度、非高斯数据更稳健。空间平稳性方面，Kim方法在发育湍流、输入输出谱不完全重合但谱形趋势相似时仍可给出可靠估计；Ritz方法仅在r=2.24 cm可靠，在r=5.76 cm不可靠。W_mis用于评估残余能量不平衡，但原文指出其值会受有限频带截断、未测场转移、向平均剖面转移和统计不确定度影响，不宜作为方法适用性的唯一判据。

研究结论（翻译）：研究人员使用Ritz和Kim方法研究了IMPED等离子体中DW和RT模之间的非线性能量转移。针对发育等离子体湍流模型的基准测试表明，Ritz方法对近高斯涨落提供可靠估计，而随着峰度增加，其性能恶化；相反，Kim方法对高斯和非高斯涨落均保持稳健。实验数据中观察到一致趋势：Kim方法在r=5.76 cm和r=2.24 cm均给出可靠估计，而Ritz方法仅在r=2.24 cm可靠。这些结果强调涨落统计性质，尤其是超额峰度，在评估能量转移估计方法适用性中的重要性。能量转移分析进一步揭示，在r=5.76 cm，10.8 kHz的RT模和8.3 kHz混合模主要向2.5 kHz的DW模转移能量；相反，在r=2.24 cm，11.3 kHz的RT模向2.5 kHz的DW模和8.8 kHz混合模转移能量。因此，分析表明在两个径向位置均存在从RT主导涨落向相对较低频DW模的净能量转移。然而，由于r=5.76 cm处非线性相互作用显著更强，整体能量转移方向主要从该径向位置结果推断。因此，主导非线性能量转移过程可解释为能量从RT模及DW–RT混合模向相对低频DW模的转移。展望部分指出，当前分析限于单场框架，未来需扩展到多场框架以描述密度、电位、温度和速度等多个耦合场之间的非线性能量转移。