scholarly journals Optimal monodomain approximations of the bidomain equations used in cardiac electrophysiology

2014 ◽  
Vol 24 (06) ◽  
pp. 1115-1140 ◽  
Author(s):  
Yves Coudière ◽  
Yves Bourgault ◽  
Myriam Rioux
Keyword(s):  
Boundary Conditions ◽  
Cardiac Electrophysiology ◽  
Propagation Speed ◽  
Neumann Boundary ◽  
Test Cases ◽  
Bidomain Model ◽  
Bidomain Equations ◽  
Monodomain Model ◽  
The Difference ◽  
Computationally Intensive

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.

scholarly journals Mathematical analysis and 2-scale convergence of a heterogeneous microscopic bidomain model

2018 ◽  
Vol 28 (05) ◽  
pp. 979-1035 ◽  
Author(s):  
Annabelle Collin ◽  
Sébastien Imperiale
Keyword(s):  
Mathematical Analysis ◽  
Cardiac Electrophysiology ◽  
Convergence Theory ◽  
Bidomain Model ◽  
Periodic Homogenization ◽  
Homogenization Procedure ◽  
Uniform Estimates ◽  
Bidomain Equations ◽  
Homogenization Process

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.

Numerical solution of the bidomain equations

2009 ◽  
Vol 367 (1895) ◽  
pp. 1931-1950 ◽  
Author(s):  
S. Linge ◽  
J. Sundnes ◽  
M. Hanslien ◽  
G.T. Lines ◽  
A. Tveito
Keyword(s):  
Mathematical Models ◽  
Cardiac Electrophysiology ◽  
Mathematical Reasoning ◽  
Full Potential ◽  
Processing Unit ◽  
Bidomain Model ◽  
Central Processing ◽  
Cellular Models ◽  
Bidomain Equations ◽  
The Past

Knowledge of cardiac electrophysiology is efficiently formulated in terms of mathematical models. However, most of these models are very complex and thus defeat direct mathematical reasoning founded on classical and analytical considerations. This is particularly so for the celebrated bidomain model that was developed almost 40 years ago for the concurrent analysis of extra- and intracellular electrical activity. Numerical simulations based on this model represent an indispensable tool for studying electrophysiology. However, complex mathematical models, steep gradients in the solutions and complicated geometries lead to extremely challenging computational problems. The greatest achievement in scientific computing over the past 50 years has been to enable the solving of linear systems of algebraic equations that arise from discretizations of partial differential equations in an optimal manner, i.e. such that the central processing unit (CPU) effort increases linearly with the number of computational nodes. Over the past decade, such optimal methods have been introduced in the simulation of electrophysiology. This development, together with the development of affordable parallel computers, has enabled the solution of the bidomain model combined with accurate cellular models, on geometries resembling a human heart. However, in spite of recent progress, the full potential of modern computational methods has yet to be exploited for the solution of the bidomain model. This paper reviews the development of numerical methods for solving the bidomain model. However, the field is huge and we thus restrict our focus to developments that have been made since the year 2000.

Mechanotransduction caused by a point force in the extracellular space

BIOMATH ◽  
2018 ◽  
Vol 7 (2) ◽  
pp. 1810197
Author(s):  
Bradley John Roth
Keyword(s):  
Extracellular Matrix ◽  
Analytical Solutions ◽  
Biological Tissue ◽  
Mathematical Description ◽  
Magnetic Nanoparticle ◽  
Biological Effects ◽  
Point Force ◽  
Bidomain Model ◽  
Bidomain Equations ◽  
The Difference

The mechanical bidomain model is a mathematical description of biological tissue that focuses on mechanotransduction. The model’s fundamental hypothesis is that differences in the intracellular and extracellular displacements activate integrins, causing a cascade of biological effects. This paper presents analytical solutions of the bidomain equations for an extracellular point force. The intra- and extracellular spaces are incompressible, isotropic, and coupled. The expressions for the intra- and extracellular displacements each contain three terms: a monodomain term that is identical in the two spaces, and two bidomain terms, one of which decays exponentially. Near the origin the intracellular displacement remains finite and the extracellular displacement diverges. Far from the origin the monodomain displacement decays in inverse proportion to the distance, the strain decays as the distance squared, and the difference between the intra- and extracellular displacements decays as the distance cubed. These predictions could be tested by applying a force to a magnetic nanoparticle embedded in the extracellular matrix and recording the mechanotransduction response.

scholarly journals Free energy and defect C-theorem in free scalar theory

2021 ◽  
Vol 2021 (5) ◽  
Author(s):  
Tatsuma Nishioka ◽  
Yoshiki Sato
Keyword(s):  
Boundary Conditions ◽  
Dirichlet Boundary ◽  
Neumann Boundary ◽  
Flat Space ◽  
Neumann Boundary Conditions ◽  
Dirichlet Boundary Conditions ◽  
Scalar Theory ◽  
Free Scalar ◽  
The Difference ◽  
Recent Classification

Abstract We describe conformal defects of p dimensions in a free scalar theory on a d-dimensional flat space as boundary conditions on the conformally flat space ℍp+1× 𝕊d−p−1. We classify two types of boundary conditions, Dirichlet type and Neumann type, on the boundary of the subspace ℍp+1 which correspond to the types of conformal defects in the free scalar theory. We find Dirichlet boundary conditions always exist while Neumann boundary conditions are allowed only for defects of lower codimensions. Our results match with a recent classification of the non-monodromy defects, showing Neumann boundary conditions are associated with non-trivial defects. We check this observation by calculating the difference of the free energies on ℍp+1× 𝕊d−p−1 between Dirichlet and Neumann boundary conditions. We also examine the defect RG flows from Neumann to Dirichlet boundary conditions and provide more support for a conjectured C-theorem in defect CFTs.

Study of fluid flow through a channel deformed by piezoelement

2018 ◽  
Vol 13 (3) ◽  
pp. 1-10 ◽  
Author(s):  
I.Sh. Nasibullayev ◽  
E.Sh Nasibullaeva ◽  
O.V. Darintsev
Keyword(s):  
Fluid Flow ◽  
Pressure Drop ◽  
Boundary Conditions ◽  
Flow Rate ◽  
Contact Surface ◽  
Piezoelectric Element ◽  
Neumann Boundary ◽  
Differential Pressure ◽  
Physical Parameters ◽  
Flow Through

The flow of a liquid through a tube deformed by a piezoelectric cell under a harmonic law is studied in this paper. Linear deformations are compared for the Dirichlet and Neumann boundary conditions on the contact surface of the tube and piezoelectric element. The flow of fluid through a deformed channel for two flow regimes is investigated: in a tube with one closed end due to deformation of the tube; for a tube with two open ends due to deformation of the tube and the differential pressure applied to the channel. The flow rate of the liquid is calculated as a function of the frequency of the deformations, the pressure drop and the physical parameters of the liquid.

scholarly journals Blow-up analyses in nonlocal reaction diffusion equations with time-dependent coefficients under Neumann boundary conditions

2020 ◽  
Vol 18 (1) ◽  
pp. 1552-1564
Author(s):  
Huimin Tian ◽  
Lingling Zhang
Keyword(s):  
Boundary Conditions ◽  
Blow Up ◽  
Sufficient Conditions ◽  
Neumann Boundary ◽  
Reaction Diffusion ◽  
Neumann Boundary Conditions ◽  
Reaction Diffusion Equations ◽  
Time Dependent ◽  
Diffusion Equations ◽  
Nonlocal Reaction

Abstract In this paper, the blow-up analyses in nonlocal reaction diffusion equations with time-dependent coefficients are investigated under Neumann boundary conditions. By constructing some suitable auxiliary functions and using differential inequality techniques, we show some sufficient conditions to ensure that the solution u ( x , t ) u(x,t) blows up at a finite time under appropriate measure sense. Furthermore, an upper and a lower bound on blow-up time are derived under some appropriate assumptions. At last, two examples are presented to illustrate the application of our main results.

scholarly journals Mode Shape-Based Damage Detection Method (MSDI): Experimental Validation

2021 ◽  
Vol 11 (10) ◽  
pp. 4589
Author(s):  
Ivan Duvnjak ◽  
Domagoj Damjanović ◽  
Marko Bartolac ◽  
Ana Skender
Keyword(s):  
Boundary Conditions ◽  
Damage Detection ◽  
Mode Shape ◽  
Experimental Validation ◽  
Detection Method ◽  
Dynamic Properties ◽  
Damage Index ◽  
Main Principle ◽  
The Difference ◽  
Different Levels

The main principle of vibration-based damage detection in structures is to interpret the changes in dynamic properties of the structure as indicators of damage. In this study, the mode shape damage index (MSDI) method was used to identify discrete damages in plate-like structures. This damage index is based on the difference between modified modal displacements in the undamaged and damaged state of the structure. In order to assess the advantages and limitations of the proposed algorithm, we performed experimental modal analysis on a reinforced concrete (RC) plate under 10 different damage cases. The MSDI values were calculated through considering single and/or multiple damage locations, different levels of damage, and boundary conditions. The experimental results confirmed that the MSDI method can be used to detect the existence of damage, identify single and/or multiple damage locations, and estimate damage severity in the case of single discrete damage.

Recovering differential operators with two constant delays under Dirichlet/Neumann boundary conditions

2020 ◽  
Vol 28 (2) ◽  
pp. 237-241
Author(s):  
Biljana M. Vojvodic ◽  
Vladimir M. Vladicic
Keyword(s):  
Boundary Conditions ◽  
Boundary Value Problems ◽  
Differential Operators ◽  
Boundary Value ◽  
Neumann Boundary ◽  
Neumann Boundary Conditions

AbstractThis paper deals with non-self-adjoint differential operators with two constant delays generated by {-y^{\prime\prime}+q_{1}(x)y(x-\tau_{1})+(-1)^{i}q_{2}(x)y(x-\tau_{2})}, where {\frac{\pi}{3}\leq\tau_{2}<\frac{\pi}{2}<2\tau_{2}\leq\tau_{1}<\pi} and potentials {q_{j}} are real-valued functions, {q_{j}\in L^{2}[0,\pi]}. We will prove that the delays and the potentials are uniquely determined from the spectra of four boundary value problems: two of them under boundary conditions {y(0)=y(\pi)=0} and the remaining two under boundary conditions {y(0)=y^{\prime}(\pi)=0}.

On Fourier Series in Eigenfunctions of Elliptic Boundary Value Problems

2003 ◽  
Vol 10 (3) ◽  
pp. 401-410
Author(s):  
M. S. Agranovich ◽  
B. A. Amosov
Keyword(s):  
Fourier Series ◽  
Boundary Conditions ◽  
Fourier Coefficients ◽  
Elliptic Boundary ◽  
A Priori ◽  
Adjoint Problem ◽  
Elliptic Boundary Value ◽  
Right Hand ◽  
The Difference ◽  
The Right

Abstract We consider a general elliptic formally self-adjoint problem in a bounded domain with homogeneous boundary conditions under the assumption that the boundary and coefficients are infinitely smooth. The operator in 𝐿2(Ω) corresponding to this problem has an orthonormal basis {𝑢𝑙} of eigenfunctions, which are infinitely smooth in . However, the system {𝑢𝑙} is not a basis in Sobolev spaces 𝐻𝑡 (Ω) of high order. We note and discuss the following possibility: for an arbitrarily large 𝑡, for each function 𝑢 ∈ 𝐻𝑡 (Ω) one can explicitly construct a function 𝑢0 ∈ 𝐻𝑡 (Ω) such that the Fourier series of the difference 𝑢 – 𝑢0 in the functions 𝑢𝑙 converges to this difference in 𝐻𝑡 (Ω). Moreover, the function 𝑢(𝑥) is viewed as a solution of the corresponding nonhomogeneous elliptic problem and is not assumed to be known a priori; only the right-hand sides of the elliptic equation and the boundary conditions for 𝑢 are assumed to be given. These data are also sufficient for the computation of the Fourier coefficients of 𝑢 – 𝑢0. The function 𝑢0 is obtained by applying some linear operator to these right-hand sides.

Sign in / Sign up

Export Citation Format

Share Document