Numerical solutions of two-dimensional nonlinear integral analysis

In this paper, the approximate solutions for two different type of two-dimensional nonlinear integral equations: two-dimensional nonlinear Volterra-Fredholm integral equations and the nonlinear mixed Volterra-Fredholm integral equations are obtained using the Laguerre wavelet method. To do this, these two-dimensional nonlinear integral equations are transformed into a system of nonlinear algebraic equations in matrix form. By solving these systems, unknown coefficients are obtained. Also, some theorems are proved for convergence analysis.Some numerical examples are presented and results are compared with the analytical solution to demonstrate the validity and applicability of the proposed method.     

A boundary element method to price time-dependent double barrier options

In this paper we propose a new method for pricing double-barrier options with moving barriers under the Black-Scholes and the CEV models. First of all, by applying a variational technique typical of the boundary element method, we derive an integral representation of the double-barrier option price in which two of the integrand functions are not given explicitly but must be obtained solving a system of Volterra integral equations of the first kind. Second, we develop an ad hoc numerical method to regularize and solve the system of integral equations obtained. Several numerical experiments are carried out showing that the overall algorithm is extraordinarily fast and accurate, even if the barriers are not differentiable functions. Moreover the numerical method presented in this paper performs significantly better than the finite difference approach.     

Numerical solution of nonlinear two-dimensional Fredholm integral equations of the second kind using Gauss product quadrature rules

The Gauss product quadrature rules and collocation method are applied to reduce the second-kind nonlinear two-dimensional Fredholm integral equations (FIE) to a nonlinear system of equations. The convergence of the proposed numerical method is proved under certain conditions on the kernel of the integral equation. An iterative method for approximating the solution of the obtained nonlinear system is provided and its convergence is proved. Also, some numerical examples are presented to show the efficiency and accuracy of the proposed method.     

An efficient cubic B-spline and bicubic B-spline collocation method for numerical solutions of multidimensional nonlinear stochastic quadratic integral equations

This paper aims to develop a novel numerical approach on the basis of B-spline collocation method to approximate the solution of one-dimensional and two-dimensional nonlinear stochastic quadratic integral equations. The proposed approach is based on the hybrid of collocation method, cubic B-spline, and bi-cubic B-spline interpolation and Itô approximation. Using this method, the problem solving turns into a nonlinear system solution of equations that is solved by a suitable numerical method. Also, the convergence analysis of this numerical approach has been discussed. In the end, examples are given to test the accuracy and the implementation of the method. The results are compared with the results obtained by other methods to verify that this method is accurate and efficient.     

Solution of nonlinear weakly singular Volterra integral equations using the fractional‐order Legendre functions and pseudospectral method

In this article, our main goal is to render an idea to convert a nonlinear weakly singular Volterra integral equation to a non‐singular one by new fractional‐order Legendre functions. The fractional‐order Legendre functions are generated by change of variable on well‐known shifted Legendre polynomials. We consider a general form of singular Volterra integral equation of the second kind. Then the fractional Legendre–Gauss–Lobatto quadratures formula eliminates the singularity of the kernel of the integral equation. Finally, the Legendre pseudospectral method reduces the solution of this problem to the solution of a system of algebraic equations. This method also can be utilized on fractional differential equations as well. The comparison of results of the presented method and other numerical solutions shows the efficiency and accuracy of this method. Also, the obtained maximum error between the results and exact solutions shows that using the present method leads to accurate results and fast convergence for solving nonlinear weakly singular Volterra integral equations. Copyright © 2015 John Wiley & Sons, Ltd.     

Modified integral current methods in electrodynamics of nonhomogeneous media

The main difficulty in numerical solution of integral equations of electrodynamics is associated with the need to solve a high-order system of linear equations with a dense matrix. It is therefore relevant to develop numerical methods that lead to linear equation systems of lower order at the cost of more complex evaluation of the coefficients. In this article we propose a method for solving linear equations of electrodynamics which is a modification of the integral current method. The main distinctive feature of the proposed method is double integration of the electric Green’s tensor in the process of algebraization of the original integral equation. The solutions of the system of linear equations are thus integral means of the electric field inside the anomaly constructed by the proposed transformation formula. We prove convergence and derive error bounds for both the solution of the integral equation and the electromagnetic field components evaluated from approximate transformation formulas.     

A note on the polynomial solution of a class of dual integral equations arising in mixed boundary value problems of elasticity

Summary The present paper is concerned with finding an effective polynomial solution to a class of dual integral equations which arise in many mixed boundary value problems in the theory of elasticity. The dual integral equations are first transformed into a Fredholm integration equation of the second kind via an auxiliary function, which is next reduced to an infinite system of linear algebraic equations by representing the unknown auxiliary function in the form of an infinite series of Jacobi polynomials. The approximate solution of this infinite system of equations can be obtained by a suitable truncation. It is shown that the unknown function involving the dual integral equations can also be expressed in the form of an infinite series of Jacobi polynomials with the same expansion coefficients with no numerical integration involved. The main advantage of the present approach is that the solution of the dual integral equations thus obtained is numerically more stable than that obtained by reducing themdirectly into an infinite system of equations, insofar as the expansion coefficients are determined essentially by solving asecond kind integral equation.     

A boundary integral equation procedure for the mixed boundary value problem of the vector helmholtz equation

A boundary integral method is developed for the mixed boundary value problem for the vector Helmholtz equation in R³. The obtained boundary integral equations for the unknown Cauchy data build a strong elliptic system of pseudodifferential equations which can therefore be used for numerical computations using Galerkin's procedure. We show existence, uniqueness and regularity of the solution of the integral equations. Especially we give the local "edge" behavior of the solution near the submanifold which divides the Dirichlet boundary from the Neumann boundary     

A numerical method for solving linear integral equations of the second kind on the non-rectangular domains based on the meshless method

The main purpose of this article is to describe a numerical scheme for solving two-dimensional linear Fredholm integral equations of the second kind on a non-rectangular domain. The method approximates the solution by the discrete collocation method based on radial basis functions (RBFs) constructed on a set of disordered data. The proposed method does not require any background mesh or cell structures, so it is meshless and consequently independent of the geometry of domain. This approach reduces the solution of the two-dimensional integral equation to the solution of a linear system of algebraic equations. The error analysis of the method is provided. The proposed scheme is also extended to linear mixed Volterra–Fredholm integral equations. Finally, some numerical examples are presented to illustrate the efficiency and accuracy of the new technique.     

Convergence of approximate solution of system of Fredholm integral equations

In this paper numerical solution of system of linear Fredholm integral equations by means of the Sinc-collocation method is considered. This approximation reduces the system of integral equations to an explicit system of algebraic equations. The exponential convergence rate of the method is proved. The method is applied to a few test examples with continuous kernels to illustrate the accuracy and the implementation of the method.     

Electrostatic problem of several arbitrarily charged unequal coaxial disks

The title problem is treated by a new method which allows a straightforward derivation of a Fredholm system of governing integral equations. An important feature of these equations is the fact that there is no need to solve the system in order to assess certain integral characteristics like, for example, the total charge at each disk. A good accuracy of this assessment is confirmed by a numerical solution of the system of integral equations. Some previously published numerical results are corrected. Several examples are considered.     

An extrapolation method for a class of boundary integral equations

Boundary value problems of the third kind are converted into boundary integral equations of the second kind with periodic logarithmic kernels by using Green's formulas. For solving the induced boundary integral equations, a Nyström scheme and its extrapolation method are derived for periodic Fredholm integral equations of the second kind with logarithmic singularity. Asymptotic expansions for the approximate solutions obtained by the Nyström scheme are developed to analyze the extrapolation method. Some computational aspects of the methods are considered, and two numerical examples are given to illustrate the acceleration of convergence.

     

On the numerical solution of loaded systems of ordinary differential equations with nonseparated multipoint and integral conditions

We propose a numerical method of solving systems of loaded linear nonautonomous ordinary differential equations with nonseparated multipoint and integral conditions. This method is based on the convolution of integral conditions to obtain local conditions. This approach allows one to reduce solving the original problem to solving a Cauchy problem for a system of ordinary differential equations and linear algebraic equations. Numerous computational experiments on several test problems with the formulas and schemes proposed for the numerical solution have been carried out. The results of the experiments show that the approach is reasonably efficient.     

关于多裂纹圆柱体的扭转*

本文在文[1]基础上,导出了含有任意分布裂纹系的圆柱扭曲函数的解析表达式,从而把问题化为以未知位错密度函数表示的奇异积分方程组.文中利用奇异积分方程的数值方法[2,7],对带有多根裂纹的圆柱的抗扭刚度和应力强度因子作了若干数值计算.此外,本文还首次将裂纹切割法[5]推广用于求解矩形柱的扭转,数值结果表明方法是成功的.     

The numerical solution of an evolution problem of second order in time on a closed smooth boundary

We consider an initial value problem for the second-order differential equation with a Dirichlet-to-Neumann operator coefficient. For the numerical solution we carry out semi-discretization by the Laguerre transformation with respect to the time variable. Then an infinite system of the stationary operator equations is obtained. By potential theory, the operator equations are reduced to boundary integral equations of the second kind with logarithmic or hypersingular kernels. The full discretization is realized by Nyström's method which is based on the trigonometric quadrature rules. Numerical tests confirm the ability of the method to solve these types of nonstationary problems.     

On the numerical solution of Cauchy type singular integral equations by the collocation method

The collocation method for the numerical solution of Fredholm integral equations of the second kind is applied, properly modified, to the numerical solution of Cauchy type singular integral equations of the first or the second kind but with constant coefficients. This direct method of numerical solution of Cauchy type singular integral equations is compared afterwards with the corresponding method resulting from applying the collocation method to the Fredholm integral equation of the second kind equivalent to the Cauchy type singular integral equation, as well as with another method, based also on the regularization procedure, for the numerical solution of the same class of equations. Finally, the convergence of the method is discussed.     

Joint availability of systems modelled by finite semi–markov processes

Consider a repairable system at the time instants t and t + x, where t, x ≥0. The joint availability of the system at these time instants is defined as the probability of the system being functional in both t and t + x. A set of integral equations is derived for the joint availability of a system modelled by a finite semi–Markov process. The result is applied to the semi–Markov model of a two–unit system with sequential preventive maintenance. The method used for the numerical solution of the resulting system of integral equations is a two–point trapezoidal rule. The system of implementation is the matrix computation package MATLAB on the Apple Macintosh SE/30. The numerical results obtained by this method are shown to be in good agreement with those from simulation.     

Boundary element methods—An overview

Variational methods for boundary integral equations deal with the weak formulations of boundary integral equations. Their numerical discretizations are known as the boundary element methods. This paper gives an overview of the method from both theoretical and numerical point of view. It summarizes the main results obtained by the author and his collaborators over the last 30 years. Fundamental theory and various applications will be illustrated through simple examples. Some numerical experiments in elasticity as well as in fluid mechanics will be included to demonstrate the efficiency of the methods.     

Newton methods for a class of nonlinear hypersingular integral equations

Different iterative schemes based on collocation methods have been well studied and widely applied to the numerical solution of nonlinear hypersingular integral equations (Capobianco et al. 2005). In this paper we apply Newton’s method and its modified version to solve the equations obtained by applying a collocation method to a nonlinear hypersingular integral equation of Prandtl’s type. The corresponding convergence results are derived in suitable Sobolev spaces. Some numerical tests are also presented to validate the theoretical results.     

Numerical Solution Based on Hat Functions for Solving Nonlinear Stochastic Itô Volterra Integral Equations Driven by Fractional Brownian Motion

This paper presents a numerical method for solving nonlinear stochastic Itô Volterra integral equations driven by fractional Brownian motion with Hurst parameter \( H \in (0,1)\) via of hat functions. Using properties of the generalized hat basis functions and fractional Brownian motion, new stochastic operational matrix of integration is achieved and the nonlinear stochastic equation is transformed into nonlinear system of algebraic equations which by solving it, an approximation solution with high accuracy is obtained. In addition, error analysis of the method is investigated, and by some examples, efficiency and accuracy of the suggested method are shown.