A meshless method based on boundary integral equations and radial basis functions for biharmonic-type problems

This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed.     

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.     

A new model for the analysis of reinforced concrete members with a coupled HdBNM/FEM

As a truly boundary-type meshless method, the hybrid boundary node method (HdBNM) does not require ‘boundary element mesh’, either for the purpose of interpolation of the solution variables or for the integration of ‘energy’. In this paper, the HdBNM is coupled with the finite element method (FEM) for predicting the mechanical behaviors of reinforced concrete. The steel bars are considered as body forces in the concrete. A bond model is presented to simulate the bond-slip between the concrete and steels using fictitious spring elements. The computational scale and cost for meshing can be further reduced. Numerical examples, in 2D and 3D cases, demonstrate the efficiency of the proposed approach.     

关于薄板的无网格局部边界积分方程方法中的友解

无网格局部边界积分方程方法是最近发展起来的一种新的数值方法,这种方法综合了伽辽金有限元、边界元和无单元伽辽金法的优点,是一种具有广阔应用前景的、真正的无网格方法．把无网格局部边界积分方程方法应用于求解薄板问题,给出了薄板无网格局部边界积分方程方法所需要的友解及其全部公式．     

A numerical method based on extended Raviart–Thomas (ER‐T) mixed finite element method for solving damped Boussinesq equation

In this paper, we present the approximate solution of damped Boussinesq equation using extended Raviart–Thomas mixed finite element method. In this method, the numerical solution of this equation is obtained using triangular meshes. Also, for discretization in time direction, we use an implicit finite difference scheme. In addition, error estimation and stability analysis of both methods are shown. Finally, some numerical examples are considered to confirm the theoretical results. Copyright © 2017 John Wiley & Sons, Ltd.     

A Burton-Miller boundary element-free method for Helmholtz problems

A Burton-Miller boundary element-free method is developed by using the Burton-Miller formulation for meshless and boundary-only analysis of Helmholtz problems. The method can produce a unique solution at all wavenumbers and is valid for Dirichlet, Neumann and mixed problems simultaneously. An efficient numerical integration procedure is presented to handle both strongly singular and hypersingular boundary integrals directly and uniformly. Numerical results reveal that this direct meshless method only involves boundary nodes and can deal with Helmholtz problems at extremely large wavenumbers.     

Wavelets collocation methods for the numerical solution of elliptic BV problems

Based on collocation with Haar and Legendre wavelets, two efficient and new numerical methods are being proposed for the numerical solution of elliptic partial differential equations having oscillatory and non-oscillatory behavior. The present methods are developed in two stages. In the initial stage, they are developed for Haar wavelets. In order to obtain higher accuracy, Haar wavelets are replaced by Legendre wavelets at the second stage. A comparative analysis of the performance of Haar wavelets collocation method and Legendre wavelets collocation method is carried out. In addition to this, comparative studies of performance of Legendre wavelets collocation method and quadratic spline collocation method, and meshless methods and Sinc–Galerkin method are also done. The analysis indicates that there is a higher accuracy obtained by Legendre wavelets decomposition, which is in the form of a multi-resolution analysis of the function. The solution is first found on the coarse grid points, and then it is refined by obtaining higher accuracy with help of increasing the level of wavelets. The accurate implementation of the classical numerical methods on Neumann’s boundary conditions has been found to involve some difficulty. It has been shown here that the present methods can be easily implemented on Neumann’s boundary conditions and the results obtained are accurate; the present methods, thus, have a clear advantage over the classical numerical methods. A distinct feature of the proposed methods is their simple applicability for a variety of boundary conditions. Numerical order of convergence of the proposed methods is calculated. The results of numerical tests show better accuracy of the proposed method based on Legendre wavelets for a variety of benchmark problems.     

One‐step boundary knot method for discontinuous coefficient elliptic equations with interface jump conditions

This study makes the first attempt to apply the boundary knot method (BKM), a meshless collocation method, to the solution of linear elliptic problems with discontinuous coefficients, also known as the elliptic interface problems. The additional jump conditions are usually required to be prescribed at the interface which is used to maintain the well‐posedness of the considered problem. To solve the problem efficiently, the original governing equation is first transformed into an equivalent inhomogeneous modified Helmholtz equation in the present numerical formulation. Then the computational domain is divided into several subdomains, and the solution on each subdomain is approximated using the BKM approach. Unlike the conventional two‐step BKM, this study presents a one‐step BKM formulation which possesses merely one influence matrix for inhomogeneous problems. Several benchmark examples with various discontinuous coefficients have been tested to demonstrate the accuracy and efficiency of the present BKM scheme. © 2016Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1509–1534, 2016     

Solving two-point boundary value problems using combined homotopy perturbation method and Green’s function method

In this paper, an algorithm is proposed for the solution of second-order boundary value problems with two-point boundary conditions. The Green’s function method is applied first to transform the ordinary differential equation into an equivalent integral one, which has already satisfied the boundary conditions. And then, the homotopy perturbation method is used to the resulting equation to construct the numerical solution for such problems. Numerical examples demonstrate the efficiency and reliability of the algorithm developed, it is quite accurate and readily implemented for both linear and nonlinear differential equations with homogeneous and nonhomogeneous boundary conditions. Furthermore, the lower order approximation is of higher accuracy for most cases. Some other extended applications of this algorithm are also exhibited.     

基于MFS的PDE-约束优化法求解热传导逆问题

在本文中,我们给出了一种有效的无网格方法来求解逆热传导问题,含有Neumann边界条件情形.所得到的PDE-约束优化法是一种在空间与时间域上的全局近似方法,其中将控制方程的基本解作为基函数.由于初始测量数据包含有噪声误差,则所得线性方程组的系数矩阵通常是病态的,文中利用广义交叉验证(GCV)的Tikhonov正则化方法来获得更加稳定的数值解.通过数值结果表明,本文给出的方法是精确、有效、鲁棒的.     

奇异杂交边界点方法求解泊松方程

As a boundary-type meshless method,the singular hybrid boundary node method(SHBNM)is based on the modified variational principle and the moving least square(MLS)approximation,so it has the advantages of both boundary element method(BEM)and meshless method.In this paper,the dual reciprocity method(DRM)is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution.The general solution is achieved by means of SHBNM,and the particular solution is approximated by using the radial basis function(RBF).Only randomly distributed nodes on the bounding surface of the domain are required and it doesn't need extra equations to compute internal parameters in the domain.The postprocess is very simple.Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.     

Spectrally accurate energy‐preserving methods for the numerical solution of the “good” Boussinesq equation

In this paper we study the geometric numerical solution of the so called “good” Boussinesq equation. This goal is achieved by using a convenient space semi‐discretization, able to preserve the corresponding Hamiltonian structure, then using energy‐conserving Runge–Kutta methods in the Hamiltonian boundary value method class for the time integration. Numerical tests are reported, confirming the effectiveness of the proposed method.     

On the numerical solution of nonlinear Burgers’-type equations using meshless method of lines

In this paper, a meshless method of lines (MMOL) is proposed for the numerical solution of nonlinear Burgers’-type equations. This technique does not require a mesh in the problem domain, and only a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions (RBFs). The scheme is tested for various examples. The results obtained by this method are compared with the exact solutions and some earlier work.     

Analysis of a meshless method for the time fractional diffusion-wave equation

In this paper a numerical technique is proposed for solving the time fractional diffusion-wave equation. We obtain a time discrete scheme based on finite difference formula. Then, we prove that the time discrete scheme is unconditionally stable and convergent using the energy method and the convergence order of the time discrete scheme is \(\mathcal {O}(\tau ^{3-\alpha })\). Firstly, we change the main problem based on Dirichlet boundary condition to a new problem based on Robin boundary condition and then, we consider a semi-discrete scheme with Robin boundary condition and show when \(\beta \rightarrow +\infty \) solution of the main semi-discrete problem with Dirichlet boundary condition is convergent to the solution of the new semi-discrete problem with Robin boundary condition. We consider the new semi-discrete problem with Robin boundary condition and use the meshless Galerkin method to approximate the spatial derivatives. Finally, we obtain an error bound for the new problem. We prove that convergence order of the numerical scheme based on Galekin meshless is \(\mathcal {O}(h)\). In the considered method the appeared integrals are approximated using Gauss Legendre quadrature formula. The main aim of the current paper is to obtain an error estimate for the meshless Galerkin method based on the radial basis functions. Numerical examples confirm the efficiency and accuracy of the proposed scheme.     

Linear B-spline finite element method for the improved Boussinesq equation

In this paper, we develop and validate a numerical procedure for solving a class of initial boundary value problems for the improved Boussinesq equation. The finite element method with linear B-spline basis functions is used to discretize the nonlinear partial differential equation in space and derive a second order system involving only ordinary derivatives. It is shown that the coefficient matrix for the second order term in this system is invertible. Consequently, for the first time, the initial boundary value problem can be reduced to an explicit initial value problem to which many accurate numerical methods are readily applicable. Various examples are presented to validate this technique and demonstrate its capacity to simulate wave splitting, wave interaction and blow-up behavior.     

Modified Tikhonov regularization method for the Cauchy problem of the Helmholtz equation

In this paper, the Cauchy problem for the Helmholtz equation is investigated. By Green’s formulation, the problem can be transformed into a moment problem. Then we propose a modified Tikhonov regularization algorithm for obtaining an approximate solution to the Neumann data on the unspecified boundary. Error estimation and convergence analysis have been given. Finally, we present numerical results for several examples and show the effectiveness of the proposed method.     

Finite volume and finite element methods for solving a one-dimensional space-fractional Boussinesq equation

In this paper, a new one-dimensional space-fractional Boussinesq equation is proposed. Two novel numerical methods with a nonlocal operator (using nodal basis functions) for the space-fractional Boussinesq equation are derived. These methods are based on the finite volume and finite element methods, respectively. Finally, some numerical results using fractional Boussinesq equation with the maximally positive skewness and the maximally negative skewness are given to demonstrate the strong potential of these approaches. The novel simulation techniques provide excellent tools for practical problems. These new numerical models can be extended to two- and three-dimensional fractional space-fractional Boussinesq equations in future research where we plan to apply these new numerical models for simulating the tidal water table fluctuations in a coastal aquifer.     

Initial boundary value problem for generalized logarithmic improved Boussinesq equation

In this paper, we consider the initial boundary value problem for generalized logarithmic improved Boussinesq equation. By using the Galerkin method, logarithmic Sobolev inequality, logarithmic Gronwall inequality, and compactness theorem, we show the existence of global weak solution to the problem. By potential well theory, we show the norm of the solution will grow up as an exponential function as time goes to infinity under some suitable conditions. Furthermore, for the generalized logarithmic improved Boussinesq equation with damped term, we obtain the decay estimate of the energy. Copyright © 2016 John Wiley & Sons, Ltd.     

A meshless method for Burgers’ equation using MQ-RBF and high-order temporal approximation

In this paper, a new numerical method is proposed to solve one-dimensional Burgers’ equation using multiquadric (MQ) radial basis function (RBF) for spatial approximation and a second-order compact finite difference scheme for temporal approximation. The numerical results obtained by this way for different Reynolds number have been compared with the existing numerical schemes to show the accuracy and efficiency of the approach. To show the superiority of this meshless method, numerical experiments with non-uniform MQ interpolation node distribution are also performed.     

The Laplace transform and polynomial Trefftz method for a class of time dependent PDEs

In this paper, we present a new method for solving 1D time dependent partial differential equations based on the Laplace transform (LT). As a result, the problem is converted into a stationary boundary value problem (BVP) which depends on the parameter of LT. The resulting BVP is solved by the polynomial Trefftz method (PTM), which can be regarded as a meshless method. In PTM, the source term is approximated by a truncated series of Chebyshev polynomials and the particular solution is obtained from a recursive procedure. Talbot’s method is employed for the numerical inversion of LT. The method is tested with the help of some numerical examples.