On the construction of solutions of terminal problems for multidimensional affine systems in quasicanonical form

We consider the construction of solutions of terminal problems for multidimensional affine systems. We show that the terminal problem for a regular system in quasicanonical form can be reduced to a boundary value problem for a system of ordinary differential equations of lower order with right-hand side depending on a vector parameter. We prove a sufficient condition for the existence of a solution of the above-mentioned boundary value problem. A method for constructing a numerical solution is developed.     

On the asymptotic convergence of a collocation method for some boundary integral equations of the first kind

In this paper we present a certain collocation method for the numerical solution of a class of boundary integral equations of the first kind with logarithmic kernel as principle part. The transformation of the boundary value problem into boundary singular integral equation of the first kind via single-layer potential is discussed. A discretization and error representation for the numerical solution of boundary integral equations has been given. Quadrature formulae have been proposed and the error arising due to the quadrature formulae used has been estimated. The convergence of the solution with respect to the proposed numerical algorithm is shown and finally some numerical results have been presented.     

Method of boundary integral equations as applied to the numerical solution of the three-dimensional Dirichlet problem for the laplace equation in a piecewise homogeneous medium

A Dirichlet problem is considered in a three-dimensional domain filled with a piecewise homogeneous medium. The uniqueness of its solution is proved. A system of Fredholm boundary integral equations of the second kind is constructed using the method of surface potentials, and a system of boundary integral equations of the first kind is derived directly from Green’s identity. A technique for the numerical solution of integral equations is proposed, and results of numerical experiments are presented.     

Numerical method of lines for parabolic functional differential equations

Initial boundary value problems for nonlinear parabolic functional differential equations are transformed by discretization in space variables into systems of ordinary functional differential equations. A comparison theorem for differential difference inequalities is proved. Sufficient conditions for the convergence of the numerical method of lines are given. An explicit Euler method is proposed for the numerical solution of systems thus obtained. This leads to difference scheme for the original problem. A complete convergence analysis for the method is given.     

On a new approach to asymptotic stabilization problems

A numerical algorithm for solving the asymptotic stabilization problem by the initial data to a fixed hyperbolic point with a given rate is proposed and justified. The stabilization problem is reduced to projecting the resolving operator of the given evolution process on a strongly stable manifold. This approach makes it possible to apply the results to a wide class of semidynamical systems including those corresponding to partial differential equations. By way of example, a numerical solution of the problem of the asymptotic stabilization of unstable trajectories of the two-dimensional Chafee-Infante equation in a circular domain by the boundary conditions is given.     

Three-level schemes of the alternating triangular method

In this paper, the schemes of the alternating triangular method are set out in the class of splitting methods used for the approximate solution of Cauchy problems for evolutionary problems. These schemes are based on splitting the problem operator into two operators that are conjugate transposes of each other. Economical schemes for the numerical solution of boundary value problems for parabolic equations are designed on the basis of an explicit-implicit splitting of the problem operator. The alternating triangular method is also of interest for the construction of numerical algorithms that solve boundary value problems for systems of partial differential equations and vector systems. The conventional schemes of the alternating triangular method used for first-order evolutionary equations are two-level ones. The approximation properties of such splitting methods can be improved by transiting to three-level schemes. Their construction is based on a general principle for improving the properties of difference schemes, namely, on the regularization principle of A.A. Samarskii. The analysis conducted in this paper is based on the general stability (or correctness) theory of operator-difference schemes.     

Boundary control problems for quasilinear systems of hyperbolic equations

For quasilinear systems of hyperbolic equations, the nonclassical boundary value problem of controlling solutions with the help of boundary conditions is considered. Previously, this problem was extensively studied in the case of the simplest hyperbolic equations, namely, the scalar wave equation and certain linear systems. The corresponding problem formulations and numerical solution algorithms are extended to nonlinear (quasilinear and conservative) systems of hyperbolic equations. Some numerical (grid-characteristic) methods are considered that were previously used to solve the above problems. They include explicit and implicit conservative difference schemes on compact stencils that are linearizations of Godunov’s method. The numerical algorithms and methods are tested as applied to well-known linear examples.     

A numerical method for a system of equations with a small parameter and a point source on an infinite interval

A method for solving a boundary-value problem on an infinite interval is considered for a linear system of second-order ordinary differential equations with a small parameter at the highest derivatives and a point source. The question is addressed of reduction of this problem to a finite interval. A mesh, condensing in the boundary layer, is used for numerical solution of a system of singularly perturbed equations on a finite interval.     

Positive solutions and dynamics of a finite difference reaction–diffusion system

We investigate the dynamics and methods of computation for some nonlinear finite difference systems that are the discretized equations of a time-dependent and a steady-state reaction–diffusion problem. The formulation of the discrete equations for the time-dependent problem is based on the implicit method for parabolic equations, and the computational algorithm is based on the method of monotone iterations using upper and lower solutions as the initial iterations. The monotone iterative method yields improved upper and lower bounds of the solution in each iteration, and the sequence of iterations converges monotonically to a solution for both the time-dependent and the steady-state problems. An important consequence of this method is that it leads to a bifurcation point that determines the dynamic behavior of the time-dependent problem in relation to the corresponding steady-state problem. This bifurcation point also determines whether the steady-state problem has one or two non-negative solutions, and is explicitly given in terms of the physical parameters of the system and the type of boundary conditions. Numerical results are presented for both the time-dependent and the steady-state problems under various boundary conditions, including a test problem with known analytical solution. These numerical results exhibit the predicted dynamic behavior of the time-dependent solution given by the theoretical analysis. Also discussed are the numerical stability of the computational algorithm and the convergence of the finite difference solution to the corresponding continuous solution of the reaction–diffusion problem. © 1993 John Wiley & Sons, Inc.     

Infinite Systems of Quasilinear Differential Difference Inequalities and Applications

The article deals with the initial boundary value problem for an infinite system of first order quasilinear functional differential equations. A comparison result concerning infinite systems of differential difference inequalities is proved. A function satisfying such inequalities is estimated by a solution of a suitable Cauchy problem for an ordinary functional differential system. The comparison result is used in an existence theorem and in the investigation of the stability of the numerical method of lines for the original problem. A theorem on the error estimate of the method is given. The infinite system of first order functional differential equations contains, as particular cases, equations with a deviated argument and integral differential equations of the Volterra type.     

Solution for a problem of linear plane elasticity with mixed boundary conditions on an ellipse by the method of boundary integrals

A numerical boundary integral scheme is proposed for the solution of the system of field equations of plane, linear elasticity in stresses for homogeneous, isotropic media in the domain bounded by an ellipse under mixed boundary conditions. The stresses are prescribed on one half of the ellipse, while the displacements are given on the other half. The method relies on previous analytical work within the Boundary Integral Method [1], [2].The considered problem with mixed boundary conditions is replaced by two subproblems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way and the problem at this stage is reduced to the solution of a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution inside the domain, and the unknown boundary values of stresses or displacements on proper parts of the boundary.On the basis of the obtained results, it is inferred that the tangential stress component on the fixed part of the boundary has a singularity at each of the two separation points, thought to be of logarithmic type. A tentative form for the singular solution is proposed to calculate the full solution in bulk directly from the given boundary conditions using the well-known Boundary Collocation Method. It is shown that this addition substantially decreases the error in satisfying the boundary conditions on some interval not containing the singular points.The obtained results are discussed and boundary curves for unknown functions are provided, as well as three-dimensional plots for quantities of practical interest. The efficiency of the used numerical schemes is discussed, in what concerns the number of boundary nodes needed to calculate the approximate solution.     

Approximate solution to a boundary value problem for the Lyapunov differential equation with a parameter

A special boundary value problem is studied for the Lyapunov differential equation which is used for investigation of the asymptotic properties of solutions to systems of periodic differential equations with a parameter. An algorithm is proposed for constructing an approximate solution to this boundary value problem, and conditions on the parameter are found under which the zero solution to the system is asymptotically stable.     

Finite difference reaction-diffusion systems with coupled boundary conditions and time delays

This paper is concerned with finite difference solutions of a coupled system of reaction-diffusion equations with nonlinear boundary conditions and time delays. The system is coupled through the reaction functions as well as the boundary conditions, and the time delays may appear in both the reaction functions and the boundary functions. The reaction-diffusion system is discretized by the finite difference method, and the investigation is devoted to the finite difference equations for both the time-dependent problem and its corresponding steady-state problem. This investigation includes the existence and uniqueness of a finite difference solution for nonquasimonotone functions, monotone convergence of the time-dependent solution to a maximal or a minimal steady-state solution for quasimonotone functions, and local and global attractors of the time-dependent system, including the convergence of the time-dependent solution to a unique steady-state solution. Also discussed are some computational algorithms for numerical solutions of the steady-state problem when the reaction function and the boundary function are quasimonotone. All the results for the coupled reaction-diffusion equations are directly applicable to systems of parabolic-ordinary equations and to reaction-diffusion systems without time delays.     

On numerical study of periodic solutions of a delay equation in biological models

Some results are presented of the numerical study of periodic solutions of a nonlinear equation with a delayed argument in connection with themathematical models having real biological prototypes. The problem is formulated as a boundary value problem for a delay equation with the conditions of periodicity and transversality. A spline-collocation finite-difference scheme of the boundary value problem using a Hermitian interpolation cubic spline of the class C ¹ with fourth order error is proposed. For the numerical study of the system of nonlinear equations of the finitedifference scheme, the parameter continuation method is used, which allows us to identify possible nonuniqueness of the solution of the boundary value problem and, hence, the nonuniqueness of periodic solutions regardless of their stability. By examples it is shown that the periodic oscillations occur for the parameter values specific to the real molecular-genetic systems of higher species, for which the principle of delay is quite easy to implement.     

Numerical Solution of the Incompressible Navier-Stokes Equations with Singular Perturbation Considerations

A numerical method for coarse grids is proposed for the numerical solution of the incompressible Navier-Stokes equations. From singular perturbation considerations, we obtain partial differential equations and boundary conditions for the outer solution and the boundary layer correction. The former problem is solved with the finite difference method and the latter with the approximate method. Numerical experiments show that accurate outer flow and boundary flux result with little computational effort.     

Determining the boundary of a two-dimensional region from the solution of the external initial boundary-value problem for the heat equation

We determine the boundary of a two-dimensional region using the solution of the external initial boundary-value problem for the nonhomogeneous heat equation. The initial values for the boundary determination include the right-hand side of the equation and the solution of the initial boundary-value problem given for finitely many points outside the region. The inverse problem is reduced to solving a system of two integral equations nonlinear in the function defining the sought boundary. An iterative procedure is proposed for numerical solution of the problem involving linearization of integral equations. The efficiency of the proposed procedure is investigated by a computer experiment.     

Generalization of a Relaxation Scheme for Systems of Forced Nonlinear Hyperbolic Conservation Laws with Spatially Dependent Flux Functions

A generalization of a finite difference method for calculating numerical solutions to systems of nonlinear hyperbolic conservation laws in one spatial variable is investigated. A previously developed numerical technique called the relaxation method is modified from its initial application to solve initial value problems for systems of nonlinear hyperbolic conservation laws. The relaxation method is generalized in three ways herein to include problems involving any combination of the following factors: systems of nonlinear hyperbolic conservation laws with spatially dependent flux functions, nonzero forcing terms, and correctly posed boundary values. An initial value problem for the forced inviscid Burgers' equation is used as an example to show excellent agreement between theoretical solutions and numerical calculations. An initial boundary value problem consisting of a system of four partial differential equations based on the two-layer shallow-water equations is solved numerically to display a more general applicability of the method than was previously known.     

On a free interface problem for linear ordinary differential equations and the one-phase Stefan problem

A numerical method for the solution of the one-phase Stefan problem is discussed. By discretizing the time variable the Stefan problem is reduced to a sequence of free boundary value problems for ordinary differential equations which are solved by conversion to initial value problems. The numerical solution is shown to converge to the solution of the Stefan problem with decreasing time increments. Sample calculations indicate that the method is stable provided the proper algorithm is chosen for integrating the initial value problems.     

Numerical analysis of a mathematical model for propagation of an electrical pulse in a neuron

This article is devoted to the study of a mathematical model arising in the mathematical modeling of pulse propagation in nerve fibers. A widely accepted model of nerve conduction is based on nonlinear parabolic partial differential equations. When considered as part of a particular initial boundary value problem the equation models the electrical activity in a neuron. A small perturbation parameter ε is introduced to the highest order derivative term. The parameter if decreased, speeds up the fast variables of the model equations whereas it does not affect the slow variables. In order to formally reduce the problem to a discussion of the moment of fronts and backs we take the limit ε → 0. This limit is singular and is therefore the solution tends to a slowly moving solution of the limiting equation. This leads to the boundary layers located in the neighborhoods of the boundary of the domain where the solution has very steep gradient. Most of the classical methods are incapable of providing helpful information about this limiting solution. To this effort a parameter robust numerical method is constructed on a piecewise uniform fitted mesh. The method consists of standard upwind finite difference operator. A rigorous analysis is carried out to obtain priori estimates on the solution of the problem and its derivatives. A parameter uniform error estimate for the numerical scheme so constructed is established in the maximum norm. It is then proven that the numerical method is unconditionally stable and provides a solution that converges to the solution of the differential equation. A set of numerical experiment is carried out in support of the predicted theory, which validates computationally the theoretical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Finite difference reaction diffusion equations with nonlinear boundary conditions

This article is concerned with numerical solutions of finite difference systems of reaction diffusion equations with nonlinear internal and boundary reaction functions. The nonlinear reaction functions are of general form and the finite difference systems are for both time-dependent and steady-state problems. For each problem a unified system of nonlinear equations is treated by the method of upper and lower solutions and its associated monotone iterations. This method leads to a monotone iterative scheme for the computation of numerical solutions as well as an existence-comparison theorem for the corresponding finite difference system. Special attention is given to the dynamical property of the time-dependent solution in relation to the steady-state solutions. Application is given to a heat-conduction problem where a nonlinear radiation boundary condition obeying the Boltzmann law of cooling is considered. This application demonstrates a bifurcation property of two steady-state solutions, and determines the dynamic behavior of the time-dependent solution. Numerical results for the heat-conduction problem, including a test problem with known analytical solution, are presented to illustrate the various theoretical conclusions. © 1995 John Wiley & Sons, Inc.