Research Article | Open Access
Nabard Habibi, Zohre Nouri, "Application of Local Fractional Homotopy Perturbation Method in Physical Problems", Advances in Mathematical Physics, vol. 2020, Article ID 2108973, 12 pages, 2020. https://doi.org/10.1155/2020/2108973
Application of Local Fractional Homotopy Perturbation Method in Physical Problems
Nonlinear phenomena have important effects on applied mathematics, physics, and issues related to engineering. Most physical phenomena are modeled according to partial differential equations. It is difficult for nonlinear models to obtain the closed form of the solution, and in many cases, only an approximation of the real solution can be obtained. The perturbation method is a wave equation solution using HPM compared with the Fourier series method, and both methods results are good agreement. The percentage of error of with , t =0.1 sec, between the present research and Yong-Ju Yang study for is less than 10. Also, the % error for in , t =0.3 sec, is less than 5, whereas for , t =0.8 and 0.7 sec, the % error for is less than 8.
The homotopy perturbation (HP) method is an analytical method that can solve linear, nonlinear, and random operator equations [1, 2]. The HP method was defined by He [3–5] as a symmetrical structure using combining HP methods in topology and was presented to be used to solve various nonlinear problems in science and engineering, including the argued nonlinear equation in fluid mechanics . He  provided a solution for the specific functional integral equations. Khan and Wu  have evolved a solution for homogenous and nonhomogeneous advection equation. Also, Xu  solved the boundary layer equation in unlimited domains. This method makes very difficult problems extremely easy, which do not need to be solved to convert nonlinear terms. The performance of nonlinear HPM has been demonstrated by many researchers. In recent years, most attention has been paid to the use of this method to solve various scientific models [6–9].
Various analytical methods are used to solve equations such as nonlinear wave equations. The variance repetition method is able to reduce the size of the calculations . This method is used in problems related to chaos and super chaos and hard equations . The HP method is a quick technique to produce convergent responses .
Recently, using HPM for local fractional solving of nonlinear problems in mathematical physics successfully achieved contains the following: diffusion equations on cantor space-time , linear and nonlinear local fractional equation of Korteweg-de Vries , fractal heat conduction problem [15, 16], wave equation in fractal strings , and the Laplace equation . Finding solutions for solving differential equations of local fractional is an interesting topic. There are some approximate and analysis methods: the local fractional functional iteration method [13–17], variationally iteration method within the Yang-Laplace transforms , Yang-Laplace transform method , local fractional Adomian decomposition and function decomposition methods , local fractional iteration of the continuously nondifferentiable functions , and Sumudu transform method .
The formulation is applied to the generalized thermoelasticity based on the fractional time derivatives under the effect of diffusion. Herzallah  criticized the research of Ref [a] and showed that their results are incorrect. New methodologies in fractional and fractal derivatives modeling in scientific and engineering were studied by Chen and Liang .
Uchaikin and Sibatov  studied fractional derivatives on cosmic scales. The numerical calculation results demonstrated in this paper do speak well for the NoRD-model as compared with the traditional one based on integer-order operators.
Chaos in a 5D hyperchaotic system with four wings in the light of nonlocal and nonsingular fractional derivatives was investigated by Bonyah . The numerical simulation results depict a new chaotic behaviors with the ABC numerical scheme.
Goufo and Toudjeu  evaluated the analysis of recent fractional evolution equations and applications. The numerical approximations of a second-order nonhomogeneous fractional Cauchy problem are performed and show regularity in the dynamics.
In the study of Turkyilmazoglu , it is theoretically shown that under a special constraint, the HPM converges to the exact solution of sought solution of nonlinear ordinary or partial differential equation. Examples (the classical Blasius flat-plate flow problem of fluid mechanics, fourth-order parabolic partial differential equation arising in the study of the transverse vibrations of a uniform flexible beam) clearly demonstrate, why and on what interval, the corresponding homotopy series generated by the HPM converges to the exact solution. The physical problems involving stronger nonlinearities can also be dealt with the homotopy method under the assurance of convergence provided by the presented theorems. This method (HPM) is a powerful device for solving a wide variety of problems. Using the HPM, it is possible to find the exact solution or an approximate solution of the problem. Some examples such as Burgers’, Schrodinger’s, and fourth-order parabolic partial differential equations are presented .
For this reason, the method of homeopathic disorder has been selected for some important physical issues that have received less attention in previous studies. A review of previous studies suggests that local fractional homotopy perturbation method is less commonly used for physical problems. The aim of this research is expansion of an HP method for the definition of local fractional derivation order that can solve a local fractional differential equation in usages of physical mathematics. Therefore, first, the HP method will be explained, and then the local damped wave equation, local fractional wave equation, and heat diffusion equation with several initial values are solved using this method. Also, the obtained results in applied examples are compared with the Fourier series method . In the present research, the homotopy perturbation method (HPM) was applied for real physical problems in engineering which has been less studied in previous studies.
2. Mathematical Fundamentals
Definition 1. Suppose that the function , local fractional derivative of of order α, , at the point is defined as follows [28–41]:where and are satisfied with the condition  , so that [28, 36] with for and
However, the partial of local fractional is defined as follows [13–18]:If there exists for any , then we denote where .
Definition 2. The diffusion equation on the Cantor sets (called local fractional diffusion equation) was recently described  aswhere denotes the fractional diffusion constant which is, in essence, a measure of the efficiency of the spreading of the underlying substance, while the local fractional wave equation is written in the following form [30, 31]:The local fractional Laplace operator is given [27, 42] as follows:We notice that the local fractional diffusion equation yieldsAnd the local fractional wave equation has the following form:where is a constant. This equation describes vibrations in a fractal medium.
The Helmholtz equation with local fractional derivative operators in the two-dimensional case was suggested in [43, 44] as follows:with the initial value conditions as follows:where is the unknown function and is a source term.
3. Analysis of HP Method
To show the ability of local fractional HP, the local fractional differential equation is written as follows:where is a local fractional differential operator. Convex nondifferentiable homotopy is defined as follows:
In nondistinguishable homotopy, these relations are called nondistinguishable deformation, and , are called nondistinguishable homotopy. By applying nondistinguishable perturbation method, relationship (15) is written as follows:
In the other words, for the calculation of the approximated solution of relation (12), relationship is expanded as a local fractional Taylor series. Therefore,
Using relation (21) and simplification, we can write
It can be written according to relation (24):
Therefore, by considering as the first approximation, it results to ; consequently,
This equation is a local fractional iteration Newton famous formula [3–11], and this relation is converged. Using relationship (28), we can calculate the relation of a local fractional iteration semi-Newton with second-order approximation parameters as follows:
When , then the approximation solution is as follows:
4. Illustrative Applied Examples
In this section, local damped wave equation, wave equation in the fractal strings, and heat diffusion equation in the case of local fractional using HPM were solved.
Example 1. The wave equation in the fractal strings.
At first, the wave equation in the fractal strings is also intended that is given the initial conditions .considering the boundary conditions:According to the local fractional HPM, the unrecognizable homotopy structure can be written as follows:On the other hand, according to equation (13), we can writeBy substituting (17) in (16), we can sayNow, for separating the sentences of relationship (18), just this relationship can be expanded and separately written to various powers ; therefore, linear fractional derivations of differential equations are adjusted as follows:Finally,:Evaluation of equations (36) to (41) for , some components of the solution of local fractional HP are as follows:Finally,Therefore, almost unrecognizable for relation (14) will be as follows:According to equation (15), is considered as boundary condition for computations as follows:With these conditions using the computer programming in the MATLAB software, can be plotted as shown in Figure 1. For Figures 1 and 2, and , respectively.
Since most references in their calculations equal to are considered, therefore, in this study, for better understanding of local fractional wave equation, is considered and for three different time versus plotted in Figure 3, and it is remarkable.
Example 2. Heat diffusion equation in the case of local fractional.
The next example, heat diffusion equation in the case of local fractional, using the HPM, is solved as follows:In this case, the initial condition is considered as follows:Now, HPM can be used, and the above equation can be written as follows:According to relationship and substitution in equation (48), consequentlyTherefore, linear fractional derivations of differential equations are as follows:and so on. By solving equations (50) to (54), it becomes clear thatWhen , then ; therefore,The equation can be rewritten asConsidering , is plotted in Figure 4.
Now, equation (48) is considered as follows:that the initial condition of the form is considered; applying the above operation is determined answering the following:By placing , is shown in Figure 5.
Example 3. Local damped wave equation.
The second example, local damped wave equation (60), using the HPM solved .In this case, the initial condition is considered as follows:Now, HPM can be used, and the above equation can be written as follows:According to relationship and substitution in equation (62), it results to the following:Therefore, linear fractional derivations of differential equations are as follows:By solving the equations (64) to (68), it becomes clear thatFinally,When , then ; therefore,By placing , versus is illustrated in Figure 6. Also, Figure 7 shows that for three different times, the variation of versus is considerable.
5. Comparison of Results of HP and Fourier Series Methods
In this section, for the validation and benchmarking of the present method, equation (60) means that the local damped wave equation by this method and the results of the Fourier series method for this equation are compared and discussed, according to the results for the following equation with initial condition:Using HPM, the solution is determined as follows:
In reference paper , the local fractional wave equation is solved:where and .
In this section, for several times and , the variations of are obtained using the Fourier series and HP methods. In these graphs, the show that the HP method is considered in the present research and illustrate the obtained results of the Fourier series method. Comparison of two graphs of for time equal to 0.1, 0.3, and 0.6 sec for two methods with is illustrated in Figures 8–13.
According to Figures 8–10 () and Figures 11–13 (), it is observed that for , the value of error for the results is very less, and by increasing time for , the values of convergence increased. In wave equation, by increasing time, due to damping, the wave domain decreased. Therefore, the best of case is equal to and .
Variations of versus with , are illustrated in Tables 1 and 2. According to Table 1, the percentage of error of between the present study and Ref  with , , and , is less than 10. Also, the % error for , in , , is less than 5, whereas for , , the % error for is less than 10. This comparison shows that by increasing from 0.5 to 1 for all of the three times, the % error between the two methods of homotopy perturbation and the Fourier series decreases sharply and converges to 1.