Convergence Analysis for the Chebyshev Collocation Methods to Volterra Integral Equations with a Weakly Singular Kernel

In this paper,a Chebyshev-collocation spectral method is developed for Volterra integral equations(VIEs)of second kind with weakly singular kernel.We first change the equation into an equivalent VIE so that the solution of the new equation possesses better regularity.The integral term in the resulting VIE is approximated by Gauss quadrature formulas using the Chebyshev collocation points.The convergence analysis of this method is based on the Lebesgue constant for the Lagrange interpolation polynomials,approximation theory for orthogonal polynomials,and the operator theory.The spectral rate of convergence for the proposed method is established in the L^(∞)-norm and weighted L^(2)-norm.Numerical results are presented to demonstrate the effectiveness of the proposed method.

THE ACCURACY COMPARISON BETWEEN CHEBYSHEV-τ METHOD AND CHEBYSHEV COLLOCATION METHOD

This paper is devoted to investigate the accuracy of the Pseudo spectral scheme with the Chebyshev tau method and Chebyshev collocation method. The computational results of the nonlinear disturbance development in plane Poiseuille flow for both methods are presented and compared in detail. It is acknowledged that the Chebyshev collocation method has higher precision than the other one, especially for near netural situation.

Simulation of electrically driven jet using Chebyshev collocation method

The model of electrically driven jet is governed by a series of quasi 1D dimensionless partial differential equations(PDEs).Following the method of lines,the Chebyshev collocation method is employed to discretize the PDEs and obtain a system of differential-algebraic equations(DAEs).By differentiating constrains in DAEs twice,the system is transformed into a set of ordinary differential equations(ODEs) with invariants.Then the implicit differential equations solver 'ddaskr' is used to solve the ODEs and post-stabilization is executed at the end of each step.Results show the distributions of radius,linear charge density,stretching ratio and also the horizontal velocity at a time point.Meanwhile,the spiral and expanding projections to X-Y plane of the jet centerline suggest the occurring of bending instability.

Laminar flow of micropolar fluid in rectangular microchannels

Compared with the classic flow on macroscale, flows in microchannels have some new phenomena such as the friction increase and the flow rate reduction. Papautsky and co-workers explained these phenomena by using a micropolar fluid model where the effects of micro-rotation of fluid molecules were taken into account. But both the curl of velocity vector and the curl of micro-rotation gyration vector were given incorrectly in the Cartesian coordinates and then the micro-rotation gyration vector had only one component in the z-direction. Besides, the gradient term of the divergence of micro-rotation gyration vector was missed improperly in the angular moment equation. In this paper, the governing equations for laminar flows of micropolar fluid in rectangular microchannels are reconstructed. The numerical results of velocity profiles and micro-rotation gyrations are obtained by a procedure based on the Chebyshev collocation method. The micropolar effects on velocity and micro-rotation gyration are discussed in detail.

Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations

An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and in- compressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangu- lar sub-domains. In each sub-domain, an improved projection scheme with second-order accuracy is used to deal with the coupling of velocity and pressure, and the Chebyshev collocation spectral method （CSM） is adopted to execute the spatial discretization. The influence matrix technique is employed to enforce the continuities of both variables and their normal derivatives between the adjacent sub-domains. The imposing of the Neu- mann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet bound- ary conditions on interface and using the first-order normal derivatives as transmission conditions to keep the continuities of variables is proposed to overcome this trouble. Three test cases are used to verify the accuracy and efficiency, and the detailed comparison be- tween the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.

Convergence analysis for delay Volterra integral equation

In this article we use Chebyshev spectral collocation method to deal with the Volterra integral equation which has two kinds of delay items. We use linear transformation to make the interval into a fixed interval [-1, 1]. Then we use the Gauss quadrature formula to approximate the solution. With the help of lemmas, we get the result that the numerical error decay exponentially in the infinity norm and the Chebyshev weighted Hilbert space norm. Some numerical experiments are given to confirm our theoretical prediction.

Instability in Three-Dimensional Magnetohydrodynamic Flows of an Electrically Conducting Fluid

The three-dimensional instability of an electrically conducting fluid between two parallel plates affected by an imposed transversal magnetic field is numerically investigated by a Chebyshev collocation method. The QZ method is utilized to obtain neutral curves of the linear instability. The details of instability are analyzed by solving the generalized Orr-Sommerfeld equation. The critical Reynolds number Rec, the stream-wise and span-wise critical wave numbers αc and βc are obtained for a wide range of Hartmann number Ha. The effects of Lorentz force and span-wise perturbation on three-dimensional instability are investigated. The results show that magnetic field would suppress the instability and critical Reynolds number tends to be larger than that for two-dimensional instability.

Collocation Methods for Hyperbolic Partial Differential Equations with Singular Sources

A numerical study is given on the spectral methods and the high order WENO finite difference scheme for the solution of linear and nonlinear hyperbolic partial differential equations with stationary and non-stationary singular sources.The singular source term is represented by theδ-function.For the approximation of theδ-function,the direct projection method is used that was proposed in[6].Theδ-function is constructed in a consistent way to the derivative operator.Nonlinear sine-Gordon equation with a stationary singular source was solved with the Chebyshev collocation method.Theδ-function with the spectral method is highly oscillatory but yields good results with small number of collocation points.The results are compared with those computed by the second order finite difference method.In modeling general hyperbolic equations with a non-stationary singular source,however,the solution of the linear scalar wave equation with the nonstationary singular source using the direct projection method yields non-physical oscillations for both the spectral method and the WENO scheme.The numerical artifacts arising when the non-stationary singular source term is considered on the discrete grids are explained.

A Hybrided Trapezoidal-Difference Scheme for Nonlinear Time-Fractional Fourth-Order Advection-Dispersion Equation Based on Chebyshev Spectral Collocation Method

In this paper,we firstly present a novel simple method based on a Picard integral type formulation for the nonlinear multi-dimensional variable coefficient fourthorder advection-dispersion equation with the time fractional derivative order a2(1,2).A new unknown function v(x,t)=■u(x,t)/■t is introduced and u(x,t)is recovered using the trapezoidal formula.As a result of the variable v(x,t)are introduced in each time step,the constraints of traditional plans considering the non-integer time situation of u(x,t)is no longer considered.The stability and solvability are proved with detailed proofs and the precise describe of error estimates is derived.Further,Chebyshev spectral collocation method supports accurate and efficient variable coefficient model with variable coefficients.Several numerical results are obtained and analyzed in multi-dimensional spatial domains and numerical convergence order are consistent with the theoretical value 3-a order for different a under infinite norm.

Mathematical modeling and numerical simulation of telephone cord buckles of elastic films

An annular sector model for the telephone cord buckles of elastic thin films on rigid substrates is established, in which the von Krman plate equations in polar coordinates are used for the elastic thin film and a discrete version of the Griffith criterion is applied to determine the shape and scale of the parameters. A numerical algorithm combining the Newmark-β scheme and the Chebyshev collocation method is designed to numerically solve the problem in a quasi-dynamic process. Numerical results are presented to show that the numerical method works well and the model agrees well with physical observations, especially successfully simulated for the first time the telephone cord buckles with two humps along the ridge of each section of a buckle.

Stability of fluid flow in a Brinkman porous medium–A numerical study

The stability of fluid flow in a horizontal layer of Brinkman porous medium with fluid viscosity different from effective viscosity is investigated. A modified Orr-Sommerfeld equation is derived and solved numerically using the Chebyshev collocation method. The critical Reynolds number Re, the critical wave number ac and the critical wave speed cc are computed for various values of porous parameter and ratio of viscosities. Based on these parameters, the stability characteristics of the system are discussed in detail. Streamlines are presented for selected values of parameters at their critical state.