Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPoisson-Nernst-Planck Systems for Narrow Tubular-like Membrane Channels
202
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We study global dynamics of the Poisson-Nernst-Planck (PNP) system for flows of two types of ions through a narrow tubular-like membrane channel. As the radius of the cross-section of the three-dimensional tubular-like membrane channel approaches zero, a one-dimensional limiting PNP system is derived. This one-dimensional limiting system differs from previous studied one-dimensional PNP systems in that it encodes the defining geometry of the three-dimensional membrane channel. To justify this limiting process, we show that the global attractors of the three-dimensional PNP systems are upper semi-continuous to that of the limiting PNP system. We then examine the dynamics of the one-dimensional limiting PNP system. For large Debye number, the steady-state of the one-dimensional limiting PNP system is completed analyzed using the geometric singular perturbation theory. For a special case, an entropy-type Lyapunov functional is constructed to show the global, asymptotic stability of the steady-state.
rate research
Read More
In this paper, we design and analyze third order positivity-preserving discontinuous Galerkin (DG) schemes for solving the time-dependent system of Poisson--Nernst--Planck (PNP) equations, which has found much use in diverse applications. Our DG method with Euler forward time discretization is shown to preserve the positivity of cell averages at all time steps. The positivity of numerical solutions is then restored by a scaling limiter in reference to positive weighted cell averages. The method is also shown to preserve steady states. Numerical examples are presented to demonstrate the third order accuracy and illustrate the positivity-preserving property in both one and two dimensions.
In this paper, we develop an adaptive finite element method for the nonlinear steady-state Poisson-Nernst-Planck equations, where the spatial adaptivity for geometrical singularities and boundary layer effects are mainly considered. As a key contribution, the steady-state Poisson-Nernst-Planck equations are studied systematically and rigorous analysis for a residual-based a posteriori error estimate of the nonlinear system is presented. With the help of Schauder fixed point theorem, we show the solution existence and uniqueness of the linearized system derived by taking $G-$derivatives of the nonlinear system, followed by the proof of the relationship between the error of solution and the a posteriori error estimator $eta$. Numerical experiments are given to validate the efficiency of the a posteriori error estimator and demonstrate the expected rate of convergence. In the further tests, adaptive mesh refinements for geometrical singularities and boundary layer effects are successfully observed.
The Poisson-Nernst-Planck equations with generalized Frumkin-Butler-Volmer boundary conditions (PNP-FBV) describe ion transport with Faradaic reactions, and have applications in a number of fields. In this article, we develop an adaptive time-stepping scheme for the solution of the PNP-FBV equations based on two time-stepping methods: a fully implicit (BDF2) method, and an implicit-explicit (SBDF2) method. We present simulations under both current and voltage boundary conditions and demonstrate the ability to simulate a large range of parameters, including any value of the singular perturbation parameter $epsilon$. When the underlying dynamics is one that would have the solutions converge to a steady-state solution, we observe that the adaptive time-stepper based on the SBDF2 method produces solutions that ``nearly converge to the steady state and that, simultaneously, the time-step sizes stabilize to a limiting size $dt_infty$. In the companion to this article cite{YPD_Part2}, we linearize the SBDF2 scheme about the steady-state solution and demonstrate that the linearized scheme is conditionally stable. This conditional stability is the cause of the adaptive time-steppers behaviour. While the adaptive time-stepper based on the fully-implicit (BDF2) method is not subject to such time-step constraints, the required nonlinear solve yields run times that are significantly longer.
We develop a modified Poisson-Nernst-Planck model which includes both the long-range Coulomb and short-range hard-sphere correlations in its free energy functional such that the model can accurately describe the ion transport in complex environment and under a nanoscale confinement. The Coulomb correlation including the dielectric polarization is treated by solving a generalized Debye-Huckel equation which is a Greens function equation with the correlation energy of a test ion described by the self Greens function. The hard-sphere correlation is modeled through the modified fundamental measure theory. The resulting model is available for problems beyond the mean-field theory such as problems with variable dielectric media, multivalent ions, and strong surface charge density. We solve the generalized Debye-Huckel equation by the Wentzel-Kramers-Brillouin approximation, and study the electrolytes between two parallel dielectric surfaces. In comparison to other modified models, the new model is shown more accurate in agreement with particle-based simulations and capturing the physical properties of ionic structures near interfaces.
Ion transport in biological and synthetic nanochannels is characterized by phenomena such as ion current fluctuations and rectification. Recently, it has been demonstrated that nanofabricated synthetic pores can mimic transport properties of biological ion channels [P. Yu. Apel, {it et al.}, Nucl. Instr. Meth. B {bf 184}, 337 (2001); Z. Siwy, {it et al.}, Europhys. Lett. {bf 60}, 349 (2002)]. Here, the ion current rectification is studied within a reduced 1D Poisson-Nernst-Planck (PNP) model of synthetic nanopores. A conical channel of a few $mathrm{nm}$ to a few hundred of nm in diameter, and of few $mu$m long is considered in the limit where the channel length considerably exceeds the Debye screening length. The rigid channel wall is assumed to be weakly charged. A one-dimensional reduction of the three-dimensional problem in terms of corresponding entropic effects is put forward. The ion transport is described by the non-equilibrium steady-state solution of the 1D Poisson-Nernst-Planck system within a singular perturbation treatment. An analytic formula for the approximate rectification current in the lowest order perturbation theory is derived. A detailed comparison between numerical results and the singular perturbation theory is presented. The crucial importance of the asymmetry in the potential jumps at the pore ends on the rectification effect is demonstrated. This so constructed 1D theory is shown to describe well the experimental data in the regime of small-to-moderate electric currents.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
