首页|一类含时Poisson-Nernst-Planck方程的虚单元计算

一类含时Poisson-Nernst-Planck方程的虚单元计算

扫码查看
针对一类含时Poisson-Nernst-Planck(PNP)方程,为避免在解决实际问题时有限元法中的网格适应性问题,构造了L2投影算子与Gummel迭代相结合的虚单元算法.该算法允许以更简单的方式设计和分析新的格式,可以灵活处理各种网格,对于多边形或多面体单元甚至非凸单元组成的网格剖分都可以很好地处理,使得虚单元法可以适应于任意多边形网格,大大降低了网格的生成难度.给出了虚单元算法在三角形网格、四边形网格、非凸网格下的数值算例.数值实验结果表明,在这3种多边形网格上,L2和H1模的收敛阶分别为二阶和一阶,均达到了最优阶.
A virtual element methods for a class of time-dependent Poisson-Nernst-Planck equations
For a class of time-dependent Poisson-Nernst-Planck(PNP)equations,in order to avoid the mesh adaptability problem in finite element method when solving practical problems,a virtual element algorithm combining L2 projection operator and Gummel iteration was proposed.This method allows new formats to be designed and analyzed in a simpler manner and flexible handling of various meshes.Virtual element method can handle polygon or polyhedron elements well,even the mesh division consisting of non-convex elements,which makes virtual element method can adapt to arbitrary polygon mesh,greatly reducing the difficulty of mesh generation.Numerical examples of virtual element algorithm under triangular mesh,quadrilateral mesh,non-convex mesh are given.Numerical experiments show that on three different grids,the convergence order of the L2 and H norms are second order and first order respectively,and them both reach the optimal order.

Poisson-Nernst-Planck equationvirtual element methodL2 projectionGummel iterationL2 normsH1 norms

刘亚、阳莺

展开 >

桂林电子科技大学数学与计算科学学院,广西桂林 541004

桂林电子科技大学广西应用数学中心,广西桂林 541004

Poisson-Nernst-Planck方程 虚单元算法 L2投影 Gummel迭代 L2模 H1模

广西自然科学基金国家自然科学基金广西科技项目桂林电子科技大学研究生优秀学位论文培育项目

2020GXNSFAA15909812161026桂科230230022020YJSPYA02

2024

桂林电子科技大学学报
桂林电子科技大学

桂林电子科技大学学报

影响因子:0.247
ISSN:1673-808X
年,卷(期):2024.44(1)