Hysteresis and Linear Stability Analysis on Multiple Steady-State Solutions to the Poisson–Nernst–Planck equations with Steric Interactions

by   Jie Ding, et al.

In this work, we numerically study linear stability of multiple steady-state solutions to a type of steric Poisson–Nernst–Planck (PNP) equations with Dirichlet boundary conditions, which are applicable to ion channels. With numerically found multiple steady-state solutions, we obtain S-shaped current-voltage and current-concentration curves, showing hysteretic response of ion conductance to voltages and boundary concentrations with memory effects. Boundary value problems are proposed to locate bifurcation points and predict the local bifurcation diagram near bifurcation points on the S-shaped curves. Numerical approaches for linear stability analysis are developed to understand the stability of the steady-state solutions that are only numerically available. Finite difference schemes are proposed to solve a derived eigenvalue problem involving differential operators. The linear stability analysis reveals that the S-shaped curves have two linearly stable branches of different conductance levels and one linearly unstable intermediate branch, exhibiting classical bistable hysteresis. As predicted in the linear stability analysis, transition dynamics, from a steady-state solution on the unstable branch to a one on the stable branches, are led by perturbations associated to the mode of the dominant eigenvalue. Further numerical tests demonstrate that the finite difference schemes proposed in the linear stability analysis are second-order accurate. Numerical approaches developed in this work can be applied to study linear stability of a class of time-dependent problems around their steady-state solutions that are computed numerically.



There are no comments yet.


page 1

page 2

page 3

page 4


Bifurcation analysis of two-dimensional Rayleigh–Bénard convection using deflation

We perform a bifurcation analysis of the steady state solutions of Rayle...

Estimating Dispersion Curves from Frequency Response Functions via Vector-Fitting

Driven by the need for describing and understanding wave propagation in ...

Stability and dynamical transition of a electrically conducting rotating fluid

In this article, we aim to study the stability and dynamic transition of...

On the relation of powerflow and Telegrapher's equations: continuous and numerical Lyapunov stability

In this contribution we analyze the exponential stability of power netwo...

On Surrogate Learning for Linear Stability Assessment of Navier-Stokes Equations with Stochastic Viscosity

We study linear stability of solutions to the Navier–Stokes equations wi...

First-order continuation method for steady-state variably saturated groundwater flow modeling

Recently, the nonlinearity continuation method has been used to numerica...

Bifurcation analysis of stationary solutions of two-dimensional coupled Gross-Pitaevskii equations using deflated continuation

Recently, a novel bifurcation technique known as the deflated continuati...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.