Fully mixed virtual element schemes for a new model of steady-state poroelastic stress-assisted diffusion in the brain

We propose a fully mixed virtual element method for the numerical approximation of the coupling between linear poroelasticity equations with strong symmetry of total poroelastic stress (using the Hellinger--Reissner principle) and stress-altered solute diffusion (where diffusive flux depends on the poroelastic stress and nonlinearly on the concentration gradient). Because of the nonlinear coupling, the function spaces associated with the nonlinear diffusion sub-problem are of Banach type. To handle this structure, the solvability of both the continuous and discrete problems is established through a decoupled fixed-point strategy. The linear poroelasticity component is analysed using the theory for perturbed saddle-point problems, whereas the nonlinear diffusion problem, relies on the classical Minty--Browder theorem for monotone global operators. The existence of solutions for the fully coupled system is rigorously proven via Schauder's fixed-point theorem. Additionally, we establish rigorous a priori error estimates for the discrete scheme, successfully handling the strongly cross-coupled nonlinearities. These findings are supported by computational evidence, demonstrating that the formulation asymptotically recovers optimal convergence rates in practice. As a key contribution, both the numerical scheme and its underlying analysis prove to be robust with respect to the poromechanical parameters. Finally, several numerical examples are presented to illustrate the properties and applicability of the proposed scheme in the study of solute transport in the context of brain multiphysics.