Vorticity‐pressure formulations for the Brinkman‐Darcy coupled problem

We introduce a new variational formulation for the Brinkman‐Darcy equations formulated in terms of the scaled Brinkman vorticity and the global pressure. The velocities in each subdomain are fully decoupled through the momentum equations, and can be later recovered from the principal unknowns. A new finite element method is also proposed, consisting in equal‐order Nédélec and piecewise continuous elements, for vorticity and pressure, respectively. The error analysis for the scheme is carried out in the natural norms, with bounds independent of the fluid viscosity. An adequate modification of the formulation and analysis permits us to specify the presentation to the case of axisymmetric configurations. We provide a set of numerical examples illustrating the robustness, accuracy, and efficiency of the proposed discretization.     

A space‐time spectral method for one‐dimensional time fractional convection diffusion equations

In this paper, we propose a space‐time spectral method for solving a class of time fractional convection diffusion equations. Because both fractional derivative and spectral method have global characteristics in bounded domains, we propose a space‐time spectral‐Galerkin method. The convergence result of the method is proved by providing a priori error estimate. Numerical results further confirm the expected convergence rate and illustrate the versatility of our method. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order numerical solution of second‐order one‐dimensional hyperbolic telegraph equation using a shifted Gegenbauer pseudospectral method

We present a high‐order shifted Gegenbauer pseudospectral method (SGPM) to solve numerically the second‐order one‐dimensional hyperbolic telegraph equation provided with some initial and Dirichlet boundary conditions. The framework of the numerical scheme involves the recast of the problem into its integral formulation followed by its discretization into a system of well‐conditioned linear algebraic equations. The integral operators are numerically approximated using some novel shifted Gegenbauer operational matrices of integration. We derive the error formula of the associated numerical quadratures. We also present a method to optimize the constructed operational matrix of integration by minimizing the associated quadrature error in some optimality sense. We study the error bounds and convergence of the optimal shifted Gegenbauer operational matrix of integration. Moreover, we construct the relation between the operational matrices of integration of the shifted Gegenbauer polynomials and standard Gegenbauer polynomials. We derive the global collocation matrix of the SGPM, and construct an efficient computational algorithm for the solution of the collocation equations. We present a study on the computational cost of the developed computational algorithm, and a rigorous convergence and error analysis of the introduced method. Four numerical test examples have been carried out to verify the effectiveness, the accuracy, and the exponential convergence of the method. The SGPM is a robust technique, which can be extended to solve a wide range of problems arising in numerous applications. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 307–349, 2016     

A local radial basis function collocation method to solve the variable‐order time fractional diffusion equation in a two‐dimensional irregular domain

The local radial basis function (RBF) method is a promising solver for variable‐order time fractional diffusion equation (TFDE), as it overcomes the computational burden of the traditional global method. Application of the local RBF method is limited to Fickian diffusion, while real‐world diffusion is usually non‐Fickian in multiple dimensions. This article is the first to extend the application of the local RBF method to two‐dimensional, variable‐order, time fractional diffusion equation in complex shaped domains. One of the main advantages of the local RBF method is that only the nodes located in the subdomain, surrounding the local point, need to be considered when calculating the numerical solution at this point. This approach can perform well with large scale problems and can also mitigate otherwise ill‐conditioned problems. The proposed numerical approach is checked against two examples with curved boundaries and known analytical solutions. Shape parameter and subdomain node number are investigated for their influence on the accuracy of the local RBF solution. Furthermore, quantitative analysis, based on root‐mean‐square error, maximum absolute error, and maximum error of the partial derivative indicates that the local RBF method is accurate and effective in approximating the variable‐order TFDE in two‐dimensional irregular domains.     

An adaptive finite element method for a time‐dependent Stokes problem

In this article, we conduct an a posteriori error analysis of the two‐dimensional time‐dependent Stokes problem with homogeneous Dirichlet boundary conditions, which can be extended to mixed boundary conditions. We present a full time–space discretization using the discontinuous Galerkin method with polynomials of any degree in time and the ? ₂ ? ?₁ Taylor–Hood finite elements in space, and propose an a posteriori residual‐type error estimator. The upper bounds involve residuals, which are global in space and local in time, and an L ²‐error term evaluated on the left‐end point of time step. From the error estimate, we compute local error indicators to develop an adaptive space/time mesh refinement strategy. Numerical experiments verify our theoretical results and the proposed adaptive strategy.     

A rational approximate solution for generalized pantograph‐delay differential equations

In this paper, a new rational approximation based on a rational interpolation and collocation method is proposed for the solutions of generalized pantograph equations. A comprehensive error analysis is provided. The first part of the error analysis gives an upper bound for the absolute error. The second part is based on residual error procedure that estimates the absolute error. Some numerical examples are given to illustrate the method. The theoretical results support the numerical results. Copyright © 2015 John Wiley & Sons, Ltd.     

Methods for the numerical solution of the Benjamin‐Bona‐Mahony‐Burgers equation

L^∞‐error estimates for finite element for Galerkin solutions for the Benjamin‐Bona‐Mahony‐Burgers (BBMB) equation are considered. A priori bound and the semidiscrete Galerkin scheme are studied using appropriate projections. For fully discrete Galerkin schemes, we consider the backward Euler method and analyze the corresponding error estimates. For a second order accuracy in time, we propose a three‐level backward method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Convergence of finite element method for linear second‐order wave equations with discontinuous coefficients

A finite element method is proposed and analyzed for hyperbolic problems with discontinuous coefficients. The main emphasize is given on the convergence of such method. Due to low global regularity of the solutions, the error analysis of the standard finite element method is difficult to adopt for such problems. For a practical finite element discretization, optimal error estimates in L^∞(L²) and L^∞(H¹) norms are established for continuous time discretization. Further, a fully discrete scheme based on a symmetric difference approximation is considered, and optimal order convergence in L^∞(H¹) norm is established. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

One‐bit sigma‐delta quantization with exponential accuracy

One‐bit quantization is a method of representing bandlimited signals by ±1 sequences that are computed from regularly spaced samples of these signals; as the sampling density λ → ∞, convolving these one‐bit sequences with appropriately chosen filters produces increasingly close approximations of the original signals. This method is widely used for analog‐to‐digital and digital‐to‐analog conversion, because it is less expensive and simpler to implement than the more familiar critical sampling followed by fine‐resolution quantization. However, unlike fine‐resolution quantization, the accuracy of one‐bit quantization is not well‐understood. A natural error lower bound that decreases like 2^?λ can easily be given using information theoretic arguments. Yet, no one‐bit quantization algorithm was known with an error decay estimate even close to exponential decay. In this paper, we construct an infinite family of one‐bit sigma‐delta quantization schemes that achieves this goal. In particular, using this family, we prove that the error signal for π‐bandlimited signals is at most O(2^?.07λ). © 2003 Wiley Periodicals, Inc.     

Weak solution to a one‐dimensional full compressible non‐Newtonian fluid

This paper is devoted to the existence of global‐in‐time weak solutions to a one‐dimensional full compressible non‐Newtonian fluid. A semi‐discrete finite element scheme is taken to generate approximate solutions, based on an exact projection technique. To enforce convergence of the approximate solutions, the uniform estimate is obtained using an iteration method and energy method, with the help of the weak compactness and convexity. Numerical simulations showing the existence of solutions are presented.     

A posteriori error estimates of edge residual type of finite element method for nonmonotone quasi‐linear elliptic problems

In this article, we study the edge residual‐based a posteriori error estimates of conforming linear finite element method for nonmonotone quasi‐linear elliptic problems. It is proven that edge residuals dominate a posteriori error estimates. Up to higher order perturbations, edge residuals can act as a posteriori error estimators. The global reliability and local efficiency bounds are established both in H ¹‐norm and L ²‐norm. Numerical experiments are provided to illustrate the performance of the proposed error estimators. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 813–837, 2014     

Finite element method for two‐dimensional time‐fractional tricomi‐type equations

In this article, we consider the finite element method (FEM) for two‐dimensional linear time‐fractional Tricomi‐type equations, which is obtained from the standard two‐dimensional linear Tricomi‐type equation by replacing the first‐order time derivative with a fractional derivative (of order α, with 1 <α< 2 ). The method is based on finite element method for space and finite difference method for time. We prove that the method is unconditionally stable, and the error estimate is presented. The comparison of the FEM results with the exact solutions is made, and numerical experiments reveal that the FEM is very effective. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Synchronization of delayed coupled reaction‐diffusion systems on networks

In this paper, the exponential synchronization problem of delayed coupled reaction‐diffusion systems on networks (DCRDSNs) is investigated. Based on graph theory, a systematic method is designed to achieve exponential synchronization between two DCRDSNs by constructing a global Lyapunov function for error system. Two different kinds of sufficient synchronization criteria are derived in the form of Lyapunov functions and coefficients of drive‐response systems, respectively. Finally, a numerical example is given to show the usefulness of the proposed criteria. Copyright © 2014 John Wiley & Sons, Ltd.     

Nonconforming finite element methods with subgrid viscosity applied to advection‐diffusion‐reaction equations

A nonconforming (Crouzeix–Raviart) finite element method with subgrid viscosity is analyzed to approximate advection‐diffusion‐reaction equations. The error estimates are quasi‐optimal in the sense that keeping the Péclet number fixed, the estimates are suboptimal of order in the mesh size for the L²‐norm and optimal for the advective derivative on quasi‐uniform meshes. The method is also reformulated as a finite volume box scheme providing a reconstruction formula for the diffusive flux with local conservation properties. Numerical results are presented to illustrate the error analysis. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

H1‐Galerkin expanded mixed finite element methods for nonlinear pseudo‐parabolic integro‐differential equations

H¹‐Galerkin mixed finite element method combined with expanded mixed element method is discussed for nonlinear pseudo‐parabolic integro‐differential equations. We conduct theoretical analysis to study the existence and uniqueness of numerical solutions to the discrete scheme. A priori error estimates are derived for the unknown function, gradient function, and flux. Numerical example is presented to illustrate the effectiveness of the proposed scheme. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

First‐Order system least squares for generalized‐Newtonian coupled Stokes‐Darcy flow

The coupled problem for a generalized Newtonian Stokes flow in one domain and a generalized Newtonian Darcy flow in a porous medium is studied in this work. Both flows are treated as a first‐order system in a stress‐velocity formulation for the Stokes problem and a volumetric flux‐hydraulic potential formulation for the Darcy problem. The coupling along an interface is done using the well‐known Beavers–Joseph–Saffman interface condition. A least squares finite element method is used for the numerical approximation of the solution. It is shown that under some assumptions on the viscosity the error is bounded from above and below by the least squares functional. An adaptive refinement strategy is examined in several numerical examples where boundary singularities are present. Due to the nonlinearity of the problem a Gauss–Newton method is used to iteratively solve the problem. It is shown that the linear variational problems arising in the Gauss–Newton method are well posed. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1150–1173, 2015     

High accuracy analysis of the characteristic‐nonconforming FEM for a convection‐dominated transport problem

A low order characteristic‐nonconforming finite element method is proposed for solving a two‐dimensional convection‐dominated transport problem. On the basis of the distinguish property of element, that is, the consistency error can be estimated as order O(h²), one order higher than that of its interpolation error, the superclose result in broken energy norm is derived for the fully discrete scheme. In the process, we use the interpolation operator instead of the so‐called elliptic projection, which is an indispensable tool in the traditional finite element analysis. Furthermore, the global superconvergence is obtained by using the interpolated postprocessing technique. Lastly, some numerical experiments are provided to verify our theoretical analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

A high‐order Godunov method for one‐dimensional convection‐diffusion‐reaction problems

In this article we develop a high‐order Godunov method for one‐dimensional convection‐diffusion‐reaction problems where convection dominates diffusion. The heart of this method comes from incorporating the diffusion term via the slope of the linear representation (recovery) of the solution on each grid cell. The method is conservative and explicit. Therefore, it is efficient in computing time. For constant coefficient linear convection, diffusion, and Lipschitz‐type reaction, the properties of the total variation stability and monotonicity preservation are proved. An error estimation is derived. Computational examples are presented and compared with the exact solutions. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 495–512, 2000     

Error expansion for a first‐order upwind difference scheme

A singularly perturbed convection‐diffusion problem is considered. The problem is discretized using a simple firstorder upwinded difference scheme on general meshes. We derive an expansion of the error of the scheme that allows us to give error bounds in the discrete maximum norm, uniformly with respect to the perturbation parameter, for both a defect correction method and the Richardson extrapolation technique. This generalizes and simplifies results obtained in earlier publications [2], [8].     

Interpolants for sixth‐order Numerov‐type methods

The classical four‐stage family of explicit sixth‐order Numerov‐type method is considered. We provide two kinds of interpolants: (a) a three‐step interpolation based on all available data at mesh points and (b) a local interpolant (ie, two steps) that is constructed after solving scaled equations of condition. These latter equations are explained and provided here. Applying these interpolants in a set of tests, we conclude that they produce global errors of the same magnitude with the underlying method.