A reduced-order extrapolated approach to solution coefficient vectors in the Crank-Nicolson finite element method for the uniform transmission line equation

This work is concerned mainly with developing and testing the reduced-order extrapolated approach to the unknown coefficient vectors in the Crank-Nicolson finite element (CNFE) solutions for the uniform transmission line equation. For this objective, the CNFE functional and matrix models and the existence, stability, and errors of the CNFE solutions of the uniform transmission line equation are first derived. Then a reduced-order extrapolated CNFE (ROECNFE) matrix model is established by means of a proper orthogonal decomposition technique, and the existence, stability, and error estimates of the ROECNFE solutions are demonstrated by matrix analysis, leading to an elegant theoretical development. Especially, our work shows that the basis functions and accuracy of the ROECNFE matrix model are the same as those of the CNFE matrix model. Finally, some numerical tests are illustrated to computationally experimentally confirm the validity and sharpness of the ROECNFE method.     

A fast numerical method for harmonic equation based on natural boundary integral

Based on the coupling of the natural boundary integral method and the finite elements method, we mainly investigate the numerical solution of Neumann problem of harmonic equation in an exterior elliptic. Using our trigonometric wavelets and Galerkin method, there obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. On the other hand, we prove that the numerical solution possesses exponential convergence rate. Especially, examples state that our method still has good accuracy for small j when the solution u ₀(θ) is almost singular.     

A Legendre spectral element method (SEM) based on the modified bases for solving neutral delay distributed‐order fractional damped diffusion‐wave equation

The main purpose of the current paper is to propose a new numerical scheme based on the spectral element procedure for simulating the neutral delay distributed‐order fractional damped diffusion‐wave equation. To this end, the temporal direction has been discretized by a finite difference formula with convergence order where 1<α<2. In the next, to obtain a full‐discrete scheme, we apply the spectral finite element method on the spatial direction. Furthermore, the unconditional stability of semidiscrete scheme and convergence order of full‐discrete scheme of new technique are discussed. Finally, 2 test problems have been considered to demonstrate the ability and efficiency of the proposed numerical technique.     

A method based on meshless approach for the numerical solution of the two‐space dimensional hyperbolic telegraph equation

In the current article, we investigate the RBF solution of second‐order two‐space dimensional linear hyperbolic telegraph equation. For this purpose, we use a combination of boundary knot method (BKM) and analog equation method (AEM). The BKM is a meshfree, boundary‐only and integration‐free technique. The BKM is an alternative to the method of fundamental solution to avoid the fictitious boundary and to deal with low accuracy, singular integration and mesh generation. Also, on the basis of the AEM, the governing operator is substituted by an equivalent nonhomogeneous linear one with known fundamental solution under the same boundary conditions. Finally, several numerical results and discussions are demonstrated to show the accuracy and efficiency of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

A Neumann–Neumann algorithm for a mortar finite element discretization of fourth‐order elliptic problems in 2D

Here we present and analyze a Neumann–Neumann algorithm for the mortar finite element discretization of elliptic fourth‐order problems with discontinuous coefficients. The fully parallel algorithm is analyzed using the abstract Schwarz framework, proving a convergence which is independent of the parameters of the problem, and depends only logarithmically on the ratio between the subdomain size and the mesh size.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

An iterative method based on equation decomposition for the fourth‐order singular perturbation problem

In this article, we propose an iterative method based on the equation decomposition technique ( 1 ) for the numerical solution of a singular perturbation problem of fourth‐order elliptic equation. At each step of the given method, we only need to solve a boundary value problem of second‐order elliptic equation and a second‐order singular perturbation problem. We prove that our approximate solution converges to the exact solution when the domain is a disc. Our numerical examples show the efficiency and accuracy of our method. Our iterative method works very well for singular perturbation problems, that is, the case of 0 < ε ? 1, and the convergence rate is very fast. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A POD‐based–optimized finite difference CN‐extrapolated implicit scheme for the 2D viscoelastic wave equation†

In this paper, we devote ourselves to establishing the unconditionally stable and absolutely convergent classical finite difference Crank‐Nicholson (CN) implicit (CFDCNI) scheme and optimized finite difference CN‐extrapolated implicit (OFDCNEI) scheme containing very few degrees of freedom but holding fully second‐order accuracy for the two‐dimensional viscoelastic wave via the proper orthogonal decomposition technique, analyzing the existence, stability, and convergence of the CFDCNI and OFDCNEI solutions, and using the numerical simulations to verify that the OFDCNEI scheme is far more superior than the CFDCNI scheme.     

Application of collocation method for solving a parabolic‐hyperbolic free boundary problem which models the growth of tumor with drug application

In this article, we want to solve a free boundary problem which models tumor growth with drug application. This problem includes five time dependent partial differential equations. The tumor considered in this model consists of three kinds of cells, proliferative cells, quiescent cells, and dead cells. Three different first‐order hyperbolic equations are given that describe the evolution of cells and other two second‐order parabolic equations describe the diffusion of nutrient and drug concentration. We solve the problem using the collocation method. Then, we prove stability and convergence of method. Also, some examples are considered to show the efficiency of method. Copyright © 2016 John Wiley & Sons, Ltd.     

A numerical method based on extended Raviart–Thomas (ER‐T) mixed finite element method for solving damped Boussinesq equation

In this paper, we present the approximate solution of damped Boussinesq equation using extended Raviart–Thomas mixed finite element method. In this method, the numerical solution of this equation is obtained using triangular meshes. Also, for discretization in time direction, we use an implicit finite difference scheme. In addition, error estimation and stability analysis of both methods are shown. Finally, some numerical examples are considered to confirm the theoretical results. Copyright © 2017 John Wiley & Sons, Ltd.     

非饱和土壤水流问题基于POD 方法的降阶有限体积元格式及外推算法实现

利用特征投影分解(proper orthogonal decomposition, 简记为POD) 方法对非饱和土壤水流问题的经典有限体积元格式做降阶处理, 建立一种具有足够高精度维数较低的降阶有限体积元格式, 并给出这种降阶有限体积元解的误差估计和外推算法的实现, 最后用数值例子说明数值结果与理论结果是相吻合的. 进一步表明了基于POD 方法的降阶有限体积元格式对求解非饱和土壤水流问题数值解是可靠和有效的.     

Bernstein Ritz‐Galerkin method for solving an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

In this article, the Ritz‐Galerkin method in Bernstein polynomial basis is implemented to give an approximate solution of a hyperbolic partial differential equation with an integral condition. We will deal here with a type of nonlocal boundary value problem, that is, the solution of a hyperbolic partial differential equation with a nonlocal boundary specification. The nonlocal conditions arise mainly when the data on the boundary cannot be measured directly. The properties of Bernstein polynomial and Ritz‐Galerkin method are first presented, then Ritz‐Galerkin method is used to reduce the given hyperbolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique presented in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A priori error estimates of an extrapolated space‐time discontinuous galerkin method for nonlinear convection‐diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection‐diffusion equation. We employ a combination of the discontinuous Galerkin finite element (DGFE) method for the space as well as time discretization. The linear diffusive and penalty terms are treated implicitly whereas the nonlinear convective term is treated by a special higher order explicit extrapolation from the previous time step, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the L^∞(L²) –norm and the L²(H¹) –seminorm with respect to the mesh size h and time step τ. Finally, we present an efficient solution strategy and numerical examples verifying the theoretical results. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1456–1482, 2010     

A spectral element method using the modal basis and its application in solving second‐order nonlinear partial differential equations

We present a high‐order spectral element method (SEM) using modal (or hierarchical) basis for modeling of some nonlinear second‐order partial differential equations in two‐dimensional spatial space. The discretization is based on the conforming spectral element technique in space and the semi‐implicit or the explicit finite difference formula in time. Unlike the nodal SEM, which is based on the Lagrange polynomials associated with the Gauss–Lobatto–Legendre or Chebyshev quadrature nodes, the Lobatto polynomials are used in this paper as modal basis. Using modal bases due to their orthogonal properties enables us to exactly obtain the elemental matrices provided that the element‐wise mapping has the constant Jacobian. The difficulty of implementation of modal approximations for nonlinear problems is treated in this paper by expanding the nonlinear terms in the weak form of differential equations in terms of the Lobatto polynomials on each element using the fast Fourier transform (FFT). Utilization of the Fourier interpolation on equidistant points in the FFT algorithm and the enough polynomial order of approximation of the nonlinear terms can lead to minimize the aliasing error. Also, this approach leads to finding numerical solution of a nonlinear differential equation through solving a system of linear algebraic equations. Numerical results for some famous nonlinear equations illustrate efficiency, stability and convergence properties of the approximation scheme, which is exponential in space and up to third‐order in time. Copyright © 2014 John Wiley & Sons, Ltd.     

A second‐order finite difference method for fractional diffusion equation with Dirichlet and fractional boundary conditions

In this paper, we develop a practical numerical method to approximate a fractional diffusion equation with Dirichlet and fractional boundary conditions. An approach based on the classical Crank–Nicolson method combined with spatial extrapolation is used to obtain temporally and spatially second‐order accurate numerical estimates. The solvability, stability, and convergence of the proposed numerical scheme are proved via the Gershgorin theorem. Numerical experiments are performed to confirm the accuracy and efficiency of our scheme.     

On the solution of an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

Numerical solution of hyperbolic partial differential equation with an integral condition continues to be a major research area with widespread applications in modern physics and technology. Many physical phenomena are modeled by nonclassical hyperbolic boundary value problems with nonlocal boundary conditions. In place of the classical specification of boundary data, we impose a nonlocal boundary condition. Partial differential equations with nonlocal boundary specifications have received much attention in last 20 years. However, most of the articles were directed to the second‐order parabolic equation, particularly to heat conduction equation. We will deal here with new type of nonlocal boundary value problem that is the solution of hyperbolic partial differential equations with nonlocal boundary specifications. These nonlocal conditions arise mainly when the data on the boundary can not be measured directly. Several finite difference methods have been proposed for the numerical solution of this one‐dimensional nonclassic boundary value problem. These computational techniques are compared using the largest error terms in the resulting modified equivalent partial differential equation. Numerical results supporting theoretical expectations are given. Restrictions on using higher order computational techniques for the studied problem are discussed. Suitable references on various physical applications and the theoretical aspects of solutions are introduced at the end of this article. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A finite element variational multiscale method for steady‐state natural convection problem based on two local gauss integrations

In this article, supposing that the velocity, pressure, and temperature are approximated by the elements , and applying the orthogonal projection technique, we introduce two Gauss integrations as a stabilizing term in the common variational multiscale (VMS) method and derive a new VMS (Two Gauss VMS) method for steady‐state natural convection problem. Comparing with the common VMS method, the Two Gauss VMS method does not need to introduce any extra variable and reduces the degrees of freedom of the discrete system a lot, but gets the same stabilized result. The effectiveness and stability of the Two Gauss VMS method are further demonstrated through two numerical examples. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 361–375, 2014     

A computational study of the one‐dimensional parabolic equation subject to nonclassical boundary specifications

Certain problems arising in engineering are modeled by nonstandard parabolic initial‐boundary value problems in one space variable, which involve an integral term over the spatial domain of a function of the desired solution. Hence, in the past few years interest has substantially increased in the solutions of these problems. As a result numerous research papers have also been devoted to the subject. Although considerable amount of work has been done in the past, there is still a lack of a completely satisfactory computational scheme. Also, there are some cases that have not been studied numerically yet. In the current article several approaches for the numerical solution of the one‐dimensional parabolic equation subject to the specification of mass, which have been considered in the literature, are reported. Finite difference methods have been proposed for the numerical solution of the new nonclassic boundary value problem. To investigate the performance of the proposed algorithm, we consider solving a test problem. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A steady barycentric Lagrange interpolation method for the 2D higher‐order time‐fractional telegraph equation with nonlocal boundary condition with error analysis

In this paper, we consider the numerical solution of the time‐fractional telegraph equation with a nonlocal boundary condition. A novel barycentric Lagrange interpolation collocation method is developed to solve this equation. Two difficulties have been sorted: the singularity of the integration and the higher accuracy. At the same, we put forward a steady barycentric Lagrange interpolation technique to overcome the new “Runge” phenomenon in computation. Error estimates of the barycentric Lagrange interpolation and the time‐fractional telegraph system for the present method are presented in Sobolev spaces. High convergence rates of the proposed method are obtained and are consisted with the numerical values. Especially in the time dimension, we get the error bound, for h‐refinement and for n_t‐density in the L₂ norms. The numerical results obtained show that the proposed numerical algorithm is accurate and computationally efficient for solving time‐fractional telegraph equation. Experiments demonstrate the high convergence rates of the proposed method are consisted with the theoretical values.     

A mixed hp FEM for the approximation of fourth‐order singularly perturbed problems on smooth domains

We consider fourth‐order singularly perturbed problems posed on smooth domains and the approximation of their solution by a mixed Finite Element Method on the so‐called Spectral Boundary Layer Mesh. We show that the method converges uniformly, with respect to the singular perturbation parameter, at an exponential rate when the error is measured in the energy norm. Numerical examples illustrate our theoretical findings.     

A backward euler orthogonal spline collocation method for the time‐fractional Fokker–Planck equation

We formulate and analyze a novel numerical method for solving a time‐fractional Fokker–Planck equation which models an anomalous subdiffusion process. In this method, orthogonal spline collocation is used for the spatial discretization and the time‐stepping is done using a backward Euler method based on the L1 approximation to the Caputo derivative. The stability and convergence of the method are considered, and the theoretical results are supported by numerical examples, which also exhibit superconvergence. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1534–1550, 2015