Related papers: Local Structure-Preserving Relaxation Method for C…
The linearized of the Poisson-Nernst-Planck (PNP) equation under closed ends around a neutral state is studied. It is reduced to a damped heat equation under non-local boundary conditions, which leads to a stochastic interpretation of the…
The Poisson-Nernst-Planck (PNP) system is a standard model for describing ion transport. In many applications, e.g., ions in biological tissues, the presence of thin boundary layers poses both modelling and computational challenges. In this…
We consider an ensemble of mass collisionless particles, which interact mutually either by an attraction of Newton's law of gravitation or by an electrostatic repulsion of Coulomb's law, under a background downward gravity in a…
The method of regularized stokeslets is extensively used in biological fluid dynamics due to its conceptual simplicity and meshlessness. This simplicity carries a degree of cost in computational expense and accuracy because the number of…
In this paper, we investigate linear first- and second-order numerical schemes for the Allen--Cahn equation with a general (possibly degenerate) mobility. Compared with existing numerical methods, our schemes employ a novel dynamic…
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in…
We consider an elastic/viscoelastic transmission problem for the Bresse system with fully Dirichlet or Dirichlet-Neumann-Neumann boundary conditions. The physical model consists of three wave equations coupled in certain pattern. The system…
The flow of ions through permeable channels causes voltage drop in physiological nanodomains such as synapses, dendrites and dendritic spines, and other protrusions. How the voltage changes around channels in these nanodomains has remained…
Compressible Mooney-Rivlin theory has been used to model hyperelastic solids, such as rubber and porous polymers, and more recently for the modeling of soft tissues for biomedical tissues, undergoing large elastic deformations. We propose a…
Efficient and stable solution of partial differential equations (PDEs) is central to scientific and engineering applications, yet existing numerical solvers rely heavily on matrix based discretizations, while learning based methods require…
We present the design and implementation of an L2-stable spectral method for the discretization of the Vlasov- Poisson model of a collisionless plasma in one space and velocity dimension. The velocity and space dependence of the Vlasov…
This work is concerned with the stationary Poisson--Nernst--Planck equation with a large parameter which describes a huge number of ions occupying an electrolytic region. Firstly, we focus on the model with a single specie of positive…
This work presents a finite element method for a modified Poisson-Nernst-Planck/Navier-Stokes (PNP/NS) model under the mechanical equilibrium, developed for compressible electrolytes. Another key contribution of this work is the reduction…
We design an arbitrary-order free energy satisfying discontinuous Galerkin (DG) method for solving time-dependent Poisson-Nernst-Planck systems. Both the semi-discrete and fully discrete DG methods are shown to satisfy the corresponding…
We consider a model problem of the scattering of linear acoustic waves in free homogeneous space by an elastic solid. The stress tensor in the solid combines the effect of a linear dependence of strains with the influence of an existing…
In this work, we propose a nonlinear stabilization technique for scalar conservation laws with implicit time stepping. The method relies on an artificial diffusion method, based on a graph-Laplacian operator. It is nonlinear, since it…
We present a method for computing locally varying nonlinear mechanical properties in particle simulations of amorphous solids. Plastic rearrangements outside a probed region are suppressed by introducing an external field that directly…
Starting with a microscopic (individual-based) Brownian dynamics model of charged particles (ions), its macroscopic description is derived as a system of partial differential equations that govern the evolution of ion concentrations in…
In this paper, we propose a robust and efficient numerical framework for simulating multicomponent gas flow in poroelastic media, with a focus on preserving fundamental thermodynamic principles and ensuring computational reliability. The…
The Poisson-Boltzmann equation (PBE) is a nonlinear elliptic PDE that arises in biomolecular modeling and is a fundamental tool for structural biology. It is used to calculate electrostatic potentials around an ensemble of fixed charges…