Numerical investigation on convergence of boundary knot method in the analysis of homogeneous Helmholtz, modified Helmholtz, and convection-diffusion problems

This paper concerns a numerical study of convergence properties of the boundary knot method (BKM) applied to the solution of 2D and 3D homogeneous Helmholtz, modified Helmholtz, and convection-diffusion problems. The BKM is a new boundary-type, meshfree radial function basis collocation technique. The method differentiates from the method of fundamental solutions (MFS) in that it does not need the controversial artificial boundary outside physical domain due to the use of non-singular general solutions instead of the singular fundamental solutions. The BKM is also generally applicable to a variety of inhomogeneous problems in conjunction with the dual reciprocity method (DRM). Therefore, when applied to inhomogeneous problems, the error of the DRM confounds the BKM accuracy in approximation of homogeneous solution, while the latter essentially distinguishes the BKM, MFS, and boundary element method. In order to avoid the interference of the DRM, this study focuses on the investigation of the convergence property of the BKM for homogeneous problems. The given numerical experiments reveal rapid convergence, high accuracy and efficiency, mathematical simplicity of the BKM.     

The solitary wave solution of the two-dimensional regularized long-wave equation in fluids and plasmas

This paper investigates the solitary wave solutions of the two-dimensional regularized long-wave equation which is arising in the investigation of the Rossby waves in rotating flows and the drift waves in plasmas. The main idea behind the numerical solution is to use a combination of boundary knot method and the analog equation method. The boundary knot method is a meshless boundary-type radial basis function collocation technique. In contrast with the method of fundamental solution, the boundary knot method uses the non-singular general solution instead of the singular fundamental solution to obtain the homogeneous solution. Similar to method of fundamental solution, the radial basis function is employed to approximate the particular solution via the dual reciprocity principle. In the current paper, we applied the idea of analog equation method. According to the analog equation method, the nonlinear governing operator is replaced by an equivalent nonhomogeneous linear one with known fundamental solution and under the same boundary conditions. Furthermore, in order to show the efficiency and accuracy of the proposed method, the present work is compared with finite difference scheme. The new method is analyzed for the local truncation error and the conservation properties. The results of several numerical experiments are given for both the single and double-soliton waves.     

The dual reciprocity boundary element formulation for nonlinear diffusion problems

This paper presents an extension of the dual reciprocity boundary element method (DRBEM) to deal with nonlinear diffusion problems in which thermal conductivity, specific heat, and density coefficients are all functions of temperature. The DRBEM, recently applied to the solution of problems governed by parabolic and hyperbolic equations, consists in the transformation of the differential equation into an integral equation involving boundary integrals only, the solution of which is achieved by employing a standard boundary element discretization coupled with a two-level finite difference time integration scheme. Contrary to previous formulations for the diffusion equation, the dual reciprocity BEM utilizes the well-known fundamental solution to Laplace's equation, which is space-dependent only. This avoids complex time integrations that normally appear in formulations employing time-dependent fundamental solutions, and permits accurate numerical solutions to be obtained in an efficient way. For nonlinear problems, the integral of conductivity is introduced as a new variable to obtain a linear diffusion equation in the Kirchhoff transform space. This equation involves a modified time variable which is itself a function of position. The problem is solved in an iterative way by using an efficient Newton-Raphson technique which is shown to be rapidly convergent.     

The numerical solution of fractional differential equations using the Volterra integral equation method based on thin plate splines

Finding the approximate solution of differential equations, including non-integer order derivatives, is one of the most important problems in numerical fractional calculus. The main idea of the current paper is to obtain a numerical scheme for solving fractional differential equations of the second order. To handle the method, we first convert these types of differential equations to linear fractional Volterra integral equations of the second kind. Afterward, the solutions of the mentioned Volterra integral equations are estimated using the discrete collocation method together with thin plate splines as a type of free-shape parameter radial basis functions. Since the scheme does not need any background meshes, it can be recognized as a meshless method. The proposed approach has a simple and computationally attractive algorithm. Error analysis is also studied for the presented method. Finally, the reliability and efficiency of the new technique are tested over several fractional differential equations and obtained results confirm the theoretical error estimates.

     

A localized approach for the method of approximate particular solutions

The method of approximate particular solutions (MAPS) has been recently developed to solve various types of partial differential equations. In the MAPS, radial basis functions play an important role in approximating the forcing term. Coupled with the concept of particular solutions and radial basis functions, a simple and effective numerical method for solving a large class of partial differential equations can be achieved. One of the difficulties of globally applying MAPS is that this method results in a large dense matrix which in turn severely restricts the number of interpolation points, thereby affecting our ability to solve large-scale science and engineering problems.In this paper we develop a localized scheme for the method of approximate particular solutions (LMAPS). The new localized approach allows the use of a small neighborhood of points to find the approximate solution of the given partial differential equation. In this paper, this local numerical scheme is used for solving large-scale problems, up to one million interpolation points. Three numerical examples in two-dimensions are used to validate the proposed numerical scheme.     

Numerical computations for singular semilinear elliptic boundary value problems

The paper studies a class of Dirichlet problems with homogeneous boundary conditions for singular semilinear elliptic equations in a bounded smooth domain in

. A numerical method is devised to construct an approximate Green's function by using radial basis functions and the method of fundamental solutions. An estimate of the error involved is also given. A weak solution of the above given problem is a solution of its corresponding nonlinear integral equation. A computational method is given to find the minimal weak solution U, and the critical index λ* (such that a weak solution U exists for λ < λ*, and U does not exist for λ > λ*).     

Applications of the complete multiple reciprocity method for solving the 1D Helmholtz equation of a semi-infinite domain

In this paper, the complete multiple reciprocity method is adopted to solve the one-dimensional (1D) Helmholtz equation for the semi-infinite domain. In order to recover the information that is missing when the conventional multiple reciprocity method is used, an appropriate complex number in the zeroth order fundamental solution is added such that the kernels derived using this proposed method are fully equivalent to those derived using the complex-valued formulation. Two examples including the Dirichlet and Neumann boundary conditions are investigated to show the validity of the proposed method analytically and numerically. The numerical results show good agreement with the analytical solutions.     

Application of the method of fundamental solutions and radial basis functions for inverse transient heat source problem

The paper considers the determination of heat sources in unsteady 2-D heat conduction problem. The determination of the strength of a heat source is achieved by using the boundary condition, initial condition and a known value of temperature in chosen points placed inside the domain. For the solution of the inverse problem of identification of the heat source the θ-method with the method of fundamental solution and radial basis functions is proposed. Due to ill conditioning of the inverse transient heat conduction problem the Tikhonov regularization method based on SVD decomposition was used. In order to determine the optimum value of the regularization parameter the L-curve criterion was used. For testing purposes of the proposed algorithm the 2-D inverse boundary-initial-value problems in square region Ω with the known analytical solutions are considered. The numerical results show that the proposed method is easy to implement and pretty accurate. Moreover the accuracy of the results does not depend on the value of the θ parameter and is greater in the case of the identification of the temperature field than in the case of the identification of the heat sources function.     

Local radial basis function–finite-difference method to simulate some models in the nonlinear wave phenomena: regularized long-wave and extended Fisher–Kolmogorov equations

In this investigation, we concentrate on solving the regularized long-wave (RLW) and extended Fisher–Kolmogorov (EFK) equations in one-, two-, and three-dimensional cases by a local meshless method called radial basis function (RBF)–finite-difference (FD) method. This method at each stencil approximates differential operators such as finite-difference method. In each stencil, it is necessary to solve a small-sized linear system with conditionally positive definite coefficient matrix. This method is relatively efficient and has low computational cost. How to choose the shape parameter is a fundamental subject in this method, since it has a palpable effect on coefficient matrix. We will employ the optimal shape parameter which results from algorithm of Sarra (Appl Math Comput 218:9853–9865, 2012). Then, we compare the approximate solutions acquired by RBF–FD method with theoretical solution and also with results obtained from other methods. The numerical results show that the RBF–FD method is suitable and robust for solving the RLW and EFK equations.

     

An operator splitting-radial basis function method for the solution of transient nonlinear Poisson problems

This paper presents an operator splitting-radial basis function (OS-RBF) method as a generic solution procedure for transient nonlinear Poisson problems by combining the concepts of operator splitting, radial basis function interpolation, particular solutions, and the method of fundamental solutions. The application of the operator splitting permits the isolation of the nonlinear part of the equation that is solved by explicit Adams-Bashforth time marching for half the time step. This leaves a nonhomogeneous, modified Helmholtz type of differential equation for the elliptic part of the operator to be solved at each time step. The resulting equation is solved by an approximate particular solution and by using the method of fundamental solution for the fitting of the boundary conditions. Radial basis functions are used to construct approximate particular solutions, and a grid-free, dimension-independent method with high computational efficiency is obtained. This method is demonstrated for some prototypical nonlinear Poisson problems in heat and mass transfer and for a problem of transient convection with diffusion. The results obtained by the OS-RBF method compare very well with those obtained by other traditional techniques that are computationally more expensive. The new OS-RBF method is useful for both general (irregular) two- and three-dimensional geometry and provides a mesh-free technique with many mathematical flexibilities, and can be used in a variety of engineering applications.     

Dual reciprocity boundary element method in Laplace domain applied to anisotropic dynamic crack problems

In this paper the dual reciprocity boundary element method in the Laplace domain for anisotropic dynamic fracture mechanic problems is presented. Crack problems are analyzed using the subregion technique. The dynamic stress intensity factors are computed using traction singular quarter-point elements positioned at the tip of the crack. Numerical inversion from the Laplace domain to the time domain is achieved by the Durbin method. Numerical examples of dynamic stress intensity factor evaluation are considered for symmetric and non-symmetric problems. The influence of the number of Laplace parameters and internal points in the solution is investigated.     

Meshless method to solve nonstationary heat conduction problems using atomic radial basis functions

The authors consider a meshless method to solve 3D nonstationary boundary-value heat conduction problems. It is implemented through an iterative scheme based on a combination of the double substitution method and the method of fundamental solutions with the use of atomic radial basis functions. The approaches to the visualization of the desired solution are considered.     

DRM multidomain mass conservative interpolation approach for the BEM solution of the two-dimensional navier-stokes equations

The dual reciprocity method (DRM) is a technique to transform the domain integrals that appear in the boundary element method into equivalent boundary integrals. In this approach, the nonlinear terms are usually approximated by an interpolation applied to the convective terms of the Navier-Stokes equations. In this paper, we introduce a radial basis function interpolation scheme for the velocity field, that satisfies the continuity equation (mass conservative). The proposed method performs better than the classical interpolation used in the DRM approach to represent such a field. The new scheme together with a subdomain variation of the dual reciprocity method allows better approximation of the nonlinear terms in the Navier-Stokes equations.     

Efficient MFS Algorithms for Problems in Thermoelasticity

We propose efficient fast Fourier transform (FFT)-based algorithms using the method of fundamental solutions (MFS) for the numerical solution of certain problems in planar thermoelasticity. In particular, we consider problems in domains possessing radial symmetry, namely disks and annuli and it is shown that the MFS matrices arising in such problems possess circulant or block-circulant structures. The solution of the resulting systems is facilitated by appropriately diagonalizing these matrices using FFTs. Numerical experiments demonstrating the applicability of these algorithms are also presented.     

Atomic radial basis functions in numerical algorithms for solving boundary-value problems for the Laplace equation

A numerical method for solution of boundary-value problems of mathematical physics is described that is based on the use of radial atomic basis functions. Atomic functions are compactly supported solutions of functional-differential equations of special form. The convergence of this numerical method is investigated for the case of using an atomic function in solving the Dirichlet boundary-value problem for the Laplace equation. __________ Translated from Kibernetika i Sistemnyi Analiz, No. 4, pp. 165–178, July–August 2008.     

A collocation and Cartesian grid methods using new radial basis function to solve class of partial differential equations

In the present paper, a class of partial differential equation represented by Poisson's type problems are solved using a proposed Cartesian grid method and a collocation technique using a new radial basis function. The advantage of using this new radial basis function represented by overcoming singularity from the diagonal elements when thin plate radial basis function is used. The new function is a combination of both multiquadric and thin plate radial basis functions. The new radial basis function contains a control parameter ?, that takes one when evaluating the singular elements and equals zero elsewhere. Collocation of the approximate solution of the potential over the governing and boundary condition equations leads to a double linear system of equations. A proposed algebraic procedure is then developed to solve the double system. Examples of Poisson and Helmholtz equations are solved and the present results are compared with the their analytical solutions. A good agreement with analytical results is achieved.     

Structural shape optimisation using boundary elements and the biological growth method

A numerical evolutionary procedure for the structural optimisation for stress reduction of two-dimensional structures is presented in this paper. The proposed procedure couples the biological growth method (BGM) with the boundary element method (BEM). The boundary-only intrinsic characteristic of BEM together with its accuracy in the boundary displacement and stress solutions make BEM especially attractive for solving shape-optimisation problems. Two formulations of BEM are used in this work: the standard for two-dimensional elastostatics for the stress analysis and the dual reciprocity method (DRM), which is used to model the swelling or shrinking of the material. Two examples are analysed to illustrate the proposed methodology and to demonstrate its versatility and robustness.     

The numerical solution of nonlinear generalized Benjamin–Bona–Mahony–Burgers and regularized long-wave equations via the meshless method of integrated radial basis functions

In this paper, a numerical technique is proposed for solving the nonlinear generalized Benjamin–Bona–Mahony–Burgers and regularized long-wave equations. The used numerical method is based on the integrated radial basis functions (IRBFs). First, the time derivative has been approximated using a finite difference scheme. Then, the IRBF method is developed to approximate the spatial derivatives. The two-dimensional version of these equations is solved using the presented method on different computational geometries such as the rectangular, triangular, circular and butterfly domains and also other irregular regions. The aim of this paper is to show that the integrated radial basis function method is also suitable for solving nonlinear partial differential equations. Numerical examples confirm the efficiency of the proposed scheme.

     

Collocation method to solve inequality-constrained optimal control problems of arbitrary order

In this paper, the generalized fractional order of the Chebyshev functions (GFCFs) based on the classical Chebyshev polynomials of the first kind is used to obtain the solution of optimal control problems governed by inequality constraints. For this purpose positive slack functions are added to inequality conditions and then the operational matrix for the fractional derivative in the Caputo sense, reduces the problems to those of solving a system of algebraic equations. It is shown that the solutions converge as the number of approximating terms increases, and the solutions approach to classical solutions as the order of the fractional derivatives approach one. The applicability and validity of the method are shown by numerical results of some examples, moreover a comparison with the existing results shows the preference of this method.

     

Algebraic multigrid methods based on element preconditioning

This paper presents a new algebraic multigrid (AMG) solution strategy for large linear systems with a sparse matrix arising from a finite element discretization of some self-adjoint, second order, scalar, elliptic partial differential equation. The AMG solver is based on Ruge/Stübens method. Ruge/Stübens algorithm is robust for M-matrices, but unfortunately the “region of robustness“ between symmetric positive definite M-matrices and general symmetric positive definite matrices is very fuzzy.

For this reason the so-called element preconditioning technique is introduced in this paper. This technique aims at the construction of an M-matrix that is spectrally equivalent to the original stiffness matrix. This is done by solving small restricted optimization problems. AMG applied to the spectrally equivalent M-matrix instead of the original stiffness matrix is then used as a preconditioner in the conjugate gradient method for solving the original problem.

The numerical experiments show the efficiency and the robustness of the new preconditioning method for a wide class of problems including problems with anisotropic elements.