Complexity
Volume 2021 (2021), Article ID 9952562, 10 pages
https://doi.org/10.1155/2021/9952562

Research Article

Numerical Solution of the Multiterm Time-Fractional Model for Heat Conductivity by Local Meshless Technique

Bander N. Almutairi,¹ Ahmed E. Abouelregal,2,3 Bandar Bin-Mohsin,¹ M. D. Alsulami,⁴ and Phatiphat Thounthong⁵

¹Department of Mathematics, College of Science, King Saud University, Riyadh 11451, Saudi Arabia
²Department of Mathematics, College of Science and Arts, Jouf University, Al-Qurayyat, Saudi Arabia
³Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
⁴University of Jeddah, College of Sciences and Arts at Alkamil, Department of Mathematics, Jeddah, Saudi Arabia
⁵Renewable Energy Research Centre, Department of Teacher Training in Electrical Engineering, Faculty of Technical Education, King Mongkut’s University of Technology North Bangkok, 1518 Pracharat 1 Road, Bangsue, Bangkok 10800, Thailand

Correspondence should be addressed to Ahmed E. Abouelregal; ahabogal@gmail.com

Received 24 March 2021; Revised 26 April 2021; Accepted 9 May 2021; Published 11 June 2021

Academic Editor: Nehad Ali Shah

Copyright © 2021 Bander N. Almutairi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Fractional partial differential equation models are frequently used to several physical phenomena. Despite the ability to express many complex phenomena in different disciplines, researchers have found that multiterm time-fractional PDEs improve the modeling accuracy for describing diffusion processes in contrast to the results of a single term. Nowadays, it attracts the attention of the active researchers. The aim of this work is concerned with the approximate numerical solutions of the three-term time-fractional Sobolev model equation using computationally attractive and reliable technique, known as a local meshless method. Because of the meshless character and the simple application in higher dimensions, there is a growing interest in meshless techniques. To assess the reliability and accuracy of the proposed method, three test problems and two types of irregular domains are taken into account.

1. Introduction

In recent years, fractional partial differential equations (FPDEs) have drawn the consideration of numerous researchers to their applications in various fields of science and technology. Partial derivatives provide a flexible model and an extraordinary tool for description of capturing the history of the variable and genetic characteristics of various dynamic systems. Extensive research has been carried out in the advancement of numerical and analytical solutions of linear and nonlinear FPDEs [1–6]. However, several researchers have not succeeded in deriving and modeling many complex phenomena utilizing linear or nonlinear PDEs with integer order [7]. Subsequently, the fractional is taken as account and is a good solution to this problem [8]. In the current work, three-term time-fractional Sobolev equation is considered which can be expressed as

With the conditions,

where is the Laplacian and denotes gradient operators, and , , and the are known constants, whereas , , and represent the Caputo derivative operator of order for the function .

In recent literature, various meshless methods have been utilized for the numerical solution of various PDE models almost in every discipline of science and engineering. In particular, the RBF-based meshless methods are the mainstream of these methods. The meshless nature is one of the main reasons behind the developing interest for such approaches. The meshless methods significantly reduce the complexity of dimensionality utilizing traditional methods such as the finite element and finite difference methods. Compared to mesh-based methods, these methods do not require mesh in the domain. The meshless methods have the ability to compute the solution in regular and irregular domain utilizing scattered or uniform nodes, which increases the priority and the advantages of meshless methods. As these facts show, these methods are really workable and useful numerical methods that can be applied to real-world challenging problems [9–17].

The RBF-based meshless methods have also some deficiencies like other numerical methods, in which the most important one is the dense ill-conditioned matrices and the selection of the optimal value of the shape parameter. To avoid these drawbacks, local meshless methods are the best alternatives, suggested by the researchers which are considered to be accurate and stable for the solution of diverse integer and fractional-order PDE models [18, 19]. The local meshless methods are less sensitive to the change in shape parameters than the global version, and it produces well-conditioned sparse matrices. Furthermore, local version of meshless methods is considered to be more effective and efficient than global ones. In recent years, the abilities of various sorts of local meshless methods in different applications have been explored [20–22].

In the current research, we have implemented the local meshless method to approximate the numerical solution of three-term time-fractional model equation (1). For this purpose, multiquadric (MQ) radial basis functions (RBFs) are used. Furthermore, two types of irregular domains are also taken in numerical examples.

2. Methodology of the Local Meshless Method

According to the local meshless method, to approximate the derivatives of at the centers by the neighborhood of , , where , we have used and for one-dimensional and two-dimensional cases, respectively.

Now, considering the following case for one-dimensional,

Substituting the multiquadric RBF in (3),

Equation (4) in matrix form is

where

for each . Equation (5) in simple form is

From (7), we obtain

(3) and (8) implies

where

The derivatives of w. r. t. and can be found as

For and , we continue as

The time derivative is discretized utilizing Caputo derivative [23], where as

Let be the time step size, and for the interval , consider , . We complete the time-fractional derivative term as

The term is approximated as follows:

Then,

Letting and , , we have

The fractional derivative of order and can be found as above.

3. Numerical Experiments

This section examines the accuracy and applicability of the proposed method for the three-term time-fractional model (1). In the test problems, we have considered regular and irregular domains. This computation is considered to be regular and scattered nodes with regular and irregular domains. In this article, we have used the Crank–Nicholson scheme and multiquadric (MQ) RBF with shape parameter value . Unless specifically stated, the spatial domain and time step size are used. Accuracy is measured as follows:

where is the exact solution, and is the approximate solution.

Problem 1. Consider the model equation:

Having the exact solution,

The proposed meshless method is implemented for generating the required numerical results for Problem 1, which are given in Table 1. Different values of a number of nodes , fractional order , and final time are used, whereas the error norms stand for and . These results revealed the fact that the recommended meshless method is capable of better results. Showing the accurate and efficient of the method, the results are compared with the exact solution for , , , , , and for various values of time step size . These results are computed using and are given in Table 2. One can observe from this table that only in few iterations, the suggested meshless method produced better results, and as the number of time iteration increases, the accuracy increase and the error norm reached up to . As the condition number, stability, and accuracy of the RBF-based meshless methods heavily depend on the value of shape parameter , a little change in shape parameter value causes instability and the results get diverge. But the suggested local meshless method is tested for Problem 1 in terms of condition number, stability, and accuracy as shown in Figure 1 for , , and . This figure revealed that the suggested meshless method is stable, accurate, and given ideal low condition number for a long range of up to 2000. Figure 2 shows the absolute error using and for and . Better accuracy of the recommended algorithm can be seen in this figure.

Table 1: Problem 1, approximate results for .

8.5869e − 08
3.9708e − 08
3.8841e − 07
1.2258e − 07
4.9568e − 06
1.7762e − 06

8.2120e − 08
3.9848e − 08
4.8466e − 07
1.4140e − 07
5.8961e − 06
1.9704e − 06

8.8415e − 08
3.9397e − 08
5.2865e − 07
1.6405e − 07
6.1379e − 06
2.1931e − 06

Table 2: Problem 1, approximate results using and .

0.2
1.1044e − 03
9.9648e − 04
7.9676e − 04
7.9063e − 04
7.1256e − 04
5.7016e − 04

0.02
1.0790e − 05
8.0577e − 06
3.9095e − 06
7.7071e − 06
5.6947e − 06
3.1952e − 06

0.002
1.0500e − 07
4.3655e − 08
3.8841e − 07
7.4936e − 08
3.0435e − 08
2.9006e − 07

(a)

(b)

Figure 1: Problem 1. (a) c and error norms, (b) c and condition number.

(a)

(b)

Figure 2: Problem 1. Absolute error for (a) and (b), where .

Problem 2. Consider the model equation:

Having the exact solution,

In Table 3, we have implemented the suggested algorithm for generating the numerical results for Problem 2 for , , , , , , and , . The results are assessed in term of and . Accurate results have been obtained in this problem as well. Showing the applicability and efficacy of the propose method, the results are compared with the exact solution for various values of and using . These results are given in Table 4. One can observe from this table that only in few iterations, the suggested meshless method produced better results, and as the number of time iteration increases, the accuracy increase and the error norm reached up to .
Just like the previous problem, the suggested method has been tested for Problem 2 in terms of condition number, stability, and accuracy as shown in Figure 3 for , , and . It can easily be seen from the figure that the suggest meshless method is stable, accurate, and given ideal low condition number for a long range of up to 2000, whereas in Figure 4, we have shown a comparison of exact and approximate solutions for various values of time and brilliant match of both the solutions can be found in this figure.
One of the principle advantages of the meshless techniques over mesh-based techniques is the implementation in the irregular domain with ease. In this article, two types of challenging irregular domains are taken into account, which are displayed in Figure 5. In Table 5, we have shown the numerical results obtained by the suggested meshless method corresponding to the irregular domains for Problem 1 and Problem 2. We have considered the various value of ’s, and the results are shown in form of and . It is observed from the table that better accuracy has been achieved in both domains.

Table 3: Problem 2, approximate results for .

1.2325e − 06
2.0404e − 07
2.3271e − 06
3.3074e − 07
1.9346e − 05
2.5361e − 06

1.5895e − 06
2.1784e − 07
3.2329e − 06
4.5136e − 07
2.6251e − 05
3.2570e − 06

1.4030e − 06
2.1123e − 07
4.2754e − 06
5.8259e − 07
3.1278e − 05
3.9590e − 06

Table 4: Problem 2, approximate results for .

0.2
2.6139e − 03
2.2018e − 03
1.5260e − 03
1.8624e − 03
1.5697e − 03
1.0945e − 03

0.02
2.3792e − 05
1.2980e − 05
4.3173e − 05
1.6857e − 05
9.2597e − 06
3.3106e − 05

0.002
2.0404e − 07
3.3074e − 07
2.5361e − 06
1.4425e − 07
2.4742e − 07
1.9042e − 06

0.0002
1.5812e − 09
1.1891e − 08
1.1014e − 07
1.1152e − 09
8.7439e − 09
8.2401e − 08

(a)

(b)

Figure 3: Problem 2. (a) c and error norms. (b) c and condition number.

Figure 4: Problem 2, approximate and exact solution for indicated time.

(a)

(b)

Figure 5: Computational domain 1 (a) and domain 2 (b).

Table 5: Approximate results corresponding to the irregular domains.

Domain

Problem 1
Domain 1
2.8188e − 07
7.0294e − 08
3.2180e − 06
8.2127e − 07
3.3601e − 05
9.0154e − 06

Domain 2
1.9344e − 06
3.4442e − 07
2.0244e − 05
3.8603e − 06
1.1285e − 04
2.6506e − 05

Problem 2
Domain 1
9.0788e − 06
1.1575e − 06
9.5402e − 05
1.0099e − 05
1.0994e − 03
1.1639e − 04

Domain 2
3.3707e − 06
6.0517e − 07
3.2748e − 05
6.3853e − 06
1.2824e − 04
3.2931e − 05

Problem 3. Consider the model equation:

Having the exact solution,

In Figure 6, we have visualized the behavior of the exact and approximate solutions for the Problem 3 using , , and , which show that the approximate solution is very compatible with the exact solution. In Figure 7, the absolute error is displayed for Problem 3.

(a)

(b)

Figure 6: Problem 3. (a) Exact solution. (b) Numerical solution.

(a)

(b)

Figure 7: Problem 3, absolute error at (a) and (b) .

4. Conclusion

In this study, our principle focused on the applicability and performance of the RBF-based local meshless method to approximate the numerical solution of three-term time-fractional Sobolev equations. The computed results show that the proposed technique can take care of these sorts of problems amazingly and accurately. The local procedure leads to a sparse system of linear equations, and the solution is approximated with good accuracy. Three test problems are taken into account to test the effectiveness and accuracy of the proposed meshless method utilizing rectangular and two irregular domains. The numerical results demonstrate the high accuracy and effectiveness of the method. Given the current research, the proposed technique is a surprisingly powerful and successful tool for solving numerical problems of multiterm time-fractional PDEs found in various fields of science and technology.

Data Availability

The data that support the findings of this study are openly available at https://hindawi.com/publish-research.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research of King Saud University for funding this work through research group no. RG-1441-327.

References

I. Ahmad, H. Ahmad, P. Thounthong, Y.-M. Chu, and C. Cesarano, “Solution of multi-term time-fractional PDE models arising in mathematical biology and physics by local meshless method,” Symmetry, vol. 12, no. 7, p. 1195, 2020.

N. A. Shah, I. Dassios, and J. D. Chung, “Numerical investigation of time-fractional equivalent width equations that describe hydromagnetic waves,” Symmetry, vol. 13, no. 3, p. 418, 2021.

H. Ahmad, T. A. Khan, I. Ahmad, P. S. Stanimirović, and Y.-M. Chu, “A new analyzing technique for nonlinear time fractional Cauchy reaction-diffusion model equations,” Results in Physics, vol. 19, Article ID 103462, 2020.

M. Inc, M. N. Khan, I. Ahmad, S.-W. Yao, H. Ahmad, and P. Thounthong, “Analysing time-fractional exotic options via efficient local meshless method,” Results in Physics, vol. 19, Article ID 103385, 2020.

N. A. Shah and J. D. Chung, “The analytical solution of fractional-order whitham–broer–kaup equations by an elzaki decomposition method,” Numerical Methods for Partial Differential Equations, 2021.

I. Ahmad, M. N. Khan, M. Inc, H. Ahmad, and K. S. Nisar, “Numerical simulation of simulate an anomalous solute transport model via local meshless method,” Alexandria Engineering Journal, vol. 59, no. 4, pp. 2827–2838, 2020.

S. T. R. Rizvi, I. Afzal, and K. Ali, “Chirped optical solitons for Triki-Biswas equation,” Modern Physics Letters B, vol. 33, no. 22, Article ID 1950264, 2019.

R. Attia, D. Lu, and M. M. A. Khater, “Chaos and relativistic energy-momentum of the nonlinear time fractional duffing equation,” Mathematical and Computational Applications, vol. 24, no. 1, p. 10, 2019.

M. N. Khan, S. U. Islam, I. Hussain, I. Ahmad, and H. Ahmad, “A local meshless method for the numerical solution of space-dependent inverse heat problems,” Mathematical Methods in the Applied Sciences, vol. 44, 2020.

M. Nawaz, I. Ahmad, and H. Ahmad, “A radial basis function collocation method for space-dependent inverse heat problems,” Journal of Applied and Computational Mechanics, vol. 6, 2020.

P. Thounthong, M. Khan, I. Hussain, I. Ahmad, and P. Kumam, “Symmetric radial basis function method for simulation of elliptic partial differential equations,” Mathematics, vol. 6, no. 12, p. 327, 2018.

N. A. Shah, I. Dassios, E. R. El-Zahar, J. D. Chung, and S. Taherifar, “The variational iteration transform method for solving the time-fractional fornberg-whitham equation and comparison with decomposition transform method,” Mathematics, vol. 9, no. 2, p. 141, 2021.

H. Ahmad, A. R. Seadawy, and T. A. Khan, “Study on numerical solution of dispersive water wave phenomena by using a reliable modification of variational iteration algorithm,” Mathematics and Computers in Simulation, vol. 177, 2020.

H. Ahmad, A. R. Seadawy, T. A. Khan, and P. Thounthong, “Analytic approximate solutions for some nonlinear parabolic dynamical wave equations,” Journal of Taibah University for Science, vol. 14, no. 1, pp. 346–358, 2020.

N. A. Shah, E. R. El-Zahar, M. D. Aljoufi, and J. D. Chung, “An efficient approach for solution of fractional-order helmholtz equations,” Advances in Difference Equations, vol. 2021, no. 1, pp. 1–15, 2021.

Y.-M. Chu, N. A. Shah, P. Agarwal, and J. D. Chung, “Analysis of fractional multi-dimensional Navier–Stokes equation,” Advances in Difference Equations, vol. 2021, no. 1, pp. 1–18, 2021.

H. Ahmad, T. A. Khan, P. S. Stanimirovic, and I. Ahmad, “Modified variational iteration technique for the numerical solution of fifth order KdV type equations,” Journal of Applied and Computational Mechanics, vol. 6, 2020.

I. Ahmad, M. Ahsan, I. Hussain, P. Kumam, and W. Kumam, “Numerical simulation of PDEs by local meshless differential quadrature collocation method,” Symmetry, vol. 11, no. 3, p. 394, 2019.

I. Ahmad, M. Ahsan, Z.-U. Din, A. Masood, and P. Kumam, “An efficient local formulation for time-dependent PDEs,” Mathematics, vol. 7, no. 3, p. 216, 2019.

I. Ahmad, S.-U. Islam, and A. Q. M. Khaliq, “Local RBF method for multi-dimensional partial differential equations,” Computers & Mathematics with Applications, vol. 74, no. 2, pp. 292–324, 2017.

C. Shu, Differential Quadrature and its Application in Engineering, Springer-Verlag, London, UK, 2000.

I. Ahmad, H. Ahmad, A. E. Abouelregal, P. Thounthong, and M. Abdel-Aty, “Numerical study of integer-order hyperbolic telegraph model arising in physical and related sciences,” The European Physical Journal Plus, vol. 135, no. 9, pp. 1–14, 2020.

M. Caputo, “Linear models of dissipation whose Q is almost frequency independent--II,” Geophysical Journal International, vol. 13, no. 5, pp. 529–539, 1967.