0rr－Sommerfeld方程数值解法中的复广义矩阵特征值问题

Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程的求解通常可以化为一个复广义矩阵特征值问题ＡＸ＝ωＢＸ。本文用酉变换分别约化Ａ，Ｂ为上Ｈｅｓｓｅｎｂｅｒｇ阵和上三角阵，然后利用Ｍｕｌｌｅｒ求根方法可以求出其全部特征值，其中特征多项式的值由Ｈｙｍａｎ方法给出。当仅需要判断有无不稳定模态时，利用一个简单的矩阵变换将其化为强特征值的求解问题，从而可使用最简单的幂迭代，Ｃｈｅｂｙｓｈｅｖ配置点法算例表明两种算法均快速有效。     

0rr－Sommerfeld方程数值解法中的复广义矩阵特征值问题

Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程的求解通常可以化为一个复广义矩阵特征值问题ＡＸ＝ωＢＸ。本文用酉变换分别约化Ａ，Ｂ为上Ｈｅｓｓｅｎｂｅｒｇ阵和上三角阵，然后利用Ｍｕｌｌｅｒ求根方法可以求出其全部特征值，其中特征多项式的值由Ｈｙｍａｎ方法给出。当仅需要判断有无不稳定模态时，利用一个简单的矩阵变换将其化为强特征值的求解问题，从而可使用最简单的幂迭代，Ｃｈｅｂｙｓｈｅｖ配置点法算例表明两种算法均快速有效。     

槽道湍流近壁区相干结构的数值模拟

从流动稳定性理论中的共振三波出发，采用类似湍流直接数值模拟中的最小单元概念，用伪谱方法，对槽道湍流近壁区的单个相干结构进行了数值模拟，得到了与实验基本一致的结果，同时也搞清楚了在谱方法中用Ｃｈｅｂｙｓｈｅｖ－τ方法和配置点法哪个更精确的问题。     

平面Poiseuille流动力系统稳定性分析

本文把平面Ｐｏｉｓｅｕｉｌｌｅ流的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程化为一非线性动力系统,用线性稳定性理论对其进行分析研究,得到一组特征方程,对离散后的特征方程施加一个矩阵变换,计算出不同波数下的临界雷诺数并与Ｏｒｓｚａｙ和Ｔｈｍｏａｓ的计算结果进行比较,最后画出了Ｐｏｉｓｅｕｉｌｌｅ流的中性稳定性曲线。     

多个共面任意分布表面裂纹的应力强度因子

采用线弹簧模型求解多个共面任意分布表面裂纹的应力强度因子。基于Ｒｅｉｓｓｎｅｒ板理论和连续分布位错思想,通过积分变换方法,将含有多个共面任意分布表面裂纹的无限平板问题归结为一组Ｃａｕｃｈｙ型奇异积分方程。利用Ｇａｕｓｓ－Ｇｈｅｂｙｓｈｅｖ笔法获得了奇异积分方程的数值解。为验证本文法的正确性,文中最后给出了有关应力强度因子或Ｐ－Ｖ曲线的数值结果并与现有的理论结果或实验结果进行了对比。结果表明了连续位     

层状介质垂直界面有限裂纹问题的积分变换解

层状弹性材料包含垂直于界面有限裂纹时，可运用富里叶变换及引用位错密度函数，导出了反映裂纹尖端奇异性的奇异积分方程组，并使用Ｌｏｂａｔｔｏ－ｃｈｅｂｙｓｈｅｖ方法解此方程组，最后得到裂纹尖端应力强度因子，为检验方法的正确性，对某两层含裂实际结构进行了计算，结果是满意的。     

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

不可压粘流N－S方程的边界积分解法

对原变量的Ｎ－Ｓ方程进行一阶时间离散，采用共轭梯度法解除压强－速度的耦合．对所得的一系列Ｌａｐｌａｃｅ方程、Ｐｏｓｓｉｏｎ方程和Ｈｅｌｍｈｏｔｚ方程均进行边界积分法求解，首次得到了粘性Ｎ－Ｓ方程的边界积分表示式．圆柱的定常、非定常尾迹计算结果表明了本文方法的有效性．     

旋转圆盘受均布压力作用时的摄动分析

本文对中心带轴匀速旋转圆盘在均布压力作用下的受力及变形问题进行了摄动分析，得到了受载条件为Ｊ＾２／ＪＨ＜＜１时，圆盘应力及形变的渐近解。其形变渐近解为Ｗｈｉｔａｋｅｒ函数形式比较表明，渐近解与Ｃｈｅｂｙｓｈｅｖ数值解有较好的一致性。     

剪切湍流大尺度相干结构的模式研究

发展了一种计算剪切湍流大尺度相干结构的新模式．该模式的基础是认为大尺度相干结构为湍流场中流体脉动能量增长最快的那部分，且包含大部分的湍流脉动能量．在此基础上。通过对湍流相干能量方程的推演。建立了描述大尺度相干结构的特征控制方程，并应用Ｃｈｅｂｙｓｈｅｖ多项式方法求得湍流相干能量的最大增长率在波数空间的分布，从而获得对应的大尺度相干结构．应用该模式研究了槽流和一自然对流中的大尺度相干结构，得到的近壁区流动结构与实验现象十分接近．     

Hydrodynamic stability,the Chebyshev tau method and spurious eigenvalues

The Chebyshev tau method is examined; a numerical technique which in recent years has been successfully applied to many hydrodynamic stability problems. The orthogonality of Chebyshev functions is used to rewrite the differential equations as a generalized eigenvalue problem. Although a very efficient technique, the occurrence of spurious eigenvalues, which are not always easy to identify, may lead one to believe that a system is unstable when it is not. Thus, the elimination of spurious eigenvalues is of great importance. Boundary conditions are included as rows in the matrices of the generalized eigenvalue problem and these have been observed to be one cause of spurious eigenvalues. Removing boundary condition rows can be difficult. This problem is addressed here, in application to the Bénard convection problem, and to the Orr-Sommerfeld equation which describes parallel flow. The procedure given here can be applied to a wide range of hydrodynamic stability problems.Received: 4 July 2002, Accepted: 13 September 2002, Published online: 27 June 2003     

On the accuracy of the calculation of transient growth in plane Poiseuille flow

This paper addresses the accuracy of numerical methods to compute the transient energy growth of plane Poiseuille flow. We show that using the Chebyshev collocation method to discretize the linearized governing equations in the wall‐normal direction can introduce numerical problems when computing the energy evolution of the flow. We demonstrate that spurious eigenmodes of the discretized linear operator and numerical integration errors are the possible sources of the numerical problems, and we also show that spurious eigenvalues with negative real parts of large magnitude can affect the calculation of energy growth. These difficulties can be avoided by using a spectral Galerkin method where the basis functions satisfy the boundary conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

MHD Falkner-Skan flow of Maxwell fluid by rational Chebyshev collocation method

The magnetohydrodynamics （MHD） Falkner-Skan flow of the Maxwell fluid is studied. Suitable transform reduces the partial differential equation into a nonlinear three order boundary value problem over a semi-infinite interval. An efficient approach based on the rational Chebyshev collocation method is performed to find the solution to the proposed boundary value problem. The rational Chebyshev collocation method is equipped with the orthogonal rational Chebyshev function which solves the problem on the semi-infinite domain without truncating it to a finite domain. The obtained results are presented through the illustrative graphs and tables which demonstrate the affectivity, stability, and convergence of the rational Chebyshev collocation method. To check the accuracy of the obtained results, a numerical method is applied for solving the problem. The variations of various embedded parameters into the problem are examined.     

弹性连接旋转柔性梁动力学分析

采用Chebyshev 谱方法对考虑根部连接弹性的平面内旋转柔性梁动力学特性进行研究. 基于Gauss-Lobatto 节点与Chebyshev 多项式方法对柔性梁变形场进行离散,通过投影矩阵法施加固定及弹性连接边界条件. 利用Chebyshev 谱方法获得了系统固有频率和模态振型数值解,通过与有限元方法及加权残余法的比较,验证了方法的有效性. 分析了弹性连接刚度、角速度比率、系统径长比及梁的长细比等参数对系统固有频率及模态振型的影响. 研究发现:由于系统弯曲模态、拉伸模态的频率随各参数的变化规律不一致,将出现频率转向与振型转换现象;随着弹性连接刚度、角速度比率及系统径长比的增大,低阶弯曲模态频率增大并超过高阶拉伸模态频率,随着梁的长细比的增大,低阶拉伸模态频率增大并超过高阶弯曲模态频率.      

Automated numerical simulation of the propagation of multiple cracks in a finite plane using the distributed dislocation method

In this paper, an automated numerical simulation of the propagation of multiple cracks in a finite elastic plane by the distributed dislocation method is developed. Firstly, a solution to the problem of a two-dimensional finite elastic plane containing multiple straight cracks and kinked cracks is presented. A serial of distributed dislocations in an infinite plane are used to model all the cracks and the boundary of the finite plane. The mixed-mode stress intensity factors of all the cracks can be calculated by solving a system of singular integral equations with the Gauss–Chebyshev quadrature method. Based on the solution, the propagation of multiple cracks is modeled according to the maximum circumferential stress criterion and Paris' law. Several numerical examples are presented to show the accuracy and efficiency of this method for the simulation of multiple cracks in a 2D finite plane.     

Variation of velocity profile of jet and its effect on interfacial stability

IntroductionSincenineteencentury ,thestabilityofgas_liquidtwo_phasejethasattractedalotofpeoplefortheoreticalstudybecauseofitswideapplicationsinindustry .Thestabilitybehavioriscloselyrelatedtotheshapeofbasicvelocityprofile,andthevelocityprofilesmeasuredfromexperimentsarenotaccurateenough ,itisworthwhiletostudythevelocitymodelinthenumericalsimulation .Insimplermodels,top_hatprofile[1,2 ]issuggestedasbasicflowfortheinviscidandincompressibleliquidandgas.Tocompareitwithrealisticflow ,Suetal.[3]ass…     

Application of the Jacobi–Davidson method to accurate analysis of singular linear hydrodynamic stability problems

The aim of this paper is to show that the Jacobi–Davidson (JD) method is an accurate and robust method for solving large generalized algebraic eigenvalue problems with a singular second matrix. Such problems are routinely encountered in linear hydrodynamic stability analysis of flows that arise in various areas of continuum mechanics. As we use the Chebyshev collocation as a discretization method, the first matrix in the pencil is nonsymmetric, full rank, and ill‐conditioned. Because of the singularity of the second matrix, QZ and Arnoldi‐type algorithms may produce spurious eigenvalues. As a systematic remedy of this situation, we use two JD methods, corresponding to real and complex situations, to compute specific parts of the spectrum of the eigenvalue problems. Both methods overcome potentially severe problems associated with spurious unstable eigenvalues and are fairly stable with respect to the order of discretization. The real JD outperforms the shift‐and‐invert Arnoldi method with respect to the CPU time for large discretizations. Three specific flows are analyzed to advocate our statements, namely a multicomponent convection–diffusion in a porous medium, a thermal convection in a variable gravity field, and the so‐called Hadley flow. Copyright © 2012 John Wiley & Sons, Ltd.     

Instability and Receptivity of Laminar Wall Jets

Results of eigenvalue analysis based on global and local eigenvalue considerations are presented. A collocation method with the Chebyshev polynomial approximation has been used for the global eigenvalue analysis. The results explain the appearance of a second unstable mode. In the case of real frequencies with Reynolds number R < 381 there is only one unstable mode. This mode coalesces at R≈ 381 with a stable mode. At R > 381 they become separated by interchannging their branches, then the second unstable mode occurs. The receptivity problem has been considered with respect to perturbations emanating from a wall. The results illustrate that high-frequency modes have a stronger response than low-frequency modes. It is shown that the method of expansion in a biorthogonal eigenfunction system and the method used by Ashpis and Reshotko are equivalent with regard to the receptivity problem solution. Received: 29 March 1996 and accepted 26 November 1996     

A Chebyshev collocation method for the Navier–Stokes equations with application to double-diffusive convection

A Chebyshev collocation method for solving the unsteady two-dimensional Navier–Stokes equations in vorticity–streamfunction variables is presented and discussed. The discretization in time is obtained through a class of semi-implicit finite difference schemes. Thus at each time cycle the problem reduces to a Stokes-type problem which is solved by means of the influence matrix technique leading to the solution of Helmholtz-type equations with Dirichlet boundary conditions. Theoretical results on the stability of the method are given. Then a matrix diagonalization procedure for solving the algebraic system resulting from the Chebyshev collocation approximation of the Helmholtz equation is developed and its accuracy is tested. Numerical results are given for the Stokes and the Navier–Stokes equations. Finally the method is applied to a double-diffusive convection problem concerning the stability of a fluid stratified by salinity and heated from below.