Russian JouRnal of numeRical analysis and mathematical modelling provides English translations of selected new original Russian papers on the theoretical aspects of numerical analysis and the application of mathematical methods to simulation and modelling. The members of the editorial board, the most prominent Russian scientists in numerical analysis and mathematical modelling, select papers on the basis of their high scientific standard, innovative approach and topical interest.

Russian JouRnal of numeRical analysis and mathematical modelling provides English translations of selected new original Russian papers on the theoretical aspects of numerical analysis and the application of mathematical methods to simulation and modelling. The members of the editorial board, the most prominent Russian scientists in numerical analysis and mathematical modelling, select papers on the basis of their high scientific standard, innovative approach and topical interest.

Abstract This paper proposes the combination of matrix exponential method with the semi-Lagrangian approach for the time integration of shallow water equations on the sphere. The second order accuracy of the developed scheme is shown. Exponential semi-Lagrangian scheme in the combination with spatial approximation on the cubed-sphere grid is verified using the standard test problems for shallow water models. The developed scheme is as good as the conventional semi-implicit semi-Lagrangian scheme in accuracy of slowly varying flow component reproduction and significantly better in the reproduction of the fast inertia-gravity waves. The accuracy of inertia-gravity waves reproduction is close to that of the explicit time-integration scheme. The computational efficiency of the proposed exponential semi-Lagrangian scheme is somewhat lower than the efficiency of semi-implicit semi-Lagrangian scheme, but significantly higher than the efficiency of explicit, semi-implicit, and exponential Eulerian schemes.

Abstract The paper discusses a stabilization of a finite element method for the equations of fluid motion in a time-dependent domain. After experimental convergence analysis, the method is applied to simulate a blood flow in the right ventricle of a post-surgery patient with the transposition of the great arteries disorder. The flow domain is reconstructed from a sequence of 4D CT images. The corresponding segmentation and triangulation algorithms are also addressed in brief.

Abstract The low-voltage cardioversion-defibrillation is a modern sparing electrotherapy method for such dangerous heart arrhythmias as paroxysmal tachycardia and fibrillation. In an excitable medium, such arrhythmias relate to appearance of spiral waves of electrical excitation, and the spiral waves are superseded to the electric boundary of the medium in the process of treatment due to high-frequency stimulation from the electrode. In this paper we consider the Aliev–Panfilov myocardial model, which provides a positive tension of three-dimensional scroll waves, and an axisymmetric model of the left ventricle of the human heart. Two relations of anisotropy are considered, namely, isotropy and physiological anisotropy. The periods of stimulation with an apical electrode are found so that the electrode successfully entrains its rhythm in the medium, the spiral wave is superseded to the base of the ventricle, and disappears. The results are compared in two-dimensional and three-dimensional media. The intervals of effective stimulation periods are sufficiently close to each other in the two-dimensional case and in the anatomical model. However, the use of the anatomical model is essential in determination of the time of superseding.

Russian JouRnal of numeRical analysis and mathematical modelling provides English translations of selected new original Russian papers on the theoretical aspects of numerical analysis and the application of mathematical methods to simulation and modelling.The members of the editorial board, the most prominent Russian scientists in numerical analysis and mathematical modelling, select papers on the basis of their high scientific standard, innovative approach and topical interest.All information regarding subscriptions, Open access, back volumes and orders is available online at www.degruyter.com/rjnamm.

Abstract For elliptic boundary value problems (the diffusion equation and elasticity theory ones) with highly varying coefficients, there are proposed iterative methods with the number of iterations independent of the coefficient jumps. In the differential case these methods take solving the boundary value problem for the Poisson equation at each step of iterations while in the finite difference (finite element) approximation it is possible to use another operator as a preconditioner.

Abstract This paper is a continuation of [38]. The analysis of the modified partial differential equation (MDE) of the constant-wind-speed linear advection equation explicit difference scheme up to the eighth-order derivatives is presented. In this paper the authors focus on the dissipative features of the Beam–Warming scheme. The modified partial differential equation is presented in the so-calledΠ-form of the first differential approximation. The most important part of this form includes the coefficientsμ(p) at the space derivatives. Analysis of these coefficients provides indications of the nature of the dissipative errors. A fragment of the stencil for determining the modified differential equation for the Beam–Warming scheme is included. The derived and presented coefficientsμ(p) as well as the analysis of the dissipative features of this scheme on the basis of these coefficients have not been published so far.

Abstract A new method is developed to calculate characteristics of contaminant transport (including non-classical regimes) in statistically homogeneous sharply contrasting media. A transport integro-differential equation in the space-time representation is formulated on the basis of the model earlier proposed by one of the authors (L. M.). Analytical expressions for transport characteristics in limiting time intervals in the one-dimensional case are derived. An interpolation form is proposed for the integral kernel of the transport equation. On a basis of this expression, an algorithm is developed for numerical modelling the contaminant transport in statistically homogeneous sharply contrasting media. Trial numerical 1D calculations are performed based on this algorithm. Good agreement was found between the numerical simulation results and the asymptotic analytical expressions.