The boundary integral equation approach for numerical solution of the one‐dimensional Sine‐Gordon equation

This article describes a numerical method based on the boundary integral equation and dual reciprocity method for solving the one‐dimensional Sine‐Gordon (SG) equation. The time derivative is approximated by the time‐stepping method and a predictor–corrector scheme is employed to deal with the nonlinearity which appears in the problem. Numerical results are presented for some problems to demonstrate the usefulness and accuracy of this approach. In addition, the conservation of energy in SG equation is investigated. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

A numerical solution for advection‐diffusion equation using dual reciprocity method

This article describes a numerical method based on the boundary integral equation and dual reciprocity method(DRM) for solving the one‐dimensional advection‐diffusion equations. The concept of DRM is used to convert the domain integral to the boundary that leads to an integration free method. The time derivative is approximated by the time‐stepping method. Numerical results are presented for some problems to demonstrate the usefulness and accuracy of the new approach. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Numerical solution of a linear elliptic partial differential equation with variable coefficients: A complex variable boundary element approach

This article presents a complex variable boundary element method for the numerical solution of a second order elliptic partial differential equation with variable coefficients. To assess the validity and accuracy of the method, it is applied to solve some specific problems with known solutions. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Determination of the source parameter in a heat equation with a non-local boundary condition

In this article we consider the inverse problem of identifying a time dependent unknown coefficient in a parabolic problem subject to initial and non-local boundary conditions along with an overspecified condition defined at a specific point in the spatial domain. Due to the non-local boundary condition, the system of linear equations resulting from the backward Euler approximation have a coefficient matrix that is a quasi-tridiagonal matrix. We consider an efficient method for solving the linear system and the predictor–corrector method for calculating the solution and updating the estimate of the unknown coefficient. Two model problems are solved to demonstrate the performance of the methods.     

A method based on meshless approach for the numerical solution of the two‐space dimensional hyperbolic telegraph equation

In the current article, we investigate the RBF solution of second‐order two‐space dimensional linear hyperbolic telegraph equation. For this purpose, we use a combination of boundary knot method (BKM) and analog equation method (AEM). The BKM is a meshfree, boundary‐only and integration‐free technique. The BKM is an alternative to the method of fundamental solution to avoid the fictitious boundary and to deal with low accuracy, singular integration and mesh generation. Also, on the basis of the AEM, the governing operator is substituted by an equivalent nonhomogeneous linear one with known fundamental solution under the same boundary conditions. Finally, several numerical results and discussions are demonstrated to show the accuracy and efficiency of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

The use of cubic B‐spline scaling functions for solving the one‐dimensional hyperbolic equation with a nonlocal conservation condition

Problems for parabolic partial differential equations with nonlocal boundary conditions have been studied in many articles, but boundary value problems for hyperbolic partial differential equations have so far remained nearly uninvestigated. In this article a numerical technique is presented for the solution of a nonclassical problem for the one‐dimensional wave equation. This method uses the cubic B‐spline scaling functions. Some numerical results are reported to support our study. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

利用边界元法求解一类重调和方程

边界元法(BEM)和多重互易法(MRM)相结合求解一类重调和方程.通过重调和基本解序列给出的MRM-方法和BEM, 推导出该类问题的MRM-边界变分方程, 用边界元法求解该变分方程, 从而得到重调和方程的近似解, 并给出了解的存在唯一性证明.通过数值算例说明了MRM-方法具有收敛速度快、计算精度高, 易编程等优点, 为使用边界元法数值求解重调和方程提供了方法和理论依据.适合于工程中的实际运算.     

一类二阶拟线性方程边值问题解的存在性

文献[1]讨论了二阶拟线性常微分方程两点边值问题的 正解的存在性,但它限制了0<α<1<β.本文则对方程 证明了正解的存在性.     

A boundary only procedure for time-dependent diffusion equations

In the recent literature, the boundary element method (BEM) is extensively used to solve time-dependent partial differential equations. However, most of these formulations yield algorithms where one has to include all interior points in the computation process if finite difference procedures are used to approximate the temporal derivative. This obviously restricts the advantages of the BEM, which is mainly considered to be a boundary only algorithm for time-independent problems. A new algorithm is demonstrated here, which extends the boundary only nature of the method to time-dependent partial differential equations. Using this procedure, one can reduce the finite difference time integration algorithm, generated in a standard manner, to a boundary only process. The proposed method is demonstrated with considerable success for diffusion problems. Results obtained in these applications are presented comparatively with analytical and other boundary element time integration procedures. The algorithm proposed may utilize several coordinate functions in the secondary reduction phase of the formulation. A summary of such functions is described here and performances of these functions are tested and compared in three applications. It is shown that some coordinate functions perform better than others under certain conditions. Using these results, we propose a general coordinate function, which may be used with satisfactory results in all parabolic partial differential equation applications.     

A FINITE DIFFERENCE METHOD FOR THE HEAT EQUATION WITH A NONLINEAR BOUNDARY CONDITION

This article is concerned with the numerical solution to a parabolic equation with a kind of nonlinear boundary conditions. A difference scheme is constructed by the method of reduction of order on uniform mesh to solve the problem. It is proved that the difference scheme is uniquely solvable and uncon-ditionaUy convergent with the convergence order 2 in both space and time in an energy norm. An effective iterative algorithm is given and a numerical example is presented to demonstrate the theoretical results.     

A Robin boundary value problem in the upper half plane for the Bitsadze equation

In this article we give the solvability conditions and an integral representation of the solution of a Robin problem for the Bitsadze equation in the upper half plane. In order to do that, we use classical results of complex analysis and carry out the composition of two Robin problems for the Cauchy Riemann operator.     

热方程非线性边界条件问题的有限差分方法

This article is concerned with the numerical solution to a parabolic equation with a kind of nonlinear boundary conditions. A difference scheme is constructed by the method of reduction of order on uniform mesh to solve the problem. It is proved that the difference scheme is uniquely solvable and uncon-ditionaUy convergent with the convergence order 2 in both space and time in an energy norm. An effective iterative algorithm is given and a numerical example is presented to demonstrate the theoretical results.     

Bernstein Ritz‐Galerkin method for solving an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

In this article, the Ritz‐Galerkin method in Bernstein polynomial basis is implemented to give an approximate solution of a hyperbolic partial differential equation with an integral condition. We will deal here with a type of nonlocal boundary value problem, that is, the solution of a hyperbolic partial differential equation with a nonlocal boundary specification. The nonlocal conditions arise mainly when the data on the boundary cannot be measured directly. The properties of Bernstein polynomial and Ritz‐Galerkin method are first presented, then Ritz‐Galerkin method is used to reduce the given hyperbolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique presented in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A stable boundary elements method for magnetohydrodynamic channel flows at high Hartmann numbers

The article is devoted to extension of boundary element method (BEM) for solving coupled equations in velocity and induced magnetic field for time dependent magnetohydrodynamic (MHD) flows through a rectangular pipe. The BEM is equipped with finite difference approach to solve MHD problem at high Hartmann numbers up to 10⁶. In fact, the finite difference approach is used to approximate partial derivatives of unknown functions at boundary points respect to outward normal vector. It yields a numerical method with no singular boundary integrals. Besides, a new approach is suggested in this article where transforms 2D singular BEM's integrals to 1D nonsingular ones. The new approach reduces computational cost, significantly. Note that the stability of the numerical scheme is proved mathematically when computational domain is discretized uniformly and Hartmann number is 40 times bigger than length of boundary elements. Numerical examples show behavior of velocity and induced magnetic field across the sections.     

Boundary element method and wave equation

In this paper the boundary integral expression for a one-dimensional wave equation with homogeneous boundary conditions is developed. This is done using the time dependent fundamental solution of the corresponding hyperbolic partial differential equation. The integral expression developed is a generalized function with the same form as the well-known D'Alembert formula. The derivatives of the solution and some useful invariants on the characteristics of the partial differential equation are also calculated.The boundary element method is applied to find the numerical solution. The results show excellent agreement with analytical solutions.A multi-step procedure for large time steps which can be used in the boundary element method is also described.In addition, the way in which boundary conditions are introduced during the time dependent process is explained in detail. In the Appendix the main properties of Dirac's delta function and the Heaviside unit step function are described.     

THE ARTIFICIAL BOUNDARY CONDITION FOR EXTERIOR OSEEN EQUATION IN 2-D SPACE

A finite element method for the solution of Oseen equation in exterior domain is proposed. In this method, a circular artificial boundary is introduced to make the computational domain finite. Then, the exact relation between the normal stress and the prescribed velocity field on the artificial boundary can be obtained analytically. This relation can serve as an boundary condition for the boundary value problem defined on the finite domain bounded by the artificial boundary. Numerical experiment is presented to demonstrate the performance of the method.     

Geometric properties of the solutions of a Hele-Shaw type equation

This article deals with the application of the methods of geometric function theory to the investigation of the free boundary problem for the equation describing flows in an unbounded simply-connected plane domain. We prove the invariance of some geometric properties of a moving boundary.

     

Exact and approximation product solutions form of heat equation with nonlocal boundary conditions using Ritz–Galerkin method with Bernoulli polynomials basis

In this article, a new method is introduced for finding the exact solution of the product form of parabolic equation with nonlocal boundary conditions. Approximation solution of the present problem is implemented by the Ritz–Galerkin method in Bernoulli polynomials basis. The properties of Bernoulli polynomials are first presented, then Ritz–Galerkin method in Bernoulli polynomials is used to reduce the given differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the techniques presented in this article for finding the exact and approximation solutions. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1143–1158, 2017     

High‐order energy‐preserving schemes for the improved Boussinesq equation

This article proposes a class of high‐order energy‐preserving schemes for the improved Boussinesq equation. To derive the energy‐preserving schemes, we first discretize the improved Boussinesq equation by Fourier pseudospectral method, which leads to a finite‐dimensional Hamiltonian system. Then, the obtained semidiscrete system is solved by Hamiltonian boundary value methods, which is a newly developed class of energy‐preserving methods. The proposed schemes can reach spectral precision in space, and in time can reach second‐order, fourth‐order, and sixth‐order accuracy, respectively. Moreover, the proposed schemes can conserve the discrete mass and energy to within machine precision. Furthermore, to show the efficiency and accuracy of the proposed methods, the proposed methods are compared with the finite difference methods and the finite volume element method. The results of several numerical experiments are given for the propagation of the single solitary wave, the interaction of two solitary waves and the wave break‐up.     

THE SINGLE-POINT QUENCHING OF NONLINEAR DEGENERATE FUNCTIONAL REACTION-DIFFUSION EQUATION

This article is concerned with the quenching phenomena of the nonlinear degenerate functional reaction-diffusion equation. Some results are obtained on the single-point quenching and the uniqueness of quenching.