Search

The transport of charged particles, which can be described by the Maxwell-Ampere Nernst-Planck (MANP) framework, is essential in various applications including ion channels and semiconductors. We propose a decoupled structure-preserving numerical scheme for the MANP model in this work. The Nernst-Planck equations are treated by the implicit exponential time differencing method associated with the Slotboom transform to preserve the positivity of the concentrations. In order to be effective with the Fast Fourier Transform, additional diffusive terms are introduced into Nernst-Planck equations. Meanwhile, the correction is introduced in the Maxwell-Ampere equation to fulfill Gauss's law. The curl-free condition for electric displacement is realized by a local curl-free relaxation algorithm whose complexity is O(N). We present sufficient restrictions on the time and spatial steps to satisfy the positivity and energy dissipation law at a discrete level. Numerical experiments are conducted to validate the expected numerical accuracy and demonstrate the structure-preserving properties of the proposed method., Comment: This paper has been accepted by Journal of Scientific Computing

In this paper we propose and analyze a finite difference numerical scheme for the Flory-Huggins-Cahn-Hilliard equation with dynamical boundary condition. The singular logarithmic potential is included in the Flory-Huggins energy expansion. Meanwhile, a dynamical evolution equation for the boundary profile corresponds to a lower-dimensional singular energy potential. In turn, a theoretical analysis for the coupled system becomes very challenging, since it contains nonlinear and singular energy potentials for both the interior region and on the boundary. In the numerical design, a convex splitting approach is applied to the chemical potential associated with the energy both at the interior region and on the boundary: implicit treatments for the singular and logarithmic terms, as well as the surface diffusion terms, combined with an explicit treatment for the concave expansive term. In addition, the discrete boundary condition for the phase variable is coupled with the evolutionary equation of the boundary profile. The resulting numerical system turns out to be highly nonlinear, singular and coupled. A careful finite difference approximation and convexity analysis reveals that such a numerical system could be represented as a minimization of a discrete numerical energy functional, which contains both the interior and boundary integrals. More importantly, all the singular terms correspond to a discrete convex functional. As a result, a unique solvability and positivity-preserving analysis could be theoretically justified, based on the subtle fact that the singular nature of the logarithmic terms around the singular limit values prevent the numerical solutions reaching these values. The total energy stability analysis could be established by a careful estimate over the finite difference inner product. Some numerical results are presented in this article., Comment: This paper has been submitted to Journal of Scientific Computing

A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment for the singular logarithmic terms and the surface diffusion term, and an explicit update for the expansive concave term. The convective term for the phase variable, as well as the coupled term in the Stokes equation, are approximated in a semi-implicit manner. In the spatial discretization, the marker and cell (MAC) difference method is applied, which evaluates the velocity components, the pressure and the phase variable at different cell locations. Such an approach ensures the divergence-free feature of the discrete velocity, and this property plays an important role in the analysis. The positivity-preserving property and the unique solvability of the proposed numerical scheme are theoretically justified, utilizing the singular nature of the logarithmic term as the phase variable approaches the singular limit values. An unconditional energy stability analysis is standard, as an outcome of the convex-concave decomposition technique. A convergence analysis with accompanying error estimate is provided for the proposed numerical scheme. In particular, a higher order consistency analysis, accomplished by supplementary functions, is performed to ensure the separation properties of numerical solution. In turn, using the approach of rough and refined error (RRE) estimates, we are able to derive an optimal rate convergence. To conclude, several numerical experiments are presented to validate the theoretical analysis., Comment: This article has been submitted to Mathematics of Computation