77 research outputs found

    Cardiac electro-mechanical activity in a deforming human cardiac tissue: modeling, existence–uniqueness, finite element computation and application to multiple ischemic disease

    No full text
    In this study, the cardiac electro-mechanical model in a deforming domain is taken with the addition of mechanical feedback and stretch-activated channel current coupled with the ten Tusscher human ventricular cell level model that results in a coupled PDE–ODE system. The existence and uniqueness of such a coupled system in a deforming domain is proved. At first, the existence of a solution is proved in the deformed domain. The local existence of the solution is proved using the regularization and the Faedo–Galerkin technique. Then, the global existence is proved using the energy estimates in appropriate Banach spaces, Gronwall lemma, and the compactness procedure. The existence of the solution in an undeformed domain is proved using the lower semi-continuity of the norms. Uniqueness is proved using Young’s inequality, Gronwall lemma, and the Cauchy–Schwartz inequality. For the application purpose, this model is applied to understand the electro-mechanical activity in ischemic cardiac tissue. It also takes care of the development of active tension, conductive, convective, and ionic feedback. The Second Piola–Kirchoff stress tensor arising in Lagrangian mapping between reference and moving frames is taken as a combination of active, passive, and volumetric components. We investigated the effect of varying strength of hyperkalemia and hypoxia, in the ischemic subregions of human cardiac tissue with local multiple ischemic subregions, on the electro-mechanical activity of healthy and ischemic zones. This system is solved numerically using the Q1 finite element method in space and the implicit–explicit Euler method in time. Discontinuities arising with the modeled multiple ischemic regions are treated to the desired order of accuracy by a simple regularization technique using the interpolating polynomials. We examined the cardiac electro-mechanical activity for several cases in multiple hyperkalemic and hypoxic human cardiac tissue. We concluded that local multiple ischemic subregions severely affect the cardiac electro-mechanical activity more, in terms of action potential (v) and mechanical parameters, intracellular calcium ion concentration [Ca2+]i, active tension (TA), stretch (λ) and stretch rate (dλdt), of a healthy cell in its vicinity, compared to a single Hyperkalemic or Hypoxic subregion. The four moderate hypoxically generated ischemic subregions affect the waveform of the stretch along the fiber and the stretch rate more than a single severe ischemic subregion

    A three-step wavelet Galerkin method for parabolic and hyperbolic partial differential equations

    No full text
    A three-step wavelet Galerkin method based on Taylor series expansion in time is proposed. The scheme is third-order accurate in time and O(2 −jp) accurate in space. Unlike Taylor–Galerkin methods, the present scheme does not contain any new higher-order derivatives which makes it suitable for solving non-linear problems. The compactly supported orthogonal wavelet bases D6 developed by Daubechies are used in the Galerkin scheme. The proposed scheme is tested with both parabolic and hyperbolic partial differential equations. The numerical results indicate the versatility and effectiveness of the proposed scheme

    Wavelet Galerkin method for fourth-order multi-dimensional elliptic partial differential equations

    No full text
    We describe a wavelet Galerkin method for numerical solutions of fourth-order linear and nonlinear partial differential equations (PDEs) in 2D and 3D based on the use of Daubechies compactly supported wavelets. Two-term connection coefficients have been used to compute higher-order derivatives accurately and economically. Localization and orthogonality properties of wavelets make the global matrix sparse. In particular, these properties reduce the computational cost significantly. Linear system of equations obtained from discretized equations have been solved using GMRES iterative solver. Quasi-linearization technique has been effectively used to handle nonlinear terms arising in nonlinear biharmonic equation. To reduce the computational cost of our method, we have proposed an efficient compression algorithm. Error and stability estimates have been derived. Accuracy of the proposed method is demonstrated through various examples. </jats:p
    corecore