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.

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.

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.

Chebyshev Collocation for Optimal Control in a Thermoconvective Flow

In this paper a Chebyshev collocation method is used for solving numerically an optimal boundary control problem in a thermoconvective fluid flow.The aim of this study is to demonstrate the capabilities of these numerical techniques for handling this kind of problems.As the problem is treated in the primitive variable formulation additional boundary conditions for the pressure and the auxiliary pressure fields are required to avoid spurious modes.A dependence of the convergence of the method on the penalizing parameter that appears in the functional cost is observed.As this parameter approaches zero some singular behaviour in the control function is observed and the order of the method decreases.These singularities are irrelevant in the problem as a regularized control function produces the same results.

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.

Multidomain pseudospectral methods for nonlinear convection-diffusion equations

Multidomain pseudospectral approximations to nonlinear convection-diffusion equations are considered. The schemes are formulated with the Legendre-Galerkin method but the nonlinear term is collocated at the Legendre/Chebyshev-Gauss-Lobatto points inside each subinterval. Appropriate base functions are introduced so that the matrix of the system is sparse, and the method can be implemented efficiently and in parallel. The stability and the optimal rate of convergence of the methods are proved. Numerical results are given for both the single domain and the multidomain methods to make a comparison.

The instability of water-mud interface in viscous two-layer flow with large viscosity contrast

The temporal instability of parallel viscous two-phase mixing layers is extended to current-fluid mud by considering a composite error function velocity profile. The influence of viscosity ratio, Reynolds number, and Froude number on the instability of the system are discussed and a new phenomenon never discussed is investigated based on our numerical results. It is shown that viscosity can enlarge the unstable wave number range, cause new instability modes, and certainly reduce the growth rate of Kelvin-Helmholtz （K-H） instability.

Probabilistic Response and Short-Term Extreme Load Estimation of Offshore Monopile Wind Turbine Towers by Probability Density Evolution Method

A new analysis framework based on probability density evolution method(PDEM)and its Chebyshev collocation solution are introduced to predict the dynamic response and short-term extreme load of offshore wind turbine(OWT)towers subjected to random sea state.With regard to the stochastic responses,random function method is employed to generate samples of sea elevation,the probability density evolution equation(PDEE)is solved to calculate time-variant probability density functions of structural responses.For the probabilistic load estimation,a FAST model of NREL 5MW offshore turbine is established to obtain samples of bending moment at the tower base.The equivalent extreme event theory is used to construct a virtual stochastic process(VSP)to assess the short-term extreme load.The results indicate that the proposed approach can predict time-variant probability density functions of the structural responses,and shows good agreement with Monte Carlo simulations.Additionally,the predicted short-term extreme load can capture the fluctuation at the tail of the extreme value distribution,thus is more rational than results from the typical distribution models.Overall,the proposed method shows good adaptation,precision and efficiency for the dynamic response analysis and load estimation of OWT towers.

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.

Floquet instability of a large density ratio liquid-gas coaxial jet with periodic fluctuation

By numerical simulation of basic flow, this paper extends Floquet stability analysis of interracial flow with periodic fluctuation into large density ratio range. Stability of a liquid aluminum jet in a coaxial nitrogen stream with velocity fluctuation is investigated by Chebyshev collocation method and the Floquet theory. Parametric resonance of the jet and the influences of different physical parameters on the instability are discussed. The results show qualitative agreement with the available experimental data.

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.

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.

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.

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.

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.

Efficient Solution of a Generalized Eigenvalue Problem Arising in a Thermoconvective Instability

The aim of this paper is to develop an efficient numerical method to compute the eigenvalues of the stability analysis of a problem describing the motion of a fluid within a cylindrical container heated non-homogeneously from below.An axisymmetric stationary motion settles in,at certain values of the external parameters appearing in the set of partial differential equations modeling the problem.This basic solution is computed by discretizing the equations with a Chebyshev collocation method.Its linear stability is formulated with a generalized eigenvalue problem.The numerical approach(generalized Arnoldi method)uses the idea of preconditioning the eigenvalue problem with a modified Cayley transformation before applying the Arnoldi method.Previous works have dealt with transformations requiring regularity to one of the submatrices.In this article we extend those results to the case in which that submatrix is singular.This method allows a fast computation of the critical eigenvalues which determine whether the steady flow is stable or unstable.The algorithm based on this method is compared to the QZ method and is found to be computationally more efficient.The reliability of the computed eigenvalues in terms of stability is confirmed via pseudospectra calculations.