- Journal Home
- Volume 18 - 2025
- Volume 17 - 2024
- Volume 16 - 2023
- Volume 15 - 2022
- Volume 14 - 2021
- Volume 13 - 2020
- Volume 12 - 2019
- Volume 11 - 2018
- Volume 10 - 2017
- Volume 9 - 2016
- Volume 8 - 2015
- Volume 7 - 2014
- Volume 6 - 2013
- Volume 5 - 2012
- Volume 4 - 2011
- Volume 3 - 2010
- Volume 2 - 2009
- Volume 1 - 2008
Numer. Math. Theor. Meth. Appl., 18 (2025), pp. 127-156.
Published online: 2025-04
Cited by
- BibTex
- RIS
- TXT
In this paper, we propose and analyze a second order accurate (in time) mass lumped mixed finite element numerical scheme for the liquid thin film coarsening model with a singular Leonard-Jones energy potential. The backward differentiation formula (BDF) stencil is applied in the temporal discretization, and a convex-concave decomposition is derived, so that the concave part corresponds to a quadratic energy. In turn, the Leonard-Jones potential term is treated implicitly and the concave part is approximated by a second order Adams-Bashforth explicit extrapolation. An artificial Douglas-Dupont regularization term is added to ensure the energy stability. Furthermore, we provide a theoretical justification that this numerical algorithm has a unique solution, such that the positivity property 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 and the mass lumped FEM approach play an essential role in the positivity-preserving property in the discrete level. In addition, an optimal rate convergence estimate in the $ℓ^∞(0, T ; H^{−1}_h )∩ ℓ^2 (0, T ; H^1_h )$ norm is presented. Finally, two numerical experiments are carried out to verify the theoretical properties.
}, issn = {2079-7338}, doi = {https://doi.org/10.4208/nmtma.OA-2024-0081}, url = {http://global-sci.org/intro/article_detail/nmtma/23944.html} }In this paper, we propose and analyze a second order accurate (in time) mass lumped mixed finite element numerical scheme for the liquid thin film coarsening model with a singular Leonard-Jones energy potential. The backward differentiation formula (BDF) stencil is applied in the temporal discretization, and a convex-concave decomposition is derived, so that the concave part corresponds to a quadratic energy. In turn, the Leonard-Jones potential term is treated implicitly and the concave part is approximated by a second order Adams-Bashforth explicit extrapolation. An artificial Douglas-Dupont regularization term is added to ensure the energy stability. Furthermore, we provide a theoretical justification that this numerical algorithm has a unique solution, such that the positivity property 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 and the mass lumped FEM approach play an essential role in the positivity-preserving property in the discrete level. In addition, an optimal rate convergence estimate in the $ℓ^∞(0, T ; H^{−1}_h )∩ ℓ^2 (0, T ; H^1_h )$ norm is presented. Finally, two numerical experiments are carried out to verify the theoretical properties.