为了对比研究有限元法(FEM)和有限体积法(FVM)2种算法以及网格结构对大气压氦气介质阻挡柱状放电仿真结果的影响,建立了2维流体模型。首先,计算了柱状放电的电流波形和电子数密度2维分布;然后,根据收敛曲线、迭代次数和运算时间对比研究了相同网格结构下FEM和FVM的计算性能;最后,根据2维电子数密度分布云图和试验结果对比研究了2种算法下的放电形态、时序的差异。结果表明:在不稳定的初始阶段,放电由均匀向柱状转化,最后达到稳定的柱状放电阶段;FVM计算方法较稳定且耗时少;网格结构变化会影响FVM计算得到的放电形态,发现方形网格结构下FVM计算获得的柱状放电与试验结果相符。
In order to investigate the influences of the finite element method(FEM), finite volume method(FVM), and the grid structures on the simulation of dielectric barrier columnar discharges in atmospheric helium, we established a twodimensional(2D) fluid model. Firstly, we calculated the current waveforms and 2D electron density distributions of the discharge; secondly, we comparatively studied the performances of the FEM and the FVM based on convergence plots, number of iterations, and computation time; finally, we comparatively studied discharge morphology and sequence obtained by the two algorithms based on the 2D electron density distributions and the experimental results. The results show that the discharge evolves from a uniform pattern into a columnar one in the unstable initial stage, and it reaches a steady columnar discharge at last. The FVM is more stable and less time consuming than the FEM. And the variation of grid structure affects the discharge morphology obtained by the FVM. Under the square grid structure, the columnar discharge obtained by the FVM coincides with the experimental results.