A comparative study of fully implicit staggered and monolithic solution methods. Part I: Coupled bidomain equations of cardiac electrophysiology

Author(s):  
Barış Cansız ◽  
Michael Kaliske
2018 ◽  
Vol 28 (05) ◽  
pp. 979-1035 ◽  
Author(s):  
Annabelle Collin ◽  
Sébastien Imperiale

The aim of this paper is to provide a complete mathematical analysis of the periodic homogenization procedure that leads to the macroscopic bidomain model in cardiac electrophysiology. We consider space-dependent and tensorial electric conductivities as well as space-dependent physiological and phenomenological nonlinear ionic models. We provide the nondimensionalization of the bidomain equations and derive uniform estimates of the solutions. The homogenization procedure is done using 2-scale convergence theory which enables us to study the behavior of the nonlinear ionic models in the homogenization process.


2005 ◽  
Vol 127 (5) ◽  
pp. 531-539 ◽  
Author(s):  
Vincent A. Mousseau

This paper will present a hybrid solution algorithm for the two-phase flow equations coupled to wall heat conduction. The partial differential equations in the physical model are the same as in RELAP5. The hybrid solution algorithm couples two solution methods, the solution method currently employed by RELAP5 and an implicitly balanced solution method. The resulting hybrid solution method is both fast and accurate. Results will be presented that show when accuracy and CPU time are considered simultaneously that there are ranges when the hybrid solution algorithm is preferred over the RELAP5 solution method.


2010 ◽  
Vol 27 (5) ◽  
pp. 1290-1312 ◽  
Author(s):  
Trygve K. Nilssen ◽  
Gunnar A. Staff ◽  
Kent-Andre Mardal

2014 ◽  
Vol 24 (06) ◽  
pp. 1115-1140 ◽  
Author(s):  
Yves Coudière ◽  
Yves Bourgault ◽  
Myriam Rioux

The bidomain model is the current most sophisticated model used in cardiac electrophysiology. The monodomain model is a simplification of the bidomain model that is less computationally intensive but only valid under equal conductivity ratio. We propose in this paper optimal monodomain approximations of the bidomain model. We first prove that the error between the bidomain and monodomain solutions is bounded by the error ‖B - A‖ between the bidomain and monodomain conductivity operators. Optimal monodomain approximations are defined by minimizing the distance ‖B - A‖, which reduces for solutions over all ℝd to minimize the Lp norm of the difference between the operator symbols. Similarly, comparing the symbols pointwise amounts to compare the propagation of planar waves in the bidomain and monodomain models. We prove that any monodomain model properly propagates at least d planar waves in ℝd. We next consider and solve the optimal problem in the L∞ and L2 norms, the former providing minimal propagation error uniformly over all directions. The quality of these optimal monodomain approximations is compared among themselves and with other published approximations using two sets of test cases. The first one uses periodic boundary conditions to mimic propagation in ℝd while the second is based on a square domain with common Neumann boundary conditions. For the first test cases, we show that the error on the propagation speed is highly correlated with the error on the symbols. The second test cases show that domain boundaries control propagation directions, with only partial impact from the conductivity operator used.


2010 ◽  
Vol 62 (11) ◽  
pp. 2678-2693
Author(s):  
Mini Mathew ◽  
M. S. Mohan Kumar

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.


2019 ◽  
Vol 37 (2) ◽  
pp. 262-302
Author(s):  
Jonathan P Whiteley

Abstract Tissue level cardiac electrophysiology is usually modelled by the bidomain equations, or the monodomain simplification of the bidomain equations. One assumption made when deriving the bidomain equations is that both the intracellular and extracellular spaces are in electrical equilibrium. This assumption neglects the disturbance of this equilibrium in thin regions close to the cell membrane known as Debye layers. We first demonstrate that the governing equations at the cell, or microscale, level may be adapted to take account of these Debye layers with little additional complexity, provided the permittivity within the Debye layers satisfies certain conditions that are believed to be satisfied for biological cells. We then homogenize the microscale equations using a technique developed for an almost periodic microstructure. Cardiac tissue is usually modelled as sheets of cardiac fibres stacked on top of one another. A common assumption is that an orthogonal coordinate system can be defined at each point of cardiac tissue, where the first axis is in the fibre direction, the second axis is orthogonal to the first axis but lies in the sheet of cardiac fibres and the third axis is orthogonal to the cardiac sheet. It is assumed further that both the intracellular and extracellular conductivity tensors are diagonal with respect to these axes and that the diagonal entries of these tensors are constant across the whole tissue. Using the homogenization technique we find that this assumption is usually valid for cardiac tissue, but highlight situations where the assumption may not be valid.


2015 ◽  
Vol 2015 ◽  
pp. 1-15 ◽  
Author(s):  
Gianmauro Cuccuru ◽  
Giorgio Fotia ◽  
Fabio Maggio ◽  
James Southern

We discuss the application of the spectral element method to the monodomain and bidomain equations describing propagation of cardiac action potential. Models of cardiac electrophysiology consist of a system of partial differential equations coupled with a system of ordinary differential equations representing cell membrane dynamics. The solution of these equations requires solving multiple length scales due to the ratio of advection to diffusion that varies among the different equations. High order approximation of spectral elements provides greater flexibility in resolving multiple length scales. Furthermore, spectral elements are extremely efficient to model propagation phenomena on complex shapes using fewer degrees of freedom than its finite element equivalent (for the same level of accuracy). We illustrate a fully unstructured all-hexahedra approach implementation of the method and we apply it to the solution of full 3D monodomain and bidomain test cases. We discuss some key elements of the proposed approach on some selected benchmarks and on an anatomically based whole heart human computational model.


Sign in / Sign up

Export Citation Format

Share Document