Resumen: Mathematical and computational models provide a powerful framework for understanding complex biological processes such as cancer cell progression. Here, we present a numerical formulation for the simulation of the evolution of glioblastoma (GBM) cancer cells in microfluidic devices which are commonly used to replicate the dynamic changes of the tumor cells in a biomimetic microenvironment. We model this physicochemical and biological complexity
with a coupled nonlinear system of transient partial differential equations involving different chemical species and cell phenotypes. In particular, we consider oxygen as the main chemical driver and the concentration of two cell phenotypes: living and dead cells. The system is solved combining a high-order continuous Galerkin finite element formulation in space with a high-order diagonally implicit Runge–Kutta (DIRK) scheme in time. This leads to a coupled nonlinear system for oxygen and living cells at each stage of the DIRK scheme that we solve using the Newton method. The same integration method is used to solve for dead cells at each mesh node. Finally, we present several examples to assess and illustrate the capabilities of the proposed formulation. The results demonstrate that this model is a valuable tool for advancing our understanding of cancer cell progression and for supporting the industrial cesign of novel devices aimed at testing new hypotheses and guiding experimental research. Idioma: Inglés DOI: 10.1016/j.finel.2026.104551 Año: 2026 Publicado en: FINITE ELEMENTS IN ANALYSIS AND DESIGN 258 (2026), 104551 [27 pp.] ISSN: 0168-874X Financiación: info:eu-repo/grantAgreement/ES/AEI/PID2021-126051OB-C41 Financiación: info:eu-repo/grantAgreement/ES/AEI/PID2021-126051OB-C42 Financiación: info:eu-repo/grantAgreement/EUR/MICINN/TED2021-129512B-I00 Tipo y forma: Article (Published version) Área (Departamento): Área Mec.Med.Cont. y Teor.Est. (Dpto. Ingeniería Mecánica)