乔恩斯（Michael E. Jones）、莱蒙斯（Don S. Lemons）、梅森（Rodney J. Mason）、托马斯（Vincent A. Thomas）和温斯克（Dan Winske）在《Journal of Computational Physics》期刊第123卷（1996年）上发表了一项名为“用于粒子云网格（PIC）代码的基于网格的库仑碰撞模型”的研究。该研究由洛斯阿拉莫斯国家实验室（Los Alamos National Laboratory）应用理论物理部完成，提出了一种在粒子云网格代码中模拟从无碰撞到库仑碰撞主导等离子体这一中间状态的新方法。

这项研究的学术背景根植于等离子体模拟领域长达数十年的方法论分化。在早期的等离子体数值模拟中，发展出了两种根本不同的方法：一种是将等离子体视为流体的流体动力学或磁流体动力学方法，该方法隐含的前提是碰撞过程保持速度分布接近麦克斯韦分布（Maxwellian），因此只需考虑玻尔兹曼方程（Boltzmann equation）的最低阶矩；另一种是动力学方法，通常使用粒子云网格代码来实现，该方法关注等离子体组分（电子和离子）的个体行为，适用于碰撞过程缺失、非麦克斯韦分布特征显著的无碰撞状态。然而，随着研究的发展，诸多新兴应用领域，例如激光产生的等离子体的相互穿透（interpenetration of laser produced plasmas）、地球等离子体层（plasmasphere）的再填充以及弱不稳定性导致的非均匀磁化等离子体中的横越磁场输运等，都迫切需要研究等离子体处于“半碰撞”（semi-collisional）或碰撞状态在时空上快速变化的区域。在这个中间区域，传统的流体或动理学方法都不完全适用。尽管当时已有混合算法（hybrid algorithms）或福克-普朗克方法（Fokker–Planck methods）等尝试弥合这一鸿沟，但无碰撞与碰撞等离子体相互作用之间的界限总体上仍未打破。正是为了填补这一空白，本研究旨在开发一种能够有效模拟多物种等离子体在广袤参数范围内相互作用的高效碰撞模型。

该研究的核心工作流程是开发并验证一套全新的碰撞算法，该算法主要分为两个部分来处理不同情况下的碰撞问题。

第一部分是用于处理不同粒子物种之间碰撞的“碰撞场”（collision field）方法。不同于传统方法中逐对处理粒子碰撞（binary Coulomb collisions），该方法通过定义在网格上的一个场量来介导短程碰撞力。这种方法的核心要求是局部（local）动量与能量守恒。研究人员通过分布函数平均，要求物种a受到物种b的碰撞力$F{ab}$满足两个条件：$kF{ab}l = -kF_{ba}l$（动量守恒）和$kva \cdot F{ab}l = -kvb \cdot F{ba}l$（能量守恒）。为满足这些条件，他们推导得到了一个碰撞力的具体函数形式，该力由流体量的网格插值决定，并引入了两个关键的碰撞频率：一个是动态摩擦（dynamic friction）频率$\nu{ab}$，另一个是与温度平衡（temperature equilibration）相关的频率$\nu^\epsilon{ab}$。这两个频率均基于经典斯必泽（Spitzer）值推导而来，能够准确反映粒子分布函数向麦克斯韦分布弛豫的速率。当粒子运动时，会将此碰撞力与洛伦兹力（Lorentz force）一同计入速度更新中。研究人员通过两个一维和二维的测试问题来验证该模型。第一个测试是让两个不同密度和初始温度、且相对静止的离子群发生热弛豫（thermal relaxation）。模拟结果显示，在$\nu^\epsilon{ab}t = 9.8$的时间点，两个分布的温度成功达到平衡，其弛豫的时间历程与分析得到的流体解高度吻合。第二个测试是让两个温度相同但存在相对漂移的离子群相互作用。模拟显示，漂移离子群在减速并最终与静止离子群达到相同速度的过程中，其分布被加热且呈现出麦克斯韦形状，最终在$\nu{ab}t = 0.48$时弛豫到相同速度，其速度时间历程也与流体计算结果一致，但偏离了动理学解。

第二部分是用于处理同一物种内部粒子间碰撞的“蒙特卡洛”式散射方法（Monte Carlo-like scattering method）。因为当(a = b)时，上述碰撞场力(F{ab})为零，模型采用了基于朗之万方程（Langevin equation）的广义形式。该力形式为(F{aa}/m_a = -\nu_a(v_a - kv_al) + a)，其中(\nu_a)是动态摩擦系数，(a)是一个为提供热化而选择的各向同性随机矢量（random, isotropic vector）。此形式在统计意义上确保动量守恒。通过将速度推进方程改写为中心有限差分格式并恰当选择随机矢量的分布函数，该算法能够保证对于任何时间步长都精确保持能量守恒。其中的动态摩擦系数(\nu_a)是通过取两个漂移麦克斯韦分布碰撞频率在相对漂移趋于零时的极限得到。该部分同样通过测试问题验证。第一个问题模拟一个初始为方形速度分布的弛豫过程，模拟显示分布迅速演变为球状并趋向麦克斯韦分布，研究人员通过计算熵函数（entropy function）来衡量弛豫速率，发现其极快便达到平衡值。第二个问题是球形诺尔（Noh）问题，模拟高密度流体撞击静止壁面并产生向外传播的激波。模拟得到的密度和温度剖面与解析解吻合良好，成功再现了正确的激波速度和密度、温度跳变幅度。

在完成两个核心算法的开发与测试后，研究人员将其整合进一个全新的多维混合代码中（该代码将离子视为粒子、电子视为无质量流体，并包含辐射效应，其基本方程在附录中给出），并对一个更复杂的现实问题——激光产生的等离子体的相互穿透进行了模拟，以展示模型的实用性。在一维计算中，两块薄金等离子体片被加热并相向膨胀，研究人员同时计算了短程的种间和种内库仑碰撞以及通过自洽电场（self-consistent electric field）作用的长程碰撞。模拟清晰地重现了整个物理过程：初始阶段等离子体等温膨胀；随后在相互作用区域，两个离子群相互穿透、加热并相对减速；最终在中心区域达到热化，且冷却由电子完成。电子密度剖面显示了预期的中心密度堆积。通过改变初始电子密度，模型还能展示从弱耦合相互穿透到强耦合并形成向外传播激波的不同物理状态。该模型还成功展示了多材料相互作用的能力，例如模拟带有CH涂层的金箔膨胀与碰撞，并显现出轻离子（如氢离子）因电荷数（Z）小而快速膨胀和弱耦合的特点。最后，研究人员还进行了二维的碰撞等离子体计算，展示了等离子体在轴向和径向同时膨胀并发生碰撞的电子密度等值线图，证明了碰撞场方法在多维模拟中比传统的粒子配对（pairing up the particles）方法效率更高。

这项研究的主要结论是，该研究所提出的新型碰撞模型为粒子代码中处理等离子体物种间的短程库仑碰撞提供了一种有效途径。该模型的核心在于，种间碰撞的碰撞场方法能够在没有大量粒子数的条件下极好地局部守恒能量和动量，而种内碰撞的朗之万方法亦能确保相关的输运系数与经典斯必泽过程一致。该方法通过结合混合代码，证明了其能够处理以前难以企及的“半碰撞”等离子体区域，成功弥合了无碰撞和强碰撞状态的模拟方法之间的鸿沟。

该研究的亮点在于其创新性、高效性和适用性。首先，它创造性地提出了“碰撞场”这一基于网格的碰撞算子概念，完全不同于传统的粒子配对散射方法，从根本上解决了动量和能量精确守恒的难题。其次，该方法计算效率极高，仅比现有PIC方法增加约50%的计算负担，却允许使用远大于碰撞频率倒数的大时间步长，大大节省了计算成本。最后，该模型被证明非常适用于多维计算，并且其开放的框架可以方便地扩展以包含更复杂的物理过程，如辐射输运和可变电离态等，为后续在激光物质相互作用、空间物理等多个领域的应用奠定了坚实的基础。