Ricardo Ruiz Baier

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.

Raimund Bürger, Ricardo Ruiz Baier, Mauricio Sepúlveda et al.

We present a fully adaptive multiresolution scheme for spatially one-dimensional quasilinear strongly degenerate parabolic equations with zero-flux and periodic boundary conditions. The numerical scheme is based on a finite volume discretization using the Engquist-Osher numerical flux and explicit time stepping. An adaptive multiresolution scheme based on cell averages is then used to speed up the CPU time and the memory requirements of the underlying finite volume scheme, whose first-order version is known to converge to an entropy solution of the problem. A particular feature of the method is the storage of the multiresolution representation of the solution in a graded tree, whose leaves are the non-uniform finite volumes on which the numerical divergence is eventually evaluated. Moreover using the $L^1$ contraction of the discrete time evolution operator we derive the optimal choice of the threshold in the adaptive multiresolution method. Numerical examples illustrate the computational efficiency together with the convergence properties.

Raimund Bürger, Ricardo Ruiz Baier, Mauricio Sepúlveda et al.

A fully adaptive finite volume multiresolution scheme for one-dimensional strongly degenerate parabolic equations with discontinuous flux is presented. The numerical scheme is based on a finite volume discretization using the Engquist--Osher approximation for the flux and explicit time--stepping. An adaptivemultiresolution scheme with cell averages is then used to speed up CPU time and meet memory requirements. A particular feature of our scheme is the storage of the multiresolution representation of the solution in a dynamic graded tree, for the sake of data compression and to facilitate navigation. Applications to traffic flow with driver reaction and a clarifier--thickener model illustrate the efficiency of this method.

Luis Miguel De Oliveira Vilaca, Michel C. Milinkovitch, Ricardo Ruiz Baier

We introduce a model for the mass transfer of molecular activators and inhibitors in two media separated by an interface, and study its interaction with the deformations exhibited by the two-layer skin tissue where they occur. The mathematical model results in a system of nonlinear advection-diffusion-reaction equations including cross-diffusion, and coupled with an interface elasticity problem. We propose a Galerkin method for the discretisation of the set of governing equations, involving also a suitable Newton linearisation, partitioned techniques, non-overlapping Schwarz alternating schemes, and high-order adaptive time stepping algorithms. The experimental accuracy and robustness of the proposed partitioned numerical methods is assessed, and some illustrating tests in 2D and 3D are provided to exemplify the coupling effects between the mechanical properties and the advection-diffusion-reaction interactions involving the two separate layers.

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.

Veronica Anaya, Afaf Bouharguane, David Mora et al.

We introduce a family of mixed methods and discontinuous Galerkin discretisations designed to numerically solve the Oseen equations written in terms of velocity, vorticity, and Bernoulli pressure. The unique solvability of the continuous problem is addressed by invoking a global inf-sup property in an adequate abstract setting for non-symmetric systems. The proposed finite element schemes, which produce exactly divergence-free discrete velocities, are shown to be well-defined and optimal convergence rates are derived in suitable norms. In addition, we establish optimal rates of convergence for a class of discontinuous Galerkin schemes, which employ stabilisation. A set of numerical examples serves to illustrate salient features of these methods.

Ricardo Ruiz Baier

A numerical method is presented to obtain approximate solutions to problems arising from sedimentation models. These processes are widely utilized in minery for recovering water from suspensions coming out of flotation processes. The main idea is to apply a multiresolution method to the existing schemes developed by Bürger et al. [2, 3, 4] and to observe the good performance of the multiresolution strategy when applied to these kind of problems. We obtain high rates of memory compression without affecting the quality of the solution.

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.

Ricardo Ruiz Baier

The main objective of this dissertation is to present an adaptation of some finite volume methods used in the resolution of problems arising in sedimentation processes of flocculated suspensions (or sedimentation with compression). This adaptation is based on the utilization of multiresolution techniques, originally designed to reduce the computational cost incurred in solving using high resolution schemes in the numerical solution of hyperbolic systems of conservation laws.