Search results

This paper aims to develop a new high accuracy numerical method based on off-step non-polynomial spline in tension approximations for the solution of Burgers-Fisher and coupled nonlinear Burgers’ equations on a graded mesh. The spline method reported here is third order accurate in space and second order accurate in time. The proposed spline method involves only two off-step points and a central point on a graded mesh. The method is two-level implicit in nature and directly derived from the continuity condition of the first order space derivative of the non-polynomial tension spline function. The linear stability analysis of the proposed method has been examined and it is shown that the proposed two-level method is unconditionally stable for a linear model problem. The method is directly applicable to problems in polar systems. To demonstrate the strength and utility of the proposed method, the authors have solved the generalized Burgers-Huxley equation, generalized Burgers-Fisher equation, coupled Burgers-equations and parabolic equation in polar coordinates. The authors show that the proposed method enables us to obtain the high accurate solution for high Reynolds number.

In this method, the authors use only two-level in time-direction, and at each time-level, the authors use three grid points for the unknown function u(x,t) and two off-step points for the known variable x in spatial direction. The methodology followed in this paper is the construction of a non-polynomial spline function and using its continuity properties to obtain consistency condition, which is third order accurate on a graded mesh and fourth order accurate on a uniform mesh. From this consistency condition, the authors derive the proposed numerical method. The proposed method, when applied to a linear equation is shown to be unconditionally stable. To assess the validity and accuracy, the method is applied to solve several benchmark problems, and numerical results are provided to demonstrate the usefulness of the proposed method.

The paper provides a third order numerical scheme on a graded mesh and fourth order spline method on a uniform mesh obtained directly from the consistency condition. In earlier methods, consistency conditions were only second order accurate. This brings an edge over other past methods. Also, the method is directly applicable to physical problems involving singular coefficients. So no modification in the method is required at singular points. This saves CPU time and computational costs.

There are no limitations. Obtaining a high accuracy spline method directly from the consistency condition is a new work. Also being an implicit method, this method is unconditionally stable.

Physical problems with singular and non-singular coefficients are directly solved by this method.

The paper develops a new method based on non-polynomial spline approximations of order two in time and three (four) in space, which is original and has lot of value because many benchmark problems of physical significance are solved in this method.

The purpose of this study is to propose a non-classical method to obtain efficient and accurate numerical solutions of the advection–diffusion–reaction equations.

Unlike conventional numerical methods, this study proposes a numerical scheme using outer Newton iteration applied to a time-dependent PDE. The linearized time dependent PDE is discretized by trapezoidal rule, which is second order in time, and by spline-based finite difference method of fourth order in space.

Using the proposed technique, even when relatively large time step sizes are used in computations, the efficiency of the proposed procedure is very clear for the numerical examples in comparison with the existing classical methods.

This study, unlike these classical methods, proposes an alternative approach based on linearizing the nonlinear problem at first, and then discretizing it by an appropriate scheme. This technique helps to avoid considering the convergence issues of Newton iteration applied to nonlinear algebraic system containing many unknowns at each time step if an implicit method is used in time discretization. The linearized PDE can be solved by implicit time integrator, which enables the use of large time step size.

The purpose of the paper is to extend the differential quadrature method (DQM) for solving time and space fractional non-linear partial differential equations on a semi-infinite domain.

The proposed method is the combination of the Legendre polynomials and differential quadrature method. The authors derived and constructed the new operational matrices for the fractional derivatives, which are used for the solutions of non-linear time and space fractional partial differential equations.

The fractional derivative of Lagrange polynomial is a big hurdle in classical DQM. To overcome this problem, the authors represent the Lagrange polynomial in terms of shifted Legendre polynomial. They construct a transformation matrix which transforms the Lagrange polynomial into shifted Legendre polynomial of arbitrary order. Then, they obtain the new weighting coefficients matrices for space fractional derivatives by shifted Legendre polynomials and use these in conversion of a non-linear fractional partial differential equation into a system of fractional ordinary differential equations. Convergence analysis for the proposed method is also discussed.

Many engineers can use the presented method for solving their time and space fractional non-linear partial differential equation models. To the best of the authors’ knowledge, the differential quadrature method has never been extended or implemented for non-linear time and space fractional partial differential equations.

The purpose of this paper is to investigate the possibility of extension to the variational iteration and the Adomian decomposition methods for solving nonlinear Huxley equation with time-fractional derivative.

Objectives achieved the main methods: the fractional derivative of f (x) in the Caputo sense is first stated. Second, the time-fractional Huxley equation is written in a differential operator form where the differential operator is in Caputo sense. After acting on both sides by the inverse operator of the fractional differential operator in Caputo sense, the Adomian's decomposition is then used to get the power series solution of the resulted time-fractional Huxley equation. Also, a second objective is achieved by applying the variational iteration method to get approximate solutions for the time-fractional Huxley equation.

There are some important findings to state and summarize here. First, the variational iteration method and the decomposition method provide the solutions in terms of convergent series with easily computable components for this considered problem. Second, it seems that the approximate solution of time-fractional Huxley equation using the decomposition method converges faster than the approximate solution using the variational iteration method. Third, the variational iteration method handles nonlinear equations without any need for the so-called Adomian polynomials. However, Adomian decomposition method provides the components of the exact solution, where these components should follow the summation given in Equation (21).

This paper presents new materials in terms of employing the variational iteration and the Adomian decomposition methods to solve the problem under consideration. It is expected that the results will give some insightful conclusions for the used techniques to handle similar fractional differential equations. This emphasizes the fact that the two methods are applicable to a broad class of nonlinear problems in fractional differential equations.

The purpose of this paper is to demonstrate how numerical solutions of the nonlinear Huxley equation are obtained by collocation‐based method using cubic B‐spline over finite elements.

For the numerical procedure, time derivative is discretized using usual finite difference scheme. Solution and its principal derivatives over the subintervals are approximated by the combination of the cubic B‐spline and unknown element parameters.

The numerical results are found to be in good agreement with the exact solution. Also the method is very accurate and conditionally stable; the results are very accurate at a small h (discretization) of x so this method can be applied for any nonlinear partial differential equations.

The paper demonstrates how numerical solutions of the nonlinear Huxley equation are obtained by collocation‐based method using cubic B‐spline over finite elements.

The purpose of this paper is to investigate the nano boundary‐layer flows over stretching surfaces with Navier boundary condition. This problem is mapped into the ordinary differential equation by presented similarity transformation. The resulting nonlinear ordinary differential equation is solved analytically by applying a newly developed method. The authors consider two types of flows: viscous flows over a two‐dimensional stretching surface; and viscous flows over an axisymmetric stretching surface.

The governing equation is solved analytically by applying a newly developed method, namely the differential transform method (DTM)‐Padé technique that is a combination of the DTM and the Padé approximation. The analytic solutions of the nonlinear ordinary differential equation are constructed in the ratio of two polynomials.

Graphical results are presented to investigate influence of the slip parameter and the suction parameter on the normal velocity and on the lateral velocity. The obtained solutions, in comparison with the numerical solutions, demonstrate remarkable accuracy. It is predicted that the DTM‐Padé can have wide application in engineering problems especially for boundary‐layer problems.

The resulting nonlinear ordinary differential equation is solved analytically by applying a newly developed method, namely the DTM‐Padé technique that is a combination of the DTM and the Padé approximation. The analytic solutions of the nonlinear ordinary differential equation are constructed in the ratio of two polynomials.

The main aim of this paper is to investigate the fractional coupled nonlinear Helmholtz equation by two new analytical methods.

This article takes an inaugural look at the fractional coupled nonlinear Helmholtz equation by using the conformable derivative. It successfully finds new fractional periodic solutions and solitary wave solutions by employing methods such as the fractional method and the fractional simple equation method. The dynamics of these fractional periodic solutions and solitary wave solutions are then graphically represented in 3D with appropriate parameters and fractal dimensions. This research contributes to a deeper comprehension and detailed exploration of the dynamics involved in high dimensional solitary wave propagation.

The proposed two mathematical approaches are simple and efficient to solve fractional evolution equations.

The fractional coupled nonlinear Helmholtz equation is described by using the conformable derivative for the first time. The obtained fractional periodic solutions and solitary wave solutions are completely new.

The purpose of this paper is to consider the thermal‐diffusion and diffusion‐thermo effects on combined heat and mass transfer of a steady magnetohydrodynamic (MHD) convective and slip flow due to a rotating disk with viscous dissipation and Ohmic heating. The main goal of the present study is to find the approximate analytic solutions by the combination of the DTM and the Padé approximants for this problem.

A new method, namely the DTM‐Padé technique, which is a combination of the differential transform method and the Padé approximation, is employed.

Graphical results for fluids of medium molecular weight (H2, air) are presented to investigate influence of the slip parameter, magnetic field parameter M, Eckert Ec, Schmidt Ec, Dufour Du and Soret Sr numbers on the profiles of the dimensionless velocity, temperature and concentration distributions. In order to show the effectiveness of the DTM‐Padé, the results obtained from the DTM‐Padé are compared with available solutions obtained using shooting method to generate the numerical solution.

This technique (DTM‐Padé) is extended to give solutions for nonlinear differential equations with boundary conditions at the infinity.

The 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).

The 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.

The 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.

To 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.