Search results

1 – 10 of over 1000
Article
Publication date: 14 September 2012

R.C. Mittal and Ram Jiwari

The purpose of this paper is to use the polynomial differential quadrature method (PDQM) to find the numerical solutions of some Burgers'‐type nonlinear partial differential…

Abstract

Purpose

The purpose of this paper is to use the polynomial differential quadrature method (PDQM) to find the numerical solutions of some Burgers'‐type nonlinear partial differential equations.

Design/methodology/approach

The PDQM changed the nonlinear partial differential equations into a system of nonlinear ordinary differential equations (ODEs). The obtained system of ODEs is solved by Runge‐Kutta fourth order method.

Findings

Numerical results for the nonlinear evolution equations such as 1D Burgers', coupled Burgers', 2D Burgers' and system of 2D Burgers' equations are obtained by applying PDQM. The numerical results are found to be in good agreement with the exact solutions.

Originality/value

A comparison is made with those which are already available in the literature and the present numerical schemes are found give better solutions. The strong point of these schemes is that they are easy to apply, even in two‐dimensional nonlinear problems.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 22 no. 7
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 16 October 2009

Dursun Irk

The purpose of this paper is to investigate the numerical solutions of the Burgers' and modified Burgers' equation using sextic B‐spline collocation method.

Abstract

Purpose

The purpose of this paper is to investigate the numerical solutions of the Burgers' and modified Burgers' equation using sextic B‐spline collocation method.

Design/methodology/approach

Crank‐Nicolson central differencing scheme has been used for the time integration and sextic B‐spline functions have been used for the space integration to the modified and time splitted modified Burgers' equation.

Findings

It has been found that the proposed method is unconditionally stable and obtained results are consistent with some earlier published studies.

Originality/value

Sextic B‐spline collocation method for the Burgers' and modified Burgers' equation is given.

Details

Kybernetes, vol. 38 no. 9
Type: Research Article
ISSN: 0368-492X

Keywords

Article
Publication date: 9 April 2020

Ranjan Kumar Mohanty and Sachin Sharma

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…

Abstract

Purpose

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

Design/methodology/approach

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.

Findings

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.

Research limitations/implications

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.

Practical implications

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

Originality/value

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.

Article
Publication date: 23 March 2012

Najeeb Alam Khan, Asmat Ara and Amir Mahmood

The purpose of this paper is to use the generalized differential transform method (GDTM) and homotopy perturbation method (HPM) for solving time‐fractional Burgers and coupled…

Abstract

Purpose

The purpose of this paper is to use the generalized differential transform method (GDTM) and homotopy perturbation method (HPM) for solving time‐fractional Burgers and coupled Burgers equations. The fractional derivatives are described in the Caputo sense.

Design/methodology/approach

In these schemes, the solutions takes the form of a convergent series. In GDTM, the differential equation and related initial conditions are transformed into a recurrence relation that finally leads to the solution of a system of algebraic equations as coefficients of a power series solution. HPM requires a homotopy with an embedding parameter which is considered as a small parameter.

Findings

The paper extends the application and numerical comparison of the GDTM and HPM to obtain analytic and approximate solutions to the time‐fractional Burgers and coupled Burgers equations.

Research limitations/implications

Burgers and coupled Burgers equations with time‐fractional derivative used.

Practical implications

The implications include traffic flow, acoustic transmission, shocks, boundary layer, the steepening of the waves and fluids, thermal radiation, chemical reaction, gas dynamics and many other phenomena.

Originality/value

The numerical results demonstrate the significant features, efficiency and reliability of the two approaches. The results show that HPM is more promising, convenient, and computationally attractive than GDTM.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 22 no. 2
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 1 June 1994

Michael M. Grigor’ev

The paper gives the description of boundary element method(BEM) with subdomains for the solution ofconvection—diffusion equations with variable coefficients and Burgersequations

Abstract

The paper gives the description of boundary element method (BEM) with subdomains for the solution of convection—diffusion equations with variable coefficients and Burgersequations. At first, the whole domain is discretized into K subdomains, in which linearization of equations by representing convective velocity by the sum of constant and variable parts is carried out. Then using fundamental solutions for convection—diffusion linear equations for each subdomain the boundary integral equation (in which the part of the convective term with the constant convective velocity is not included into the pseudo‐body force) is formulated. Only part of the convective term with the variable velocity, which is, as a rule, more than one order less than convective velocity constant part contribution, is left as the pseudo‐source. On the one hand, this does not disturb the numerical BEM—algorithm stability and, on the other hand, this leads to significant improvement in the accuracy of solution. The global matrix, similar to the case of finite element method, has block band structure whereas its width depends only on the numeration order of nodes and subdomains. It is noted, that in comparison with the direct boundary element method the number of global matrix non‐zero elements is not proportional to the square of the number of nodes, but only to the total number of nodal points. This allows us to use the BEM for the solution of problems with very fine space discretization. The proposed BEM with subdomains technique has been used for the numerical solution of one‐dimensional linear steady‐state convective—diffusion problem with variable coefficients and one‐dimensional non‐linear Burgersequation for which exact analytical solutions are available. It made it possible to find out the BEM correctness according to both time and space. High precision of the numerical method is noted. The good point of the BEM is the high iteration convergence, which is disturbed neither by high Reynolds numbers nor by the presence of negative velocity zones.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 4 no. 6
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 1 January 2013

R.C. Mittal, Ram Jiwari and Kapil K. Sharma

The purpose of this paper is to propose a numerical method to solve time dependent Burgers' equation with appropriate initial and boundary conditions.

Abstract

Purpose

The purpose of this paper is to propose a numerical method to solve time dependent Burgers' equation with appropriate initial and boundary conditions.

Design/methodology/approach

The presence of the nonlinearity in the problem leads to severe difficulties in the solution approximation. In construction of the numerical scheme, quasilinearization is used to tackle the nonlinearity of the problem which is followed by semi discretization for spatial direction using differential quadrature method (DQM). Semi discretization of the problem leads to a system of first order initial value problems which are followed by fully discretization using RK4 scheme. The method is analyzed for stability and convergence.

Findings

The method is illustrated and compared with existing methods via numerical experiments and it is found that the proposed method gives better accuracy and is quite easy to implement.

Originality/value

The new scheme is developed by using some numerical schemes. The scheme is analyzed for stability and convergence. In support of predicted theory some test examples are solved using the presented method.

Article
Publication date: 29 March 2013

Alper Korkmaz and İdris Dağ

The purpose of this paper is to simulate numerical solutions of nonlinear Burgers' equation with two well‐known problems in order to verify the accuracy of the cubic B‐spline…

Abstract

Purpose

The purpose of this paper is to simulate numerical solutions of nonlinear Burgers' equation with two well‐known problems in order to verify the accuracy of the cubic B‐spline differential quadrature methods.

Design/methodology/approach

Cubic B‐spline differential quadrature methods have been used to discretize the Burgers' equation in space and the resultant ordinary equation system is integrated via Runge‐Kutta method of order four in time. Numerical results are compared with each other and some former results by calculating discrete root mean square and maximum error norms in each case. A matrix stability analysis is also performed by determining eigenvalues of the coefficient matrices numerically.

Findings

Numerical results show that differential quadrature methods based on cubic B‐splines generate acceptable solutions of nonlinear Burgers' equation. Constructing hybrid algorithms containing various basis to determine the weighting coefficients for higher order derivative approximations is also possible.

Originality/value

Nonlinear Burgers' equation is solved by cubic B‐spline differential quadrature methods.

Article
Publication date: 31 December 2020

J.I. Ramos

The purpose of this paper is to determine both analytically and numerically the kink solutions to a new one-dimensional, viscoelastic generalization of Burgersequation, which…

Abstract

Purpose

The purpose of this paper is to determine both analytically and numerically the kink solutions to a new one-dimensional, viscoelastic generalization of Burgersequation, which includes a non-linear constitutive law, and the number of kinks as functions of the non-linearity and relaxation parameters.

Design/methodology/approach

An analytical procedure and two explicit finite difference methods based on first-order accurate approximations to the first-order derivatives are used to determine the single- and double-kink solutions.

Findings

It is shown that only two parameters characterize the solution and that the existence of a shock wave requires that the (semi-positive) relaxation parameter be less than unity and the non-linearity parameter be less than two. It is also shown that negative values of the non-linearity parameter result in kinks with a single inflection point and strain and dissipation rates with a single relative minimum and a single, relative maximum, respectively. For non-linearity parameters between one and two, it is shown that the kink has three inflection points that merge into a single one as this parameter approaches one and that the strain and dissipation rates exhibit relative maxima and minima whose magnitudes decrease and increase as the relaxation and nonlinearity coefficients, respectively, are increased. It is also shown that the viscoelastic generalization of the Burgers equation presented here is related to an ϕ8−scalar field.

Originality/value

A new, one-dimensional, viscoelastic generalization of Burgersequation, which includes a non-linear constitutive law and relaxation is proposed, and its kink solutions are determined both analytically and numerically. The equation and its solutions are connected with scalar field theories and may be used to both studies the effects of the non-linearity and relaxation and assess the accuracy of numerical methods for first-order, non-linear partial differential equations.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 31 no. 9
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 13 August 2019

Mohamed El-Beltagy

The paper aims to compare and clarify the differences and between the two well-known decomposition spectral techniques; the Winer–Chaos expansion (WCE) and the Winer–Hermite…

Abstract

Purpose

The paper aims to compare and clarify the differences and between the two well-known decomposition spectral techniques; the Winer–Chaos expansion (WCE) and the Winer–Hermite expansion (WHE). The details of the two decompositions are outlined. The difficulties arise when using the two techniques are also mentioned along with the convergence orders. The reader can also find a collection of references to understand the two decompositions with their origins. The geometrical Brownian motion is considered as an example for an important process with exact solution for the sake of comparison. The two decompositions are found practical in analysing the SDEs. The WCE is, in general, simpler, while WHE is more efficient as it is the limit of WCE when using infinite number of random variables. The Burgers turbulence is considered as a nonlinear example and WHE is shown to be more efficient in detecting the turbulence. In general, WHE is more efficient especially in case of nonlinear and/or non-Gaussian processes.

Design/methodology/approach

The paper outlined the technical and literature review of the WCE and WHE techniques. Linear and nonlinear processes are compared to outline the comparison along with the convergence of both techniques.

Findings

The paper shows that both decompositions are practical in solving the stochastic differential equations. The WCE is found simpler and WHE is the limit when using infinite number of random variables in WCE. The WHE is more efficient especially in case of nonlinear problems.

Research limitations/implications

Applicable for SDEs with square integrable processes and coefficients satisfying Lipschitz conditions.

Originality/value

This paper fulfils a comparison required by the researchers in the stochastic analysis area. It also introduces a simple efficient technique to model the flow turbulence in the physical domain.

Details

Engineering Computations, vol. 36 no. 7
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 6 November 2017

Fahimeh Saberi Zafarghandi, Maryam Mohammadi, Esmail Babolian and Shahnam Javadi

The purpose of this paper is to introduce a local Newton basis functions collocation method for solving the 2D nonlinear coupled Burgersequations. It needs less computer storage…

Abstract

Purpose

The purpose of this paper is to introduce a local Newton basis functions collocation method for solving the 2D nonlinear coupled Burgersequations. It needs less computer storage and flops than the usual global radial basis functions collocation method and also stabilizes the numerical solutions of the convection-dominated equations by using the Newton basis functions.

Design/methodology/approach

A meshless method based on spatial trial space spanned by the local Newton basis functions in the “native” Hilbert space of the reproducing kernel is presented. With the selected local sub-clusters of domain nodes, an approximation function is introduced as a sum of weighted local Newton basis functions. Then the collocation approach is used to determine weights. The method leads to a system of ordinary differential equations (ODEs) for the time-dependent partial differential equations (PDEs).

Findings

The method is successfully used for solving the 2D nonlinear coupled Burgersequations for reasonably high values of Reynolds number (Re). It is a well-known issue in the analysis of the convection-diffusion problems that the solution becomes oscillatory when the problem becomes convection-dominated if the standard methods are followed without special treatments. In the proposed method, the authors do not detect any instability near the front, hence no technique is needed. The numerical results show that the proposed method is efficient, accurate and stable for flow with reasonably high values of Re.

Originality/value

The authors used more stable basis functions than the standard basis of translated kernels for representing of kernel-based approximants for the numerical solution of partial differential equations (PDEs). The local character of the method, having a well-structured implementation including enforcing the Dirichlet and Neuman boundary conditions, and producing accurate and stable results for flow with reasonably high values of Re for the numerical solution of the 2D nonlinear coupled Burgersequations without any special technique are the main values of the paper.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 27 no. 11
Type: Research Article
ISSN: 0961-5539

Keywords

1 – 10 of over 1000