We study the dispersion and dissipation of the numerical scheme obtained by taking a weighted averaging of the consistent (finite element) mass matrix and lumped (spectral element) mass matrix for the small wavenumber limit. We find and prove that for the optimum blending the resulting scheme (a) provides 2p+4 order accuracy for p-th order method (two orders more accurate compared with finite a...

In this article, we discuss the numerical approximation of transport phenomena occurring at material interfaces between physical subdomains with heterogenous properties. The model in each subdomain consists of a partial differential equation with diffusive, convective and reactive terms, the coupling between each subdomain being realized through an interface transmission condition of Robin type...

In this paper, the semi-numerical techniques known as the optimal homotopy analysis method (HAM) and Differential Transform Method (DTM) are applied to study the magneto-hemodynamic laminar viscous flow of a conducting physiological fluid in a semi-porous channel under a transverse magnetic field. The two-dimensional momentum conservation partial differential equations are reduced to ordinary f...

A high-order accurate numerical quadrature algorithm is presented for the evaluation of integrals over curved surfaces and volumes which are defined implicitly via a fixed isosurface of a given function restricted to a given hyperrectangle. By converting the implicitly defined geometry into the graph of an implicitly defined height function, the approach leads to a recursive algorithm on the nu...

A finite element method with numerical quadrature is considered for the solution of a class of second-order quasilinear elliptic problems of nonmonotone type. Optimal a priori error estimates for the H and the L norms are derived. The uniqueness of the finite element solution is established for a sufficiently fine mesh. Application to numerical homogenization methods is considered.

Keywords: Operator-splitting method Finite element method Parabolic equations High-dimensional problems a b s t r a c t An operator-splitting finite element method for solving high-dimensional parabolic equations is presented. The stability and the error estimates are derived for the proposed numerical scheme. Furthermore, two variants of fully-practical operator-splitting finite element algori...

Numerical solutions of the generalized Burgers-Huxley equation are obtained using a polynomial differential quadrature method with minimal computational effort. To achieve this, a combination of a polynomial-based differential quadrature method in space and a low-storage third-order total variation diminishing Runge-Kutta scheme in time has been used. The computed results with the use of this t...

The RK5GL3 method is a numerical method for solving initial value problems in ordinary differential equations, and is based on a combination of a fifth-order Runge-Kutta method and 3-point Gauss-Legendre quadrature. In this paper we describe an effective local error control algorithm for RK5GL3, which uses local extrapolation with an eighth-order Runge-Kutta method in tandem with RK5GL3, and a ...