Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

Linearized and rational approximation method for solving non‐linear Burgers' equation

A new numerical method called linearized and rational approximation method is presented to solve non‐linear evolution equations. The utility of the method is demonstrated for the case of differentiation of functions involving steep gradients. The solution of Burgers' equation is presented to illustrate the effectiveness of the technique for the solution of non‐linear evolution equations exhibiting nearly discontinuous solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

A six-dimensional quadrature procedure

An account is given of the use of Gaussian quadrature product formulae in the evaluation of certain six-dimensional, two-centre integrals involving one-electron Green's functions. These integrals occur in a new molecular variational principle recently proposed by Hall, Hyslop and Rees [1] from which an approximate energy may be derived which can be shown to be at least as good as that obtained from the Rayleigh-Ritz principle. Reductions in computing time are realized by removing certain singularities using a subtraction technique and also by using an empirically determined Richardson-type extrapolation formula.This paper was presented during the session on numerical integration methods for molecules of the 1970 Quantum Theory Conference in Nottingham. It has been revised in the light of the interesting discussion which followed.     

QUADRATURE METHODS FOR HIGHLY OSCILLATORY SINGULAR INTEGRALS

We address the evaluation of highly oscillatory integrals,with power-law and logarithmic singularities.Such problems arise in numerical methods in engineering.Notably,the evaluation of oscillatory integrals dominates the run-time for wave-enriched boundary integral formulations for wave scattering,and many of these exhibit singularities.We show that the asymptotic behaviour of the integral depends on the integrand and its derivatives at the singular point of the integrand,the stationary points and the endpoints of the integral.A truncated asymptotic expansion achieves an error that decays faster for increasing frequency.Based on the asymptotic analysis,a Filon-type method is constructed to approximate the integral.Unlike an asymptotic expansion,the Filon method achieves high accuracy for both small and large frequency.Complex-valued quadrature involves interpolation at the zeros of polynomials orthogonal to a complex weight function.Numerical results indicate that the complex-valued Gaussian quadrature achieves the highest accuracy when the three methods are compared.However,while it achieves higher accuracy for the same number of function evaluations,it requires signi cant additional cost of computation of orthogonal polynomials and their zeros.     

Chebyshev spectral variational integrator and applications

The Chebyshev spectral variational integrator(CSVI) is presented in this paper. Spectral methods have aroused great interest in approximating numerically a smooth problem for their attractive geometric convergence rates. The geometric numerical methods are praised for their excellent long-time geometric structure-preserving properties.According to the generalized Galerkin framework, we combine two methods together to construct a variational integrator, which captures the merits of both methods. Since the interpolating points of the variational integrator are chosen as the Chebyshev points,the integration of Lagrangian can be approximated by the Clenshaw-Curtis quadrature rule, and the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of configuration variables and the corresponding derivatives. The numerical float errors of the first-order spectral differentiation matrix can be alleviated by using a trigonometric identity especially when the number of Chebyshev points is large. Furthermore, the spectral variational integrator(SVI) constructed by the Gauss-Legendre quadrature rule and the multi-interval spectral method are carried out to compare with the CSVI, and the interesting kink phenomena for the Clenshaw-Curtis quadrature rule are discovered. The numerical results reveal that the CSVI has an advantage on the computing time over the whole progress and a higher accuracy than the SVI before the kink position. The effectiveness of the proposed method is demonstrated and verified perfectly through the numerical simulations for several classical mechanics examples and the orbital propagation for the planet systems and the Solar system.     

Torsional vibration of functionally graded nano-rod under magnetic field supported by a generalized torsional foundation based on nonlocal elasticity theory

Abstract

This article contains the nonlocal elasticity theory to capture size effects in functionally graded (FG) nano-rod under magnetic field supported by a torsional foundation. Torque effect of an axial magnetic field on an FG nano-rod has been defined using Maxwell’s relation. The material properties were assumed to vary according to the power law in radial direction. The Navier equation and boundary conditions of the size-dependent FG nano-rod were derived by the Hamilton’s principle. These equations were solved by employing the generalized differential quadrature method (GDQM). Presented model has the ability to turn into the classical model if the material length scale parameter is taken to be zero. The effects of some parameters, such as inhomogeneity constant, magnetic field and small-scale parameter, were studied. As an important result of this study can be stated that an FG nano-rod model based on the nonlocal elasticity theory behaves softer and has smaller natural frequency.     

On Gauss-Type Quadrature Rules

Some Gauss-type Quadrature rules over [0, 1], which involve values and/or the derivative of the integrand at 0 and/or 1, are investigated. Our work is based on the orthogonal polynomials with respect to linear weight function ω(t): = 1 ? t over [0, 1]. These polynomials are also linked with a class of recently developed “identity-type functions”. Along the lines of Golub's work, the nodes and weights of the quadrature rules are computed from Jacobi-type matrices with simple rational entries. Computational procedures for the derived rules are tested on different integrands. The proposed methods have some advantage over the respective Gauss-type rules with respect to the Gauss weight function ω(t): = 1 over [0, 1].     

Mass and momentum conservation of the least-squares spectral collocation method for the Navier-Stokes equations

From the literature, it is known that the Least-Squares Spectral Element Method (LSSEM) for the stationary Stokes equations performs poorly with respect to mass conservation but compensates this lack by a superior conservation of momentum. Furthermore, it is known that the Least-Squares Spectral Collocation Method (LSSCM) leads to superior conservation of mass and momentum for the stationary Stokes equations. In the present paper, we consider mass and momentum conservation of the LSSCM for time-dependent Stokes and Navier-Stokes equations. We observe that the LSSCM leads to improved conservation of mass (and momentum) for these problems. Furthermore, the LSSCM leads to the well-known time-dependent profiles for the velocity and the pressure profiles. To obtain these results, we use only a few elements, each with high polynomial degree, avoid normal equations for solving the overdetermined linear systems of equations and introduce the Clenshaw-Curtis quadrature rule for imposing the average pressure to be zero. Furthermore, we combined the transformation of Gordon and Hall (transfinite mapping) with the least-squares spectral collocation scheme to discretize the internal flow problems.     

Recursion relations for the extended Krylov subspace method

The evaluation of matrix functions of the form f(A)v, where A is a large sparse or structured symmetric matrix, f is a nonlinear function, and v is a vector, is frequently subdivided into two steps: first an orthonormal basis of an extended Krylov subspace of fairly small dimension is determined, and then a projection onto this subspace is evaluated by a method designed for small problems. This paper derives short recursion relations for orthonormal bases of extended Krylov subspaces of the type K^m,mi+1(A)=span{A^-m+1v,…,A^-1v,v,Av,…,A^miv}, m=1,2,3,…, with i a positive integer, and describes applications to the evaluation of matrix functions and the computation of rational Gauss quadrature rules.