Adaptive and iterative methods for simulations of nanopores with the PNP-Stokes equations
From MaRDI portal
Publication:2424486
Abstract: We present a 3D finite element solver for the nonlinear Poisson-Nernst-Planck (PNP) equations for electrodiffusion, coupled to the Stokes system of fluid dynamics. The model serves as a building block for the simulation of macromolecule dynamics inside nanopore sensors. We add to existing numerical approaches by deploying goal-oriented adaptive mesh refinement. To reduce the computation overhead of mesh adaptivity, our error estimator uses the much cheaper Poisson-Boltzmann equation as a simplified model, which is justified on heuristic grounds but shown to work well in practice. To address the nonlinearity in the full PNP-Stokes system, three different linearization schemes are proposed and investigated, with two segregated iterative approaches both outperforming a naive application of Newton's method. Numerical experiments are reported on a real-world nanopore sensor geometry. We also investigate two different models for the interaction of target molecules with the nanopore sensor through the PNP-Stokes equations. In one model, the molecule is of finite size and is explicitly built into the geometry; while in the other, the molecule is located at a single point and only modeled implicitly -- after solution of the system -- which is computationally favorable. We compare the resulting force profiles of the electric and velocity fields acting on the molecule, and conclude that the point-size model fails to capture important physical effects such as the dependence of charge selectivity of the sensor on the molecule radius.
Recommendations
- A stabilized finite element method for modified Poisson-Nernst-Planck equations to determine ion flow through a nanopore
- Poisson-Nernst-Planck equations for simulating biomolecular diffusion-reaction processes. I: Finite element solutions
- An Effective Finite Element Iterative Solver for a Poisson--Nernst--Planck Ion Channel Model with Periodic Boundary Conditions
- A stabilized finite element method for the Poisson-Nernst-Planck equations in three-dimensional ion channel simulations
- A conservative discretization of the Poisson-Nernst-Planck equations on adaptive Cartesian grids
Cites work
- A Convergent Adaptive Algorithm for Poisson’s Equation
- A new finite element formulation for computational fluid dynamics. VII. The Stokes problem with various well-posed boundary conditions: Symmetric formulations that converge for all velocity/pressure spaces
- A transport equation for confined structures derived from the Boltzmann equation
- A Transpose-Free Quasi-Minimal Residual Algorithm for Non-Hermitian Linear Systems
- An abstract analysis of optimal goal-oriented adaptivity
- An optimal control approach to a posteriori error estimation in finite element methods
- Automated goal-oriented error control. I: Stationary variational problems
- Axioms of adaptivity
- Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems
- Convergence of adaptive BEM and adaptive FEM-BEM coupling for estimators without \(h\)-weighting factor
- Convergent finite element discretizations of the Navier-Stokes-Nernst-Planck-Poisson system
- Energetically stable discretizations for charge transport and electrokinetic models
- Gmsh: a 3-D finite element mesh generator with built-in pre- and post-processing facilities
- scientific article; zbMATH DE number 2088244 (Why is no real title available?)
- scientific article; zbMATH DE number 1878659 (Why is no real title available?)
- On the Effectiveness of Gummel’s Method
- Poisson-Nernst-Planck equations for simulating biomolecular diffusion-reaction processes. I: Finite element solutions
- Preconditioning discretizations of systems of partial differential equations.
- Stochastic Problems in Physics and Astronomy
- Variational multiscale models for charge transport
Cited in
(9)- Use and comparative assessment of the CVFEM method for Poisson-Boltzmann and Poisson-Nernst-Planck three dimensional simulations of impedimetric nano-biosensors operated in the DC and AC small signal regimes
- Existence and local uniqueness for the Stokes-Nernst-Planck-drift-diffusion-Poisson system modeling nanopore and nanowire sensors
- Modelling of resistive pulse sensing: flexible methods for submicron particles
- A stabilized finite element method for modified Poisson-Nernst-Planck equations to determine ion flow through a nanopore
- New mixed finite element methods for the coupled Stokes and Poisson–Nernst–Planck equations in Banach spaces
- Banach spaces-based mixed finite element methods for the coupled Navier-Stokes and Poisson-Nernst-Planck equations
- Analysis of a finite element method for the Stokes-Poisson-Boltzmann equations
- The small-convection limit in a Nernst-Planck-Navier-Stokes system involving Dirichlet potential boundary condition
- Transient electrohydrodynamic flow with concentration-dependent fluid properties: modelling and energy-stable numerical schemes
This page was built for publication: Adaptive and iterative methods for simulations of nanopores with the PNP-Stokes equations
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2424486)
