Inge CuC

Método sin malla para el estudio de la terapia fototermal en tumores cancerígenos bidimensionales

Juan Felipe

Universidad Pontificia Bolivariana

Carlos Andrés Bustamante Chaverra

Universidad Pontificia Bolivariana

Raúl Adolfo

Universidad Pontificia Bolivariana

Whady Felipe

Universidad Pontificia Bolivariana

DOI: https://doi.org/10.17981/ingecuc.21.2.2025.12

Palabras clave: Ecuación de biocalor, Terapia fototermal plasmónica, Nanoparticulas de oro, Funciones de base radial, LMAPS, Muerte celular

---

Resumen

Introducción: La terapia fototermal plasmónica (PPTT) con nanopartículas metálicas ha ganado relevancia como alternativa menos invasiva en el tratamiento de tumores cancerígenos del tipo adenocarcinoma.  La PPTT utiliza radiación láser para generar efectos plasmónicos en las nanopartículas distribuidas en el tejido cancerígeno, lo que resulta en hipertermia y muerte celular por apoptosis.

Objetivo: Obtener distribuciones de temperatura y a partir de esta, distribuciones bidimiensionales de daño celular en tejidos tumorales sometidos a PPTT.

Metodología: Se implementa una metodología numérica basada en funciones de base radial (RBF) para la solución de la ecuación de Pennes o biocalor y los modelos de Arrhenius y de tres estados para la estimación de la muerte celular. La metodología numérica desarrollada, validada y aplicada está basada en el método de soluciones particulares aproximadas en formulación local (LMAPS).

Resultados: Se demuestra la capacidad del método de solucionar problemas con fuentes variables, regiones múltiples y diferentes tipos de condiciones de fronteras mediante la comparación con la herramienta computacional OpenFOAM basada en volúmenes finitos y resultados numéricos reportados por otros autores. Esto se realiza a partir de la solución de situaciones hipotéticas de transferencia de calor en tejidos incluyendo dominios 1d y 2d con fuentes metabólica, término de perfusión y conversión térmica de la radiación del láser.

Conclusiones: Mediante una situación de PPTT en tejido superficial, se aplica la metodología numérica desarrollada desde la descripción de distribución de energía aportada por el láser hasta la estimación de los porcentajes locales de muerte celular. Esta metodología numérica es la base para el análisis, optimización y diseño de procesos de aplicación de PPTT en ambientes clínicos considerando su potencial para resolver geometrías complejas, condiciones de frontera y parámetros variables en el tiempo y dominios con múltiples regiones.

---

Biografía del autor/a

Juan Felipe, Universidad Pontificia Bolivariana

Received his degree in Mechanical Engineering from Universidad Pontificia Bolivariana in 2022, during which he actively participated in projects and research groups working on subjects related to thermofluids and energy transition. He is interested in fields related to renewable energy.

Raúl Adolfo, Universidad Pontificia Bolivariana

Textile Engineer, Specialization in Mechanical Design, Master in Engineering, PhD degree in Engineering. Research professor of the Faculty of Nanotechnology Engineering of the Universidad Pontificia Bolivariana, Director of the Research Group of Automatics and Design A+D. Knowledge in mathematical modeling and Computational Mechanics.

Whady Felipe, Universidad Pontificia Bolivariana

Ingeniero Mecánico de la Universidad Pontificia Bolivariana, doctorado en mecánica computacional, simulación y modelamiento matemático y numérico del Wessex Institute Of Technology, University of Wales (Reino Unido). con intereses profesionales orientados a la generación de conocimiento, investigación, educación y aplicación del conocimiento y la ciencia. Experiencia en más de 26 años en el sector educativo, par evaluador en procesos de acreditación y certificación, investigación en las áreas de simulación, matemática aplicada, fluidos, fenómenos de transporte y energía.

---

Citas

S. Jain, D. G. Hirst, and J. M. O’Sullivan, “Gold nanoparticles as novel agents for cancer therapy,” Br. J. Radiol., vol. 85, no. 1010, pp. 101–113, 2012, doi: 10.1259/bjr/59448833.

L. Jauffred, A. Samadi, H. Klingberg, P. M. Bendix, and L. B. Oddershede, “Plasmonic Heating of Nanostructures,” Chem. Rev., vol. 119, no. 13, pp. 8087–8130, 2019, doi: 10.1021/acs.chemrev.8b00738.

C. D. Kaddi, J. H. Phan, and M. D. Wang, “Computational nanomedicine: Modeling of nanoparticle-mediated hyperthermal cancer therapy,” Nanomedicine, vol. 8, no. 8, pp. 1323–1333, 2013, doi: 10.2217/nnm.13.117.

I. Raouf, S. Khalid, A. Khan, J. Lee, H. S. Kim, and M. H. Kim, “A review on numerical modeling for magnetic nanoparticle hyperthermia: Progress and challenges,” J. Therm. Biol., vol. 91, p. 102644, 2020, doi: 10.1016/j.jtherbio.2020.102644.

H. Pennes, “Analysis of tissue and arterial blood temperatures in the resting human forearm,” J Appl Physiol, vol. 1, pp. 93–148, 1962, doi: 10.5005/jp/books/12678_10.

E. H. Wissler, “50 years of JAP: Pennes’ 1948 paper revisited,” J. Appl. Physiol., vol. 85, pp. 35–41, 1998.

K. Das and S. C. Mishra, “Study of thermal behavior of a biological tissue: An equivalence of Pennes bioheat equation and Wulff continuum model,” J. Therm. Biol., vol. 45, pp. 103–109, 2014, doi: 10.1016/j.jtherbio.2014.08.007.

J. Ghazanfarian, R. Saghatchi, and D. V. Patil, “Implementation of Smoothed-Particle Hydrodynamics for non-linear Pennes’ bioheat transfer equation,” Appl. Math. Comput., vol. 259, pp. 21–31, 2015, doi: 10.1016/j.amc.2015.02.036.

B. Mochnacki and E. Majchrzak, “Numerical model of thermal interactions between cylindrical cryoprobe and biological tissue using the dual-phase lag equation,” Int. J. Heat Mass Transf., vol. 108, pp. 1–10, 2017, doi: 10.1016/j.ijheatmasstransfer.2016.11.103.

I. K. Tjahjono and Y. Bayazitoglu, “Near-infrared light heating of a slab by embedded nanoparticles,” Int. J. Heat Mass Transf., vol. 51, no. 7–8, pp. 1505–1515, 2008, doi: 10.1016/j.ijheatmasstransfer.2007.07.047.

J. Vera and Y. Bayazitoglu, “Gold nanoshell density variation with laser power for induced hyperthermia,” Int. J. Heat Mass Transf., vol. 52, no. 3–4, pp. 564–573, 2009, doi: 10.1016/j.ijheatmasstransfer.2008.06.036.

Y. Feng et al., “Nanoshell-mediated laser surgery simulation for prostate cancer treatment,” Eng. Comput., vol. 25, no. 1, pp. 3–13, 2009, doi: 10.1007/s00366-008-0109-y.

S. Motaei, M. Ghazavi, and G. Rezazadeh, “Incorporating temperature-dependent properties into the modeling of photo-thermo-mechanical interactions in cancer tissues”, Thermal Science and Engineering Progress, vol. 47, pp. 102351, 2024, doi: 10.1016/j.tsep.2023.102351.

X. Xu, A. Meade, and Y. Bayazitoglu, “Feasibility of selective nanoparticle-assisted photothermal treatment for an embedded liver tumor,” Lasers Med. Sci., vol. 28, no. 4, pp. 1159–1168, 2013, doi: 10.1007/s10103-012-1195-z.

H. C. Huang, K. Rege, and J. J. Heys, “Spatiotemporal temperature distribution and cancer cell death in response to extracellular hyperthermia induced by gold nanorods,” ACS Nano, vol. 4, no. 5, pp. 2892–2900, 2010, doi: 10.1021/nn901884d.

S. Soni and R. K. Sinha, “Controlling parameters for plasmonic photothermal ablation of a tumor,” IEEE J. Sel. Top. Quantum Electron., vol. 22, no. 4, pp. 21–28, 2016, doi: 10.1109/JSTQE.2016.2514359.

A. K. Shaw, D. Khurana, and S. Soni, “Thermal damage analysis of sub-surface soft tissue sarcoma for indocyanine green mediated photothermal cancer therapy", Thermal Science and Engineering Progress, vol. 46, pp. 102168, 2023, doi: 10.1016/j.tsep.2023.102168.

E. J. Kansa, “Multiquadrics-A scattered data approximation scheme with applications to computational fluid-dynamics-II solutions to parabolic, hyperbolic and elliptic partial differential equations,” Comput. Math. with Appl., vol. 19, no. 8–9, pp. 147–161, 1990, doi: 10.1016/0898-1221(90)90271-K.

R. Schaback, “Basis Function Interpolation,” vol. 3, pp. 251–263, 1995.

C. K. Lee, X. Liu, and S. C. Fan, “Local multiquadric approximation for solving boundary value problems,” Comput. Mech., vol. 30, no. 5–6, pp. 396–409, 2003, doi: 10.1007/s00466-003-0416-5.

B. Šarler and R. Vertnik, “Meshfree explicit local radial basis function collocation method for diffusion problems,” Comput. Math. with Appl., vol. 51, no. 8 SPEC. ISS., pp. 1269–1282, 2006, doi: 10.1016/j.camwa.2006.04.013.

T. J. Moroney and I. W. Turner, “A finite volume method based on radial basis functions for two-dimensional nonlinear diffusion equations,” Appl. Math. Model., vol. 30, no. 10, pp. 1118–1133, 2006, doi: 10.1016/j.apm.2005.07.007.

P. H. Chen, C. S., Fan, C. M., & Wen, “The method of approximate particular solutions for solving certain partial differential equations,” Numer. Methods Partial Differ. Equ., vol. 28, no. 2, pp. 506–522, 2012, doi: 10.1002/num.20631.

C. S. Chen, C. M. Fan, and P. H. Wen, “The method of approximate particular solutions for solving elliptic problems with variable coefficients,” Int. J. Comput. Methods, vol. 8, no. 3, pp. 545–559, 2011, doi: 10.1142/S0219876211002484.

K. E. Atkinson, “The numerical evaluation of particular solutions for poisson’s equation,” IMA J. Numer. Anal., vol. 5, no. 3, pp. 319–338, 1985, doi: 10.1093/imanum/5.3.319.

S. Chantasiriwan, “Cartesian grid methods using radial basis functions for solving Poisson, Helmholtz, and diffusion-convection equations,” Eng. Anal. Bound. Elem., vol. 28, no. 12, pp. 1417–1425, 2004, doi: 10.1016/j.enganabound.2004.08.004.

C. A. Bustamante, H. Power, and W. F. Florez, “A global meshless collocation particular solution method for solving the two-dimensional Navier-Stokes system of equations,” Comput. Math. with Appl., vol. 65, no. 12, pp. 1939–1955, 2013, doi: 10.1016/j.camwa.2013.04.014.

G. Yao, J. Kolibal, and C. S. Chen, “A localized approach for the method of approximate particular solutions,” Comput. Math. with Appl., vol. 61, no. 9, pp. 2376–2387, 2011, doi: 10.1016/j.camwa.2011.02.007.

C. A. Bustamante, H. Power, and W. F. Florez, “An efficient accurate Local Method of Approximate Particular Solutions for solving convection-diffusion problems,” Eng. Anal. Bound. Elem., vol. 47, no. 1, pp. 32–37, 2014, doi: 10.1016/j.enganabound.2014.06.004.

X. Zhang, and Yangjiong Wu, “High-resolution strategy for localized method of approximate particular solutions to solve unsteady Navier–Stokes problems," Eng. Anal. Bound. Elem., vol. 159, no. 1, pp. 11–16, 2024, doi: 10.1016/j.enganabound.2023.11.018.

H. Zhang, “Lattice Boltzmann method for solving the bioheat equation,” Phys. Med. Biol., vol. 53, no. 3, 2008, doi: 10.1088/0031-9155/53/3/N01.

S. Soni, H. Tyagi, R. A. Taylor, and A. Kumar, “Investigation on nanoparticle distribution for thermal ablation of a tumour subjected to nanoparticle assisted thermal therapy,” J. Therm. Biol., vol. 43, no. 1, pp. 70–80, 2014, doi: 10.1016/j.jtherbio.2014.05.003.

M. J. C. van Gemert and A. J. Welch, “Time constants in thermal laser medicine,” Lasers Surg. Med., vol. 9, no. 4, pp. 405–421, 1989, doi: 10.1002/lsm.1900090414.

N. Sahoo, S. Ghosh, A. Narasimhan, and S. K. Das, “Investigation of non-Fourier effects in bio-tissues during laser assisted photothermal therapy,” Int. J. Therm. Sci., vol. 76, pp. 208–220, 2014, doi: 10.1016/j.ijthermalsci.2013.08.014.

N. M. Dimitriou, A. Pavlopoulou, I. Tremi, V. Kouloulias, G. Tsigaridas, and A. G. Georgakilas, “Prediction of gold nanoparticle and microwave-induced hyperthermia effects on tumor control via a simulation approach,” Nanomaterials, vol. 9, no. 2, 2019, doi: 10.3390/nano9020167.

C. J. Greenshields, “The Open Source CFD Toolbox - Programmer’s Guide,” no. December. OpenCFD Ltd., United Kingdom, 2015.

U. K. OpenCFD Ltd., “The Open Source CFD Toolbox - User’s Guide,” no. May. OpenFOAM Foundation, 2012