Subscribe Now Subscribe Today
Research Article
 

A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE



Sayan Kaennakham and Nissaya Chuathong
 
Facebook Twitter Digg Reddit Linkedin StumbleUpon E-mail
ABSTRACT

Background and Objective: The effect of shape parameter is known to play a crucial role in determining the final results of a collocation-based numerical method. This is also the case for boundary element method where radial basis functions are used to collocate the non-homogeneous term. However, finding an optimal shape parameter is known not to be simple particularly when dealing with complex PDE problems. This investigation was carried out focusing on three purposes. Firstly, it is to propose a new form of shape parameter contained in the inverse-multiquadric RBF that behaves both linearly and exponentially. Secondly, it is to integrate the effect of the local phenomena of the problem at hand into the mechanism of the proposed shape via. the local Reynolds number (Re). Thirdly, the methodology of dual reciprocity boundary element method is studied, applied and computationally implemented to one of the most challenging types of PDEs, ‘Burgers equations’, famous for its rich in transient, couple and nonlinear phenomena. Materials and Methods: The study began with gathering mostly-used and proposed forms of shape parameter and analyzing the general aspects in terms of their affectiveness. A new form of shape parameter that was hoped to alleviate the drawbacks commonly found when using those previously proposed shape parameter. The investigation then moved on applying the methodology of the dual reciprocity boundary element method in conjunction with the newly formulated shape parameter to one of the most challenging forms of PDE namely Burgers’ equations. The overall effectiveness of the was evaluated by comparing the results against both their analytical solutions and other numerical works when available in literature. Results: Main findings of this work are as follows. Firstly, the method has successfully been applied to Burgers’ equations using inverse multiquadric radial basis function at relatively high Reynolds number. Secondly, it is found from all the results obtained in this work that the proposed shape parameter can outperform the fixed ones and certainly deserves further investigation. Lastly, when compared with other numerical works, the accuracy lies in an acceptable level and moreover, gets better when the problem become advective dominated, at high Reynolds number. Conclusion: It is found in this work that with its ability to adapt itself locally, the proposed choice of variable shape provides reasonable solutions while requiring only the upper and lower bounding values. This makes choosing the suitable shape much more effective and simpler, particularly when the flow reaches the stage of instability, high Reynolds number.

Services
Related Articles in ASCI
Similar Articles in this Journal
Search in Google Scholar
View Citation
Report Citation

 
  How to cite this article:

Sayan Kaennakham and Nissaya Chuathong, 2017. A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE. Journal of Applied Sciences, 17: 491-501.

DOI: 10.3923/jas.2017.491.501

URL: https://scialert.net/abstract/?doi=jas.2017.491.501
 
Received: August 10, 2017; Accepted: December 14, 2017; Published: February 09, 2018


Copyright: © 2017. This is an open access article distributed under the terms of the creative commons attribution License, which permits unrestricted use, distribution and reproduction in any medium, provided the original author and source are credited.

INTRODUCTION

The radial basis functions (RBFs), ϕ are commonly found as multivariate functions whose values are dependent only on the distance from the origin and commonly assumed to be strictly positive definite. This means that ϕ(x) = ϕ(r)∈R, the set of real numbers, with x∈Rn and r∈R or in other words, on the distance from a point of a given set {xj} and ϕ(x-xj) = ϕ(rj)∈R where can normally defined as expressed in Eq. 1:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(1)

for some fixed points x∈Rn and rj = ||x-xj||2 is the Euclidean distance. Two classes of numerical tools that involve the use of RBFs are the so-called ‘meshfree method’ and ‘boundary element method’. There are many forms of RBFs invented, proposed, tested and documented nowadays under the context of finding solution to different type partial differential equations (PDEs). Below is the list of some well-known and most widely used RBFs:

•  Linear (LR): 1+r
•  Gaussian (GU): Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
•  Cubic (CU) : r3
•  Polyharmonic (PY) in R3: r2n-1, n∈ℕ
•  Polyharmonic (PY) in R2: r2n ln(r), n∈ℕ
•  Multiquadric (MQ): Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
•  Inverse multiquadric (IMQ): Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
•  Thin-plate spline (TPS): r2 ln(r)
•  Matern/sobolev (MS): *[KV (r) rV]
•  *[KV is an order Bessel function]

In this study, nevertheless, the RBF type being focused on is the ‘Inverse-multiquadric (IMQ)’ form expressed as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(2)

where, β = ..., -3/2, -1/2, 1/2, 3/2,... and ε is the so-called ‘shape parameter’ and is normally given in an ‘ad hoc’ manner. Different values of ε, illustrated in Fig. 1, lead to different impact on the final results as well acknowledged. This kind of radial basis function has been proven to be strictly positive definite as nicely documented by Buhmann1, meaning that the distance matrix of the interpolation problem is invertible. Throughout this study, it focused on the most popular form of the inverse multiquadric, i.e., β = -1/2, only.

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
Fig. 1:
Infinitely smooth inverse-multiquadric RBF measured numerically at different shape values, ε

Appearing as an alternative numerical method, finite element, finite volume and finite difference, over the last two decades, the boundary element method (BEM) is now known to be another important tool for solving a wide range of applied sciences and engineering that involve linear as well as certain types of nonlinear partial differential equations (PDEs). Like other numerical schemes, the method is not without difficulties or challenges when it faces problems with nonlinear, transient, coupled and nonhomogeneous form of PDEs.

Brebbia and Butterfield2 proposed an improved version of the scheme and they named it as ‘dual reciprocity boundary element method (DRBEM)’. In the process of DRBEM, the solution is divided into two parts: Complementary solutions of its homogeneous form and the particular solutions of the inhomogeneous counterpart. Since the particular solutions are not always available especially in complex problems, the inhomogeneous term of the PDE is approximated by a series of simple functions and transformed to the boundary integrals employing particular solutions of considered problem. The most widely used approximating functions in DRBEM are radial basis functions (RBFs) for which particular solutions can be easily determined3.

When it comes to using a radial basis function, particularly those forms with the parameter, it is then down to the user’s decision to choose what they believe is best and very often this is done under rather ‘ad hoc’ manner. For multiquadric (MQ) type of RBF, there are some investigations done on providing information regarding choosing optimal shape value. Hardy4 suggests that by fixing the shape at ε = 1/(0.815d), where, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDEand di is the distance from the node to its nearest neighbor, good results should be anticipated. Also, in the work of Franke and Schaback5 where the choice of a fixed shape of the form Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE where, D is the diameter of the smallest circle containing all data nodes, can also be a good alternative.

Some recent attempts to pinpoint the optimal value of ε involve the work of Zhang et al.6 where they demonstrated and concluded that the optimal shape parameter is problem dependent. In 2002, Wang and Lui7 pointed out that by analyzing the condition number of the collocation matrix, a suitable range of derivable values of ε can be found. Later in 2003, Lee et al.8 suggested that the final numerical solutions obtained are found to be less affected by the method when the approximation is applied locally rather than globally.

At this point, nevertheless, there are several facts to be noticed at. Firstly, most of the proposed forms are concerned with either numerical aspect or fixed values only, meaning that no real physical aspects included or being under consideration. Secondly, all the works were done under the context of pure RBF-collocation methods, nothing is done under DRBEM. Thirdly, the shape parameter contained in the inverse-multiquadric type of RBF has not, by far, been taken a look at. This prompts our initiative idea of this study.

One of the classical non-linear parabolic equation system widely acknowledged as the Burgers’ equations, named after the great Physicist9. The equations retain the nonlinear of the governing equation in a number of applications, flow through a shock wave traveling in viscous fluid, the phenomena of turbulence, air/water pollution or chemical compound carried by fluids, sedimentation of two kinds of particles in fluid suspensions under the effect of gravity (Burgers9, Nee and Duan10).

For large values of the Reynolds number (Re), one of the major difficulties is due to in-viscid boundary layers produced by the steepening effect of the nonlinear advection term in Burgers’ equations. This difficulty is also encountered in an inviscid Navier–Stokes equation for a convection dominated flow and in fact, the famous coupled-Burgers’ equations are one of the principle equations used to evaluate any newly proposed numerical methods11.

This all leads to our main objectives of this investigation which are, firstly, it proposed a new form of shape parameter contained in the inverse-multiquadric RBF. Secondly, the proposed shape was designed to take into consideration the local physical of the problem at hand governed by the dimensionless Reynolds number (Re). Thirdly, the methodology of DRBEM was studied, applied and computationally implemented to one of the most challenging types of PDEs known as ‘Burgers’ equation’, famous for its rich in transient, couple and nonlinear phenomena. The results obtained from this investigation were validated against other alternative numerical works in literature when available.

MATERIALS AND METHODS

DRBEM for Burgers’ equations: Let us consider the system of two-dimensional unsteady, nonlinear and coupled Burgers’ equations, expressed as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(3)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(4)

Subject to the initial conditions:

u(x, y, 0) = β1(x, y) and v(x, y, 0) = β2 (x, y)

and the boundary conditions:

u(x, y, t) = γ1(x, y, t) and v(x, y, t) = γ2 (x, y, t)

for (x, y)∈∂Ω, t>0, u(x, y, t) and v(x, y, t) are the velocity components, β1, β2, γ1 and γ2 are known functions and Re is the Reynolds number, described by Keannakham et al.12

The mathematical construction of the dual reciprocity boundary element method can start with the Poisson equation as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(5)

Which as its equivalent integral form, given by Brebbia and Butterfield2 as:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(6)

where, u* the fundamental solution and the term is defined as Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE, where n is the unit outward normal to Γ and can be written as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(7)

Next, the method was applied with N and L being the number of boundary and internal nodes, respectively, b can be now approximated by:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Or, can be written as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(8)

Or:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

where, the function f is the radial basis function which is, in this work, the inverse-multiquadric type. With this radial basis function, then the following form was obtained:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(9)

for some particular solution, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE Applying Green’s theorem, the boundary element approximation to Eq. 9, then it becomes, at a node ith:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(10)

For i = 1, 2,..., N.

After introducing the interpolation function and integrating over each boundary elements, the above Eq. 10 can be re-written in terms of nodal values as:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(11)

where, the definition of the terms Hik and Gik are described by Toutip13. The index k is used for the boundary nodes which are the field points. After application to all boundary nodes, using a collocation technique, Eq. 11 can be compactly expressed in matrix form as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(12)

By substituting α = F–1b from Eq. 8, into Eq. 12 making the right hand side of Eq. 12 a known vector. Therefore, it can be rewritten as:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(13)

where, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE Applying boundary conditions to Eq. 13, then it can be seen as the simple form as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(14)

where, x contains N unknown boundary values of u’s and q’s.

After Eq. 14 is solved using standard techniques such as Gaussian elimination, the values at any internal node can be calculated from Eq. 11, i.e., ci = 1 as expressed in Eq. 15 where each one involving a separate multiplication of known vectors and matrices:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(15)

Firstly, substituting Eq. 12 into Eq. 11 to get the equation system matrix which expressed as:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(16)

It is then possible to approximate the terms, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE and Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE as below:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(17)

From Eq. 16 and Eq. 17, the following was obtained:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(18)

and:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(19)

where, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE and:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(20)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(21)

In the case, if Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE were considered to be generated by using the same redial basis function, then:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(22)

For the time derivatives, the forward difference method was used to approximate time derivative Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE and Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDENow, substituting Eq. 32-35, 38 and 39 in Eq. 43 and 44, then the following was obtained:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(23)

and

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(24)

by setting Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE and Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE then the final forms of DRBEM for this type of equations A = πr2 as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(25)

and

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(26)

Note that the elements of matrices H, G and A depend only on geometrical data. Thus, they can all be computed once and stored.

Table 1:
Some choices proposed in literature
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Proposed variable IMQ shape: Regarding the studies in the search for optimal choice of the shape parameter, many outstanding and well-known forms proposed in the past are listed in Table 1. In this table, it can be seen that one of the first attempts in searching for variable form of shape parameter was done in 1990 by Kansa14 which the experimental manner was looked at. This form was then modified by Kansa and Carlson15. In their work, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE and Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE were observed that the better quality of results can be achieved if these two values vary by several orders of magnitude. Sarra16 and Sarra and Sturgill17 later invented a linear form of variable shape parameter and applied to both interpolation and some benchmark partial differential equations. The command ‘rand’ is the MATLAB function that returns uniformly distributed pseudo-random numbers on the unit interval. In their study, the variable form of shape parameter was shown to be superior to the fixed form particularly when the system includes some information about the minimum distance of a center to its nearest neighbor hn, with also a user input value μ. It is, nevertheless, to be noticed that all the attempts mentioned so far are under the context of the collocation meshfree method where there is no, by far, work carried out under the DRBEM. This leads to one of our main objectives of this investigation.

In this study, it proposed a new form of variable shape parameter where both linear and exponential manners are taken into consideration, expressed as in Eq. 27-30:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(27)

Where:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(28)

and

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(29)

and ζ is set to correspond to local Reynolds number (Reij) and is defined as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(30)

and the dimensionless number called The Reynolds Number, (Re), which is the ratio of the inertia to the viscous or friction force and is defined, in a global sense, as:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(31)

where, L is the characteristic length scale, U is the velocity and v is the kinematic viscosity. In this investigation, this number is locally defined as:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(32)

where, vt is defined in a local sense as:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(33)

where, N+L is the total number of nodes (both boundary and internal) involved in the system, the superscript t indicates the time level. The subscript i and j mean at ith and jth node, respectively.

Numerical setup: The computation process is done following the below algorithm:

•  Initializing the matrix C by setting it to zero for both U and V
•  Setting all the values appearing in Eq. 25 and 26 and then applying the initial conditions to the right hand side and boundary condition to the left hand side, reducing to the form Ax = y
•  Using the iterative scheme as A(k)x(k) = y(k), where, k = 1, 2, 3.... The solutions of for both U and Vcan then be achieved when the algorithm stops i.e. when it reaches the following condition:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(34)

with an user’s given value of a tolerance ω.

•  The solutions obtained from step 3 will be used as the initial conditions in next time-step and the computing process continues until the time level given is reached

All numerical simulations under this investigation were carried out on an Intel(R) Core(TM) i7-5500U CPU @ 2.4 GHz with 8.00 GB of RAM.

RESULTS

To demonstrate the effectiveness of the proposed shape parameter, the methodology was applied to two benchmark examples governed by the Burgers’ form of equations. The results obtained were validated against both the corresponding exact solutions and those published in literature. In some cases, the root mean square (RMS) error norm was used and it was defined as below:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(35)

It was noted here that the Reynolds number written as ‘Re’ from now on is the global Reynolds number described and used in Eq. 31-33.

Experiment 1: The first example was the well-known form where its exact solutions for validation were provided by using a Hopf-Cole transformation nicely done by Fletcher18 and was expressed as follows:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(36)

and:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(37)

where, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE and both the initial and boundary conditions to be imposed to the equation system was generated directly from the above exact form over the domain Ω = {(x, y):0<x<1, 0<y<1}.

Table 2-9 showed the computed values of both U- and V-values for Re = 1 at t = 0.01 and Re = 10 at t = 0.5, respectively, obtained from this work compared against both its’ exact solutions and those provided by Biazar and Aminikhah19. In this case, only 13 internal nodes were chosen and it can be seen that the method constructed in this study performed well and results agreed nicely with both references.

Table 2:
Comparison of U-velocity for Re = 1, t = 0.01 with boundary node N = 80 and internal nodes L = 13, for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 3:
Comparison of U-velocity for Re = 1, t = 0.01 with boundary nodes N = 80 and internal nodes L = 13, for (εmin, εmax) = (0.1, 10)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 4:
Comparison of V-velocity for Re = 1, t = 0.01 with boundary nodes N = 80 and internal nodes L = 13, for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 5:
Comparison of V-velocity for Re = 1, t = 0.01 with boundary nodes N = 80 and internal nodes L = 13, for (εmin, εmax) = (0.1, 10)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 6:
Comparison of U-velocity for Re = 10, t = 0.5 with boundary nodes N = 80 and internal nodes L = 13 for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 7:
Comparison of U-velocity for Re = 10, t = 0.5 with boundary nodes N = 80 and internal nodes L = 13, for (εmin, εmax) = (0.1, 10)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 8:
Comparison of V-values for Re = 10, t = 0.5 with boundary nodes N = 80 and internal nodes L = 13, for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 9:
Comparison of V-values for Re = 10, t = 0.5 with boundary nodes N = 80 and internal nodes L = 13, for (εmin, εmax) = (0.1, 10)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
Fig. 2(a-b):
Comparison of V-values for (εmin, εmax) = (0.1, 10) at t = 0.1 with Re = 80 between the exact (left) and the computed ones (right)

Solution profiles of both U- and V-velocity, at moderate Reynolds number, were illustrated in Fig. 2-4 showed the root mean square (RMS) error.

Experiment 2: For this example, it was adopted from the work of Aminikhah20 for solution validation. This problem is on the domain Ω = {(x, y): 0<x<1, 0<y<1)} with the following initial and boundary conditions:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(38)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(39)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(40)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(41)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
Fig. 3(a-b):
Comparison of U-values for (εmin, εmax) = (0.1, 10) at t = 0.1 with Re = 80 between the exact (left) and the computed ones (right)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
Fig. 4(a-b):
Root mean square (RMS) error comparison at two intervals of (εmin, εmax)and that obtained using a fixed value of 3, computed with boundary nodes N = 120 and internal nodes L = 24 at t = 1.0

This leads to the exact solutions as the following forms:

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(42)

Table 10:
Comparison of U-velocity for Re = 100, t = 0.5 with boundary nodes N = 80 and internal nodes L = 12, for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 11:
Comparison of V-values for Re = 100, t = 0.5 with boundary nodes N = 80 and internal nodes L = 12, for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 12:
Comparison of U-velocity for Re = 500, t = 0 with boundary nodes N = 80 and internal nodes L = 12, for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

and

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
(43)

where, Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE = exp (-5π2 t/Re).

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
Fig. 5(a-b):
Comparison of V-velocity for (εmin, εmax) = (0.1, 20) at t = 1 with Re = 1,000 between the exact (left) and the computed ones (right)

Table 13:
Comparison of V-velocity for Re = 500, t = 0.5 with boundary nodes N = 80 and internal nodes L = 12, for (εmin, εmax) = (0.1, 1)
Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE

Table 10-13 provided results comparison amongst those produced in this study and both the exact solutions and other numerical investigation, Aminikhah20. The solution profiles obtained at high Reynolds number were shown in Fig. 5-6 where Fig. 7 provided the comparison of RMS at different values of both Reynolds numbers and ranges of (εmin, εmax).

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
Fig. 6(a-b):
Comparison of U-velocity for (εmin, εmax) = (0.1, 20) at t = 1 with Re = 1,000 between the exact (left) and the computed ones (right)

Image for - A Proposed Adaptive Inverse Multiquadric Shape Parameter Applied with the Dual Reciprocity BEM to Nonlinear and Coupled PDE
Fig. 7(a-b):
Root mean square (RMS) error comparison at two intervals of (εmin, εmax) and that obtained using a fixed value of 1.5, computed with boundary nodes N = 120 and internal nodes L = 24 at t = 0.75

DISCUSSION

Regarding the use of inverse multiquadric (IMQ) type of radial basis functions (RBFs), it is not surprising to realize that only a few works available in literature. Most of these works, moreover, were carried out under the context of meshfree/meshless method. Some recent works that adopted this type of RBF are Chantawara21 and Chuathong22, where some fixed values of shape parameters were adopted. To validate results obtained in this paper, their suggested values of shape parameter were adopted.

The computational results produced in this work for the two benchmark examples are promising in several aspects. Firstly, by providing only two thresholds or bounding values of (εmin, εmax), searching for optimal choice becomes simpler. This feature can be seen more clearly when compared with those obtained by using a fixed value of ε, as shown in all Tables. The effect of the wide of the range (εmin, εmax) was seen not to have any significant effect on the final results quality ensuring the simplicity of the proposed shape form.

Another aspect that is desirable about the main findings of this study is the capability of the method to handle the instability of the problem when the Reynolds number increases. It is commonly found and is known to be a great challenge for any numerical procedure to tackle Burgers’ equations at high Reynolds number (>500), the flow becomes more turbulent. One of our previous works and many other works, Keannakham et al.12 and references herein, have shown that beyond this value of Reynolds number, the error norm is growing dramatically. It is shown clearly in this work that the proposed adaptive parameter is able to take into consideration the local physics that effects both the computing collocation matrix and the numerical algorithm itself.

CONCLUSION

In this study, the effectiveness of shape parameter contained in the inverse multiquadric type of RBF in conjunction with the method of DRBEM was investigated numerically. The investigation began with applying DRBEM to one of the most complicated PDEs namely Burgers’ equations. Then a new form of shape parameter that behaves locally-adaptive, i.e., it varies accordingly to the local change of the physics of the problem which, here, is the local Reynolds number(Re) was proposed. The proposed variable shape form also contains both linear and exponential aspects based on the distance between the center node i and pointed node j in the RBF-collocation numerical method adopted. Some important conclusions can now be withdrawn from the investigation and they are as follows:

•  DRBEM has successfully been applied to Burgers’ equations using inverse multiquadric radial basis function at relatively high Reynolds number (Re)
•  It is found from all the results obtained in this work that the proposed shape parameter can outperform the fixed ones and certainly deserves further investigation
•  When compared with other numerical works, the accuracy lies in an acceptable level and moreover, gets better when the problem become advective dominated, at high Reynolds number

SIGNIFICANCE STATEMENT

This study discovers the benefit of using the local physical aspect of the problem as the indicator for the radial basis function’s shape parameter which is designed to be locally adaptive. This study will help the researcher to uncover the critical area of choosing the optimal shape parameter when dealing with this particular form of equations. Thus, a new theory on these radial basis function’s adaptive shape parameters and possibly other related fields, may be arrived at.

ACKNOWLEDGMENTS

The corresponding author would like to express his sincere gratitude and appreciation to the Center of Excellence in Mathematics, Thailand for their financial support.

REFERENCES
1:  Buhmann, M.D., 1990. Multivariate cardinal interpolation with radial-basis functions. Construct. Approximation, 6: 225-225.
Direct Link  |  

2:  Brebbia, C.A. and R. Butterfield, 1978. Formal equivalence of direct and indirect boundary element methods. Applied Math. Model., 2: 132-134.
CrossRef  |  Direct Link  |  

3:  Carlson, R.E. and T.A. Foley, 1991. The parameter R2 in multiquadric interpolation. Comput. Math. Applic., 21: 29-42.
CrossRef  |  Direct Link  |  

4:  Hardy, R.L., 1971. Multiquadric equations of topography and other irregular surfaces. J. Geophys. Res., 76: 1905-1915.
CrossRef  |  Direct Link  |  

5:  Franke, C. and R. Schaback, 1998. Convergence order estimates of meshless collocation methods using radial basis functions. Adv. Comput. Math., 93: 381-392.
CrossRef  |  Direct Link  |  

6:  Zhang, X., K.Z. Song, M.W. Lu and X. Liu, 2000. Meshless methods based on collocation with radial basis functions. Comput. Mech., 26: 333-343.
CrossRef  |  Direct Link  |  

7:  Wang, J.G. and G.R. Liu, 2002. On the optimal shape parameters of radial basis functions used for 2-D meshless methods. Comput. Meth. Applied Mech. Eng., 191: 2611-2630.
CrossRef  |  Direct Link  |  

8:  Lee, C.K., X. Liu and S.C. Fan, 2003. Local multiquadric approximation for solving boundary value problems. Comput. Mech., 30: 396-409.
CrossRef  |  Direct Link  |  

9:  Burgers, J.M., 1948. A mathematical model illustrating the theory of turbulence. Adv. Applied Mech., 11: 171-199.
CrossRef  |  Direct Link  |  

10:  Nee, J. and J. Duan, 1998. Limit set of trajectories of the coupled viscous Burgers' equations. Applied Math. Lett., 11: 57-61.
CrossRef  |  Direct Link  |  

11:  Zhang, D.S., G.W. Wei, D.J. Kouri and D.K. Hoffman, 1997. Burgers' equation with high Reynolds number. Phys. Fluids, Vol. 9. 10.1063/1.869305

12:  Keannakham, S., K. Chantawara and W. Toutip, 2014. [Optimal Radial Basis Function (RBF) for Dual Reciprocity Boundary Element Method (DRBEM) applied to coupled Burgers' equations with increasing Reynolds number]. Aust. J. Basic Applied Sci., 8: 462-476.

13:  Toutip, W., 2001. The dual reciprocity boundary element method for linear and nonlinear problems. Ph.D. Thesis, University of Hertfordshire, Hertfordshire, United Kingdom.

14:  Kansa, E.J., 1990. Multiquadrics: A scattered data approximation scheme with applications to computational fluid-dynamics-I surface approximations and partial derivative estimates. Comput. Math. Applied, 19: 127-145.
CrossRef  |  Direct Link  |  

15:  Kansa, E.J. and R.E. Carlson, 1992. Improved accuracy of multiquadric interpolation using variable shape parameters. Comput. Math. Applic., 24: 99-120.
CrossRef  |  Direct Link  |  

16:  Sarra, S.A., 2005. Adaptive radial basis function methods for time dependent partial differential equations. Applied Numer. Math., 54: 79-94.
CrossRef  |  Direct Link  |  

17:  Sarra, S.A. and D. Sturgill, 2009. A random variable shape parameter strategy for radial basis function approximation methods. Eng. Anal. Bound. Elements, 33: 1239-1245.
CrossRef  |  Direct Link  |  

18:  Fletcher, C.A., 1983. Generating exact solutions of the two‚Äźdimensional Burgers' equations. Int. J. Numer. Meth. Fluids, 3: 213-216.
CrossRef  |  Direct Link  |  

19:  Biazar, J. and H. Aminikhah, 2009. Exact and numerical solutions for non-linear Burger's equation by VIM. Math. Comput. Model., 49: 1394-1400.
CrossRef  |  Direct Link  |  

20:  Aminikhah, H., 2012. A new efficient method for solving two-dimensional Burgers' equation. ISRN Comput. Math., 10.5402/2012/603280

21:  Chantawara, K., 2016. A numerical comparison of radial basis functions in applications of the dual reciprocity boundary element method to convection-diffusion problems. Ph.D. Thesis, Khon Kaen University, Thailand.

22:  Chuathong, N., 2016. An automatic Centroid-Node Adaptive Meshless (CNAM) method for solving convection-diffusion problems. Ph.D. Thesis, Khon Kaen University, Thailand.

©  2021 Science Alert. All Rights Reserved