B-spline Collocation Method Research Articles

Numerical solution and analysis of extended Fisher–Kolmogorov equation using an improved collocation algorithm

The paper investigated the approximate solution of the extended Fisher–Kolmogorov. For this, an advanced approach: the improved quintic B-spline collocation technique has been employed, which is an enhancement over the conventional B-spline collocation method. The B-spline interpolant of degree five has been refined through the posteriori corrections, leading to the development of the improved B-spline solution. This proposed technique demonstrates superior convergence compared to the standard B-spline method, as assessed both theoretically and numerically. In tackling the extended Fisher–Kolmogorov equation, the space discretization is conducted using the improved quintic B-spline collocation methodology (IQSCM), and for the temporal domain discretization, the Crank-Nicolson scheme is applied. The error bounds and convergence analysis is established using the Green's functions. Von-Neumann stability analysis is carried out to discuss the stability of the technique. A few examples are solved and are represented graphically which helps in determining the nature of the solution. Also, L 2 , L ∞ error norms, and order of convergence are calculated to demonstrate the contribution of the new improved technique over the standard spline collocation technique.

Application of the Cubic Exponential B-spline Collocation Operator Splitting Method for Numerical solutions of Burgers Equation

In this study, the cubic exponential B-spline collocation method has been proposed for the numerical solutions of the Burgers equation with the operator splitting. To apply the operator splitting method, the Burgers' equation has decomposed into two sub-equations based on the time term: the linear part (diffusion) and the nonlinear part (convection). Subsequently, for each sub-equation, Crank-Nicolson finite difference schemes in the temporal direction and cubic exponential B-spline functions and their derivatives, have applied at the x_m nodal points in the spatial direction. The algebraic equation systems obtained have been solved numerically using the Lie-Trotter and Strang splitting schemes to get the solutions of the main equation. Some advantages of the splitting methods include preserving the physical characteristics of the solution, yielding more convergent results over long time intervals, enabling simpler algorithms, and facilitating the storage of solution vectors on computer. To assess the accuracy of the computed numerical results the L_2 and L_∞ error norms have been used. Additionally, the obtained results have been compared with some studies in the literature. The stability analysis of the applied method has been investigated using the von Neumann Fourier series method.

Extended B-Spline Collocation Method for Singularly Perturbed Time Delay Parabolic Differential-Difference Problems

The paper presents a robust collocation B-spline numerical scheme for solving singularly perturbed time lag parabolic differential–difference problems. The method uses Taylor’s series expansion for the approximation of retarded terms, the [Formula: see text]-method (when [Formula: see text]) for time discretization, and the extended third-degree B-spline collocation method for spatial discretization in a piecewise uniform Shishkin mesh. The advantage of using an extended cubic B-spline rather than the classical third-degree B-spline method is that it introduces one additional free parameter to control the global shape parameter. The stability and uniform convergence of the proposed scheme are investigated. The method is first-order accurate in [Formula: see text] and almost second-order accurate in [Formula: see text], surpassing existing literature methods.

Cubic B-spline based elastic and viscoelastic wave propagation method

The requirement for high computational efficiency in inversion imaging poses challenges to further enhancing the accuracy of imaging large elastic and viscoelastic models. Considering that forward modeling serves as a prerequisite for inversion imaging, it is essential to introduce a wave propagation simulation method that can enhance simulation accuracy without significantly compromising computational speed. In our paper, cubic B-spline is introduced into the forward modeling of the wave equation. The B-spline collocation method offers a linear complexity for computation and exhibits fourth-order convergence in terms of computational errors. We also find that it can perfectly adapt to the boundary conditions of the perfectly matched layer in the forward modeling solution of the wave simulation. Several typical numerical examples, including 2-D acoustic wave equation, 2-D elastic wave propagation, and 2-D viscoelastic propagation in homogeneous media, are provided to validate its convergence, accuracy, and computational efficiency.

Simulation of coupled groundwater flow and contaminant transport using quintic B-spline collocation method

PurposeThe purpose of this study is to introduce a new numerical scheme with no stability condition and high-order accuracy for the solution of two-dimensional coupled groundwater flow and transport simulation problems with regular and irregular geometries and compare the results with widely acceptable programs such as Modular Three-Dimensional Finite-Difference Ground-Water Flow Model (MODFLOW) and Modular Three-Dimensional Multispecies Transport Model (MT3DMS).Design/methodology/approachThe newly proposed numerical scheme is based on the method of lines (MOL) approach and uses high-order approximations both in space and time. Quintic B-spline (QBS) functions are used in space to transform partial differential equations, representing the relevant physical phenomena in the system of ordinary differential equations. Then this system is solved with the DOPRI5 algorithm that requires no stability condition. The obtained results are compared with the results of the MODFLOW and MT3DMS programs to verify the accuracy of the proposed scheme.FindingsThe results indicate that the proposed numerical scheme can successfully simulate the two-dimensional coupled groundwater flow and transport problems with complex geometry and parameter structures. All the results are in good agreement with the reference solutions.Originality/valueTo the best of the authors' knowledge, the QBS-DOPRI5 method is used for the first time for solving two-dimensional coupled groundwater flow and transport problems with complex geometries and can be extended to high-dimensional problems. In the future, considering the success of the proposed numerical scheme, it can be used successfully for the identification of groundwater contaminant source characteristics.

Numerical Solution of Third-Order Rosenau–Hyman and Fornberg–Whitham Equations via B-Spline Interpolation Approach

This study aims to find the numerical solution of the Rosenau–Hyman and Fornberg–Whitham equations via the quintic B-spline collocation method. Quintic B-spline, along with finite difference and theta-weighted schemes, is used for the discretization and approximation purposes. The effectiveness and robustness of the procedure is assessed by comparing the computed results with the exact and available results in the literature using absolute and relative error norms. The stability of the proposed scheme is studied using von Neumann stability analysis. Graphical representations are drawn to analyze the behavior of the solution.

Cubic B-spline based numerical schemes for delayed time-fractional advection-diffusion equations involving mild singularities

This article presents two efficient layer-adaptive numerical schemes for a class of time-fractional advection-diffusion equations with a large time delay. The fractional derivative of order α with α ∈ (0, 1) is taken in the Caputo sense. The solution to this type of problem generally has a layer due to the mild singularity near the time t = 0. Consequently, the polynomial interpolation discretizing scheme degrades the convergence rate in the case of uniform meshes. In the presence of a singularity, the temporal fractional operator is discretized by employing the L1 technique on a layer-resolving mesh. In contrast, the cubic B-spline collocation method is used in the spatial direction. The convergence analysis and estimation of error are presented for the proposed scheme under reasonable regularity assumptions on the coefficients. The scheme achieves its optimal convergence rate (2 − α) for suitable choice of grading parameter (γ ≥ (2 − α)/α). Furthermore, we modified the proposed scheme by discretizing the fractional operator with the help of the L1-2 technique. The modified scheme gets a quadratic order convergence for γ ≥ 2/α. In addition, we extend the proposed schemes to solve the corresponding semilinear problem. Numerical examples demonstrate the efficiency and applicability of the proposed techniques.

Numerical Solution to the Time-Fractional Burgers–Huxley Equation Involving the Mittag-Leffler Function

Fractional differential equations play a significant role in various scientific and engineering disciplines, offering a more sophisticated framework for modeling complex behaviors and phenomena that involve multiple independent variables and non-integer-order derivatives. In the current research, an effective cubic B-spline collocation method is used to obtain the numerical solution of the nonlinear inhomogeneous time-fractional Burgers–Huxley equation. It is implemented with the help of a θ-weighted scheme to solve the proposed problem. The spatial derivative is interpolated using cubic B-spline functions, whereas the temporal derivative is discretized by the Atangana–Baleanu operator and finite difference scheme. The proposed approach is stable across each temporal direction as well as second-order convergent. The study investigates the convergence order, error norms, and graphical visualization of the solution for various values of the non-integer parameter. The efficacy of the technique is assessed by implementing it on three test examples and we find that it is more efficient than some existing methods in the literature. To our knowledge, no prior application of this approach has been made for the numerical solution of the given problem, making it a first in this regard.

Numerical solutions of the EW and MEW equations using a fourth-order improvised B-spline collocation method

Numerical solutions of the EW and MEW equations using a fourth-order improvised B-spline collocation method

Numerical Solution of Heat Equation using Modified Cubic B-spline Collocation Method

In this paper, a collocation method is presented based on the Modified Cubic B-spline Method (MCBSM) for the numerical solution of the heat equation. The PDE is fully discretized by using the Modified Cubic B-spline basis collocation for spatial discretization and the finite difference method is used for the time discretization. A numerical example from PDE is used to evaluate the accuracy of the proposed method. The numerical results are evaluated in comparison to the exact solutions. The findings consistently indicate that the suggested technique provides good error estimates. We also discovered that our proposed method was unconditionally stable. Hence, based on the results and the efficiency of the method, the method is suitable for solving heat equation.

A novel perspective for simulations of the Modified Equal-Width Wave equation by cubic Hermite B-spline collocation method

In the current study, the Modified Equal-Width (MEW) equation will be handled numerically by a novel technique using collocation finite element method where cubic Hermite B-splines are used as trial functions. To test the accuracy and efficiency of the method, four different experimental problems; namely, the motion of a single solitary wave, interaction of two solitary waves, interaction of three solitary waves and the birth of solitons with the Maxwellian initial condition will be investigated. In order to verify, the validity and reliability of the proposed method, the newly obtained error norms L2 and L∞ as well as three conservation constants have been compared with some of the other numerical results given in the literature at the same parameters. Furthermore, some wave profiles of the newly obtained numerical results have been given to demonstrate that each test problem exhibits accurate physical simulations. The advantage of the proposed method over other methods is the usage of inner points at Legendre and Chebyshev polynomial roots. This advantage results in better accuracy with less number of elements in spatial direction. The results of the numerical experiments clearly reveal that the presented scheme produces more accurate results even with comparatively coarser grids.

Exact Solution of Couple Burgers’ Equation using Cubic B-Spline Collocation Method for Fluid Suspension/Colloid under the Influence of Gravity

In this research, Cubic B-Spline plane method was used to solve numerically the one-dimensional Burger’s Equation with initial condition , boundary conditions; . The cubic trigonometric B-spline was used for interpolating the solutions at each time and using the Von-Neumann method to check the stability. The obtained numerical result showed that the method was efficient, robust and reliable for solving Burgers’ Equation accurately even involving high Reynolds numbers for which the exact solutions have failed.

Utilizing Cubic B-Spline Collocation Technique for Solving Linear and Nonlinear Fractional Integro-Differential Equations of Volterra and Fredholm Types

Fractional integro-differential equations (FIDEs) of both Volterra and Fredholm types present considerable challenges in numerical analysis and scientific computing due to their complex structures. This paper introduces a novel approach to address such equations by employing a Cubic B-spline collocation method. This method offers a robust and systematic framework for approximating solutions to the FIDEs, facilitating precise representations of complex phenomena. Within this research, we establish the mathematical foundations of the proposed scheme, elucidate its advantages over existing methods, and demonstrate its practical utility through numerical examples. We adopt the Caputo definition for fractional derivatives and conduct a stability analysis to validate the accuracy of the method. The findings showcase the precision and efficiency of the scheme in solving FIDEs, highlighting its potential as a valuable tool for addressing a wide array of practical problems.

Numerical Analysis of the Discrete MRLW Equation for a Nonlinear System Using the Cubic B-Spline Collocation Method

By employing the cubic B-spline functions, a collocation approach was devised in this study to address the Modified Regularized Long Wave (MRLW) equation. Then, we derived the corresponding nonlinear system and easily solved it using Newton’s iterative approach. It was established that the cubic B-spline collocation technique exhibits unconditional stability. The dynamics of solitary waves, including their pairwise and triadic interactions, were meticulously investigated utilizing the proposed numerical method. Additionally, the transformation of the Maxwellian initial condition into solitary wave formations is presented. To validate the current work, three distinct scenarios were compared against the analytical solution and outcomes from alternative methods under both L2- and L∞-error norms. Primarily, the key strength of the suggested scheme lies in its capacity to yield enhanced numerical resolutions when employed to solve the MRLW equation, and these conservation laws show that the solitary waves have time and space translational symmetry in the propagation process. Finally, this paper concludes with a summary of our findings.

B-spline Collocation Methods for Solving Linear Two-Point Boundary Value Problems

B-spline Collocation Methods for Solving Linear Two-Point Boundary Value Problems

A high-order B-spline collocation method for solving a class of nonlinear singular boundary value problems

A high-order B-spline collocation method for solving a class of nonlinear singular boundary value problems

Hybrid cubic and hyperbolic b-spline collocation methods for solving fractional Painlevé and Bagley-Torvik equations in the Conformable, Caputo and Caputo-Fabrizio fractional derivatives

In this paper, we approximate the solution of fractional Painlevé and Bagley-Torvik equations in the Conformable (Co), Caputo (C), and Caputo-Fabrizio (CF) fractional derivatives using hybrid hyperbolic and cubic B-spline collocation methods, which is an extension of the third-degree B-spline function with more smoothness. The hybrid B-spline function is flexible and produces a system of band matrices that can be solved with little computational effort. In this method, three parameters m, η, and λ play an important role in producing accurate results. The proposed methods reduce to the system of linear or nonlinear algebraic equations. The stability and convergence analysis of the methods have been discussed. The numerical examples are presented to illustrate the applications of the methods and compare the computed results with those obtained using other methods.

Efficient numerical schemes based on the cubic B-spline collocation method for time-fractional partial integro-differential equations of Volterra type

Efficient numerical schemes based on the cubic B-spline collocation method for time-fractional partial integro-differential equations of Volterra type

Third-degree B-spline collocation method for singularly perturbed time delay parabolic problem with two parameters

This study deals with a fitted third-degree B-spline collocation method for two parametric singularly perturbed parabolic problems with a time lag. The proposed method comprises the Cranck-Nicolson method for time discretization and the third-degree B-spline method spatial variable discretization. Rigorous numerical experimentations were carried out on some test examples. The obtained numerical results depict that the proposed scheme is more accurate than some methods existing in the literature. Parameter convergence analysis of the scheme is carried out and shows the present scheme is (ε−μ)−uniform convergent with the order of convergence ((Δt)2 + ℓ2).

Numerical Solution of the Heat Equation by Cubic B-spline Collocation Method

Numerical Solution of the Heat Equation by Cubic B-spline Collocation Method