本文探讨了等离子体模拟的两种主要方法——流体模拟和粒子模拟，并详细介绍了使用OpenFOAM进行Particle In Cell（PIC）方法的尝试，以处理复杂结构下的等离子体行为。针对非结构网格的挑战，文中提到了有限元法粒子模拟（FEM-PIC）和OpenFOAM的潜力，同时分享了在COMSOL和OpenFOAM中实现PIC方法遇到的问题和经验。
等离子体模拟主要有两种方法，流体模拟和粒子模拟。
等离子体的流体模拟其实和CFD通用的数值方法一样，主要是有限差分方法、有限元方法和有限体积方法。有限元商用软件comsol，有专门的等离子体模块，采用流体模拟，易用性非常不错。开源的低温等离子体流体模拟程序好像非常少，磁约束核聚变领域（即高温等离子体）的流体模拟程序源代码基本都公开，比如SOLPS、BOUT++等。因为磁约束聚变装置，比如托卡马克和仿星器，磁场结构特殊，模拟程序通常采用特殊的与磁场有关的坐标系，因此专业性强，通用性差，基本只流行与聚变领域。
粒子模拟，英文 应该是Particle Simulation，包括PIC、MCC、DSMC等。PIC（Particle In Cell）模拟粒子在电磁场中的运动和电磁场的随时间变化，即带电粒子和电磁场存在相互作用；MCC，蒙特卡洛碰撞，模拟带电粒子之间，以及与中性粒子的碰撞，我觉得库伦碰撞的二体碰撞方法也可以归为MCC；DSMC，直接蒙特卡洛模拟，模拟中性原子分子之间的碰撞，这个貌似通常是一个单独的方向，但是可以和等离子体结合模拟。我接触的粒子模拟程序基本都是均匀的结构化网格，商业粒子程序做的比较好的是VSIM，PEGASUS差一些。VSIM也是均匀网格，采用cut-cell技术模拟不规则边界。我基于法国一个团队开发的Smilei程序（Home - Smilei v4.2 documentation）开发了模拟静电问题的粒子程序，采用类似cut-cell的技术处理不规则边界，求解泊松方程涉及求解线性方程组，我采用SuperLU、PETSc等库。采用均匀网格的好处是&#x
最低0.47元/天 解锁文章
扫一扫
3020
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
