Abstract:In this paper, two finite difference numerical schemes are proposed and analyzed for the droplet liquid film model, with a singular Leonard-Jones energy potential involved. Both first and second order accurate temporal algorithms are considered. In the first order scheme, the convex potential and the surface diffusion terms are implicitly, while the concave potential term is updated explicitly. Furthermore, we provide a theoretical justification that this numerical algorithm has a unique solution, such that the positivity is always preserved for the phase variable at a point-wise level, so that a singularity is avoided in the scheme. In fact, the singular nature of the Leonard-Jones potential term around the value of 0 prevents the numerical solution reaching such singular value, so that the positivity structure is always preserved. Moreover, an unconditional energy stability of the numerical scheme is derived, without any restriction for the time step size. In the second order numerical scheme, the BDF temporal stencil is applied, and an alternate convex-concave decomposition is derived, so that the concave part corresponds to a quadratic energy. In turn, the combined Leonard-Jones potential term is treated implicitly, and the concave part the is approximated by a second order Adams-Bashforth explicit extrapolation, and an artificial Douglas-Dupont regularization term is added to ensure the energy stability. The unique solvability and the positivity-preserving property for the second order scheme could be similarly established. In addition, optimal rate convergence analysis is provided for both the first and second order accurate schemes. A few numerical simulation results are also presented, which demonstrate the robustness of the numerical schemes.

An asymptotically stable scheme for diffusive coagulation-fragmentation models

Numerical study of finite volume scheme for coagulation–fragmentation equations with singular rates

Numerical schemes for a fully nonlinear coagulation-fragmentation model coming from wave kinetic theory

Discrete Coagulation-Fragmentation equations with multiplicative coagulation kernel and constant fragmentation kernel

Mass-conserving solutions to coagulation-fragmentation equations with balanced growth

Mass-conserving self-similar solutions to coagulation-fragmentation equations

A finite volume scheme for boundary-driven convection–diffusion equations with relative entropy structure

High Resolution Finite Difference Schemes for a Size Structured Coagulation-Fragmentation Model in the Space of Radon Measures

On the Courant-Friedrichs-Lewy condition for numerical solvers of the coagulation equation

Simple Positivity Preserving Nonlinear Finite Volume Scheme for Subdiffusion Equations on General Non-conforming Distorted Meshes

Boundary value for a nonlinear transport equation emerging from a stochastic coagulation-fragmentation type model

Applications of improved duality lemmas to the discrete coagulation-fragmentation equations with diffusion

A stable and accurate scheme for solving the Stefan problem coupled with natural convection using the Immersed Boundary Smooth Extension method

Structure-preserving, energy stable numerical schemes for a liquid thin film coarsening model

Asymptotic preserving and uniformly unconditionally stable finite difference schemes for kinetic transport equations

A positive and bounded convergent scheme for general space-fractional diffusion-reaction systems with inertial times

Structure-preserving scheme for fractional nonlinear diffusion equations

A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources

Smoothness of moments of the solutions of discrete coagulation equations with diffusion

Numerical analysis and physical simulations for the time fractional radial diffusion equation

Existence of mass-conserving weak solutions to the singular coagulation equation with multiple fragmentation