In this paper we propose a novel coupled poroelasticity-diffusion model for the formation of extracellular oedema and infectious myocarditis valid in large deformations, manifested as an interaction between interstitial flow and the immune-driven dynamics between leukocytes and pathogens. The governing partial differential equations are formulated in terms of skeleton displacement, fluid pressure, Lagrangian porosity, and the concentrations of pathogens and leukocytes. A five-field mixed-primal finite element scheme is proposed for the numerical approximation of the problem, and we provide the stability analysis for a simplified system emanating from linearisation. We also discuss the construction of an adequate, Schur complement based, nested preconditioner. The produced computational tests exemplify the properties of the new model and of the mixed schemes.

Sequential block-partitioned solvers have in the recent past been quite popular for multi-physics and in particular poroelasticity models. Such enable tailored solver technology for the respective single-physics problems via iterative coupling, as well as suggest suitable block-preconditioners for monolithic solvers.In this talk, we focus on a thermodynamically consistent poroelasticity model recently proposed. It extends the classical quasi-static Biot equations by incoporating inertia contributions in both solid and fluid equations, aiming at biomedical applications; for instance, the perfusion of the heart.Following ideas and techniques from previous works, we present block-partitioned solvers for the fully dynamic poroelasticity model supported by theoretical convergence analysis.