Open Access Open Access  Restricted Access Access granted  Restricted Access Subscription Access

Vol 28, No 3 (2017)

I. Mathematical Modeling

Numerical Method of Determining the Excitation Source for the Fitzhugh–Nagumo Mathematical Model

Tuikina S.R., Solov’eva S.I.

Abstract

We consider the inverse problem for the two-dimensional modified FitzHugh–Nagumo model in the presence of an infarct. The inverse problem determines the myocardium excitation source function (a function of space variables and time) from a system of partial differential equations. Additional dynamic measurements of the potential are carried out on the entire inside boundary of the region representing a cross-section of the heart and its ventricles by a horizontal plane, which fits the real heart geometry. A numerical method is proposed for the solution of this inverse problem with the discrepancy functional; numerical results are reported.

Computational Mathematics and Modeling. 2017;28(3):301-309
pages 301-309 views

Article

Spectral Method to Investigate the Dependence of MHD Stability of an Electrolysis Bath on the Shape of the Skull

Savenkova N.P., Mokin A.Y., Il’yutko V.P.

Abstract

We consider the dependence of the MHD stability of an electrolysis bath on the shape of the work space. As the optimal work-space shape we choose the one that achieves the best separation of the eigenvalues in the spectrum of the multidimensional problem posed for the kinematic equation for the electrolytealuminum interface in a particular electrolysis bath.

Computational Mathematics and Modeling. 2017;28(3):310-315
pages 310-315 views

Optimal Resource Allocation in a Two-Sector Economic Model with an Integral Functional

Kiselev Y.N., Avvakumov S.N., Orlov M.V., Orlov S.M.

Abstract

We investigate the resource allocation problem in a two-sector economic model with a Cobb-Douglas production function with different depreciation rates. The problem is considered on a finite time horizon with an integral type functional. Optimality of the extremum solution constructed by the Pontryagin maximum principle is established. When the planning horizon is sufficiently long, the optimal control has two or three switching points, contains one singular section, and vanishes on the terminal section. A transitional “calibration” regime exists between the singular section, where the motion is along a singular ray, and the terminal section. The solution of the maximum-principle boundary-value problem is presented in explicit form, accompanied by graphs based on numerical results.

Computational Mathematics and Modeling. 2017;28(3):316-338
pages 316-338 views

Analysis of Pulse Wave Similarity in Photoplethysmograms

Golovina A.M., D’yakonov A.G., Kharatsidi O.A.

Abstract

The article focuses on the analysis of photoplethysmograms. The quality of a photoplethysmogram is assessed by solving the classification problem, i.e., identifying the patient from the photoplethysmogram. Efficient patient identification methods are proposed, including methods based on DTW and TWED metrics. Identification accuracy reaches 70%.

Computational Mathematics and Modeling. 2017;28(3):339-349
pages 339-349 views

Numerical Investigation of Transient Free Convective Flow in Vertical Channel Filled with Porous Material in the Presence of Thermal Dispersion

Jha B.K., Aina B.

Abstract

The present work consists of a numerical investigation of transient free convective flow in vertical channel formed by two infinite vertical parallel plates filled with porous material in the presence of thermal dispersion. The governing coupled-nonlinear equations of momentum and energy transport are solved numerically using the implicit finite difference method, while the approximate analytical solution is also presented to find the expression for velocity, temperature, skin friction, and rate of heat transfer for the steady fully developed flow using the perturbation technique. The main objective is to investigate the effects of the dimensionless time, Darcy number, thermal dispersion, and Prandtl number on the fluid flow and heat transfer characteristics. Solutions are presented in graphical form and given in terms of fluid velocity, fluid temperature, skin friction, and rate of heat transfer for various parametric values. The significant result from this study is that velocity and temperature is enhanced with increase in thermal dispersion parameter and time. Furthermore, excellent agreement is found between the steady-state solution and the transient solution at large values of time.

Computational Mathematics and Modeling. 2017;28(3):350-367
pages 350-367 views

Mathematical Models of Investment Cycles

Kurkina E.S.

Abstract

We propose and investigate three mathematical models describing production cycles. They incorporate various mechanisms of endogenous fluctuations in economic systems. The models are based on ODE systems. The first model is a Keynesian IS-LM model of business cycles. The interest rate is determined by the money market and influences the relationship between savings and investments, allowing funds to flow from one to the other and vice versa. In the second case, the fluctuation mechanism is associated with time lags between investment growth, capital growth, and rate of return on capital. As a result, the economy periodically “overheats”, as rapid growth of capital suppresses the return rates, production becomes unprofitable, and investments sharply decline. Two models realizing this mechanism are proposed. One is a minimalist model based on a system of three ODEs. The other is an augmented model that sufficiently fully describes modern economic systems of developed countries and consists of nine ODEs and nine algebraic equations. It encompasses all the principal markets: labor market, capital market, financial market, and commodity market. Bifurcation analysis of the three models is carried out, oscillation regions are determined, and oscillation mechanisms are examined in detail. The model parameters are chosen so that the cycle periods are 12–17 years long.

Computational Mathematics and Modeling. 2017;28(3):377-399
pages 377-399 views

Methods to Increase Fault Tolerance of Combinational Integrated Microcircuits by Redundancy Coding

Gavrilov S.V., Gurov S.I., Zhukova T.D., Rukhlov V.S., Ryzhova D.I., Tel’pukhov D.V.

Abstract

Increasing the operating reliability of integrated microcircuits (IMC) remains, on the whole, an unsolved design problem. An important aspect of this problem is the stability of the circuits under transient faults (malfunctions) in large integrated circuits. Faults appear due to various disturbances: radiation, supply voltage jumps, signal degradation over time, etc. Investigations show that the probability of an error due to these factors may vary between very wide limits: from less than 0.1% for large circuits and up to 30% for very small circuits. In this article, we consider various methods of enhancing the fault tolerance of combinational circuits and also assess the effect of a single fault and a stuck-at fault on circuit operation for the case of combinational circuits from the ISCAS’85 set.

Computational Mathematics and Modeling. 2017;28(3):400-406
pages 400-406 views

Convexity/Concavity and Stability Aspects of Rational Cubic Fractal Interpolation Surfaces

Chand A.K., Vijender N., Navascués M.A.

Abstract

Fractal interpolation is more general than the classical piecewise interpolation due to the presence of the scaling factors that describe smooth or non-smooth shape of a fractal curve/surface. We develop the rational cubic fractal interpolation surfaces (FISs) by using the blending functions and rational cubic fractal interpolation functions (FIFs) with two shape parameters in each sub-interval along the grid lines of the interpolation domain. The properties of blending functions and C1-smoothness of rational cubic FIFs render C1-smoothness to our rational cubic FISs. We study the stability aspects of the rational cubic FIS with respect to its independent variables, dependent variable, and first order partial derivatives at the grids. The scaling factors and shape parameters seeded in the rational cubic FIFs are constrained so that these rational cubic FIFs are convex/concave whenever the univariate data sets along the grid lines are convex/concave. For these constrained scaling factors and shape parameters, our rational cubic FIS is convex/concave to given convex/concave surface data.

Computational Mathematics and Modeling. 2017;28(3):407-430
pages 407-430 views

Numerical Simulation of Singularly Perturbed Reaction-Diffusion Equation Using Finite Element Method

Srivastava A.

Abstract

This article deals with the study of sign-changing solutions of the nonlinear singularly perturbed reaction-diffusion equation. Sign changing solutions of the nonlinear problem do not appear to have been previously studied in detail computationally, and it is hoped that this paper will help to provide a new idea in this direction. A variant of Newton’s method having tenth order of convergence has been established to linearize the nonlinear system of equations. Examples of the nonlinear problem having nonlinearities in homogeneous/nonhomogeneous form are considered to show the existence of solutions.

Computational Mathematics and Modeling. 2017;28(3):431-447
pages 431-447 views

II. Informatics

A Bound on the Probability of Ruin in Merton’s Model

Morozov V.V., Babin V.A.

Abstract

We consider a modified Merton’s model of optimal consumption that allows for the utility of continuous and terminal consumption. An explicit solution of the Hamilton-Jacoby-Bellman equation is found. An upper bound is constructed on the probability of an event involving either investor ruin or negative consumption.

Computational Mathematics and Modeling. 2017;28(3):368-376
pages 368-376 views