Fully mixed finite element methods for the coupling of viscoelasticity and reaction-diffusion models
Abstract
We propose and analyse a fully mixed finite element method for a two-way coupled mechanochemical model of calcium signalling in viscoelastic tissue. The mechanical response follows a Kelvin--Voigt law, with the elastic and viscous stresses, velocity, and spin tensor as unknowns, and the symmetry of the total stress imposed weakly. The reaction-diffusion system for calcium and inactivated inositol triphosphate is also formulated in mixed form, introducing the diffusive flux as an additional unknown. The coupling is bidirectional: calcium generates an active stress, while local dilation feeds back into the calcium dynamics. We establish global existence and uniqueness by reducing the coupled system to an integro-differential problem and subsequently recovering the mixed variables. The semidiscrete scheme is formulated for general conforming finite element spaces and shown to be well posed under a single inf-sup condition. An interpolation operator that commutes with the divergence and preserves the weak symmetry pairing yields optimal a priori error estimates. We specialise the analysis to three families of weakly symmetric finite elements, deriving explicit convergence rates and, for one family, a superconvergent local post-processing of the displacement. Numerical tests in two and three dimensions confirm the predicted convergence rates and illustrate the method through viscoelastic benchmarks and mechanically driven calcium wavetrains.
Related
- A Multiscale Finite Element Method for reaction-diffusion eigenproblems arising from neutronics
- Subdivision-based collocation methods for time fractional convection diffusion reaction equation with environmental and biomedical applications
- A structure-aware construction of quantum BPX preconditioners for Q1 finite element discretizations of reaction-diffusion problems
- Finite difference methods for three kinds of reaction-diffusion equations with free boundaries
- A virtual element method for Kelvin--Voigt viscoelasticity
BibTeX
@misc{bustos2026fully,
title = {{Fully mixed finite element methods for the coupling of viscoelasticity and reaction-diffusion models}},
author = {Alonso J. Bustos and Jeonghun J. Lee and Ricardo Ruiz-Baier},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.03317},
eprint = {2610.03317},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.03317},
}
RIS
TY - UNPB
TI - Fully mixed finite element methods for the coupling of viscoelasticity and reaction-diffusion models
AU - Alonso J. Bustos
AU - Jeonghun J. Lee
AU - Ricardo Ruiz-Baier
PY - 2026
DA - 2026/10/02
PB - arXiv
DO - 10.48550/arxiv.2610.03317
UR - https://arxiv.org/abs/2610.03317
AB - We propose and analyse a fully mixed finite element method for a two-way coupled mechanochemical model of calcium signalling in viscoelastic tissue. The mechanical response follows a Kelvin--Voigt law, with the elastic and viscous stresses, velocity, and spin tensor as unknowns, and the symmetry of the total stress imposed weakly. The reaction-diffusion system for calcium and inactivated inositol triphosphate is also formulated in mixed form, introducing the diffusive flux as an additional unknown. The coupling is bidirectional: calcium generates an active stress, while local dilation feeds back into the calcium dynamics. We establish global existence and uniqueness by reducing the coupled system to an integro-differential problem and subsequently recovering the mixed variables. The semidiscrete scheme is formulated for general conforming finite element spaces and shown to be well posed under a single inf-sup condition. An interpolation operator that commutes with the divergence and preserves the weak symmetry pairing yields optimal a priori error estimates. We specialise the analysis to three families of weakly symmetric finite elements, deriving explicit convergence rates and, for one family, a superconvergent local post-processing of the displacement. Numerical tests in two and three dimensions confirm the predicted convergence rates and illustrate the method through viscoelastic benchmarks and mechanically driven calcium wavetrains.
ER -