Abstract
In this paper, the two-dimensional Bernoulli wavelets (BWs) with Ritz–Galerkin method are applied for the numerical solution of the time fractional diffusion-wave equation. In this way, a satisfier function which satisfies all the initial and boundary conditions is derived. The two-dimensional BWs and Ritz–Galerkin method with satisfier function are used to transform the problem under consideration into a linear system of algebraic equations. The proposed scheme is applied for numerical solution of some examples. It has high accuracy in computation that leads to obtaining the exact solutions in some cases.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
Many phenomena in various field of the science, can be modeled very successfully by time-fractional differential equations. In this paper we focus on the following fractional diffusion-wave equation (FDWE) with damping [1]:
with the initial conditions:
and the boundary conditions:
where \(L>0, T>0, 1<q \le 2\) is the order of the fractional derivative in the Caputo sense, \(f_0, f_1, g_0\) and \(g_1\) are known and sufficiently smooth functions, while the function u is to be determined. In the case \(q=2\), this equation is named telegraph equation.
Recently, considerable amount of papers have been proposed methods for solving the FDWE [2–13]. Chen et al. [1] obtained the analytical solution by the method of separation of variables and proposed the numerical solution with finite difference method. In [2], Bhrawy et al. applied a spectral tau method based on the Jacobi operational matrix to solve the problem. Liu et al. [3] proposed the fractional predictor–corrector method to solve this problem. In [4], Mainardi derived the fundamental solutions for the FDWE. The combination of the compact difference method and alternating direction implicit method are used for solving two-dimensional fractional Cattaneo equation in [5]. A fully discrete difference scheme is recommended for a diffusion-wave system by Wess [6]. Heydari et al. [7] applied fractional operational matrix (FOM) of integration for the Legendre wavelets (LWs) to solve the problem. In [8] a compact finite-difference scheme is used for the fourth-order fractional diffusion-wave system. In [9] finite difference schemes of second-order are proposed for the time-fractional diffusion-wave equation. In [10], Hu and Zhang used finite-difference methods for fourth-order fractional diffusion-wave. Sumudu transform method for solving fractional differential equations and fractional diffusion-wave equation applied by Darzi et al. [11]. Hosseini et al. [12] employed the meshless local radial point interpolation method which is based on the Galerkin weak form and radial point interpolation approximation for solving FDWE.
The Ritz–Galerkin method is the method to transform a continuous problem to a discrete problem. Several partial differential equations are numerically solved by Ritz–Galerkin method, but using of the appropriate satisfier function in the Ritz–Galerkin method is taken into consideration recently, see for instance [14–20]. The satisfier function fulfills all the problem conditions. In conclusion, employing of it in Ritz–Galerkin method provides the facility to satisfy the problem conditions, also leads to a system of algebraic equations of smaller size and hence reduces the computation time.
In mathematical research, wavelet theory is a relatively new and growing area. It has been used in a wide range of engineering; for instance, wavelets are very successfully applied in signal analysis for waveform representations and segmentations [21]. Wavelets allow the accurate representation of many types of functions and operators [22, 23]. Furthermore, wavelets make a connection with fast numerical algorithms [24].
The shifted Legendre polynomials \(p_{n}(x), n=0,1,2,...\), that \(0 \le x \le 1\), are more efficient in approximation theory, among other orthogonal polynomials [25, 26]. The Bernoulli polynomials are not orthogonal functions. In [27], the superiority of Bernoulli polynomials \(\beta _n(x),\, n=0,1,2,...,\) so that \(0 \le x \le 1,\) to shifted Legendre polynomials in approximation of functions is proposed.
In this paper, we define two-dimensional BWs for the first time. Moreover, this is the first time the Ritz–Galerkin method in the two-dimensional BWs basis and with utilizing the satisfier function is employed to give an approximate solution of FDWE. We also compare our results with those results obtained by [3] and [7]. Comparison for the numerical examples shows the more accuracy and less computations of our scheme in comparison to other published methods.
This paper is separated in to the following sections: In Sect. 2, we introduce basic formulation of wavelets and the Bernoulli wavelets. In Sect. 3, we construct two different satisfier functions and apply the Ritz–Galerkin method in the two-dimensional BWs basis for numerically discretize the problem. Section 4 presents and discusses the numerical results for two test examples, whilst Sect. 5 includes the conclusions of this paper.
Properties of Bernoulli wavelets
Wavelets and Bernoulli wavelets
Dilation and translation of a function (mother wavelet) construct a family of functions called wavelets and is defined as follows [28]
where \(a,b \in {\mathbb {R}}\) are dilation and translation parameters that vary continuously. The family of discrete wavelets that form a basis for \(L^{2}({\mathbb {R}})\) is defined as follows
where n and k are positive integers and \(a_{0}>1, b_{0}>1\).
Bernoulli wavelets are obtained, when we choose Bernoulli polynomial \(\beta _{m}(t)\) as mother wavelet. BWs \(\varphi _{n,m}(t)=\varphi (k,{\hat{n}},m,t)\) have four arguments; \({\hat{n}}=n-1, n=1,2,4,...,2^{k-1}, m\) is the order for Bernoulli polynomials and k can assume any positive integer. They are defined on the interval [0, 1) for \(m=1,2,...,M-1\), that \(M>0\) is a fixed integer by [29]
and
where \(n=1,2,...,2^{k-1}\). The coefficient \(\frac{1}{\sqrt{\frac{(-1)^{m-1}(m!)^{2}}{(2m)!}\alpha _{2m}}}\) is applied to normalize the Bernoulli wavelets. Bernoulli polynomials form a complete basis over the interval [0, 1] [30] and are defined by [31]
that \(\alpha _i,\,\,i=0,1,...,m\) are Bernoulli numbers.
For Bernoulli polynomials we have [32]
Now, let
We define, for the first time, the two-dimensional Bernoulli wavelets \(\varphi _{n_1m_1n_2m_2}(x,y)\) as
where \(m_1=0,1,2,...,M_1-1, m_2=0,1,2,...,M_2-1\). Here, \({\hat{n}}_1\) and \({\hat{n}}_2\) are defined similarly to \({\hat{n}}, k_1\) and \(k_2\) can be any positive integer, \({\tilde{\beta }}_{m_1}\) and \({\tilde{\beta }}_{m_2}\) are defined similarly to \({\tilde{\beta }}_{m}\) of order \(m_1\) and \(m_2\), respectively.
Function approximation
A function \(f(x,y)\in L^2([0,1)\times [0,1))\) may be expanded as in terms of two-dimensional Bernoulli wavelets as
If the infinite series for f is truncated, then Eq. (2.4) can be written as
where \(\Phi _1(x)\) and \(\Phi _2(y)\) are \(2^{k_1-1}M_1\times 1\) and \(2^{k_2-1}M_2\times 1\) matrices, respectively, given by
In Eq. (2.5), F is a \(2^{k_1-1}M_1\times 2^{k_2-1}M_2\) matrix that can be calculated from [33]
where \(\langle . \rangle\) denotes the inner product, \(D_1=\langle \Phi _1(x),\Phi _1(x) \rangle =\int _{0}^{1}\Phi _1(x)\Phi ^T_1(x)\mathrm{d}x\) and \(D_2=\langle \Phi _2(y),\Phi _2(y)\rangle =\int _{0}^{1}\Phi _2(y)\Phi ^T_2(y)\mathrm{d}y\) are \(2^{k_1-1}M_1\times 2^{k_1-1}M_1\) and \(2^{k_2-1}M_2\times 2^{k_2-1}M_2\) matrices, respectively.
Satisfier function
In the Ritz–Galerkin method with the two-dimensional BWs basis, the approximation \({\widetilde{u}}(x,t)\) of the solution u(x, t) in (1.1) is sought in the form of the truncated series
where \(\sigma _{nilj}(x,t)=x(x-L)t^2\varphi _{ni}(x)\varphi _{lj}(t)\) and w(x, t) is a satisfier function. The most important point in using the Ritz–Galerkin method is finding the satisfier function, which satisfies all the problem conditions [16]. Interpolation is one of the methods that is usually used to derive Satisfier functions. On the other hand, we have seen from experience, when in constructing of the satisfier function we use only the problem’s data, we get a satisfier function that is closer to the exact solution. Therefore, we obtain the cost-effective computational results [14, 16]. Here, we construct two different satisfier functions for the initial and boundary conditions:
It is worth pointing out that \(f_0(x), f_1(x), g_0(t)\) and \(g_1(t)\) satisfy the following compatibility conditions:
The first technique for obtaining the satisfier function is as follows:
We set
then we construct \(k_{1}(x)\) and \(k_{2}(x)\), such that (3.8) fulfils the conditions (3.2)–(3.5).
Clearly, if \(k_{1}(x)\) and \(k_{2}(x)\) satisfy the following conditions:
then (3.8) can be a satisfier function for (3.2)–(3.5).
Equations (3.10) form a system of linear equations, which can be solved for \(k_{1}(x)\) and \(k_{2}(x)\) when \(g_0(0)\acute{g_1}(0)-g_1(0)\acute{g_0}(0)\ne 0\). By solving this system, we obtain
From compatibility conditions (3.6) and (3.7), it is easy to see that \(k_{1}(x)\) and \(k_{2}(x)\) have properties (3.9). Therefore, we introduce the satisfier function w(x, t) which satisfies the initial conditions (3.2) and (3.3), and the boundary conditions (3.4) and (3.5) when
as:
where \(k_{1}(x)\) and \(k_{2}(x)\) are obtained from (3.11) and (3.12).
The second satisfier function is constructed as follows:
We firstly transform the nonhomogeneous boundary condition into a homogeneous boundary condition. Let
where
The function v(x, t) then satisfies the problem with homogeneous boundary conditions:
that \(F_0(x)=f_0(x)-\phi (x,0),\) and \(F_1(x)=f_1(x)-\phi _{t}(x,0).\) From conditions (3.14)–(3.17), we drive compatibility conditions:
Therefore, \(W(x,t)=F_0(x)+t F_1(x)\) satisfies conditions (3.14)–(3.17) and eventually, we introduce the satisfier function for (3.2)–(3.5) as:
It is worth to mention that if \(g_0(0)\acute{g_1}(0)-g_1(0)\acute{g_0}(0)\ne 0\), we prefer the first attained satisfier function in (3.13), since in constructing of (3.13) we used only the problem’s data.
Returning now to the Ritz–Galerkin approximation (3.1), the expansion coefficients \(c_{nilj}\) are determined by the Galerkin equations:
where
and
where \(\varphi _{ni}(x),\varphi _{lj}(t)\) are BWs. Equations (3.19) form a linear system of equations which can be solved for the elements of \(c_{nilj}, n=1,...,2^{k_1-1}, i=0,...,M_1-1, l=1,...,2^{k_2-1}, j=0,...,M_2-1\) using mathematical softwares.
Numerical results and comparisons
In this section, we apply the numerical scheme in the previous section for finding the approximation solutions of two examples of FDWE. We compare our results with obtained results in [3] and [7]. In all examples the package of Mathematica ver. 10.4 has been used. The approximate norm-2 of absolute error is given as
Example 1
Notice the following FDWE [7]:
with the homogenous initial and boundary conditions, and let
Now, we apply the numerical method presented in this paper for \(k_1=k_2=1\) and \(M_1=M_2=3\). From Eq. (3.18) we have \(w(x,t)=0\) and from Eq. (3.19), for \(q=1.1, 1.3, 1.5, 1.7, 1.9\) we obtain
Thus, from (3.1) we have
which is the exact solution.
In [7], Heydari et al. used fractional operational matrix of integration for the LWs to solve this problem for \(q=1.1, 1.3, 1.5, 1.7, 1.9\). The best absolute error of the approximate solutions at some different points, in [7], with \(k_1=k_2=3\) and \(M_1=M_2=3\) is \(1.6695 \times 10^{-6}\).
Example 2
Consider the following FDWE [3, 7]:
where
The initial and boundary conditions are determined correspondingly to the exact solution \(u(x, t) =e^xt^3.\)
From (3.18), we obtain \(w(x,t)=t^3(1-x+ex)\), then apply the numerical method presented in this paper for \(k_1=k_2=1\) and \(M_1=M_2=3\). Tables 1 and 2 present, respectively, the absolute error and the \(L^2\) norm error for \(u(x,t)-{\widetilde{u}}(x,t)\) with different values of q. In Table 3, the absolute error of u(x, t) and CPU times for \(q = 1.5\) with \(k_1=k_2=1\) and different values of \(M_1=M_2\) are given. It is seen from Table 3 that, with increase in the number of the two-dimensional Bernoulli wavelets basis, the approximate values of u(x, t) converge to the exact solutions. In Figs. 1 and 2, the exact and approximate solutions, and also the absolute difference between exact and approximate solutions of u(x, t) with \(q=1.1\) and \(q=1.9\) are plotted, respectively. Also, the graphs of the absolute errors for \(q = 1.5\) at \(t = 0.5\) and \(x = 0.5\) are shown in Fig. 3.
Liu et al. [3] employed the fractional predictor–corrector method and solved this problem with \(q=1.85\) and different values of time and space step sizes. They obtained \(1.6341 \times 10^{-3}\) for the best maximum absolute error. Moreover, in Table 1 we compare our results with obtained results in [7]. These comparisons show the more accuracy and less computations of our technique in comparison to other published methods.
Conclusion
In this paper, the two-dimensional BWs was defined. Then the satisfier function in Ritz–Galerkin method with the two-dimensional BWs basis was successfully applied to solve the second-order time FDWE. Using of satisfier function in the Ritz–Galerkin method is an efficient tool to put on the initial and boundary conditions. Furthermore, a small number of basis elements were sufficient to derive accurate numerical solutions. Also, our results were compared with obtained results in [3] and [7]. Comparison for the numerical examples shows the more accuracy and less computations of the proposed method in comparison to other published methods.
References
Chen, J., Liu, F., Anh, V., Shen, S., Liu, Q., Liao, C.: The analytical solution and numerical solution of the fractional diffusion-wave equation with damping. Appl. Math. Comput. 219, 1737–1748 (2012)
Bhrawy, A.H., Dohac, E.H., Baleanud, D., Ezz-Eldien, S.S.: A spectral tau algorithm based on Jacobi operational matrix for numerical solution of time fractional diffusion-wave equations. J. Comput. Phys. 293, 142–156 (2015)
Liu, F., Meerschaert, M.M., McGough, R.J., Zhuang, P., Liu, Q.: Numerical methods for solving the multi-term time-fractional wave-diffusion equation. Fract. Calc. Appl. Anal. 16, 9–25 (2013)
Mainardi, F.: The fundamental solutions for the fractional diffusion-wave equation. Appl. Math. Lett. 9, 23–28 (1996)
Ren, J., Gao, G.: Efficient and stable numerical methods for the two-dimensional fractional Cattaneo equation. Numer. Algorithms 1(24), 876–895 (2014)
Wess, W.: The fractional diffusion equation. J. Math. Phys. 27, 2782–2785 (1996)
Heydari, M., Hooshmandasl, M.R., Ghaini, F.M., Cattani, C.: Wavelets method for the time fractional diffusion-wave equation. Phys. Lett. A 379, 71–76 (2015)
Hu, X., Zhang, L.: A compact finite difference scheme for the fourth-order fractional diffusion-wave system. Comput. Phys. Commun. 182, 1645–1650 (2011)
Zeng, F.: Second-order stable finite difference schemes for the time-fractional diffusion-wave equation. J. Sci. Comput. 1–20 (2014)
Hu, X., Zhang, L.: On finite difference methods for fourth-order fractional diffusion-wave and subdiffusion systems. Appl. Math. Comput. 218, 5019–5034 (2012)
Darzi, R., Mohammadzade, B., Mousavi, S., Beheshti, R.: Sumudu transform method for solving fractional differential equations and fractional diffusion-wave equation. J. Math. Comput. Sci. 6, 79–84 (2013)
Hosseini, V., Shivanian, E., Chen, W.: Local radial point interpolation (MLRPI) method for solving time fractional diffusion-wave equation with damping. J. Comput. Phys. 312, 307–332 (2016)
Cui, M.: Convergence analysis of high-order compact alternating direction implicit schemes for the two-dimensional time fractional diffusion equation. Numer. Algorithms 62, 383–409 (2013)
Yousefi, S.A., Barikbin, Z.: Ritz Legendre multiwavelet method for the damped generalized regularized long-wave equation. J. Comput. Nonlinear Dyn. 7, 1–4 (2011)
Yousefi, S.A., Barikbin, Z., Dehghan, M.: Ritz–Galerkin method with Bernstein polynomial basis for finding the product solution form of heat equation with non-classic boundary conditions. Int. J. Numer. Methods Heat Fluid Flow 22, 39–48 (2012)
Yousefi, S.A., Lesnic, D., Barikbin, Z.: Satisfier function in Ritz–Galerkin method for the identification of a time-dependent diffusivity. J. Inverse Ill Posed Probl. 20, 701–722 (2012)
Barikbin, Z., Ellahi, R., Abbasbandy, S.: The Ritz–Galerkin method for MHD Couette flow of non-Newtonian fluid. Int. J. Ind. Math. 6, 235–243 (2014)
Lesnic, D., Yousefi, S.A., Ivanchov, M.: Determination of a time-dependent diffusivity from nonlocal conditions. J. Appl. Math. Comput. 41, 301–320 (2013)
Rashedi, K., Adibi, H., Dehghan, M.: Application of the Ritz–Galerkin method for recovering the spacewise-coefficients in the wave equation. Comput. Math. Appl. 65, 1990–2008 (2013)
Rashedi, K., Adibi, H., Dehghan, M.: Determination of space-time dependent heat source in a parabolic inverse problem via the Ritz–Galerkin technique. Inverse Probl. Sci. Eng. 22, 1077–1108 (2014)
Chui, C.K.: Wavelets: A Mathematical Tool for Signal Analysis. SIAM, Philadelphia (1997)
Shamsi, M., Razzaghi, M.: Solution of Hallen’s integral equation using multiwavelets. Comput. Phys. Commun. 168, 187–197 (2005)
Lakestani, M., Razzaghi, M., Dehghan, M.: Semi orthogonal spline wavelets approximation for Fredholm integro-differential equations. Math. Probl. Eng. 1–12 (2006)
Beylkin, G., Coifman, R., Rokhlin, V.: Fast wavelet transforms and numerical algorithms I. Commun. Pure Appl. Math. 44, 141–183 (1991)
Marzban, H., Razzaghi, M.: Hybrid functions approach for linearly constrained quadratic optimal control problems. Appl. Math. Model. 27, 471–485 (2003)
Razzaghi, M., Elnagar, G.: Linear quadratic optimal control problems via shifted Legendre state parametrization. Int. J. Syst. Sci. 25, 393–399 (1994)
Mashayekhi, S., Ordokhani, Y., Razzaghi, M.: Hybrid functions approach for optimal control of systems described by integro-differential equations. Appl. Math. Model. 37, 3355–3368 (2013)
Guf, J.S., Jiang, W.S.: The Haar wavelets operational matrix of integration. Int. J. Syst. Sci. 27, 623–628 (1996)
Keshavarz, E., Ordokhani, Y., Razzaghi, M.: Bernoulli wavelet operational matrix of fractional order integration and its applications in solving the fractional order differential equations. Appl. Math. Model. 38(24), 6038–6051 (2014)
Kreyszig, E.: Introductory Functional Analysis with Applications. Wiley, New York (1978)
Costabile, F., Dellaccio, F., Gualtieri, M.I.: A new approach to Bernoulli polynomials. Rend. Mat. 26, 1–12 (2006)
Arfken, G.: Mathematical Methods for Physicists, 3rd edn. Academic press, San Diego (1985)
Maleknejad, K., Hashemizadeh, E., Basirat, B.: Computational method based on Bernstein operational matrices for nonlinear Volterra–Fredholm–Hammerstein integral equations. Commun. Nonlinear Sci. Numer. Simulat. 17, 52–61 (2012)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Barikbin, Z. Two-dimensional Bernoulli wavelets with satisfier function in the Ritz–Galerkin method for the time fractional diffusion-wave equation with damping. Math Sci 11, 195–202 (2017). https://doi.org/10.1007/s40096-017-0214-4
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s40096-017-0214-4