Numerical Solutions with Linearization Techniques of the Fractional Harry Dym Equation

Open access

Abstract

In this study, numerical solutions of the fractional Harry Dym equation are investigated. Linearization techniques are utilized for non-linear terms existing in the fractional Harry Dym equation. The error norms L2 and L are computed. Stability of the finite difference method is studied with the aid of Von Neumann stabity analysis.

Abstract

In this study, numerical solutions of the fractional Harry Dym equation are investigated. Linearization techniques are utilized for non-linear terms existing in the fractional Harry Dym equation. The error norms L2 and L are computed. Stability of the finite difference method is studied with the aid of Von Neumann stabity analysis.

1 Introduction

Fractional analysis; is the generalization of the classic analysis of integration and differentiation of process (noninteger) order. This issue is an old issue as much as differential calculus. G. W. Leibniz and Marquis de L' Hospital correspondence in 1695, is known as the first exit point of fractional calculus. Leibniz expressed fractional order derivatives of noninteger order f (x) = emx, m ∈ R, as follows:

dnemxdxn=mnemx,

here, n is a value, where i is the noninteger.

Later, many scientists, such as Liouville, Riemann,Weyl, Lacroix, Leibniz, Grunward and Letnikov (cf. [4]), expanded the range of this derivative.

Since the beginning of the definition of fractional order derivatives first handled by Leibniz, fractional partial differential equations have attracted the attention of many scientists and have also shown a progressive development.(cf. [1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 13, 14, 15, 16, 17, 18, 19]).

Beside these facts, the third order fractional Harry Dym partial differential equation is studied in mathematics and especially in the theory of solitons. This equation is given as

γutγ=u33ux3.

Here, 0 < γ ≤ 1 is the order of fractional derivative and u(x, t) is a function of x and t.

The Harry Dym equation first appeared in a study by Kruskal [10]. Harry Dym equation represents a system in which dispersion and nonlinearity were coupled. Furthermore, the Harry Dym equation is a completely integrable nonlinear evolution that may be solved by means of the inverse scattering transform. It does not possess the Painlevé property.

This paper is organized as follows:

In the second section, some basic facts dealing with the finite difference method are mentioned and three Linearization techniques are presented. In the third section, stability analysis of the proposed method is investigated and it is shown that the Harry Dym equation is stable under which conditions. Also, numerical examples are given. In the fourth section, conclusions obtained throughout the paper are discussed.

2 Finite Difference Methods

In this section, we first need to define a set of grid points in a domain D to obtain a numerical solution to Eq. (1:1) using finite difference methods as follows:

Let Δx(h)=baN(N is an integer) denotes a state step size and Δt denotes a time step size. Draw a set of horizontal and vertical line across D, and get all intersection points (xj, tn) or simply ( j,n) where x j = a+ jΔx, j = 0,1,2, · · · ,N, and tn = nΔt, n = 0,1, . . . ,M. If we write D = [a,b]×[0,T] then we may choose Δt=TM(Mis an integer) and tn = nΔt, n = 0,1, . . . ,M.

Then, an appropriate finite difference approximation is given in Eq. (1.1) instead of derivatives and its variable. In this case, the solution problem of Eq. (1.1) is reduced to solution problem of algebraic differential systems of linear and nonlinear equations consisting of finite difference equation. But when applied to non-linear problems, it normally leads to nonlinear system of equations and they cannot be solved directly. Therefore, we use three linearization techniques for a nonlinear term as given in Eq. (1.1).

2.1 Linearization 1

First, we use the Caputo fractional derivative approximation for γutγdefined by

γutγ1h2ΔtγΓ2γk=0numn+1kumnkk+11γk1γn1ΔtγΓ2γum1um0n=0,

and Crank–Nicolson derivative approximation given by

u=un+1+un2,

in Eq. (1.1) at the nodal point (m,n+1) [12]. Then, if we discretize time derivative of the fractional Harry Dym equation by using Caputo fractional derivative formula between two successive time levels n and n +1, Crank Nicolson derivative formula and usual finite difference formula between two successive time levels n and n +1,respectively, we obtain

(Δt)γΓ(2γ)k=0n[umn+1kumnk][(k+1)1γk1γ]=[(u3)mn+1+(u3)mn2][um+2n2um+1n+2um1num2n2h3].

The nonlinear term in the above mentioned equation is linearized by using the following equation:

(u3)mn+1=3(u2)nun+12(u3)n.

If the necessary arrangements are made in Eq. (2.1), we have the following equation:

(Δt)γΓ(2γ)k=0n[umn+1kumnk][(k+1)1γk1γ]=14h3[3(u2)mnumn+1(u3)n][um+2n2um+1n+2um1num2n].

2.1.1 Linearization 2

Let us use the Caputo fractional derivative approximation for γutγ.,Crank–Nicolson derivative approximation and usual finite difference approximation for Uxxx given by

Uxxx12h3(um+2n2um+1n+2um1num2n),

in Eq. (1.1) at the nodal point (m,n+1), respectively [12]:

If we use the following linearization technique for the non-linear term U 3Uxxx:

U3Uxxx(Umn)312h3(um+2n2um+1n+2um1num2n),

then we have the following system of algebraic equation:

(Δt)γΓ(2γ)k=0n[umn+1kumnk][(k+1)1γk1γ]=(Umn)312h3(um+2n+12um+1n+1+2um1n+1um2n+1).

Eq. (2.3) can be solved using an approximate algorithm.

2.1.2 Linearization 3

Let us use the Caputo fractional derivative approximation for γutγand usual finite difference approximation for Uxxx in Eq. (1.1) at the nodal point (m,n+1) respectively [12]: If we use the following linearization technique for the nonlinear term U 3Uxxx, then we have

U3UxxxUmn+Um+1n2312h3um+2n2um+1n+2um1num2n.

Thus, we get the following system of algebraic equation:.

(Δt)γΓ(2γ)k=0n1[umnkumnk1][(k+1)1γk1γ]=((Umn+Um+1n)2)312h3(um+2n2um+1n+2um1num2n).

Eq. (2.4) can be solved using an approximate algorithm.

3 Stability analysis

In this section, we investigate whether this method is stable based on von-Neumann analysis. If the Fourier method analyzes the stability, then the growth factor of a typical Fourier mode is defined as:

umn=ξnemφ,=1,

where ξ n is considered as the amplification factor. First, by substituting the Fourier mode (3.1) into the recurrence relationship (2.3), one can obtain

4h3(Δt)γΓ(2γ)k=0n1[(k+1)1γk1γ][ξnk+1ξnk]emφ=(ξnemφ)3(ξn+1e(m+2)φ2ξn+1e(m+1)φ+2ξn+1e(m1)φξn+1e(m2)φ).

Next, if we assume that ξ n+1 = ζ ξ n and ζ = ζ (φ) are independent of time, we can easily obtain the following expression:

4h3ΔtγΓ2γk=0n1k+11γk1γζk3n+1ξk3n=ξe3mφe2φ2eφ+2eφe2φ.

Hence, we get

ξ=X1+iX2Y1+iY2,
X1=2(Δt)γΓ(2γ)k=0n1[(k+1)1γk1γ][ζ13nkζ3nk],X2=0,Y1=24cos2[mγ]sin[mφ]sin2[φ2]sin[φ]+8sin3[mφ]sin2[φ2]sin[φ],Y2=24cos[mγ]sin2[mφ]sin2[φ2]sin[φ]+8cos3[mφ]sin2[φ2]sin[φ],

which shows that

|ζ|=|X1+X2Y1+Y2|.

For the Fourier stability definition and for the examined scheme to be stable, the condition |ξ | ≤ 1 must be gratified. Therefore, if the following inequality is provided, the schema is unconditionally stable.

X12Y12+Y221.

Other schemes can be studied by following a similar way.

3.1 L2 and L error norms

The equation of the numerical results was obtained for the test problem used in this study and all computations have been run on using double precision arithmetic. To show how accurate the results, both the error norm L2

L2=||UexactUN||2=hJ=0N|Ujexact(UN)j|2,

and L

L=||UexactUN||=maxj|Ujexact(UN)j|,

are going to be computed and presented.

3.2 Test Problem

The analytical solution of the fractional Harry Dym equation is given as follows [11].

u(x,t)=(432(x+t))23,tt0,0x1.

In our computations, for the numerical solution of the test problem three different linearization techniques have been applied. The values of the error norms L2 and L have been computed at t = 1 for different values of Δt, h. The comparison of the error norms L2 and L obtained by the linearization techniques is summarized in Table 1. As summarized in Table 1, it is obvious that the obtained results using linearization 1 are better than the obtained results using other linearizations. As shown in Figure 1, one can compare the exact and approximate solution of the collocation method using four radial basis functions for h = 0.01. Similarly, as shown in Figures 2 and 3 the exact and approximate solutions can be compared for linearization 2 and linearization 3 successively.

Fig. 1

Download Figure

Fig. 1

Comparison of the exact and approximate solution for linearization 2 at h = 0.01

Citation: Applied Mathematics and Nonlinear Sciences 4, 1; 10.2478/AMNS.2019.1.00004

Fig. 2

Download Figure

Fig. 2

Comparison of the exact and approximate solution for linearization 3 at h=0.01

Citation: Applied Mathematics and Nonlinear Sciences 4, 1; 10.2478/AMNS.2019.1.00004

Fig. 3

Download Figure

Fig. 3

Comparison of the exact and approximate solution for linearization 1, linearization 2, and linearization 3 at h = 0.01

Citation: Applied Mathematics and Nonlinear Sciences 4, 1; 10.2478/AMNS.2019.1.00004

Table 1

Comparison of the error norms L2 and L that are obtained using the linearization techniques at t = 1 = 0.9 for different values of Δt, h.


Lin. I
Lin. II
Lin. III
L2×102L×102L2×102L×102L2×102L×102
Δt = h = 0.0010.000020.000020.005660.056640.005650.05659
Δt = h = 0.050.052540.083440.645560.941720.585160.84499
Δt = h = 0.040.029570.051280.506900.821060.472240.76100
Δt = h = 0.030.014180.027730.368000.683910.351130.65086
Δt = h = 0.020.005060.011860.231240.523130.225020.50855
Δt = h = 0.010.000880.002850.101610.323180.100460.31945
Δt = h = 0.22.445680.002853.127563.788420.807561.29370
Δt = h = 0.13.276200.416451.312111.412870.940270.94563

4 Conclusions

Finite difference methods based on using three different linearization techniques have been proposed for the numerical solutions of the fractional Harry Dym equation. In addition, numerical results were obtained by using three different linearization techniques. The proposed method has been tested on a problem and demonstrated how effective it is. The error norms L2 and L have been calculated and given. The third linearization technique, as shown in Figure 3, yielded better results. Because the third linearization technique gives better results, this technique can be suggested in the next problems and the situation of the problem should be considered. The obtained results show that the error norms are sufficiently small during all computer runs. It has been observed that the considered method is a power numerical scheme to solve the fractional Harry Dym equation.

Communicated by Juan Luis García Guirao

References

  • [1]

    Al-Khaled, K. and Alquran, M. (2014). An approximate solution for a fractional model of generalized Harry Dym equation. Mathematical Sciences, 8(4), 125-130.

  • [2]

    Biswas, A., Mirzazadeh, M., Eslami, M., Milovic, D., & Belic, M. (2014). Solitons in optical metamaterials by functional variable method and first integral approach. Frequenz, 68(11-12), 525-530.

  • [3]

    Baskonus, H., Mekkaoui, T., Hammouch, Z., and Bulut, H. (2015). Active control of a chaotic fractional order economic system. Entropy, 17(8), 5771-5783.

  • [4]

    Das, S. (2011). Functional fractional calculus. Springer Science & Business Media.

  • [5]

    Esen, A., Ucar, Y., Yagmurlu, N. and Tasbozan, O. (2013). A Galerkin finite element method to solve fractional diffusion and fractional diffusion-wave equations. Mathematical Modelling and Analysis, 18(2), 260-273.

  • [6]

    Ravichandran, C., Jothimani, K., Baskonus, H. M., and Valliammal, N. (2018). New results on nondensely characterized integrodifferential equations with fractional order. The European Physical Journal Plus, 133(3), 109.

  • [7]

    Kaya, D., Gülbahar, S., Yokuş, A. and Gülbahar, M. (2018). Solutions of the fractional combined KdV–mKdV equation with collocation method using radial basis function and their geometrical obstructions. Advances in Difference Equations, 2018(1), 77.

  • [8]

    Kumar, D., Singh, J. and Kılıçman, A. (2013). An efficient approach for fractional Harry Dym Equation by Using Sumudu Transform, Abstr. Appl. Anal., 2013, Article ID 608943, 8.

  • [9]

    Kumar, S., Tripathi, M.P. and Singh, O.P. (2013). A fractional model of Harry Dym equation and its approximate solution, Ain Shams Eng. J., 4 , 111-115.

  • [10]

    Kruskal, M. (1975). Nonlinear Wave Equations. In Jürgen Moser, editor, Dynamical Systems, Theory and Applications,Heidelberg. Springer. 38, 310–354.

  • [11]

    Rawashdeh, M., S. (2014). New approach to solve the fractional Harry Dym equation using the FRDTM, Int. J. Pure Appl. Math., 95, 553-566.

  • [12]

    Kaya, D., Yokuş, A. ( 2014). Stability Analysis and Numerical Solutions for Time Fractional KdVB Equation, International Comference on Computational Experimental Science and Engineering, Antalya.

  • [13]

    Esen, A., Sulaiman, T. A., Bulut, H., and Baskonus, H. M. (2018). Optical solitons to the space-time fractional (1+ 1)-dimensional coupled nonlinear Schrdinger equation. Optik, 167, 150-156.

  • [14]

    Yang, X. J., Gao, F., Machado, J. A., and Baleanu, D. (2017). A new fractional derivative involving the normalized sinc function without singular kernel. arXiv preprint arXiv:1701.05590.

  • [15]

    Yavuz, M., Ozdemir, N. and Baskonus, H. M. (2018). Solutions of partial differential equations using the fractional operator involving Mittag-Leffler kernel. The European Physical Journal Plus, 133(6), 215.

  • [16]

    Yavuz, M. (2019). Characterizations of two different fractional operators without singular kernel. Mathematical Modelling of Natural Phenomena, 14(3), 302.

  • [17]

    Bulut, H., Kumar, D., Singh, J., Swroop, R., and Baskonus, H. M. (2018). Analytic study for a fractional model of HIV infection of CD4+ TCD4+ T lymphocyte cells. Math. Nat. Sci, 2(1), 33-43.

  • [18]

    Baskonus, H. M., Hammouch, Z., Mekkaoui, T., and Bulut, H. (2016, June). Chaos in the fractional order logistic delay system: Circuit realization and synchronization. In AIP Conference Proceedings (Vol. 1738, No. 1, p. 290005). AIP Publishing.

  • [19]

    Bulut, H., Yel, G., Bakonu, H. M. (2016). An application of improved Bernoulli sub-equation function method to the nonlinear time-fractional burgers equation. Turkish Journal of Mathematics and Computer Science, 5, 1-7.

[1]

Al-Khaled, K. and Alquran, M. (2014). An approximate solution for a fractional model of generalized Harry Dym equation. Mathematical Sciences, 8(4), 125-130.

[2]

Biswas, A., Mirzazadeh, M., Eslami, M., Milovic, D., & Belic, M. (2014). Solitons in optical metamaterials by functional variable method and first integral approach. Frequenz, 68(11-12), 525-530.

[3]

Baskonus, H., Mekkaoui, T., Hammouch, Z., and Bulut, H. (2015). Active control of a chaotic fractional order economic system. Entropy, 17(8), 5771-5783.

[4]

Das, S. (2011). Functional fractional calculus. Springer Science & Business Media.

[5]

Esen, A., Ucar, Y., Yagmurlu, N. and Tasbozan, O. (2013). A Galerkin finite element method to solve fractional diffusion and fractional diffusion-wave equations. Mathematical Modelling and Analysis, 18(2), 260-273.

[6]

Ravichandran, C., Jothimani, K., Baskonus, H. M., and Valliammal, N. (2018). New results on nondensely characterized integrodifferential equations with fractional order. The European Physical Journal Plus, 133(3), 109.

[7]

Kaya, D., Gülbahar, S., Yokuş, A. and Gülbahar, M. (2018). Solutions of the fractional combined KdV–mKdV equation with collocation method using radial basis function and their geometrical obstructions. Advances in Difference Equations, 2018(1), 77.

[8]

Kumar, D., Singh, J. and Kılıçman, A. (2013). An efficient approach for fractional Harry Dym Equation by Using Sumudu Transform, Abstr. Appl. Anal., 2013, Article ID 608943, 8.

[9]

Kumar, S., Tripathi, M.P. and Singh, O.P. (2013). A fractional model of Harry Dym equation and its approximate solution, Ain Shams Eng. J., 4 , 111-115.

[10]

Kruskal, M. (1975). Nonlinear Wave Equations. In Jürgen Moser, editor, Dynamical Systems, Theory and Applications,Heidelberg. Springer. 38, 310–354.

[11]

Rawashdeh, M., S. (2014). New approach to solve the fractional Harry Dym equation using the FRDTM, Int. J. Pure Appl. Math., 95, 553-566.

[12]

Kaya, D., Yokuş, A. ( 2014). Stability Analysis and Numerical Solutions for Time Fractional KdVB Equation, International Comference on Computational Experimental Science and Engineering, Antalya.

[13]

Esen, A., Sulaiman, T. A., Bulut, H., and Baskonus, H. M. (2018). Optical solitons to the space-time fractional (1+ 1)-dimensional coupled nonlinear Schrdinger equation. Optik, 167, 150-156.

[14]

Yang, X. J., Gao, F., Machado, J. A., and Baleanu, D. (2017). A new fractional derivative involving the normalized sinc function without singular kernel. arXiv preprint arXiv:1701.05590.

[15]

Yavuz, M., Ozdemir, N. and Baskonus, H. M. (2018). Solutions of partial differential equations using the fractional operator involving Mittag-Leffler kernel. The European Physical Journal Plus, 133(6), 215.

[16]

Yavuz, M. (2019). Characterizations of two different fractional operators without singular kernel. Mathematical Modelling of Natural Phenomena, 14(3), 302.

[17]

Bulut, H., Kumar, D., Singh, J., Swroop, R., and Baskonus, H. M. (2018). Analytic study for a fractional model of HIV infection of CD4+ TCD4+ T lymphocyte cells. Math. Nat. Sci, 2(1), 33-43.

[18]

Baskonus, H. M., Hammouch, Z., Mekkaoui, T., and Bulut, H. (2016, June). Chaos in the fractional order logistic delay system: Circuit realization and synchronization. In AIP Conference Proceedings (Vol. 1738, No. 1, p. 290005). AIP Publishing.

[19]

Bulut, H., Yel, G., Bakonu, H. M. (2016). An application of improved Bernoulli sub-equation function method to the nonlinear time-fractional burgers equation. Turkish Journal of Mathematics and Computer Science, 5, 1-7.

Journal Information

Figures

  • View in gallery

    Comparison of the exact and approximate solution for linearization 2 at h = 0.01

  • View in gallery

    Comparison of the exact and approximate solution for linearization 3 at h=0.01

  • View in gallery

    Comparison of the exact and approximate solution for linearization 1, linearization 2, and linearization 3 at h = 0.01

Metrics

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 163 163 67
PDF Downloads 92 92 26