A numerical algorithm based on modified cubic trigonometric B-spline functions for computational modelling of hyperbolic-type wave equations

2017 ◽  
Vol 34 (4) ◽  
pp. 1257-1276 ◽  
Author(s):  
Ali Saleh Alshomrani ◽  
Sapna Pandit ◽  
Abdullah K. Alzahrani ◽  
Metib Said Alghamdi ◽  
Ram Jiwari

Purpose The main purpose of this work is the development of a numerical algorithm based on modified cubic trigonometric B-spline functions for computational modelling of hyperbolic-type wave equations. These types of equations describe a variety of physical models in the vibrations of structures, nonlinear optics, quantum field theory and solid-state physics, etc. Design/methodology/approach Dirichlet boundary conditions cannot be handled easily by cubic trigonometric B-spline functions. Then, a modification is made in cubic trigonometric B-spline functions to handle the Dirichlet boundary conditions and a numerical algorithm is developed. The proposed algorithm reduced the hyperbolic-type wave equations into a system of first-order ordinary differential equations (ODEs) in time variable. Then, stability-preserving SSP-RK54 scheme and the Thomas algorithm are used to solve the obtained system. The stability of the algorithm is also discussed. Findings A different technique based on modified cubic trigonometric B-spline functions is proposed which is quite different from the schemes developed (Abbas et al., 2014; Nazir et al., 2016) and depicts the computational modelling of hyperbolic-type wave equations. Originality/value To the best of the authors’ knowledge, this technique is novel for solving hyperbolic-type wave equations and the developed algorithm is free from quasi-linearization process and finite difference operators for time derivatives. This algorithm gives better results than the results discussed in literature (Dehghan and Shokri, 2008; Batiha et al., 2007; Mittal and Bhatia, 2013; Jiwari, 2015).

2017 ◽  
Vol 27 (8) ◽  
pp. 1638-1661 ◽  
Author(s):  
Ram Jiwari ◽  
Ali Saleh Alshomrani

Purpose The main aim of the paper is to develop a new B-splines collocation algorithm based on modified cubic trigonometric B-spline functions to find approximate solutions of nonlinear parabolic Burgers’-type equations with Dirichlet boundary conditions. Design/methodology/approach A modification is made in cubic trigonometric B-spline functions to handle the Dirichlet boundary conditions and an algorithm is developed with the help of modified cubic trigonometric B-spline functions. The proposed algorithm reduced the Burgers’ equations into a system of first-order nonlinear ordinary differential equations in time variable. Then, strong stability preserving Runge-Kutta 3rd order (SSP-RK3) scheme is used to solve the obtained system. Findings A different technique based on modified cubic trigonometric B-spline functions is proposed which is quite different from to the schemes developed in Abbas et al. (2014) and Nazir et al. (2016), and the developed algorithms are free from linearization process and finite difference operators. Originality/value To the best knowledge of the authors, this technique is novel for solving nonlinear partial differential equations, and the new proposed technique gives better results than the results discussed in Ozis et al. (2003), Kutluay et al. (1999), Khater et al. (2008), Korkmaz and Dag (2011), Kutluay et al. (2004), Rashidi et al. (2009), Mittal and Jain (2012), Mittal and Jiwari (2012), Mittal and Tripathi (2014), Xie et al. (2008) and Kadalbajoo et al. (2005).


2017 ◽  
Vol 34 (8) ◽  
pp. 2793-2814 ◽  
Author(s):  
Sapna Pandit ◽  
Ram Jiwari ◽  
Karan Bedi ◽  
Mehmet Emir Koksal

Purpose The purpose of this study is to develop an algorithm for approximate solutions of nonlinear hyperbolic partial differential equations. Design/methodology/approach In this paper, an algorithm based on the Haar wavelets operational matrix for computational modelling of nonlinear hyperbolic type wave equations has been developed. These types of equations describe a variety of physical models in nonlinear optics, relativistic quantum mechanics, solitons and condensed matter physics, interaction of solitons in collision-less plasma and solid-state physics, etc. The algorithm reduces the equations into a system of algebraic equations and then the system is solved by the Gauss-elimination procedure. Some well-known hyperbolic-type wave problems are considered as numerical problems to check the accuracy and efficiency of the proposed algorithm. The numerical results are shown in figures and Linf, RMS and L2 error forms. Findings The developed algorithm is used to find the computational modelling of nonlinear hyperbolic-type wave equations. The algorithm is well suited for some well-known wave equations. Originality/value This paper extends the idea of one dimensional Haar wavelets algorithms (Jiwari, 2015, 2012; Pandit et al., 2015; Kumar and Pandit, 2014, 2015) for two-dimensional hyperbolic problems and the idea of this algorithm is quite different from the idea for elliptic problems (Lepik, 2011; Shi et al., 2012). Second, the algorithm and error analysis are new for two-dimensional hyperbolic-type problems.


2018 ◽  
Vol 28 (11) ◽  
pp. 2620-2649 ◽  
Author(s):  
Rajni Rohila ◽  
R.C. Mittal

Purpose This paper aims to develop a novel numerical method based on bi-cubic B-spline functions and alternating direction (ADI) scheme to study numerical solutions of advection diffusion equation. The method captures important properties in the advection of fluids very efficiently. C.P.U. time has been shown to be very less as compared with other numerical schemes. Problems of great practical importance have been simulated through the proposed numerical scheme to test the efficiency and applicability of method. Design/methodology/approach A bi-cubic B-spline ADI method has been proposed to capture many complex properties in the advection of fluids. Findings Bi-cubic B-spline ADI technique to investigate numerical solutions of partial differential equations has been studied. Presented numerical procedure has been applied to important two-dimensional advection diffusion equations. Computed results are efficient and reliable, have been depicted by graphs and several contour forms and confirm the accuracy of the applied technique. Stability analysis has been performed by von Neumann method and the proposed method is shown to satisfy stability criteria unconditionally. In future, the authors aim to extend this study by applying more complex partial differential equations. Though the structure of the method seems to be little complex, the method has the advantage of using small processing time. Consequently, the method may be used to find solutions at higher time levels also. Originality/value ADI technique has never been applied with bi-cubic B-spline functions for numerical solutions of partial differential equations.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Neeraj Dhiman ◽  
M.J. Huntul ◽  
Mohammad Tamsir

Purpose The purpose of this paper is to present a stable and efficient numerical technique based on modified trigonometric cubic B-spline functions for solving the time-fractional diffusion equation (TFDE). The TFDE has numerous applications to model many real objects and processes. Design/methodology/approach The time-fractional derivative is used in the Caputo sense. A modification is made in trigonometric cubic B-spline (TCB) functions for handling the Dirichlet boundary conditions. The modified TCB functions have been used to discretize the space derivatives. The stability of the technique is also discussed. Findings The obtained results are compared with those reported earlier showing that the present technique gives highly accurate results. The stability analysis shows that the method is unconditionally stable. Furthermore, this technique is efficient and requires less storage. Originality/value The current work is novel for solving TFDE. This technique is unconditionally stable and gives better results than existing results (Ford et al., 2011; Sayevand et al., 2016; Ghanbari and Atangana, 2020).


2012 ◽  
Vol 2012 ◽  
pp. 1-15 ◽  
Author(s):  
Gang Li ◽  
Linghui Hong ◽  
Wenjun Liu

We consider viscoelastic wave equations of the Kirchhoff typeutt-M(∥∇u∥22)Δu+∫0tg(t-s)Δu(s)ds+ut=|u|p-1uwith Dirichlet boundary conditions, where∥⋅∥pdenotes the norm in the Lebesgue spaceLp. Under some suitable assumptions ongand the initial data, we establish a global nonexistence result for certain solutions with arbitrarily high energy, in the sense thatlim⁡t→T*-(∥u(t)∥22+∫0t∥u(s)∥22ds)=∞for some0<T*<+∞.


2015 ◽  
Vol 32 (5) ◽  
pp. 1275-1306 ◽  
Author(s):  
R C Mittal ◽  
Amit Tripathi

Purpose – The purpose of this paper is to develop an efficient numerical scheme for non-linear two-dimensional (2D) parabolic partial differential equations using modified bi-cubic B-spline functions. As a test case, method has been applied successfully to 2D Burgers equations. Design/methodology/approach – The scheme is based on collocation of modified bi-cubic B-Spline functions. The authors used these functions for space variable and for its derivatives. Collocation form of the partial differential equation results into system of first-order ordinary differential equations (ODEs). The obtained system of ODEs has been solved by strong stability preserving Runge-Kutta method. The computational complexity of the method is O(p log(p)), where p denotes total number of mesh points. Findings – Obtained numerical solutions are better than those available in literature. Ease of implementation and very small size of computational work are two major advantages of the present method. Moreover, this method provides approximate solutions not only at the grid points but also at any point in the solution domain. Originality/value – First time, modified bi-cubic B-spline functions have been applied to non-linear 2D parabolic partial differential equations. Efficiency of the proposed method has been confirmed with numerical experiments. The authors conclude that the method provides convergent approximations and handles the equations very well in different cases.


2018 ◽  
Vol 35 (1) ◽  
pp. 315-333
Author(s):  
Ramesh Chand Mittal ◽  
Sumita Dahiya

Purpose In this study, a second-order standard wave equation extended to a two-dimensional viscous wave equation with timely differentiated advection-diffusion terms has been solved by differential quadrature methods (DQM) using a modification of cubic B-spline functions. Two numerical schemes are proposed and compared to achieve numerical approximations for the solutions of nonlinear viscous wave equations. Design/methodology/approach Two schemes are adopted to reduce the given system into two systems of nonlinear first-order partial differential equations (PDE). For each scheme, modified cubic B-spline (MCB)-DQM is used for calculating the spatial variables and their derivatives that reduces the system of PDEs into a system of nonlinear ODEs. The solutions of these systems of ODEs are determined by SSP-RK43 scheme. The CPU time is also calculated and compared. Matrix stability analysis has been performed for each scheme and both are found to be unconditionally stable. The results of both schemes have been extensively discussed and compared. The accuracy and reliability of the methods have been successfully tested on several examples. Findings A comparative study has been carried out for two different schemes. Results from both schemes are also compared with analytical solutions and the results available in literature. Experiments show that MCB-DQM with Scheme II yield more accurate and reliable results in solving viscous wave equations. But Scheme I is comparatively less expensive in terms of CPU time. For MCB-DQM, less depository requirements lead to less aggregation of approximation errors which in turn enhances the correctness and readiness of the numerical techniques. Approximate solutions to the two-dimensional nonlinear viscous wave equation have been found without linearizing the equation. Ease of implementation and low computation cost are the strengths of the method. Originality/value For the first time, a comparative study has been carried out for the solution of nonlinear viscous wave equation. Comparisons are done in terms of accuracy and CPU time. It is concluded that Scheme II is more suitable.


Sign in / Sign up

Export Citation Format

Share Document