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

09/17/2020
by   Jie Ding, et al.
0

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.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/21/2021

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

We perform a bifurcation analysis of the steady state solutions of Rayle...
research
05/13/2023

Bifurcation analysis of waning-boosting epidemiological models with repeat infections and varying immunity periods

We consider the SIRWJS epidemiological model that includes the waning an...
research
12/27/2019

Estimating Dispersion Curves from Frequency Response Functions via Vector-Fitting

Driven by the need for describing and understanding wave propagation in ...
research
01/25/2023

Lyapunov Stability of First and Second Order GeCo and gBBKS Schemes

In this paper we investigate the stability properties of fixed points of...
research
10/15/2021

Stability and dynamical transition of a electrically conducting rotating fluid

In this article, we aim to study the stability and dynamic transition of...
research
09/10/2022

A constrained gentlest ascent dynamics and its applications to find excited states of Bose-Einstein condensates

In this paper, the gentlest ascent dynamics (GAD) developed in [W. E and...
research
01/13/2022

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

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

Please sign up or login with your details

Forgot password? Click here to reset