Numerical solution of convection-diffusion equation using cubic B-splines collocation methods with neumann boundary conditions

International Journal of Applied Mathematics and Computation

Ramesh C Mittal¹  Rakesh Kumar Jain² 

[1] Professor, Department of Mathematics, Indian Institute of Technology Roorkee$$;Research Scholar, Department of Mathematics, Indian Institute of Technology, Roorkee$$

关键词: : convection-diffusion partial differential equation;    Neumann boundary conditions;    cubic B-splines basis functions;    SSP-RK54 scheme;    Thomas algorithm.;   

DOI  :  10.0000/ijamc.2012.4.2.470

来源: PSIT Kanpur

【 摘 要 】

In this paper, two numerical methods are proposed to approximate the solutions of the convection-diffusion partial differential equations with Neumann boundary conditions. The methods are based on collocation of cubic B-splines over finite elements so that we have continuity of the dependent variable and its first two derivatives throughout the solution range. In Method-I, we discretize the time derivative with Crank Nicolson scheme and handle spatial derivatives with cubic B-splines. Stability of this method has been discussed and shown that it is unconditionally stable. In Method-II, we apply cubic B-splines for spatial variable and derivatives which produce a system of first order ordinary differential equations. We solve this system by using SSP-RK54 scheme. These methods needs less storage space that causes to less accumulation of numerical errors. In numerical test problems, the performance of these methods is shown by computing for different time levels. Illustrative five examples are included to demonstrate the validity and applicability of these methods. Results shown by these methods are found to be in good agreement with the exact solutions. Easy and economical implementation is the strength of these methods.

【 授权许可】

Unknown