In this work, we introduce an inverse averaging finite element method (IAFEM) for solving the size-modified Poisson-Nernst-Planck (SMPNP) equations. Comparing with the classical Poisson-Nernst-Planck (PNP) equations, the SMPNP equations add a nonlinear term to each of the Nernst-Planck (NP) fluxes to describe the steric repulsion which can treat multiple nonuniform particle sizes in simulations. Since the new terms include sums and gradients of ion concentrations, the nonlinear coupling of SMPNP equations is much stronger than that of PNP equations. By introducing a generalized Slotboom transform, each of the size-modified NP equation is transformed into a self-adjoint equation with exponentially behaved coefficient, which has similar simple form to the standard NP equation with the Slotboom transformation. This treatment enables employing our recently developed inverse averaging technique to deal with the exponential coefficients of the reformulated formulations, featured with advantages of numerical stability and flux conservation especially in strong nonlinear and convection-dominated cases. Comparing with previous stabilization methods, the IAFEM proposed in this paper can still possess the numerical stability when dealing with convection-dominated problems. And it is more concise and easier to be numerically implemented. Numerical experiments about a model problem with analytic solutions are presented to verify the accuracy and order of IAFEM for SMPNP equations. Studies about the size-effects of a sphere model and an ion channel system are presented to show that our IAFEM is more effective and robust than the traditional finite element method (FEM) when solving SMPNP equations in simulations of biological systems.
翻译:在这项工作中,我们引入了一种反平均限量元素法(IAFEM),用于解决大小修改后的 Poisson-Nernst-Planck(SMPNP)方程式。与古典Poisson-Nernst-Planck(PNP)方程式相比,SMPNP方程式在Nernst-Planck(NPNP)每种通量中都增加了一个非线性词,用来描述在模拟中处理多种非单式粒子大小。由于新术语包括了离子浓度的数值和梯度,SMPNPP方程式的非线性联变异性比PNPNP的等式要强得多。通过引入通用的Slotboopmus(SMPNPP)方程式变异性变异性变异性(SMPMNP)等式模型可以使用我们最近开发的逆性平均技术来处理重新配制的指数系数,其特点是数字稳定性的优点和变异性变式变异性变异性变异性变异性变法,在坚固的IMFIMA的系统中,这个变型变变的SMA的变型变型变变变变变型中, 和变型的SMA变型的系统则则将显示的变的变的变现的变式系统将比前的变的变的变的变的变型号都的变的变的变型法, 的变型法是更易性变式的变式的变型式的变型式的变型式的变型式, 。