Monotone-iterative method for mixed boundary value problems for generalized difference equations with “maxima”

In this paper the authors investigate special type of difference equations which involve both delays and the maximum value of the unknown function over a past time interval. This type of equations is used to model a real process which present state depends significantly on its maximal value over a past time interval. An appropriate mixed boundary value problem for the given nonlinear difference equation is set up. An algorithm, namely, the monotone iterative technique is suggested to solve this problem approximately. An important feature of our algorithm is that each successive approximation of the unknown solution is equal to the unique solution of an appropriately constructed initial value problem for a linear difference equation with “maxima”, and a formula for its explicit form is given. Also, each approximation is a lower/upper solution of the given nonlinear boundary value problem. Several numerical examples are considered to illustrate the practical application of the suggested algorithm.     

Approximation of Solution Components for Ill-Posed Problems by the Tikhonov Method with Total Variation

An ill-posed problem in the form of a linear operator equation given on a pair of Banach spaces is considered. Its solution is representable as a sum of a smooth and a discontinuous component. A stable approximation of the solution is obtained using a modified Tikhonov method in which the stabilizer is constructed as a sum of the Lebesgue norm and total variation. Each of the functionals involved in the stabilizer depends only on one component and takes into account its properties. Theorems on the componentwise convergence of the regularization method are stated, and a general scheme for the finite-difference approximation of the regularized family of approximate solutions is substantiated in the n-dimensional case.     

A numerical solution of the elasticity problem of settled of the wronkler ground with variable coefficients

In this paper, we have given numerical solution of the elasticity problem of settled on the wronkler ground with variable coefficient. The approximation solution of boundary value problem which is pertinent to this has been converted to integral equations, and then by using the successive approximation methods, has been reached. In addition to this, the approximation solution of the problem was put into Padé series form. We applied these methods to an example which is the elasticity problem of unit length homogeny beam, which is a special form of boundary value problem. First we calculate the successive approximation of the given boundary value problem then transform it into Padé series form, which give an arbitrary order for solving differential equation numerically.     

Two types of interface defects

The solution of a plane problem in the theory of elasticity for a two-component body with an interface, a finite part of which is either weakly distorted or is a weakly curved crack is constructed using the perturbation method. In the first case, it is assumed that the discontinuities in the forces and displacements at the interface are known, and, in the second case, the non-equilibrium nature of the load in the crack is taken into account. General quadrature formulae are derived for the complex potentials, which enable any approximation to be obtained in terms of elementary functions in many important practical cases. An algorithm is indicated for calculating each approximation. Families of defects are studied, the form of which is determined by power functions. The effect of the amplitude of the distortion and the shape of the interface crack on the Cherepanov–Rice integral as well as the shape of the distorted part of the interface on the stress concentration is investigated in the first approximation. An analysis of the applicability of the oscillating solution for a distorted interface crack is carried out. The results of the calculations are shown in the form of graphical relations.     

Analytical and numerical treatment of a fractional integro-differential equation of volterra-type

An analytical treatment of a fractional generalisation of the free electron laser (FEL) equation is given, whereby a closed form solution is obtained in terms of Kummer functions Φ(a, c; z). Subsequently, Tau method approximation of these functions facilitates the numerical evaluation of the solution of the FEL equation.     

Negative norm error control for second-kind convolution Volterra equations

Summary. We consider a piecewise constant finite element approximation to the convolution Volterra equation problem of the second kind: find such that in a time interval . An a posteriori estimate of the error measured in the norm is developed and used to provide a time step selection criterion for an adaptive solution algorithm. Numerical examples are given for problems in which is of a form typical in viscoelasticity theory. Received March 5, 1998 / Revised version received November 30, 1998 / Published online December 6, 1999     

The anti-symmetric ortho-symmetric solution of a linear matrix equation and its optimal approximation

This paper discusses the anti-symmetric ortho-symmetric solution of a linear matrix equation and its optimal approximation. By the generalized singular value decomposition of the matrices, the necessary and sufficient conditions for the solvability of the matrix equation and the general form of the anti-symmetric ortho-symmetric solution are given. In addition, the existence and uniqueness of the optimal approximation are proved. Numerical methods of the optimal approximation to a given matrix and numerical experiments are described.     

On the convergence of a numerical method for a hypersingular integral equation on a closed surface

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises in the solution of the Neumann boundary value problem for the Laplace equation with a representation of a solution in the form of a double-layer potential. We consider the case in which the interior or exterior boundary value problem is solved in a domain; whose boundary is a smooth closed surface, and an integral equation is written out on that surface. For the integral operator in that equation, we suggest quadrature formulas like the method of vortical frames with a regularization, which provides its approximation on the entire surface for the use of a nonstructured partition. We construct a numerical scheme for the integral equation on the basis of suggested quadrature formulas, prove an estimate for the norm of the inverse matrix of the related system of linear equations and the uniform convergence of numerical solutions to the exact solution of the hypersingular integral equation on the grid.     

An exponentially convergent algorithm for nonlinear differential equations in Banach spaces

An exponentially convergent approximation to the solution of a nonlinear first order differential equation with an operator coefficient in Banach space is proposed. The algorithm is based on an equivalent Volterra integral equation including the operator exponential generated by the operator coefficient. The operator exponential is represented by a Dunford-Cauchy integral along a hyperbola enveloping the spectrum of the operator coefficient, and then the integrals involved are approximated using the Chebyshev interpolation and an appropriate Sinc quadrature. Numerical examples are given which confirm theoretical results.

     

Résolution mathématique d'équations décrivant l'évolution d'une surface d'un matériau contraint

We consider a cristal structure, constituted by an elastic substrate and a film with a small thickness. The lattice parameters between the film and the substrate are not the same; consequently, a strain appears in the structure. This strain generates morphologies(see [1,2]).The difficulty consists in finding the profile of the film-vapor surface at any time, which depends on the elastic displacement of the structure. To this end, a physical model, detailed in [2], consists in solving a coupled system of partial derivative equations. The unknowns are the elastic displacement of the structure and the profile of the evolution surface. The elastic displacement solves the linearized elasticity equations posed over the domain occupied by the structure. The boundary of this domain depends on the evolution surface. The second equation is the evolution equation, depending on the elastic displacement by a term of the surface energy. This model is greatly simplified in order to obtain a decoupled two-dimensional model: the map of the film-vapor surface solves a non-linear partial derivatives equation, which is independent of the displacement of the structure.In this Note, we give some results of the existence and uniqueness of a solution for this model under some assumptions about the first derivative of the map.     

Pade Approximants for Linear Integral Equations?

Padé approximation is often used to find a matrix elementof the solution of a Fredholm integral equation of the secondkind. An example is given in which this process is necessarilydivergent, even though the kernel of the integral equation iscompact and the elements lie in a Hilbert space.     

Testing the validity of the effective rate constant approximation for surface reaction with transport

When one incorporates transport effects into a surface-volume reaction, an integrodifferential equation for the bound state concentration occurs. Such a form is inconvenient for data analysis. An effective rate constant approximation for the solution is correct to O(Da²) as the Damköhler number Da → 0. A numerical simulation of the integrodifferential equation is performed which shows that the effective rate constant approximation is useful even outside this regime.     

Algorithms for solving matrix polynomial equations of special form

In this paper we consider a series of algorithms for calculating radicals of matrix polynomial equations. A particular aspect of this problem arise in author’s work, concerning parameter identification of linear dynamic stochastic system. Special attention is given to searching the solution of an equation in a neighbourhood of some initial approximation. The offered approaches and algorithms allow us to receive fast and quite exact solution. We give some recommendations for application of given algorithms.     

Convergence of a numerical scheme of the type of the vortex loop method on a closed surface with an approximation to the surface shape

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises when solving the Neumann boundary value problem for the Laplace equation with the use of the representation of the solution in the form of a double layer potential. We study the case in which an exterior or interior boundary value problem is solved in a domain whose boundary is a smooth closed surface and the integral equation is written out on that surface. For the numerical solution of the integral equation, the surface is approximated by spatial polygons whose vertices lie on the surface. We construct a numerical scheme for solving the integral equation on the basis of such an approximation to the surface with the use of quadrature formulas of the type of the method of discrete singularities with regularization. We prove that the numerical solutions converge to the exact solution of the hypersingular integral equation uniformly on the grid.     

The dynamics of the motion of a flat punch on the boundary of an elastic half-plane

The dynamic contact problem of the motion of a flat punch on the boundary of an elastic half-plane is considered. During motion, the punch deforms the elastic half-plane, penetrating it in such a manner that its base remains parallel to the boundary of the half-plane at each instant of time. In movable coordinates connected to the moving punch, the contact problem reduces to solving a two-dimensional integral equation, whose two-dimensional kernel depends on the difference between the arguments for each of the variables. An approximate solution of the integral equation of the problem is constructed in the form of a Neumann series, whose zeroth term is represented in the form of the superposition of the solutions of two-dimensional integral equations on the coordinate semiaxis minus the solution of the integral equation on the entire axis. This approach provides a way to construct the solution of the two-dimensional integral equation of the problem in four velocity ranges of motion of the punch, which cover the entire spectrum of its velocities, as well as to perform a detailed analysis of the special features of the contact stresses and vertical displacements of the free surface on the boundary of the contract area. An approximate method for solving the integral equation, which is based on a special approximation of the integrand of the kernel of the integral equation in the complex plane, is proposed for obtaining effective solutions of the problem that do not contain singular quadratures.     

A remark on finding the coefficient of the dissipative boundary condition via the enclosure method in the time domain

An inverse problem for the wave equation outside an obstacle with a dissipative boundary condition is considered. The observed data are given by a single solution of the wave equation generated by an initial data supported on an open ball. An explicit analytical formula for the computation of the coefficient at a point on the surface of the obstacle, which is nearest to the center of the support of the initial data, is given. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical solutions of the Laplace’s equation

This paper uses the sinc methods to construct a solution of the Laplace’s equation using two solutions of the heat equation. A numerical approximation is obtained with an exponential accuracy. We also present a reliable algorithm of Adomian decomposition method to construct a numerical solution of the Laplace’s equation in the form a rapidly convergence series and not at grid points. Numerical examples are given and comparisons are made to the sinc solution with the Adomian decomposition method. The comparison shows that the Adomian decomposition method is efficient and easy to use.     

Convergence of the Neumann Series in BEM for the Neumann Problem of the Stokes System

A weak solution of the Neumann problem for the Stokes system in Sobolev space is studied in a bounded Lipschitz domain with connected boundary. A solution is looked for in the form of a hydrodynamical single layer potential. It leads to an integral equation on the boundary of the domain. Necessary and sufficient conditions for the solvability of the problem are given. Moreover, it is shown that we can obtain a solution of this integral equation using the successive approximation method. Then the consequences for the direct boundary integral equation method are treated. A solution of the Neumann problem for the Stokes system is the sum of the hydrodynamical single layer potential corresponding to the boundary condition and the hydrodynamical double layer potential corresponding to the trace of the velocity part of the solution. Using boundary behavior of potentials we get an integral equation on the boundary of the domain where the trace of the velocity part of the solution is unknown. It is shown that we can obtain a solution of this integral equation using the successive approximation method.     

An efficient meshfree point collocation moving least squares method to solve the interface problems with nonhomogeneous jump conditions

We are going to study a simple and effective method for the numerical solution of the closed interface boundary value problem with both discontinuities in the solution and its derivatives. It uses a strong‐form meshfree method based on the moving least squares (MLS) approximation. In this method, for the solution of elliptic equation, the second‐order derivatives of the shape functions are needed in constructing the global stiffness matrix. It is well‐known that the calculation of full derivatives of the MLS approximation, especially in high dimensions, is quite costly. In the current work, we apply the diffuse derivatives using an efficient technique. In this technique, we calculate the higher‐order derivatives using the approximation of lower‐order derivatives, instead of calculating directly derivatives. This technique can improve the accuracy of meshfree point collocation method for interface problems with nonhomogeneous jump conditions and can efficiently estimate diffuse derivatives of second‐ and higher‐orders using only linear basis functions. To introduce the appropriate discontinuous shape functions in the vicinity of interface, we choose the visibility criterion method that modifies the support of weight function in MLS approximation and leads to an efficient computational procedure for the solution of closed interface problems. The proposed method is applied for elliptic and biharmonic interface problems. For the biharmonic equation, we use a mixed scheme, which replaces this equation by a coupled elliptic system. Also the application of the present method to elasticity equation with discontinuities in the coefficients across a closed interface has been provided. Representative numerical examples demonstrate the accuracy and robustness of the proposed methodology for the closed interface problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1031–1053, 2015     

An iteration method for the symmetric solutions and the optimal approximation solution of the matrix equation AXB=C

An iteration method is constructed to solve the linear matrix equation AXB=C over symmetric X. By this iteration method, the solvability of the equation AXB=C over symmetric X can be determined automatically, when the equation AXB=C is consistent over symmetric X, its solution can be obtained within finite iteration steps, and its least-norm symmetric solution can be obtained by choosing a special kind of initial iteration matrix, furthermore, its optimal approximation solution to a given matrix can be derived by finding the least-norm symmetric solution of a new matrix equation . Finally, numerical examples are given for finding the symmetric solution and the optimal approximation symmetric solution of the matrix equation AXB=C.