Related papers: A structure-preserving, operator splitting scheme …
We present and analyze finite difference numerical schemes for the Allen Cahn/Cahn-Hilliard equation with a logarithmic Flory Huggins energy potential. Both the first order and second order accurate temporal algorithms are considered. In…
This work proposes a suite of numerical techniques to facilitate the design of structure-preserving integrators for nonlinear dynamics. The celebrated LaBudde-Greenspan integrator and various energy-momentum schemes adopt a difference…
We propose fully discrete, implicit-in-time finite-volume schemes for a general family of non-linear and non-local Fokker-Planck equations with a gradient-flow structure, usually known as aggregation-diffusion equations, in any dimension.…
This work is devoted to study unconditionally energy stable and mass-conservative numerical schemes for the following repulsive-productive chemotaxis model: Find $u \geq 0$, the cell density, and $v \geq 0$, the chemical concentration, such…
We propose a nonlinear Discrete Duality Finite Volume scheme to approximate the solutions of drift diffusion equations. The scheme is built to preserve at the discrete level even on severely distorted meshes the energy / energy dissipation…
Stochastic chemical systems with diffusion are modeled with a reaction-diffusion master equation. On a macroscopic level, the governing equation is a reaction-diffusion equation for the averages of the chemical species. On a mesoscopic…
We consider unsteady poroelasticity problem in fractured porous medium within the classical Barenblatt double-porosity model. For numerical solution of double-porosity poroelasticity problems we construct splitting schemes with respect to…
In this work, the MMC-TDGL equation, a stochastic Cahn-Hilliard equation is solved numerically by using the finite difference method in combination with a convex splitting technique of the energy functional. For the non-stochastic case, we…
In this paper, we present a first-order finite element scheme for the viscoelastic electrohydrodynamic model. The model incorporates the Poisson-Nernst-Planck equations to describe the transport of ions and the Oldroyd-B constitutive model…
Different relaxation approximations to partial differential equations, including conservation laws, Hamilton-Jacobi equations, convection-diffusion problems, gas dynamics problems, have been recently proposed. The present paper focuses onto…
We present a novel energy-based numerical analysis of semilinear diffusion-reaction boundary value problems. Based on a suitable variational setting, the proposed computational scheme can be seen as an energy minimisation approach. More…
Split form schemes for Euler and Navier-Stokes equations are useful for computation of turbulent flows due to their better robustness. This is because they satisfy additional conservation properties of the governing equations like kinetic…
In this study, a family of second order process based modified Patankar Runge-Kutta schemes is proposed with both the mass and mole maintained in balance while preserving the positivity of density and pressure with the time step determined…
This paper develops a two-level fourth-order scheme for solving time-fractional convection-diffusion-reaction equation with variable coefficients subjected to suitable initial and boundary conditions. The basis properties of the new…
In this paper we propose and analyze a finite difference numerical scheme for the Flory-Huggins-Cahn-Hilliard equation with dynamical boundary condition. The singular logarithmic potential is included in the Flory-Huggins energy expansion.…
We develop a set of numerical schemes for the Poisson--Nernst--Planck equations. We prove that our schemes are mass conservative, uniquely solvable and keep positivity unconditionally. Furthermore, the first-order scheme is proven to be…
We study numerical methods for dissipative particle dynamics (DPD), which is a system of stochastic differential equations and a popular stochastic momentum-conserving thermostat for simulating complex hydrodynamic behavior at mesoscales.…
This paper concerns the construction and analysis of a numerical scheme for a mixed discrete-continuous fragmentation equation. A finite volume scheme is developed, based on a conservative formulation of a truncated version of the…
Many applications involve partial differential equations which admits nontrivial steady state solutions. The design of schemes which are able to describe correctly these equilibrium states may be challenging for numerical methods, in…
The energy balance of a partially molten rocky planet can be expressed as a non-linear diffusion equation using mixing length theory to quantify heat transport by both convection and mixing of the melt and solid phases. In this formulation…