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.     

Determination of periodic orbits of nonlinear oscillators by means of piecewise-linearization methods

An approximate method based on piecewise linearization is developed for the determination of periodic orbits of nonlinear oscillators. The method is based on Taylor series expansions, provides piecewise analytical solutions in three-point intervals which are continuous everywhere and explicit three-point difference equations which are P-stable and have an infinite interval of periodicity. It is shown that the method presented here reduces to the well-known Störmer technique, is second-order accurate, and yields, upon applying Taylor series expansion and a Padé approximation, another P-stable technique whenever the Jacobian is different from zero. The method is generalized for single degree-of-freedom problems that contain the velocity, and (approximate) analytical solutions are presented. Finally, by introducing the inverse of a vector and the vector product and quotient, and using Taylor series expansions and a Padé approximation, the method has been generalized to multiple degree-of-freedom problems and results in explicit three-point finite difference equations which only involve vector multiplications.     

Improving the accuracy of solutions of the linear second kind volterra integral equations system by using the taylor expansion method

This paper presents a new and an efficient method for determining solutions of the linear second kind Volterra integral equations system. In this method, the linear Volterra integral equations system using the Taylor series expansion of the unknown functions transformed to a linear system of ordinary differential equations. For determining boundary conditions we use a new method. This method is effective to approximate solutions of integral equations system with a smooth kernel, and a convolution kernel. An error analysis for the proposed method is provided. And illustrative examples are given to represent the efficiency and the accuracy of the proposed method.     

Numerical solution of integral equations of the first kind

An approximate method for solving integral equations of the first kind is considered, with the approximate solution represented as a finite expansion in some basis. Solution examples for a number of model problems are given. The dependence of the approximation error on the accuracy of the initial data is analyzed numerically.Translated from Metody Matematicheskogo Modelirovaniya, Avtomatizatsiya Obrabotki Nablyudenii i Ikh Primeneniya, pp. 72–79, 1986.     

非线性积分方程迭代配置法的渐近展开及其外推

非线性积分方程迭代配置法的渐近展开及其外推韩国强（华南理工大学计算机工程与科学系）ＡＳＹＭＰＴＯＴＩＣＥＲＲＯＲＥＸＭＮＳＩＯＮＳＡＮＤＥＸＴＲＡＰＯＬＡＴＩＯＮＦＯＲＴＨＥＩＴＥＲＡＴＥＤＣＯＬＬＯＣＡＴＩＯＮＭＥＴＨＯＤＳＯＦＮＯＮＬＩＮＥＡＲＩ...     

On the separation of isolated solutions of nonlinear integral equations

We present the theoretical justification and a method for practical realization of the process of separation of solutions isolated in a bounded domain for some classes of nonlinear integral equations. We study the problem of construction of a sequence of approximation equations by the method of mechanical quadratures and the problem of existence of solutions of these equations. We also present methods for approximate solution of these equations and obtaina posteriori error estimates.     

Application of wavelet transforms to the solution of boundary value problems for linear parabolic equations

A method based on wavelet transforms is proposed for finding weak solutions to initial-boundary value problems for linear parabolic equations with discontinuous coefficients and inexact data. In the framework of multiresolution analysis, the general scheme for finite-dimensional approximation in the regularization method is combined with the discrepancy principle. An error estimate is obtained for the stable approximate solution obtained by solving a set of linear algebraic equations for the wavelet coefficients of the desired solution.     

MULTILEVEL AUGMENTATION METHODS FOR SOLVING OPERATOR EQUATIONS

We introduce multilevel augmentation methods for solving operator equations based on direct sum decompositions of the range space of the operator and the solution space of the operator equation and a matrix splitting scheme. We establish a general setting for the analysis of these methods, showing that the methods yield approximate solutions of the same convergence order as the best approximation from the subspace. These augmentation methods allow us to develop fast, accurate and stable nonconventional numerical algorithms for solving operator equations. In particular, for second kind equations, special splitting techniques are proposed to develop such algorithms. These algorithms are then applied to solve the linear systems resulting from matrix compression schemes using wavelet-like functions for solving Fredholm integral equations of the second kind. For this special case, a complete analysis for computational complexity and convergence order is presented. Numerical examples are included to demonstra     

Branching of solutions for the problem of mean-square approximation of a real finite function of two variables by the modulus of double Fourier transformation

We investigate the branching of solutions of a Hammerstein-type nonlinear two-dimensional integral equation that arises in the problems of mean-square approximation of a real finite nonnegative function of two variables by the modulus of double Fourier integral, depending on two parameters [Mat. Metody Fiz.-Mekh. Polya, 51, No. 1, 53–64; No. 4, 80–85 (2008)]. We have derived analytical expressions for eigenfunctions of the corresponding linear homogeneous integral equation, necessary for the construction of branched-off solutions, and obtained systems of transcendental equations for finding the points of their branching. We also present, in the first approximation, the analytical representations of complex solutions branched-off from the real solution for the two-dimensional case of branching.     

The use of Chebyshev series for approximate analytic solution of ordinary differential equations

Application of Chebyshev series to solve ordinary differential equations is described. This approach is based on the approximation of the solution to a given Cauchy problem and its derivatives by partial sums of shifted Chebyshev series. The coefficients of the series are determined by an iterative process using Markov quadrature formulas. It is shown that the proposed approach can be applied to formulate an approximate analytical method for solving Cauchy problems. A number of examples are considered to illustrate the obtaining of approximate analytical solutions in the form of partial sums of shifted Chebyshev series.     

Optimal homotopy analysis and control of error for solutions to the non-local Whitham equation

The Whitham equation is a non-local model for nonlinear dispersive water waves. Since this equation is both nonlinear and non-local, exact or analytical solutions are rare except for in a few special cases. As such, an analytical method which results in minimal error is highly desirable for general forms of the Whitham equation. We obtain approximate analytical solutions to the non-local Whitham equation for general initial data by way of the optimal homotopy analysis method, through the use of a partial differential auxiliary linear operator. A method to control the residual error of these approximate solutions, through the use of the embedded convergence control parameter, is discussed in the context of optimal homotopy analysis. We obtain residual error minimizing solutions, using relatively few terms in the solution series, in the case of several different kernels and associated initial data. Interestingly, we find that for a specific class of initial data, there exists an exact solution given by the first term in the homotopy expansion. A specific example of initial data which satisfies the condition producing an exact solution is included. These results demonstrate the applicability of optimal homotopy analysis to equations which are simultaneously nonlinear and non-local.     

Volterra积分方程的蒙特卡罗数值求解方法

本文讨论了第一类、第二类以及具有奇异核的Volterra积分方程的数值解问题.利用重要抽样蒙特卡罗随机模拟方法获得积分方程解的近似计算结果.通过对文献中算例的实现表明文中所提方法扩展了Volterra型积分方程的数值求解方法,     

An approximation theory for conjugate surfaces and solutions of elliptic multiple integral problems: Application to numerical solutions of generalized Laplace's equation

An approximation theory is given for a class of elliptic quadratic forms which include the study of conjugate surfaces for elliptic multiple integral problems. These ideas follow from the quadratic form theory of Hestenes, applied to multiple integral problems by Dennemeyer, and extended with applications for approximation problems by Gregory.The application of this theory to a variety of approximation problem areas in this setting is given. These include conjugate surfaces and conjugate solutions in the calculus of variations, oscillation problems for elliptic partial differential equations, eigenvalue problems for compact operators, numerical approximation problems, and, finally, the intersection of these problem areas.In the final part of this paper the ideas are specifically applied to the construction and counting of negative vectors in order to obtain new numerical methods for solving Laplace-type equations and to obtain the “Euler-Lagrange equations” for symmetric-banded tridiagonal matrices. In this new result (which will allow the reexamination of both the theory and applications of symmetricbanded matrices) one can construct, in a meaningful way, negative vectors, oscillation vectors, eigenvectors, and extremal solutions of classical problems as well as efficient algorithms for the numerical solution of partial differential equations. Numerical examples (test runs) are given.     

An approximation method using approximate approximations

The aim of this article is to extend the method of approximate approximations to boundary value problems. This method was introduced by V. Maz'ya in 1991 and has been used until now for the approximation of smooth functions defined on the whole space and for the approximation of volume potentials. In the present article we develop an approximation procedure for the solution of the interior Dirichlet problem for the Laplace equation in two dimensions using approximate approximations. The procedure is based on potential theoretical considerations in connection with a boundary integral equations method and consists of three approximation steps as follows. In the first step, the unknown source density in the potential representation of the solution is replaced by approximate approximations. In the second, the decay behavior of the generating functions is used to gain a suitable approximation for the potential kernel, and in the third, Nyström's method leads to a linear algebraic system for the approximate source density. For every step a convergence analysis is established and corresponding error estimates are given.     

Modified method of splitting with respect to physical processes for solving radiation gas dynamics equations

An approach based on a modified splitting method is proposed for solving the radiation gas dynamics equations in the multigroup kinetic approximation. The idea of the approach is that the original system of equations is split using the thermal radiation transfer equation rather than the energy equation. As a result, analytical methods can be used to solve integrodifferential equations and problems can be computed in the multigroup kinetic approximation without iteration with respect to the collision integral or matrix inversion. Moreover, the approach can naturally be extended to multidimensional problems. A high-order accurate difference scheme is constructed using an approximate Godunov solver for the Riemann problem in two-temperature gas dynamics.     

Reduced differential transform method for nonlinear integral member of Kadomtsev–Petviashvili hierarchy differential equations

The main objective of this paper is to use the reduced differential to transform method (RDTM) for finding the analytical approximate solutions of two integral members of nonlinear Kadomtsev–Petviashvili (KP) hierarchy equations. Comparing the approximate solutions which obtained by RDTM with the exact solutions to show that the RDTM is quite accurate, reliable and can be applied for many other nonlinear partial differential equations. The RDTM produces a solution with few and easy computation. This method is a simple and efficient method for solving the nonlinear partial differential equations. The analysis shows that our analytical approximate solutions converge very rapidly to the exact solutions.     

The formally variable separation approach for the dust-acoustic solitary waves with dust charge variation

The formally variable separation approach is used for handling the dust-acoustic solitary waves in a dusty plasma, including consideration of dust charge variation. New analytical solutions of nonlinear waves are formally derived for the governing equation of the system. We have triumphantly derived the exact analytical expressions and some approximate expressions of the nonlinear dust-acoustic waves in a dusty plasma under some special cases. The work introduces entirely new solutions and emphasizes the power of the newly developed method that can be used in problems with identical nonlinearities.     

Contact problems for functionally graded materials of complicated structure

An approximate analytical method allowing one to efficiently solve, to a preassigned accuracy, contact problems for materials with properties arbitrarily varying in depth is developed. Its possibilities are illustrated with the example of torsion of an elastic half-space, having a coating inhomogeneous across its thickness, by a circular stamp. All the results obtained are rigorously substantiated. For the approximate solutions constructed, their error is analyzed. The asymptotic properties of the solutions are investigated. The cases of a nonmonotonic change in the elastic properties are considered. In particular, the analytical solutions are examined in the case where the variation gradient of the elastic properties changes its sign many times. The results derived allow one to solve the inverse problems of elasticity theory of inhomogeneous media (e.g., the problem on controlling the variation in the elastic properties of a covering across its thickness).     

Application of homotopy-perturbation method to Klein–Gordon and sine-Gordon equations

In this paper, the homotopy-perturbation method (HPM) is employed to obtain approximate analytical solutions of the Klein–Gordon and sine-Gordon equations. An efficient way of choosing the initial approximation is presented. Comparisons with the exact solutions, the solutions obtained by the Adomian decomposition method (ADM) and the variational iteration method (VIM) show the potential of HPM in solving nonlinear partial differential equations.     

A recurrence method for a special class of continuous time linear programming problems

This article studies a numerical solution method for a special class of continuous time linear programming problems denoted by (SP). We will present an efficient method for finding numerical solutions of (SP). The presented method is a discrete approximation algorithm, however, the main work of computing a numerical solution in our method is only to solve finite linear programming problems by using recurrence relations. By our constructive manner, we provide a computational procedure which would yield an error bound introduced by the numerical approximation. We also demonstrate that the searched approximate solutions weakly converge to an optimal solution. Some numerical examples are given to illustrate the provided procedure.