A C0‐weak Galerkin finite element method for fourth‐order elliptic problems

This article proposes and analyzes a C⁰‐weak Galerkin (WG) finite element method for solving the biharmonic equation in two‐dimensional and three‐dimensional. The new WG method uses continuous piecewise‐polynomial approximations of degree for the unknown u and discontinuous piecewise‐polynomial approximations of degree k for the trace of on the interelement boundaries. Optimal error estimates are obtained in H², H¹, and L² norms. Numerical experiments illustrate and confirm the theoretical results. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1090–1104, 2016     

A Jacobi spectral Galerkin method for the integrated forms of fourth‐order elliptic differential equations

This article analyzes the solution of the integrated forms of fourth‐order elliptic differential equations on a rectilinear domain using a spectral Galerkin method. The spatial approximation is based on Jacobi polynomials P (x), with α, β ∈ (?1, ∞) and n the polynomial degree. For α = β, one recovers the ultraspherical polynomials (symmetric Jacobi polynomials) and for α = β = ?½, α = β = 0, the Chebyshev of the first and second kinds and Legendre polynomials respectively; and for the nonsymmetric Jacobi polynomials, the two important special cases α = ?β = ±½ (Chebyshev polynomials of the third and fourth kinds) are also recovered. The two‐dimensional version of the approximations is obtained by tensor products of the one‐dimensional bases. The various matrix systems resulting from these discretizations are carefully investigated, especially their condition number. An algebraic preconditioning yields a condition number of O(N), N being the polynomial degree of approximation, which is an improvement with respect to the well‐known condition number O(N⁸) of spectral methods for biharmonic elliptic operators. The numerical complexity of the solver is proportional to N^d+1 for a d‐dimensional problem. This operational count is the best one can achieve with a spectral method. The numerical results illustrate the theory and constitute a convincing argument for the feasibility of the method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

B-splines methods with redefined basis functions for solving fourth order parabolic partial differential equations

In this work, we discuss two methods for solving a fourth order parabolic partial differential equation. In Method-I, we decompose the given equation into a system of second order equations and solve them by using cubic B-spline method with redefined basis functions. In Method-II, the equation is solved directly by applying quintic B-spline method with redefined basis functions. Stability of these methods have been discussed. Both methods are unconditionally stable. These methods are tested on four examples. The computed results are compared wherever possible with those already available in literature. We have developed Method-I for fourth order non homogeneous parabolic partial differential equation from which we can obtain displacement and bending moment both simultaneously, while Method-II gives only displacement. The results show that the derived methods are easily implemented and approximate the exact solution very well.     

Minimum time‐step criteria for the Galerkin finite element methods applied to one‐dimensional parabolic partial differential equations

The finite element method has been well established for numerically solving parabolic partial differential equations (PDEs). Also it is well known that a too large time step should not be chosen in order to obtain a stable and accurate numerical solution. In this article, accuracy analysis shows that a too small time step should not be chosen either for some time‐stepping schemes. Otherwise, the accuracy of the numerical solution cannot be improved or can even be worsened in some cases. Furthermore, the so‐called minimum time step criteria are established for the Crank‐Nicolson scheme, the Galerkin‐time scheme, and the backward‐difference scheme used in the temporal discretization. For the forward‐difference scheme, no minimum time step exists as far as the accuracy is concerned. In the accuracy analysis, no specific initial and boundary conditions are invoked so that such established criteria can be applied to the parabolic PDEs subject to any initial and boundary conditions. These minimum time step criteria are verified in a series of numerical experiments for a one‐dimensional transient field problem with a known analytical solution. The minimum time step criteria developed in this study are useful for choosing appropriate time steps in numerical simulations of practical engineering problems. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Higher‐order discontinuous Galerkin method for pyramidal elements using orthogonal bases

We study finite elements of arbitrarily high‐order defined on pyramids for discontinuous Galerkin methods. We propose a new family of high‐order pyramidal finite elements using orthogonal basis functions which can be used in hybrid meshes including hexahedra, tetrahedra, wedges, and pyramids. We perform a comparison between these orthogonal functions and nodal functions for affine and non‐affine elements. Different strategies for the inversion of the mass matrix are also considered and discussed. Numerical experiments are conducted for the three dimensional Maxwell's equations. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Weak Galerkin finite element methods for linear parabolic integro‐differential equations

The semidiscrete and fully discrete weak Galerkin finite element schemes for the linear parabolic integro‐differential equations are proposed. Optimal order error estimates are established for the corresponding numerical approximations in both and norms. Numerical experiments illustrating the error behaviors are provided.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1357–1377, 2016     

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     

Oscillation criteria for fourth‐order differential equations with middle term

Oscillation criteria for self‐adjoint fourth‐order differential equations were established for various conditions on the coefficients r(x) > 0, q(x) and p(x). However, most of these results deal with the case when lim_{x → ∞}∫^x₁q(s) ds < +∞. In this note we give a new oscillation criterion in the case when this condition is not fulfilled, in particular when q(x)↗ + ∞ (even with exponential growth).     

Applications of cubic B‐splines collocation method for solving nonlinear inverse parabolic partial differential equations

In this article, we discuss a numerical method for solving some nonlinear inverse parabolic partial differential equations with Dirichlet's boundary conditions. The approach used, is based on collocation of cubic B‐splines over finite elements, so we have continuity of the dependent variable and its first two derivatives throughout the solution range. We apply cubic B‐splines for spatial variable and derivatives, which produce an ill‐posed system. We solve this system using the Tikhonov regularization method. The accuracy of the proposed method is demonstrated by applying it on two test problems. The figures and comparisons have been presented for clarity. Also the stability of this method has been discussed. The main advantage of the resulting scheme is that the algorithm is very simple, so it is very easy to implement. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 88–104, 2017     

Continuous Galerkin finite element methods for a forward‐backward heat equation

A space‐time finite element method is introduced to solve a model forward‐backward heat equation. The scheme uses the continuous Galerkin method for the time discretization. An error analysis for the method is presented. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 257–265, 1999     

A fourth‐order H1‐Galerkin mixed finite element method for Kuramoto–Sivashinsky equation

A H¹‐Galerkin mixed finite element method is applied to the Kuramoto–Sivashinsky equation by using a splitting technique, which results in a coupled system. The method described in this article may also be considered as a Petrov–Galerkin method with cubic spline space as trial space and piecewise linear space as test space, since the second derivative of a cubic spline is a linear spline. Optimal‐order error estimates are obtained without any restriction on the mesh for both semi‐discrete and fully discrete schemes. The advantage of this method over that presented in Manickam et al., Comput. Math. Appl. vol. 35(6) (1998) pp. 5–25; for the same problem is that the size (i.e., (n + 1) × (n + 1)) of each resulting linear system is less than half of the size of the linear system of the earlier method, where n is the number of subintervals in the partition. Further, there is a requirement of less regularity on exact solution in this method. The results are validated with numerical examples. Finally, instability behavior of the solution is numerically captured with this method.     

An efficient numerical approach to solve a class of variable‐order fractional integro‐partial differential equations

The main purpose of this work is to investigate an initial boundary value problem related to a suitable class of variable order fractional integro‐partial differential equations with a weakly singular kernel. To discretize the problem in the time direction, a finite difference method will be used. Then, the Sinc‐collocation approach combined with the double exponential transformation is employed to solve the problem in each time level. The proposed numerical algorithm is completely described and the convergence analysis of the numerical solution is presented. Finally, some illustrative examples are given to demonstrate the pertinent features of the proposed algorithm.     

Multiquadric quasi‐interpolation methods for solving partial differential algebraic equations

In this article, we propose two meshless collocation approaches for solving time dependent partial differential algebraic equations (PDAEs) in terms of the multiquadric quasi‐interpolation schemes. In presenting the process of the solution, the error is estimated. Furthermore, the comparisons on condition numbers of the collocation matrices using different methods and the sensitivity of the shape parameter c are given. With the use of the appropriate collocation points, the method for PDAEs with index‐2 is improved. The results show that the methods have some advantages over some known methods, such as the smaller condition numbers or more accurate solutions for PDAEs which has an modal index‐2 or an impulse solution with index‐2. Copyright © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 95–119, 2014     

Efficient numerical methods for boundary data and right‐hand side reconstructions in elliptic partial differential equations

In this article, we discuss the application of two important numerical methods, Ritz–Galerkin and Method of Fundamental Solutions (MFS), for solving some inverse problems, arising in the context of two‐dimensional elliptic equations. The main incentive for studying the considered problems is their wide applications in engineering fields. In the previous literature, the use of these methods, particularly MFS for right hand side reconstruction has been limited, partly due to stability concerns. We demonstrate that these diculties may be surmounted if the aforementioned methods are combined with techniques such as dual reciprocity method(DRM). Moreover, we incorporate some iterative regularization techniques. This fact is especially veried by taking into account the noisy data with boundary conditions. In addition, parts of this article are dedicated to the problem of boundary data approximation and the issue of numerical stability, ending with a general discussion on the advantages and disadvantages of various methods. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1995–2026, 2015     

Nonintrusive reduced‐order modeling of parametrized time‐dependent partial differential equations

We propose a nonintrusive reduced‐order modeling method based on the notion of space‐time‐parameter proper orthogonal decomposition (POD) for approximating the solution of nonlinear parametrized time‐dependent partial differential equations. A two‐level POD method is introduced for constructing spatial and temporal basis functions with special properties such that the reduced‐order model satisfies the boundary and initial conditions by construction. A radial basis function approximation method is used to estimate the undetermined coefficients in the reduced‐order model without resorting to Galerkin projection. This nonintrusive approach enables the application of our approach to general problems with complicated nonlinearity terms. Numerical studies are presented for the parametrized Burgers' equation and a parametrized convection‐reaction‐diffusion problem. We demonstrate that our approach leads to reduced‐order models that accurately capture the behavior of the field variables as a function of the spatial coordinates, the parameter vector and time. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Exact solutions for nonlinear partial differential equations via Exp‐function method

In this article, the Exp‐function method is applied to nonlinear Burgers equation and special fifth‐order partial differential equation. Using this method, we obtain exact solutions for these equations. The method is straightforward and concise, and its applications are promising. This method can be used as an alternative to obtain analytical and approximate solutions of different types of nonlinear differential equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Adaptive pseudo‐transient‐continuation‐Galerkin methods for semilinear elliptic partial differential equations

In this article, we investigate the application of pseudo‐transient‐continuation (PTC) schemes for the numerical solution of semilinear elliptic partial differential equations, with possible singular perturbations. We will outline a residual reduction analysis within the framework of general Hilbert spaces, and, subsequently, use the PTC‐methodology in the context of finite element discretizations of semilinear boundary value problems. Our approach combines both a prediction‐type PTC‐method (for infinite dimensional problems) and an adaptive finite element discretization (based on a robust a posteriori residual analysis), thereby leading to a fully adaptive PTC ‐Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for different examples.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2005–2022, 2017     

Optimal convergence analysis of mixed finite element methods for fourth‐order elliptic and parabolic problems

By using a special interpolation operator developed by Girault and Raviart (finite element methods for Navier‐Stokes Equations, Springer‐Verlag, Berlin, 1986), we prove that optimal error bounds can be obtained for a fourth‐order elliptic problem and a fourth‐order parabolic problem solved by mixed finite element methods on quasi‐uniform rectangular meshes. Optimal convergence is proved for all continuous tensor product elements of order k ≥ 1. A numerical example is provided for solving the fourth‐order elliptic problem using the bilinear element. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

A basis function for approximation and the solutions of partial differential equations

In this article, we introduce a type of basis functions to approximate a set of scattered data. Each of the basis functions is in the form of a truncated series over some orthogonal system of eigenfunctions. In particular, the trigonometric eigenfunctions are used. We test our basis functions on recovering the well‐known Franke's and Peaks functions given by scattered data, and on the extension of a singular function from an irregular domain onto a square. These basis functions are further used in Kansa's method for solving Helmholtz‐type equations on arbitrary domains. Proper one level and two level approximation techniques are discussed. A comparison of numerical with analytic solutions is given. The numerical results show that our approach is accurate and efficient. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

The dual reciprocity boundary elements method for the linear and nonlinear two‐dimensional time‐fractional partial differential equations

In this paper, we apply the dual reciprocity boundary elements method for the numerical solution of two‐dimensional linear and nonlinear time‐fractional modified anomalous subdiffusion equations and time‐fractional convection–diffusion equation. The fractional derivative of problems is described in the Riemann–Liouville and Caputo senses. We employ the linear radial basis function for interpolation of the nonlinear, inhomogeneous and time derivative terms. This method is improved by using a predictor–corrector scheme to overcome the nonlinearity which appears in the nonlinear problems under consideration. The accuracy and efficiency of the proposed schemes are checked by five test problems. The proposed method is employed for solving some examples in two dimensions on unit square and also in complex regions to demonstrate the efficiency of the new technique. Copyright © 2016 John Wiley & Sons, Ltd.