Mostafa Bendahmane

Boris Andreianov, Mostafa Bendahmane, Kenneth H. Karlsen

We consider a class of doubly nonlinear degenerate hyperbolic-parabolic equations with homogeneous Dirichlet boundary conditions, for which we first establish the existence and uniqueness of entropy solutions. We then turn to the construction and analysis of discrete duality finite volume schemes (in the spirit of Domelevo and Omnès \cite{DomOmnes}) for these problems in two and three spatial dimensions. We derive a series of discrete duality formulas and entropy dissipation inequalities for the schemes. We establish the existence of solutions to the discrete problems, and prove that sequences of approximate solutions generated by the discrete duality finite volume schemes converge strongly to the entropy solution of the continuous problem. The proof revolves around some basic a priori estimates, the discrete duality features, Minty-Browder type arguments, and "hyperbolic" $L^\infty$ weak-$\star$ compactness arguments (i.e., propagation of compactness along the lines of Tartar, DiPerna, ...). Our results cover the case of non-Lipschitz nonlinearities.

Boris Andreianov, Mostafa Bendahmane, Kenneth H. Karlsen et al.

We prove convergence of discrete duality finite volume (DDFV) schemes on distorted meshes for a class of simplified macroscopic bidomain models of the electrical activity in the heart. Both time-implicit and linearised time-implicit schemes are treated. A short description is given of the 3D DDFV meshes and of some of the associated discrete calculus tools. Several numerical tests are presented.

Mostafa Bendahmane, Raimund Bürger, Ricardo Ruiz Baier et al.

We present a fully adaptive multiresolution scheme for spatially two-dimensional, possibly degenerate reaction-diffusion systems, focusing on combustion models and models of pattern formation and chemotaxis in mathematical biology. Solutions of these equations in these applications exhibit steep gradients, and in the degenerate case, sharp fronts and discontinuities. The multiresolution scheme is based on finite volume discretizations with explicit time stepping. The multiresolution representation of the solution is stored in a graded tree. By a thresholding procedure, namely the elimination of leaves that are smaller than a threshold value, substantial data compression and CPU time reduction is attained. The threshold value is chosen optimally, in the sense that the total error of the adaptive scheme is of the same slope as that of the reference finite volume scheme. Since chemical reactions involve a large range of temporal scales, but are spatially well localized (especially in the combustion model), a locally varying adaptive time stepping strategy is applied. It turns out that local time stepping accelerates the adaptive multiresolution method by a factor of two, while the error remains controlled.

Mostafa Bendahmane, Raimund Bürger, Ricardo Ruiz Baier

This work deals with the numerical solution of the monodomain and bidomain models of electrical activity of myocardial tissue. The bidomain model is a system consisting of a possibly degenerate parabolic PDE coupled with an elliptic PDE for the transmembrane and extracellular potentials, respectively. This system of two scalar PDEs is supplemented by a time-dependent ODE modeling the evolution of the so-called gating variable. In the simpler sub-case of the monodomain model, the elliptic PDE reduces to an algebraic equation. Two simple models for the membrane and ionic currents are considered, the Mitchell-Schaeffer model and the simpler FitzHugh-Nagumo model. Since typical solutions of the bidomain and monodomain models exhibit wavefronts with steep gradients, we propose a finite volume scheme enriched by a fully adaptive multiresolution method, whose basic purpose is to concentrate computational effort on zones of strong variation of the solution. Time adaptivity is achieved by two alternative devices, namely locally varying time stepping and a Runge-Kutta-Fehlberg-type adaptive time integration. A series of numerical examples demonstrates thatthese methods are efficient and sufficiently accurate to simulate the electrical activity in myocardial tissue with affordable effort. In addition, an optimalthreshold for discarding non-significant information in the multiresolution representation of the solution is derived, and the numerical efficiency and accuracy of the method is measured in terms of CPU time speed-up, memory compression, and errors in different norms.

Mostafa Bendahmane, Ziad Khalil, Mazen Saad

A classical model for water-gas flows in porous media is considered. The degenerate coupled system of equations obtained by mass conservation is usually approximated by finite volume schemes in the oil reservoir simulations. The convergence properties of these schemes are only known for incompressible fluids. This chapter deals with construction and convergence analysis of a finite volume scheme for compressible and immiscible flow in porous media. In comparison with incompressible fluid, compressible fluids requires more powerful techniques. We present a new result of convergence in a two or three dimensional porous medium and under the only modification that the density of gas depends on global pressure.

Verónica Anaya, Mostafa Bendahmane, David Mora et al.

We present a Virtual Element Method (VEM) for a nonlocal reaction-diffusion system of the cardiac electric field. To this system, we analyze an $H^1(Ω)$-conforming discretization by means of VEM which can make use of general polygonal meshes. Under standard assumptions on the computational domain, we establish the convergence of the discrete solution by considering a series of a priori estimates and by using a general $L^p$ compactness criterion. Moreover, we obtain optimal order space-time error estimates in the $L^2$ norm. Finally, we report some numerical tests supporting the theoretical results.

Mostafa Bendahmane, Mauricio Sepulveda

We analyze a finite volume scheme for nonlocal SIR model, which is a nonlocal reaction-diffusion system modeling an epidemic disease. We establish existence solutions to the finite volume scheme, and show that it converges to a weak solution. The convergence proof is based on deriving series of a priori estimates and using a general $L^p$ compactness criterion.

Mostafa Bendahmane, Felipe Wallison Chaves-Silva

This paper is devoted to the analysis of the uniform null controllability for a family of nonlinear reaction-diffusion systems approximating a parabolic-elliptic system which models the electrical activity of the heart. The uniform, with respect to the degenerating parameter, null controllability of the approximating system by means of a single control is shown. The proof is based on the combination of Carleman estimates and weighted energy inequalities.

Mostafa Bendahmane, Raimund Bürger, Ricardo Ruiz Baier

This paper presents a finite-volume method, together with fully adaptive multi-resolution scheme to obtain spatial adaptation, and a Runge-Kutta-Fehlberg scheme with a local time-varying step to obtain temporal adaptation, to solve numerically the known "bidominio" equations that model the electrical activity of the tissue in the myocardium. Two simple models are considered for membrane flows and ionic currents. First we define an approximate solution and we verify its convergence to the corresponding weak solution of the continuum problem, obtaining in this way an alternative demonstration that the continuum problem is well-posed. Next we introduce the multiresolution technique and derive an optimal noise reduction threshold. The efficiency and precision of our method is seen in the reduction of machine time, memory usage, and errors in comparison to other methods. ----- En este trabajo se presenta un metodo de volumenes finitos enriquecido con un esquema de multiresolucion completamente adaptativo para obtener adaptatividad espacial, y un esquema Runge-Kutta-Fehlberg con paso temporal de variacion local para obtener adaptatividad temporal, para resolver numericamente las conocidas ecuaciones "bidominio" que modelan la actividad electrica del tejido en el miocardio. Se consideran dos modelos simples para las corrientes de membrana y corrientes ionicas. En primer lugar definimos una solucion aproximada y nos referimos a su convergencia a la correspondiente solucion debil del problema continuo, obteniendo de este modo una demostracion alternativa de que el problema continuo es bien puesto. Luego de introducir la tecnica de multiresolucion, se deriva un umbral optimo para descartar la informacion no significativa, y tanto la eficiencia como la precision de nuestro metodo es vista en terminos de la aceleracion de tiempo de maquina, compresion de memoria computacional y errores en diferentes normas.