The Modified Picard Iteration Method for Solving the Moving Boundary Problems with a Dissolution Term Encountered in the Drug Release Systems

Document Type : Research Article

Authors

Department of Chemical Engineering, Yasouj University, Yasouj 75918-74831, Iran

Abstract

The primary goal of this work is to develop an accurate analytical solution for the moving boundary problems with a dissolution term encountered in drug release from nanoporous structures. To achieve this, we propose applying the modified Picard iteration method for the first time. Using this approach, the moving boundary problem with Dirichlet boundary conditions is transformed into a single integral and then solved iteratively by selecting a suitable starting function. To validate the proposed method, we first analyze the behavior of the equations in specific cases. For these examples, the modified Picard iteration method provides more accurate results than the Picard iteration method, the variational iteration method, and the Adomian decomposition method. Furthermore, the solutions obtained from the proposed method satisfy both boundary conditions, whereas other iterative approaches do not fulfill one of the boundary conditions. The proposed method is also verified by comparing its results with experimental data on drug release from halloysite tubules reported in the literature. The findings indicate that simultaneous dissolution from the moving interface and the perimeter of the pores may be a possible mechanism for drug release from the halloysite tubules.

Keywords


[1] Garshasbi, M., Dastour, H., 2016. A mollified marching solution of an inverse ablation-type moving boundary problem. Computational and Applied Mathematics, 35(1), 61-73. https://doi.org/10.1007/s40314-014-0180-5
[2] Singh, A.K., Kumar, A., Rajeev. 2019. A Stefan problem with variable thermal coefficients and moving phase change material. Journal of King Saud University - Science, 31(4), 1064-1069. https://doi.org/10.1016/j.jksus.2018.09.009
[3] Beiranvand, A., Ivaz, K., 2014. Solving the Stefan Problem with Kinetics. Computational Methods for Differential Equations, 2(1), 37-49.
[4] Crank, J., 1987. Free and Moving Boundary Problems, Clarendon Press, Oxford shire: New York, p. 424.
[5] Kim, S.H., 2014. Two Simple Numerical Methods for the Free Boundary in One-Phase Stefan Problem. Journal of Applied Mathematics, 2014, 1-10. https://doi.org/10.1155/2014/764532
[6] Sadoun, N., Si-Ahmed, E.-K., Colinet, P., Legrand, J., 2012. On the boundary immobilization and variable space grid methods for tra nsient heat conduction problems with phase change: Discussion and refinement. Comptes Rendus Mécanique, 340(7), 501-511. https://doi.org/10.1016/j.crme.2012.03.003
[7] Karabenli, H., Uçar, Y., Aksan, N., 2016. A collocation finite element solution for Stefan problems with periodic boundary conditions. Filomat, 30(3), 699-709. https://doi.org/10.2298/FIL1603699K
[8] Ivanovic, M., Svicevic, M., Savovic, S., 2017. Numerical solution of Stefan problem with variable space grid method based on mixed finite element/finite difference approach. International Journal of Numerical Methods for Heat & Fluid Flow, 27(12), 2682-2695. https://doi.org/10.1108/HFF-11-2016-0443
[9] Singh, J., Jitendra, Rai, K.N., 2020. Legendre wavelet based numerical solution of variable latent heat moving boundary problem. Mathematics and Computers in Simulation, 178, 485-500. https://doi.org/10.1016/j.matcom.2020.06.020
[10] Chaurasiya, V., Chaudhary, R.K., Awad, M.M., Singh, J., 2022. A numerical study of a moving boundary problem with variable thermal c onductivity and temperature-dependent moving PCM under periodic boundary condition. The European Physical Journal Plus, 137(6), 714. https://doi.org/10.1140/epjp/s13360-022-02927-w
[11] Lee, P.I., 2011. Modeling of drug release from matrix systems involving moving boundaries: Approximate analytical solutions. International Journal of Pharmaceutics, 418(1), 18-27. https://doi.org/10.1016/j.ijpharm.2011.01.019
[12] Fazli-Abukheyli, R., Rahimi, M.R., Ghaedi, M., 2020. Experimental study and modeling of in vitro agrochemicals release from nanoporous anodic alumina. Chemical Papers. https://doi.org/10.1007/s11696-019-01045-9
[13] Yu, Z.-T., Fan, L.-W., Hu, Y.-C., Cen, K.-F., 2010. Perturbation solution to heat conduction in melting or solidification with heat generation. Heat and Mass Transfer, 46(4), 479-483. https://doi.org/10.1007/s00231-010-0596-4
[14] Parhizi, M., Jain, A., 2019. Solution of the Phase Change Stefan Problem with Time-Dependent Heat F lux Using Perturbation Method. Journal of Heat Transfer, 141(2), 024503. https://doi.org/10.1115/1.4041956
[15] Bollati, J., Natale, M.F., Semitiel, J.A., Tarzia, D.A., 2018. Integral balance methods applied to a non-classical Stefan problem. https://doi.org/10.2298/TSCI180901310B
[16] Ribera, H., Myers, T.G., MacDevette, M.M., 2019. Optimising the heat balance integral method in spherical and cylindric al Stefan problems. Applied Mathematics and Computation, 354, 216-231. https://doi.org/10.1016/j.amc.2019.02.039
[17] Wood, A.S., 2001. A new look at the heat balance integral method. Applied Mathematical Modelling, 25(10), 815-824. https://doi.org/10.1016/S0307-904X(01)00016-6
[18] Sadoun, N., Si-Ahmed, E.-K., Colinet, P., 2006. On the refined integral method for the one-phase Stefan problem with t ime-dependent boundary conditions. Applied Mathematical Modelling, 30(6), 531-544. https://doi.org/10.1016/j.apm.2005.06.003
[19] Myers, T.G., 2010. Optimal exponent heat balance and refined integral methods applied to Stefan problems. International Journal of Heat and Mass Transfer, 53(5), 1119-1127. https://doi.org/10.1016/j.ijheatmasstransfer.2009.10.045
[20] Mennig, J., Özişik, M.N., 1985. Coupled integral equation approach for solving melting or solidification. International Journal of Heat and Mass Transfer, 28(8), 1481-1485. https://doi.org/10.1016/0017-9310(85)90250-9
[21] Sudhakar, B., 1992. On integral iterative formulations in classical Stefan problems. Chemical Engineering Science, 47(12), 3158-3162. https://doi.org/10.1016/0009-2509(92)87020-Q
[22] Sudhakar, B., 1992. An integral method for non-linear moving boundary problems. Chemical Engineering Science, 47(2), 475-479. https://doi.org/10.1016/0009-2509(92)80035-B
[23] Rajeev, Kushwaha, M.S., 2013. Homotopy perturbation method for a limit case Stefan problem governed by fractional diffusion equation. Applied Mathematical Modelling, 37(5), 3589-3599. https://doi.org/10.1016/j.apm.2012.07.047
[24] Slota, D., 2010. The application of the homotopy perturbation method to one-phase inverse Stefan problem. International Communications in Heat and Mass Transfer, 37(6), 587-592. https://doi.org/10.1016/j.icheatmasstransfer.2010.03.009
 [25] Słota, D., Chmielowska, A., Brociek, R., Szczygieł, M., 2020. Application of the Homotopy Method for Fractional Inverse Stefan Probl em. Energies, 13(20), 5474. https://doi.org/10.3390/en13205474
[26] Hetmaniok, E., Słota, D., Wituła, R., Zielonka, A., 2011. Comparison of the Adomian decomposition method and the variational ite ration method in solving the moving boundary problem. Computers & Mathematics with Applications, 61(8), 1931-1934. https://doi.org/10.1016/j.camwa.2010.07.050
[27] Qin, X.-Y., Duan, Y.-X., Yin, M.-R., 2014. Approximate Analytic Solutions for the Two-Phase Stefan Problem Using the Adomian Decomposition Method. Journal of Applied Mathematics, 2014, 1-6. https://doi.org/10.1155/2014/391606
[28] Słota, D., 2007. Direct and inverse one-phase Stefan problem solved by the variational iteration method. Computers & Mathematics with Applications, 54(7-8), 1139-1146. https://doi.org/10.1016/j.camwa.2006.12.061
[29] Witula, R., Hetmaniok, E., Slota, D., Zielonka, A., 2010. Application of the Picard's Iterative Method for the Solution of One-P hase Stefan Problem. Archives of Foundry Engineering, 10 (sp. is. ~4), 83-88.
[30] Witula, R., Hetmaniok, E., Slota, D., Zielonka, A., 2011. Solution of the two-phase Stefan problem by using the Picard's iterati ve method. Thermal Science, 15(suppl. 1), 21-26. https://doi.org/10.2298/TSCI11S1021W
[31] Higuchi, T., 1961. Rate of release of medicaments from ointment bases containing drugs in suspension. Journal of Pharmaceutical Sciences, 50(10), 874-875. https://doi.org/10.1002/jps.2600501018
[32] Paul, D.R., McSpadden, S.K., 1976. Diffusional release of a solute from a polymer matrix. Journal of Membrane Science, 1, 33-48. https://doi.org/10.1016/s0376-7388(00)82256-5
[33] Lee, P.I., 1980. Diffusional release of a solute from a polymeric matrix — approximate analytical solutions. Journal of Membrane Science, 7(3), 255-275. https://doi.org/10.1016/s0376-7388(00)80472-x
[34] Abdekhodaie, M.J., Cheng, Y.L., 1997. Diffusional release of a dispersed solute from planar and spherical matrices into finite external volume. Journal of Controlled Release, 43(2), 175-182. https://doi.org/10.1016/s0168-3659(96)01482-4
[35] Cohen, D.S., Erneux, T., 1998. Controlled Drug Release Asymptotics. SIAM Journal on Applied Mathematics, 58(4), 1193-1204. https://doi.org/10.1137/S0036139995293269
[36] Jain, A., McGinty, S., Pontrelli, G., 2022. Drug diffusion and release from a bioerodible spherical capsule. International Journal of Pharmaceutics, 616, 121442. https://doi.org/10.1016/j.ijpharm.2021.121442
[37] Garshasbi, M., Malek Bagomghaleh, S., 2022. On a moving boundary problem associated with the swelling drug release platforms. International Journal of Computer Mathematics, 99(12), 2499-2523. https://doi.org/10.1080/00207160.2022.2069466
[38] Robin, W.A., 2010. Solving differential equations using modified Picard iteration. International Journal of Mathematical Education in Science and Technol ogy, 41(5), 649-665. https://doi.org/10.1080/00207391003675182
[39] Mani, G., Macias, C.E., Feldman, M.D., Marton, D., Oh, S., Mauli Agrawal, C., 2010. Delivery of paclitaxel from cobalt–chromium alloy surfaces without polymeric carriers. Biomaterials, 31(20), 5372-5384. https://doi.org/10.1016/j.biomaterials.2010.03.043
[40] Zhang, Y., He, F., Sun, Z., Li, L., Huang, Y., 2014. Controlled delivery of dexamethasone from TiO2 film with nanoporous structure on Ti–25Nb–3Mo–2Sn–3Zr biomedical alloy without polymeric carrier. Materials Letters, 128, 384-387. https://doi.org/10.1016/j.matlet.2014.04.149
[41] Veerabadran, N.G., Price, R.R., Lvov, Y.M., 2007. Clay nanotubes for encapsulation and sustained release of drugs. Nano, 02(02), 115-120. https://doi.org/10.1142/s1793292007000441
[42] Elumalai, D.N., Lvov, Y., Derosa, P., 2015. Implementation of a Simulation Model of the Controlled Release of Molecular Species from Halloysite Nanotubes. Journal of Encapsulation and Adsorption Sciences, 05(01), 74-92. https://doi.org/10.4236/jeas.2015.51006