Accurate numerical simulation of electrodiffusion and water movement in brain tissue

02/04/2021
by   Ada J. Ellingsrud, et al.
0

Mathematical modelling of ionic electrodiffusion and water movement is emerging as a powerful avenue of investigation to provide new physiological insight into brain homeostasis. However, in order to provide solid answers and resolve controversies, the accuracy of the predictions is essential. Ionic electrodiffusion models typically comprise non-trivial systems of non-linear and highly coupled partial and ordinary differential equations that govern phenomena on disparate time scales. Here, we study numerical challenges related to approximating these systems. We consider a homogenized model for electrodiffusion and osmosis in brain tissue and present and evaluate different associated finite element-based splitting schemes in terms of their numerical properties, including accuracy, convergence, and computational efficiency for both idealized scenarios and for the physiologically relevant setting of cortical spreading depression (CSD). We find that the schemes display optimal convergence rates in space for problems with smooth manufactured solutions. However, the physiological CSD setting is challenging: we find that the accurate computation of CSD wave characteristics (wave speed and wave width) requires a very fine spatial and fine temporal resolution.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
11/14/2019

A mixture theory-based finite element formulation for the study of biodegradation of poroelastic scaffolds

We derive a mixture theory-based mathematical model of the degradation o...
research
03/20/2020

Uniform-in-time convergence of numerical methods for non-linear degenerate parabolic equations

Gradient schemes is a framework that enables the unified convergence ana...
research
04/21/2023

Experimental Convergence Rate Study for Three Shock-Capturing Schemes and Development of Highly Accurate Combined Schemes

We study experimental convergence rates of three shock-capturing schemes...
research
07/01/2019

On the development of symmetry-preserving finite element schemes for ordinary differential equations

In this paper we introduce a procedure, based on the method of equivaria...
research
12/13/2021

A flux-vector splitting scheme for the shallow water equations extended to high-order on unstructured meshes

We present an advection-pressure flux-vector splitting method for the on...
research
03/08/2022

Computational Benchmarks with Optimal Multilevel Argyris FEM

The main drawback for the application of the conforming Argyris FEM is t...

Please sign up or login with your details

Forgot password? Click here to reset