Discrete mass conservation for porous media saturated flow

Global and local mass conservation for velocity fields associated with saturated porous media flow have long been recognized as integral components of any numerical scheme attempting to simulate these flows. In this work, we study finite element discretizations for saturated porous media flow that use Taylor–Hood (TH) and Scott–Vogelius (SV) finite elements. The governing equations are modified to include a stabilization term when using the TH elements, and we provide a theoretical result that shows convergence (with respect to the stabilization parameter) to pointwise mass‐conservative solutions. We also provide results using the SV approximation pair. These elements are pointwise divergence free, leading to optimal convergence rates and numerical solutions. We give numerical results to verify our theory and a comparison with standard mixed methods for saturated flow problems. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 625–640, 2014     

The transport dynamics in complex systems governing by anomalous diffusion modelled with Riesz fractional partial differential equations

In this paper, numerical solutions of fractional Fokker–Planck equations with Riesz space fractional derivatives have been developed. Here, the fractional Fokker–Planck equations have been considered in a finite domain. In order to deal with the Riesz fractional derivative operator, shifted Grünwald approximation and fractional centred difference approaches have been used. The explicit finite difference method and Crank–Nicolson implicit method have been applied to obtain the numerical solutions of fractional diffusion equation and fractional Fokker–Planck equations, respectively. Numerical results are presented to demonstrate the accuracy and effectiveness of the proposed numerical solution techniques. Copyright © 2016 John Wiley & Sons, Ltd.     

Schwarz alternating domain decomposition approach for the solution of two‐dimensional Navier–Stokes flow problems by the method of approximate particular solutions

The method of approximate particular solutions (MAPS) is used to solve the two‐dimensional Navier–Stokes equations. This method uses particular solutions of a nonhomogeneous Stokes problem, with the multiquadric radial basis function as a nonhomogeneous term, to approximate the velocity and pressure fields. The continuity equation is not explicitly imposed since the used particular solutions are mass conservative. To improve the computational efficiency of the global MAPS, the domain is split into overlapped subdomains where the Schwarz Alternating Algorithm is employed using velocity or traction values from neighboring subdomains as boundary conditions. When imposing only velocity boundary conditions, an extra step is required to find a reference value for the pressure at each subdomain to guarantee continuity of pressure across subdomains. The Stokes lid‐driven cavity flow problem is solved to assess the performance of the Schwarz algorithm in comparison to a finite‐difference‐type localized MAPS. The Kovasznay flow problem is used to validate the proposed numerical scheme. Despite the use of relative coarse nodal distributions, numerical results show excellent agreement with respect to results reported in literature when solving the lid‐driven cavity (up to Re = 10,000) and the backward facing step (at Re = 800) problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 777–797, 2015     

An h–p version of the Chebyshev spectral collocation method for nonlinear delay differential equations

We develop and analyze a spectral collocation method based on the Chebyshev–Gauss–Lobatto points for nonlinear delay differential equations with vanishing delays. We derive an a priori error estimate in the H¹‐norm that is completely explicit with respect to the local time steps and the local polynomial degrees. Several numerical examples are provided to illustrate the theoretical results.     

Numerical solution of the Rosenau–KdV–RLW equation by operator splitting techniques based on B‐spline collocation method

In the present study, the operator splitting techniques based on the quintic B‐spline collocation finite element method are presented for calculating the numerical solutions of the Rosenau–KdV–RLW equation. Two test problems having exact solutions have been considered. To demonstrate the efficiency and accuracy of the present methods, the error norms L₂ and L_∞ with the discrete mass Q and energy E conservative properties have been calculated. The results obtained by the method have been compared with the exact solution of each problem and other numerical results in the literature, and also found to be in good agreement with each other. A Fourier stability analysis of each presented method is also investigated.     

Numerical and theoretical discussions for solving nonlinear generalized Benjamin–Bona–Mahony–Burgers equation based on the Legendre spectral element method

This article is devoted to solving numerically the nonlinear generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation that has several applications in physics and applied sciences. First, the time derivative is approximated by using a finite difference formula. Afterward, the stability and convergence analyses of the obtained time semi‐discrete are proven by applying the energy method. Also, it has been demonstrated that the convergence order in the temporal direction is O(dt) . Second, a fully discrete formula is acquired by approximating the spatial derivatives via Legendre spectral element method. This method uses Lagrange polynomial based on Gauss–Legendre–Lobatto points. An error estimation is also given in detail for full discretization scheme. Ultimately, the GBBMB equation in the one‐ and two‐dimension is solved by using the proposed method. Also, the calculated solutions are compared with theoretical solutions and results obtained from other techniques in the literature. The accuracy and efficiency of the mentioned procedure are revealed by numerical samples.     

An analytical and numerical study of the Stefan problem with convection by means of an enthalpy method

In this paper, we consider a theoretical and numerical study of the Stefan problem with convection, described by the Navier–Stokes equations with no‐slip boundary conditions. The mathematical formulation adopted is based on the enthalpy method. The existence of a weak solution is proved in the bidimensional case. The numerical effectiveness of the model considered is confirmed by some numerical results. Copyright © 2001 John Wiley & Sons, Ltd.     

A finite element approach for analysis and computational modelling of coupled reaction diffusion models

In this article, we establish the existence and uniqueness of solutions to the coupled reaction–diffusion models using Banach fixed point theorem. The Galerkin finite element method is used for the approximation of solutions, and an a priori error estimate is derived for such approximations. A scheme is proposed by combining the Crank–Nicolson and the predictor–corrector methods for the time discretization. Some numerical examples are considered to illustrate the accuracy and efficiency of the proposed scheme. It is found that the scheme is second‐order convergent. In addition, nonuniform grids are used in some cases to enhance the accuracy of the scheme.     

Nonstandard finite difference schemes for a class of generalized convection–diffusion–reaction equations

In this work, a class of nonstandard finite difference (NSFD) schemes are proposed to approximate the solutions of a class of generalized convection–diffusion–reaction equations. First, in the case of no diffusion, two exact finite difference schemes are presented using the method of characteristics. Based on these two exact schemes, a class of exact schemes are presented by introducing a parameter α. Second, since the forms of these exact schemes are so complicated that they are not convenient to use, a class of NSFD schemes are derived from the exact schemes using numerical approximations. It follows that, under certain conditions about denominator function of time‐step sizes, these NSFD schemes are elementary stable and the solutions are positive and bounded. Third, by means of the Mickens' technique of subequations, a new class of implicit NSFD schemes are constructed for the full convection–diffusion–reaction equations. It is shown that, under certain parameters set, these NSFD schemes are capable of preserving the non‐negativity and boundedness of the analytical solutions. Finally, some numerical simulations are provided to verify the validity of our analytical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1288–1309, 2015     

A characteristic domain decomposition and space–time local refinement method for first‐order linear hyperbolic equations with interfaces

We develop a characteristic‐based domain decomposition and space–time local refinement method for first‐order linear hyperbolic equations. The method naturally incorporates various physical and numerical interfaces into its formulation and generates accurate numerical solutions even if large time‐steps are used. The method fully utilizes the transient and strongly local behavior of the solutions of hyperbolic equations and provides solutions with significantly improved accuracy and efficiency. Several numerical experiments are presented to illustrate the performance of the method and for comparison with other domain decomposition and local refinement schemes. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 1–28, 1999     

Long-time behaviour for solutions of the Vlasov–Poisson–Fokker–Planck equation

We study the long-time behaviour of solutions of the Vlasov–Poisson–Fokker–Planck equation for initial data small enough and satisfying some suitable integrability conditions. Our analysis relies on the study of the linearized problems with bounded potentials decaying fast enough for large times. We obtain global bounds in time for the fundamental solutions of such problems and their derivatives. This allows to get sharp bounds for the decay of the difference between the solutions of the Vlasov–Poisson–Fokker–Planck equation and the solution of the free equation with the same initial data. Thanks to these bounds, we get an explicit form for the second term in the asymptotic expansion of the solutions for large times. © 1998 B. G. Teubner Stuttgart—John Wiley & Sons, Ltd.     

Mobile and immobile gaseous transport: Embedded analytical solutions to finite volume methods

The motivation is driven by deposition processes based on chemical vapor problems. The underlying model problem is based on coupled transport–reaction equations with mobile and immobile areas. We deal with systems of ordinary and partial differential equations. Such equation systems are delicate to solve and we introduce a novel solver method, that takes into account ways to solve analytically parts of the transport and reaction equations. The main idea is to embed the analytical and semianalytical solutions, which can then be explicitly given to standard numerical schemes of higher order. The numerical scheme is based on flux‐based characteristic methods, which is a finite volume method. Such a method is an attractive alternative to the standard numerical schemes, which fully discretize the full equations. We instead reduce the computational time while embedding fast computable analytical parts. Here, we can accelerate the solver process, with a priori explicitly given solutions. We will focus on the derivation of the analytical solutions for general and special solutions of the characteristic methods that are embedded into a finite volume method. In the numerical examples, we illustrate the higher‐order method for different benchmark problems. Finally, the method is verified with realistic results. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

The Crank–Nicolson–Galerkin finite element method for a nonlocal parabolic equation with moving boundaries

The aim of this article is to establish the convergence and error bounds for the fully discrete solutions of a class of nonlinear equations of reaction–diffusion nonlocal type with moving boundaries, using a linearized Crank–Nicolson–Galerkin finite element method with polynomial approximations of any degree. A coordinate transformation which fixes the boundaries is used. Some numerical tests to compare our Matlab code with some existing moving finite element methods are investigated. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1515–1533, 2015     

Chebyshev Wavelet collocation method for solving generalized Burgers–Huxley equation

In this paper, new and efficient numerical method, called as Chebyshev wavelet collocation method, is proposed for the solutions of generalized Burgers–Huxley equation. This method is based on the approximation by the truncated Chebyshev wavelet series. By using the Chebyshev collocation points, algebraic equation system has been obtained and solved. Approximate solutions of the generalized Burgers–Huxley equation are compared with exact solutions. These calculations demonstrate that the accuracy of the Chebyshev wavelet collocation solutions is quite high even in the case of a small number of grid points. Copyright © 2015 John Wiley & Sons, Ltd.     

The comparison of two reliable methods for accurate solution of time‐fractional Kaup–Kupershmidt equation arising in capillary gravity waves

In this paper, we compared two different methods, one numerical technique, viz Legendre multiwavelet method, and the other analytical technique, viz optimal homotopy asymptotic method (OHAM), for solving fractional‐order Kaup–Kupershmidt (KK) equation. Two‐dimensional Legendre multiwavelet expansion together with operational matrices of fractional integration and derivative of wavelet functions is used to compute the numerical solution of nonlinear time‐fractional KK equation. The approximate solutions of time fractional Kaup–Kupershmidt equation thus obtained by Legendre multiwavelet method are compared with the exact solutions as well as with OHAM. The present numerical scheme is quite simple, effective, and expedient for obtaining numerical solution of fractional KK equation in comparison to analytical approach of OHAM. Copyright © 2016 John Wiley & Sons, Ltd.     

Legendre–Gauss–Lobatto spectral collocation method for nonlinear delay differential equations

A Legendre–Gauss–Lobatto spectral collocation method is introduced for the numerical solutions of a class of nonlinear delay differential equations. An efficient algorithm is designed for the single‐step scheme and applied to the multiple‐domain case. As a theoretical result, we obtain a general convergence theorem for the single‐step case. Numerical results show that the suggested algorithm enjoys high‐order accuracy both in time and in the delayed argument and can be implemented in a robust and efficient manner. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical scheme for the solution of viscous Cahn–Hilliard equation

In this paper, we present a numerical scheme for the solution of viscous Cahn–Hilliard equation. The scheme is based on Adomian's decomposition approach and the solutions are calculated in the form of a convergent series with easily computable components. Some numerical examples are presented. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

A compact finite difference method for the solution of the generalized Burgers–Fisher equation

In this article, numerical solutions of the generalized Burgers–Fisher equation are obtained using a compact finite difference method with minimal computational effort. To verify this, a combination of a sixth‐order compact finite difference scheme in space and a low‐storage third‐order total variation diminishing Runge–Kutta scheme in time have been used. The computed results with the use of this technique have been compared with the exact solution to show the accuracy of it. The approximate solutions to the equation have been computed without transforming the equation and without using linearization. Comparisons indicate that there is a very good agreement between the numerical solutions and the exact solutions in terms of accuracy. The present method is seen to be a very good alternative to some existing techniques for realistic problems. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

A Stable Numerical Approach for Implicit Non-Linear Neutral Delay Differential Equations

In this paper we consider implicit non-linear neutral delay differential equations to derive efficient numerical schemes with good stability properties. The basic idea is to reformulate the original problem eliminating the dependence on the derivative of the solution in the past values. Our hypothesis on the original equation allow us to study the boundedness and asymptotic stability of the true and numerical solutions by the theory of stability with respect to the forcing term.     

Exact Asymptotics of the Density of the Transition Probability for Discontinuous Markov Processes

In this paper we consider some Kolmogorov–Feller equations with a small parameter h. We present a method for constructing the exact (exponential) asymptotics of the fundamental solution of these equations for finite time intervals uniformly with respect to h. This means that we construct an asymptotics of the density of the transition probability for discontinuous Markov processes. We justify the asymptotic solutions constructed. We also present an algorithm for constructing all terms of the asymptotics of the logarithmic limit (logarithmic asymptotics) of the fundamental solution as t → +0 uniformly with respect to h. We write formulas of the asymptotics of the logarithmic limit for some special cases as t → +0. The method presented in this paper also allows us to construct exact asymptotics of solutions of initial–boundary value problems that are of probability meaning.