Related papers: An Efficient Finite Element Iterative Method for S…
We present a multiscale atomistic-to-continuum method for ionic crystals with defects. Defects often play a central role in ionic and electronic solids, not only to limit reliability, but more importantly to enable the functionalities that…
High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on…
We study the geometric particle-in-cell methods for an electrostatic hybrid plasma model. In this model, ions are described by the fully kinetic equations, electron density is determined by the Boltzmann relation, and space-charge effects…
In our recent work [H. Zhang, F.X. Trias, A. Oliva, D. Yang, Y. Tan, Y. Sheng. PIBM: Particulate immersed boundary method for fluid-particle interaction problems. Powder Technology. 272(2015), 1-13.], a particulate immersed boundary method…
A finite element methodology for large classes of variational boundary value problems is defined which involves discretizing two linear operators: (1) the differential operator defining the spatial boundary value problem; and (2) a Riesz…
We consider the mathematical analysis and numerical approximation of a system of nonlinear partial differential equations that arises in models that have relevance to steady isochoric flows of colloidal suspensions. The symmetric velocity…
The gas of the interacted electrons is usually described within Kohn-Sham approximation by the set of Poisson and Schr\"{o}dinger equations with an effective potential for the single-particle wave functions. The solution of these equations…
We consider the integrable family of symmetric boundary-driven interacting particle systems that arise from the non-compact XXX Heisenberg model in one dimension with open boundaries. In contrast to the well-known symmetric exclusion…
We present a novel framework for Finite Element Particle-in-Cell methods based on the discretization of the underlying Hamiltonian structure of the Vlasov-Maxwell system. We derive a semi-discrete Poisson bracket, which retains the defining…
The Poisson-Boltzmann equation offers an efficient way to study electrostatics in molecular settings. Its numerical solution with the boundary element method is widely used, as the complicated molecular surface is accurately represented by…
We present a spectral element model for general-purpose simulation of non-overturning nonlinear water waves using the incompressible Navier-Stokes equations (INSE) with a free surface. The numerical implementation of the spectral element…
Most organic and inorganic surfaces (e.g., glass, nucleic acids or lipid membranes) become charged in aqueous solutions. The resulting ionic distribution induces effective interactions between the charged surfaces. Stacks of like-charged…
The steady-state Poisson-Nernst-Planck (ssPNP) equations are an effective model for the description of ionic transport in ion channels. It is observed that an ion channel exhibits voltage-dependent switching between open and closed states.…
We show strong uniform convergence of monotone P1 finite element methods to the viscosity solution of isotropic parabolic Hamilton-Jacobi-Bellman equations with mixed boundary conditions on unstructured meshes and for possibly degenerate…
A numerical method using implicit surface representations is proposed to solve the linearized Poisson-Boltzmann equations that arise in mathematical models for the electrostatics of molecules in solvent. The proposed method used an implicit…
We have developed efficient numerical algorithms for solving 3D steady-state Poisson-Nernst-Planck (PNP) equations with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are…
In situations involving large potentials or surface charges, the Poisson Boltzman(PB) equation has shortcomings because it neglects ion-ion interactions and steric effects. This has been widely recognized by the electrochemistry community,…
An efficient nonlinear multigrid method for a mixed finite element method of the Darcy-Forchheimer model is constructed in this paper. A Peaceman-Rachford type iteration is used as a smoother to decouple the nonlinearity from the divergence…
In this paper, we consider the numerical solution of some nonlinear poroelasticity problems that are of Biot type and develop a general algorithm for solving nonlinear coupled systems. We discuss the difficulties associated with flow and…
We present a linear scaling formulation for the solution of the all-electron Coulomb problem in crystalline solids. The resulting method is systematically improvable and well suited to large-scale quantum mechanical calculations in which…