A second order method for a drug release process defined by a differential Maxwell-Wichert stress-strain relation
DOI: https://doi.org/10.3846/mma.2026.22960Abstract
Polymeric drug delivery platforms offer promising capabilities for controlled drug release thanks to their ability to be custom-designed with specific properties. In this paper we present a model to simulate the complex interplay between solvent absorption, polymer swelling, drug release and stress development within these platforms. A system of nonlinear partial differential equations coupled with an ordinary differential equation is introduced to avoid drawbacks from other models found in the literature. These incorporated a memory effect but from a numerical standpoint, required storing all previous time steps, making them computationally expensive. This paper proposes a new numerical method to simulate such devices based on nonuniform grids and an implicit midpoint time discretization. Our main results are the second order convergence of the method for nonsmooth solutions and the scheme’s stability under the assumption of quasiuniform grids and a small enough timestep.
Keywords:
drug delivery, polymeric devices, Maxwell-Wichert, numerical method, second-order convergenceHow to Cite
Share
License
Copyright (c) 2026 The Author(s). Published by Vilnius Gediminas Technical University.
This work is licensed under a Creative Commons Attribution 4.0 International License.
References
L. Angermann and S. Wang. A super-convergent unsymmetric finite volume method for convection-diffusion equations. Journal of Computational and Applied Mathematics, 358:179–189, 2019. https://doi.org/10.1016/j.cam.2019.03.017
S. Barbeiro, S. Bardeji, J.A. Ferreira and L. Pinto. Non-Fickian diffusion coupled models in porous media. Numerische Mathematik, 138:869–904, 2018. https://doi.org/10.1007/s00211-017-0922-6
S. Barbeiro, J.A. Ferreira and R. Grigorieff. Supraconvergence of a finite difference scheme for solutions in Hs(0,L). IMA Journal of Numerical Analysis, 25:797–811, 2005. https://doi.org/10.1093/imanum/dri018
J.S. Borges, G.C.M. Campos, J.A. Ferreira and G. Romanazzi. Drug release from polymeric platforms for nonsmooth solutions. Applied Numerical Mathematics, 213:12–37, 2025. https://doi.org/10.1016/j.apnum.2025.02.016
J.S. Borges, J.A. Ferreira, G. Romanazzi and E. Abreu. Drug release from viscoelastic polymeric matrices - a stable and supraconvergent FDM. Computers & Mathematics with Applications, 99:257–269, 2021. https://doi.org/10.1016/j.camwa.2021.08.007
J. Bramble and S. Hilbert. Estimation of linear functionals on Sobolev spaces with application to Fourier transforms and spline interpolation. SIAM Journal on Numerical Analysis, 7:112–124, 1970. https://doi.org/10.1137/0707006
H.F. Brinson and L.C. Brinson. Polymer Engineering Science and Viscoelasticity. Springer US, 2008. https://doi.org/10.1007/978-0-387-73861-1
F. De Hoog and D. Jackett. On the rate of convergence of finite difference schemes on nonuniform grids. The ANZIAM Journal, 26(3):247–256, 1985. https://doi.org/10.1017/S0334270000004495
X. Feng, R. He and Z. Chen. H1-superconvergence of finite difference method base on Q1-element on quasi-uniform mesh for the 3D Poisson equation. Numerical Methods for Partial Differential Equations, 36:29–48, 2020. https://doi.org/10.1002/num.22386
J.A. Ferreira, P. de Oliveira, M. Grassi and G. Romanazzi. Drug Release from Viscoelastic Swelling Polymeric Platforms. SIAM Journal on Applied Mathematics, 78(3):1378–1401, 2018. https://doi.org/10.1137/16M1100344
J.A. Ferreira, E. Gudin˜o and P. de Oliveira. A Second Order Approximation for Quasilinear Non-Fickian Diffusion Models. Computational Methods in Applied Mathematics, 13(4):471–493, 2013. https://doi.org/10.1515/cmam-2013-0017
J.A. Ferreira, H.P. Go´mez and L. Pinto. A Mathematical Model for NIR Light Protocol Optimization in Controlled Transdermal Drug Delivery. Applied Mathematical Modelling, 112:1–17, 2022. https://doi.org/10.1016/j.apm.2022.07.014
J.A. Ferreira and G. Pena. FDM/FEM for nonlinear convection-diffusionreaction equations with Neumann boundary conditions – Convergence analysis for smooth and nonsmooth solutions. Journal of Computational and Applied Mathematics, 446:115866, 2024. https://doi.org/10.1016/j.cam.2024.115866
R. Grigorieff. Some stability inequalities for compact finite difference operators. Mathematische Nachrichten, 135:93–102, 1988. https://doi.org/10.1002/mana.19881350110
H.O. Kreiss, T.A. Manteuffel, B. Swartz, B. Wendroffand and Jr A.B. White. Supraconvergent schemes on irregular grids. Mathematics Of Computation, 47:537–554, 1986. https://doi.org/10.1090/S0025-5718-1986-0856701-5
N. Kumar, J. Tushar and J. Yuan. Convergence analysis of virtual element method for electric interface model on polygonal meshes with small edges. Computers and Mathematics with Applications, 159:254–266, 2024. https://doi.org/10.1016/j.camwa.2024.02.029
D. Li and J. Wang. Unconditionally optimal errors analysis of Crank-Nicolson Galerkin FEMs for strongly nonlinear parabolic systems. Journal of Scientific Computing, 72:892–915, 2017. https://doi.org/10.1007/s10915-017-0381-3
J. Lo´pez-Marcos and J. Sanz-Serna. A definition of stability for nonlinear problems. In Numerical Treatment of Differential Equations, Proceedings of NUMDIFF-4, pp. 216–226, 1988.
J.C. Lo´pez-Marcos and J.M. Sanz-Serna. Stability and Convergence in Numerical Analysis III: Linear Investigation of Nonlinear Stability. IMA Journal of Numerical Analysis, 8(1):71–84, 1988. https://doi.org/10.1093/imanum/8.1.71
T.A. Manteuffel and A.B. White. The Numerical Solution of Second-Order Boundary Value Problems on Nonuniform Meshes. Mathematics of Computation, 47(176):511–535, 1986. https://doi.org/10.1090/S0025-5718-1986-0856700-3
T. Ortega and J. Sanz-Serna. Nonlinear stability and convergence of finitedifference methods for the ”good” Boussineq equation. Numerische Mathematik, 58:215–229, 1990. https://doi.org/10.1007/BF01385620
E. Su¨li. Convergence and nonlinear stability of the LagrangeGalerkin method. Numerich Mathematik, 53:459–483, 1988. https://doi.org/10.1007/BF01396329
J. Wang. Superconvergence analysis of an energy stable scheme for nonlinear reaction-diffusion equation with BDF mixed FEM. Applied Numerical Mathematics, 153:457–472, 2020. https://doi.org/10.1016/j.apnum.2020.03.007
S. Wang and G. Kreiss. Convergence of summation-by-parts finite difference methods for the wave equation. Journal of Scientific Computing, 71:219–245, 2017. https://doi.org/10.1007/s10915-016-0297-3
Y. Wang, H. Yi, X. Fang and G. Li. Unconditionally optimal error estimates of linearized Crank-Nicolson virtual elements for quasilinear parabolic problems on general polygonal meshes. ESAIM: Mathematical Modelling and Numerical Analysis, 58:881–926, 2024. https://doi.org/10.1051/m2an/2024017
M. Wheeler. A priori l2 error estimates for Galerkin approximations to parabolic partial differential equations. SIAM Journal on Numerical Analysis, 10:723–759, 1973. https://doi.org/10.1137/0710062
J. Yan, Z. Yin and A. Zhu. H1-Galerkin mixed finite element method the vibration equation of beam with structural damping. Computational and Applied Mathematics, 43:308, 2024. https://doi.org/10.1007/s40314-024-02831-2
J. Zhang and H. Rui. A coupling of Galerkin and mixed finite element methods for the quasi-static thermo-poroelasticity with nonlinear convective transport. Journal of Computational and Applied Mathematics, 441:115672, 2024. https://doi.org/10.1016/j.cam.2023.115672
D. Zhu, X. Feng and L. Qian. Error analysis of Crank-Nicolson-Leapfrog scheme for the two-phase Cahn-Hilliard-Navier-Stokes incompressible flow. Computers and Mathematics with Applications, 172:78–93, 2024. https://doi.org/10.1016/j.camwa.2024.07.026
View article in other formats
CrossMark check
Published
Issue
Section
Copyright
Copyright (c) 2026 The Author(s). Published by Vilnius Gediminas Technical University.
License
This work is licensed under a Creative Commons Attribution 4.0 International License.