The new numerical method for solving the system of two-dimensional Burgers’ equations

In this paper, the system of two-dimensional Burgers’ equations are solved by local discontinuous Galerkin (LDG) finite element method. The new method is based on the two-dimensional Hopf–Cole transformations, which transform the system of two-dimensional Burgers’ equations into a linear heat equation. Then the linear heat equation is solved by the LDG finite element method. The numerical solution of the heat equation is used to derive the numerical solutions of Burgers’ equations directly. Such a LDG method can also be used to find the numerical solution of the two-dimensional Burgers’ equation by rewriting Burgers’ equation as a system of the two-dimensional Burgers’ equations. Three numerical examples are used to demonstrate the efficiency and accuracy of the method.     

Sensitivity analysis of shock wave Burgers’ equation via a novel algorithm based on scale-3 Haar wavelets

In this paper, a novel technique is being formulated for the numerical solutions of Shock wave Burgers' equations for planar and non-planar geometry. It is well known that Burgers' equation is sensitive to the perturbations in the diffusion term. Thus we use robustness of wavelets generated by dilation and translation of Haar wavelets on third scale to capture the sensitivity information. The present approach is an improved form of the scale-2 Haar wavelet method. The scheme is based on the forward finite difference scheme for time integration, scale-3 Haar wavelets for space integration and the nonlinearity has been tackled via quasilinearzation technique. Through scale-3 Haar wavelet analysis once the wavelet coefficient is calculated then we can compute the solutions at near the perturbation point. The computation cost of the present scheme is negligible. The proposed method is tested on six test problems to check its computational efficiency where the convergence analysis of scale-3 Haar wavelet method is the proof of our computational arguments.     

An efficient numerical algorithm for the study of time fractional Tricomi and Keldysh type equations

Engineering with Computers - This work addresses a hybrid scheme for the numerical solutions of time fractional Tricomi and Keldysh type equations. In proposed methodology, Haar wavelets are used...     

A new approach for numerical solution of integro-differential equations via Haar wavelets

A new method is proposed for numerical solution of Fredholm and Volterra integro-differential equations of second kind. The proposed method is based on Haar wavelets approximation. Special characteristics of Haar wavelets approximation has been used in the derivation of this method. The new method is the extension of the recent work [Aziz and Siraj-ul-Islam, New algorithms for numerical solution of nonlinear Fredholm and Volterra integral equations using Haar wavelets, J. Comput. Appl. Math. 239 (2013), pp. 333–345] from integral equations to integro-differential equations. The method is specifically derived for nonlinear problems. Two new algorithms are also proposed based on this new method, one each for numerical solution of Fredholm and Volterra integro-differential equations. The proposed algorithms are generic and are applicable to all types of both nonlinear Fredholm and Volterra integro-differential equations of second kind. The cost of the new algorithms is considerably reduced by using the Broyden's method instead of Newton's method for solution of system of nonlinear equations. Most of the numerical methods designed for solution of integro-differential equations rely on some other technique for numerical integration. The advantage of our method is that it does not use numerical integration. The integrand is approximated using Haar wavelets approximation and then exact integration is performed. The method is tested on number of problems and numerical results are compared with existing methods in the literature. The numerical results indicate that accuracy of the obtained solutions is reasonably high even when the number of collocation points is small.     

A Haar wavelet quasilinearization approach for numerical simulation of Burgers’ equation

In this paper, an efficient numerical scheme based on uniform Haar wavelets and the quasilinearization process is proposed for the numerical simulation of time dependent nonlinear Burgers’ equation. The equation has great importance in many physical problems such as fluid dynamics, turbulence, sound waves in a viscous medium etc. The Haar wavelet basis permits to enlarge the class of functions used so far in the collocation framework. More accurate solutions are obtained by wavelet decomposition in the form of a multi-resolution analysis of the function which represents a solution of boundary value problems. The accuracy of the proposed method is demonstrated by three test problems. The numerical results are compared with existing numerical solutions found in the literature. The use of the uniform Haar wavelet is found to be accurate, simple, fast, flexible, convenient and has small computation costs.     

A comparative study of numerical integration based on Haar wavelets and hybrid functions

A quadrature rule based on uniform Haar wavelets and hybrid functions is proposed to find approximate values of definite integrals. The wavelet-based algorithm can be easily extended to find numerical approximations for double, triple and improper integrals. The main advantage of this method is its efficiency and simple applicability. Error estimates of the proposed method alongside numerical examples are given to test the convergence and accuracy of the method.     

A novel numerical method on the resolution of the time-invariant systems response

An algorithm for evaluation of time response to a time-invariant system based on Haar wavelets is proposed. The system is assumed to be governed by a high order linear differential equation with constant coefficients. The main objective of this study is to convert a differential equation into an algebraic-form equation and provide an elegant approach for computer programming. By the method, the computation complexity can be greatly reduced or much simplified in calculate the time response of a time-invariant system. The accurate and fastness of the method have been demonstrated by the examples of integration of the stiff systems     

An efficient numerical algorithm for multi-dimensional time dependent partial differential equations

An efficient and robust numerical scheme based on Haar wavelets and finite differences is suggested for the solution of two-dimensional time dependent linear and nonlinear partial differential equations (PDEs). Excellent feature of the scheme is the conversion of linear and non-linear PDEs to algebraic equations which are comparatively easy to handle. Convergence of the scheme, which guarantees small error norm as the resolution level increases, is also an important part of this work. Different error norms are computed to check efficiency of the technique. Computations verify accuracy, flexibility and low computational cost of the method.     

A non-uniform Haar wavelet method for numerically solving two-dimensional convection-dominated equations and two-dimensional near singular elliptic equations

A non-uniform Haar wavelet based collocation method has been developed in this paper for two-dimensional convection dominated equations and two-dimensional near singular elliptic partial differential equations, in which traditional Haar wavelet method produces oscillatory solutions or low accurate solutions. The main idea behind the proposed method is to transform the computation of numerical solution of considered partial differential equations to computation of solution of a linear system of equations. This process is done by discretizing space variables with non-uniform Haar wavelets. To confirm efficiency of the proposed method seven benchmark problems are solved and the obtained results are compared with exact solutions and with local meshless methods, finite element method, finite difference method and polynomial collocation method. Numerical experiments show that the proposed method gives convincing results even in less number of collocation nodes.     

Identification of linear time varying systems by Haar wavelet

This study addresses the identification of linear time varying systems. The identification is based on the expansion of all time functions in the state equations by Haar wavelets. The unknown time function can thus be identified in terms of Haar wavelets. A Haar wavelet is a set of complete, orthogonal basis and is easy to use computations. Several good properties of Haar wavelets are utilized in the algorithm. Both numerical and experimental results verify the analysis.     

Numerical solutions of two dimensional Sobolev and generalized Benjamin–Bona–Mahony–Burgers equations via Haar wavelets

This work addresses, numerical method based on Haar wavelets and finite differences to solve two dimensional linear, nonlinear Sobolev and non-linear generalized Benjamin–Bona–Mahony–Burgers (NGBBMB) equations. The temporal part is discretized using finite differences while spatial part is approximated by two dimensional Haar wavelets. With this strategy, computing solution of two dimensional PDEs reduces to computing solution of linear system of algebraic equations. Collocation approach is then applied to determine the wavelet coefficients from linear system. This paper shows that two dimensional Haar wavelets are suitable and effective for two dimensional linear and non-linear PDEs. For validation of the proposed scheme different problems have been solved and error norms $L_{\infty },L_2$ are computed. Computation verifies that current scheme has good outcome.     

一维对流扩散反应方程的高精度数值方法

通过将原方程变换为对流扩散方程,将所得方程的对流项采用四阶组合紧致迎风格式离散,扩散项采用四阶对称紧致格式离散之后,对得到的空间半离散格式采用四阶龙格库塔方法进行时间推进,得到了一种求解非定常对流扩散反应问题的高精度方法,其收敛阶为O(h4+τ4).经数值实验并与文献结果进行对比,表明该格式适用于对流占优问题的数值模拟,验证了格式的良好性能.     

A haar wavelet approximation for two-dimensional time fractional reaction–subdiffusion equation

Engineering with Computers - In this study, we established a wavelet method, based on Haar wavelets and finite difference scheme for two-dimensional time fractional reaction–subdiffusion...     

Numerical solution of stiff differential equations via Haar wavelets

A simple and effective algorithm based on Haar wavelets is proposed to the solution of stiff problems in this article. It can integrate the stiff equation with very accurate results for any length of time. The simulation result shows that the single-term Haar wavelet method is better than the improved Runge–Kutta–Fehlberg method, while the terms of both expansions are the same.     

基于提升方案的高维形态小波构造

研究基于提升方案的高维形态小波构造方法.结合数学形态学和高维小波分析来构造 高维形态小波.提出高维多通道形态小波构造的一般性框架--高维多通道提升方案.采用最大 值(最小值)形态算子构造一维多通道Haar形态小波和高维多通道Haar形态小波.最后通过两 个例子说明形态小波变换的优点和可行性.     

Quadrature rules for numerical integration based on Haar wavelets and hybrid functions

In this paper Haar wavelets and hybrid functions have been applied for numerical solution of double and triple integrals with variable limits of integration. This approach is the generalization and improvement of the methods (Siraj-ul-Islam et al. (2010) [9]) where the numerical methods are only applicable to the integrals with constant limits. Apart from generalization of the methods [9], the new approach has two major advantages over the classical methods based on quadrature rule: (i) No need of finding optimum weights as the wavelet and hybrid coefficients serve the purpose of optimal weights automatically (ii) Mesh points of the wavelets algorithm are used as nodal values instead of considering the n nodes as unknown roots of polynomial of degree n. The new methods are more efficient. The novel methods are compared with existing methods and applied to a number of benchmark problems. Accuracy of the methods are measured in terms of absolute errors.     

Wavelet based iterative methods for a class of 2D-partial integro differential equations

In this paper, an iterative method based on quasilinearization is presented to solve a class of two dimensional partial integro differential equations that arise in nuclear reactor models and population models. Two different approaches based on Haar and Legendre wavelets are studied to develop numerical methods. In the first approach, time domain is approximated with the help of forward finite difference approach. In the second approach, both time as well as space domains are approximated by wavelets. Appropriate examples are solved using these methods and the obtained results are compared with the methods available in the recent literature.     

On the numerical solution of space–time fractional diffusion models

A flexible numerical scheme for the discretization of the space–time fractional diffusion equation is presented. The model solution is discretized in time with a pseudo-spectral expansion of Mittag–Leffler functions. For the space discretization, the proposed scheme can accommodate either low-order finite-difference and finite-element discretizations or high-order pseudo-spectral discretizations. A number of examples of numerical solutions of the space–time fractional diffusion equation are presented with various combinations of the time and space derivatives. The proposed numerical scheme is shown to be both efficient and flexible.     

Sufficient conditions for the preservation of the boundedness in a numerical method for a physical model with transport memory and nonlinear damping

Departing from a finite-difference scheme to approximate solutions of a nonlinear, hyperbolic partial differential equation which generalizes the Burgers–Huxley equation from fluid dynamics, we investigate conditions on the model coefficients and the computational parameters under which positive and bounded initial data evolve into positive and bounded new approximations. The model under investigation includes nonlinear coefficients of damping and advection, and the reaction term extends the reaction law of the classical Fisher–Kolmogorov–Petrovsky–Piscounov equation. The method can be expressed in vector form in terms of a multiplicative matrix which, under certain parametric conditions, becomes an M-matrix. Using the fact that every M-matrix is non-singular and that the entries of its inverse are positive, real numbers, we establish sufficient conditions under which the method provides new, positive and bounded approximations from previous, positive and bounded data and boundary conditions. The numerical results confirm the fact that the conditions derived here are sufficient for the positivity and the boundedness of the approximations; moreover, computational experiments evidence the fact that the method still preserves these properties for values of the model and the numerical parameters outside of the analytic regions of positivity and boundedness. We point out that our simulations show a good agreement between the numerical approximations computed through our method and the corresponding, analytical solutions.