A spectral collocation method based on integrated Chebyshev polynomials for two-dimensional biharmonic boundary-value problems

This paper reports a new spectral collocation method for numerically solving two-dimensional biharmonic boundary-value problems. The construction of the Chebyshev approximations is based on integration rather than conventional differentiation. This use of integration allows: (i) the imposition of the governing equation at the whole set of grid points including the boundary points and (ii) the straightforward implementation of multiple boundary conditions. The performance of the proposed method is investigated by considering several biharmonic problems of first and second kinds; more accurate results and higher convergence rates are achieved than with conventional differential methods.     

A second-order continuity domain decomposition technique based on integrated Chebyshev polynomials for two-dimensional elliptic problems

This paper presents a second-order continuity non-overlapping domain decomposition (DD) technique for numerically solving second-order elliptic problems in two-dimensional space. The proposed DD technique uses integrated Chebyshev polynomials to represent the solution in subdomains. The constants of integration are utilized to impose continuity of the second-order normal derivative of the solution at the interior points of subdomain interfaces. To also achieve a C²$C^2$ function at the intersection of interfaces, two additional unknowns are introduced at each intersection point. Numerical results show that the present DD method yields a higher level of accuracy than conventional DD techniques based on differentiated Chebyshev polynomials.     

Numerical technique based on generalized Laguerre and shifted Chebyshev polynomials

In this study, we present a numerical scheme for solving a class of fractional partial differential equations. First, we introduce psi -Laguerre polynomials like psi-shifted Chebyshev polynomials and employ these newly introduced polynomials for the solution of space-time fractional differential equations. In our approach, we project these polynomials to develop operational matrices of fractional integration. The use of these orthogonal polynomials converts the problem under consideration into a system of algebraic equations. The solution of this system provide us the desired results. The convergence of the proposed method is analyzed. Finally, some illustrative examples are included to observe the validity and applicability of the proposed method.     

Numerical solutions based on Chebyshev collocation method for singularly perturbed delay parabolic PDEs

The main purpose of this work is to provide a numerical approach for the delay partial differential equations based on a spectral collocation approach. In this research, a rigorous error analysis for the proposed method is provided. The effectiveness of this approach is illustrated by numerical experiments on two delay partial differential equations. Copyright © 2013 John Wiley & Sons, Ltd.     

Some problems of A. Kroó on multiple Chebyshev polynomials

Three problems of A. Kroó on multiple Chebyshev polynomials are solved using the Borsuk–Ulam antipodal theorem.     

A collocation approach for solving linear complex differential equations in rectangular domains

In this paper, a collocation method is presented to find the approximate solution of high‐order linear complex differential equations in rectangular domain. By using collocation points defined in a rectangular domain and the Bessel polynomials, this method transforms the linear complex differential equations into a matrix equation. The matrix equation corresponds to a system of linear equations with the unknown Bessel coefficients. The proposed method gives the analytic solution when the exact solutions are polynomials. Numerical examples are included to demonstrate the validity and applicability of the technique and the comparisons are made with existing results. The results show the efficiency and accuracy of the present work. All of the numerical computations have been performed on a computer using a program written in MATLAB v7.6.0 (R2008a). Copyright © 2012 John Wiley & Sons, Ltd.     

A fast finite difference method for biharmonic equations on irregular domains and its application to an incompressible Stokes flow

Biharmonic equations have many applications, especially in fluid and solid mechanics, but is difficult to solve due to the fourth order derivatives in the differential equation. In this paper a fast second order accurate algorithm based on a finite difference discretization and a Cartesian grid is developed for two dimensional biharmonic equations on irregular domains with essential boundary conditions. The irregular domain is embedded into a rectangular region and the biharmonic equation is decoupled to two Poisson equations. An auxiliary unknown quantity Δu along the boundary is introduced so that fast Poisson solvers on irregular domains can be used. Non-trivial numerical examples show the efficiency of the proposed method. The number of iterations of the method is independent of the mesh size. Another key to the method is a new interpolation scheme to evaluate the residual of the Schur complement system. The new biharmonic solver has been applied to solve the incompressible Stokes flow on an irregular domain.      

A new method for two dimensional hyperbolic problems in semi‐unbounded irregular domains

Anewmethod called the full rational differential quadrature method is presented to deal with two dimensional linear and nonlinear hyperbolic problems in semi‐unbounded irregular domains. The spacial and temporal discretizations are both implemented by the rational differential quadrature method (RDQM). The RDQM, proven to be A‐stable (Chen and Tanaka, Comput Mech 28 (2002) 331–338) in the temporal discretization, is much more efficient than the finite difference schemes widely used in earlier works. In addition, the irregular boundary conditions are treated by the direct expansion method(DEM). Numerical experiments show that the present method is of high efficiency and easy to implement. © 2011Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Discrete weighted least-squares method for the Poisson and biharmonic problems on domains with smooth boundary

In this article a discrete weighted least-squares method for the numerical solution of elliptic partial differential equations exhibiting smooth solution is presented. It is shown how to create well-conditioned matrices of the resulting system of linear equations using algebraic polynomials, carefully selected matching points and weight factors. Two simple algorithms generating suitable matching points, the Chebyshev matching points for standard two-dimensional domains and the approximate Fekete points of Sommariva and Vianello for general domains, are described. The efficiency of the presented method is demonstrated by solving the Poisson and biharmonic problems with the homogeneous Dirichlet boundary conditions defined on circular and annular domains using basis functions in the form satisfying and in the form not satisfying the prescribed boundary conditions.     

On a rational differential quadrature method in irregular domains for problems with boundary layers

A rational differential quadrature method in irregular domains (RDQMID) is investigated to deal with a kind of singularly perturbed problems with boundary layers. Through a transformation, the boundary layer, which may be not straight, is transformed into a segment of a line parallel to one of the Cartesian axes. The rational differential quadrature method (RDQM) is applied to discretize the governing equation. Finally, a direct expansion method of the boundary conditions (DEMBC) is raised to deal with the boundary conditions. Numerical experiments show that RDQMID is of high accuracy, efficiency and easy to programme.     

The collocation method for mixed boundary value problems on domains with curved polygonal boundaries

Summary. We consider an indirect boundary integral equation formulation for the mixed Dirichlet-Neumann boundary value problem for the Laplace equation on a plane domain with a polygonal boundary. The resulting system of integral equations is solved by a collocation method which uses a mesh grading transformation and a cosine approximating space. The mesh grading transformation method yields fast convergence of the collocation solution by smoothing the singularities of the exact solution. A complete stability and solvability analysis of the transformed integral equations is given by use of a Mellin transform technique, in a setting in which each arc of the polygon has associated with it a periodic Sobolev space. Received April 15, 1995 / Revised version received April 10, 1996     

A Cartesian‐grid collocation method based on radial‐basis‐function networks for solving PDEs in irregular domains

This paper reports a new Cartesian‐grid collocation method based on radial‐basis‐function networks (RBFNs) for numerically solving elliptic partial differential equations in irregular domains. The domain of interest is embedded in a Cartesian grid, and the governing equation is discretized by using a collocation approach. The new features here are (a) one‐dimensional integrated RBFNs are employed to represent the variable along each line of the grid, resulting in a significant improvement of computational efficiency, (b) the present method does not require complicated interpolation techniques for the treatment of Dirichlet boundary conditions in order to achieve a high level of accuracy, and (c) normal derivative boundary conditions are imposed by means of integration constants. The method is verified through the solution of second‐ and fourth‐order PDEs; accurate results and fast convergence rates are obtained. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Some results on the asymptotic behavior for hyperbolic problems in cylindrical domains becoming unbounded

We study here the asymptotic behavior of the solution of a hyperbolic problem defined on a cylindrical domain [0,T]×p(−?,?)×ω⊂[0,T]×Rn when ?→∞. We show that, under very general assumptions, the solution of this problem converges faster than any power of towards the solution of another hyperbolic problem, defined on [0,T]×ω, in any bounded subdomain. We give both necessary and sufficient conditions for this convergence to occur.     

A mixed hybrid formulation based on oscillated finite element polynomials for solving Helmholtz problems

A mixed-hybrid-type formulation is proposed for solving Helmholtz problems. This method is based on (a) a local approximation of the solution by oscillated finite element polynomials and (b) the use of Lagrange multipliers to “weakly” enforce the continuity across element boundaries. The computational complexity of the proposed discretization method is determined mainly by the total number of Lagrange multiplier degrees of freedom introduced at the interior edges of the finite element mesh, and the sparsity pattern of the corresponding system matrix. Preliminary numerical results are reported to illustrate the potential of the proposed solution methodology for solving efficiently Helmholtz problems in the mid- and high-frequency regimes.     

A rational approximation based on Bernstein polynomials for high order initial and boundary values problems

We introduce a new method to solve high order linear differential equations with initial and boundary conditions numerically. In this method, the approximate solution is based on rational interpolation and collocation method. Since controlling the occurrence of poles in rational interpolation is difficult, a construction which is found by Floater and Hormann [1] is used with no poles in real numbers. We use the Bernstein series solution instead of the interpolation polynomials in their construction. We find that our approximate solution has better convergence rate than the one found by using collocation method. The error of the approximate solution is given in the case of the exact solution f ∈ C^d+2[a, b].     

A reduced MFE formulation based on POD for the non-stationary conduction-convection problems

In this article, a reduced mixed finite element (MFE) formulation based on proper orthogonal decomposition (POD) for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.     

A method of Fourier series for solution of problems in piecewise inhomogeneous domains with rectilinear crack (screen)

In the framework of the theory of harmonic functions, potentials of steady state processes (heat conduction, filtration, or electrostatics) in the piecewise inhomogeneous plane separated by a rectilinear strongly permeable crack or by a weakly permeable screen into two half-planes with quadratic permeability functions are constructed. The motion is induced by given singular points of the potential (sources, sinks, etc.). Compact formulas that directly express potentials in these domains in terms of harmonic functions are obtained; the resulting functions map the set of harmonic functions to the set of potentials conserving the type of singularities.     

A layer potential analysis for transmission problems for Brinkman‐type systems in Lipschitz domains in

The aim of this paper is to establish a well‐posedness result for a boundary value problem of transmission‐type for the standard and generalized Brinkman systems in two Lipschitz domains in , the former being bounded, and the latter, its complement in . As a first step, we establish a well‐posedness result for a transmission problem for the standard Brinkman systems on complementary Lipschitz domains in by making use of the Potential theory developed for such a system. As a second step, we prove our desired result (in L²‐based Sobolev spaces) by using a method based on Fredholm operator theory and the well‐posedness result from the previous step.     

A collocation method based on the Bessel functions of the first kind for singular perturbated differential equations and residual correction

In this paper, a collocation method is given to solve singularly perturbated two‐point boundary value problems. By using the collocation points, matrix operations and the matrix relations of the Bessel functions of the first kind and their derivatives, the boundary value problem is converted to a system of the matrix equations. By solving this system, the approximate solution is obtained. Also, an error problem is constructed by the residual function, and it is solved by the presented method. Thus, the error function is estimated, and the approximate solutions are improved. Finally, numerical examples are given to show the applicability of the method, and also, our results are compared by existing results. Copyright © 2014 JohnWiley & Sons, Ltd.