谱元法求解不可压缩流体流动和传热问题及其柱坐标系下线性稳定性分析.pdf

Spectral Element Method for Incompressible Fluid Flow and Heat Transfer Problems and Linear Stability Analysis in Cylindrical Coordinates A Thesis Submitted to Chongqing University in Partial Fulfillment of the Requirement for the Doctor ’s Degree of Engineering By Zhouhua Qiu Supervised by Prof. Zeng Zhong Specialty: Fluid Mechanics College of Resource and Environmental Sciences of Chongqing University, Chongqing, China September 20 13 中文摘要 摘 要 流体流动的稳定性一直是流体力学的中心问题之一。研究流动的稳定性，就 是研究在什么条件下流动会从一个流动状态转换成另外一个流动状态。线性稳定 性分析是研究流动第一次对流失稳比较有效的方法。目前，线性稳定性分析的数 值方法主要集中于有限体积法、有限差分法、有限单元法和谱方法。谱方法因为 其高精度在线性稳定性分析中有广泛应用。利用谱元法进行线性稳定性分析的研 究则较少。谱元法不仅具有谱方法的精度，而且能够很好地通过区域划分技术避 免谱方法在柱坐标和极坐标系统下其极点附近的节点聚集问题。本文研究了不同 坐标系下的谱元法，并最终将其应用到柱坐标系下的线性稳定性分析。本文的主 要工作有： （1）介绍谱离散的理论基础，通过求解直角坐标系的 Poisson-型方程介绍谱 元法的具体算法。 （2 ）提出了极坐标系下的Fourier-Legendre 谱元方法，此方法在周向上采用 Fourier 展开离散，在极点所在单元的径向采用Legendre-Gauss-Radau 积分点离散， 其它单元的径向采用Legendre-Gauss-Lobatto 积分点，有效地消除了极点处的坐标 奇异性，并通过区域分解技术控制节点的分布，避免节点在圆心附近聚集。 （3 ）使用谱元方法作为空间离散结合时间分裂法给出了求解非定常不可压缩 Navier-Stokes(N-S)方程的具体过程，将其推广到极坐标系和柱坐标系下。 （4 ）提出了同位网格下直接求解定常不可压缩N-S 方程的P ×P 谱元法， N N–2 基本思想是消去方程中的时间项，采用 Picard 线性化迭代，给出一种求解二维定 常不可压缩N-S 方程的P ×P 谱元迭代法，方法中压力展开阶数比速度低两阶， N N–2 不仅满足了LBB 条件避免了压力振荡现象，并且不会因使用交错网格导致的插值 计算而影响精度。 （5 ）将谱元法应用到柱坐标系下线性稳定性分析，研究流体流动的第一次失 稳。使用基于Picard 线性化的P ×P 谱元法求解得到不可压缩N-S 方程的基态解， N N–2 基态解加入扰动代入三维方程推导得到柱坐标系下基于谱元法的广义特征值问 题，最后利用ARPACK 程序求解得到实部最大的特征值，从而寻找出流动的相关 控制参数的失稳临界值，利用能量分析探索流动线性失稳的物理机制。 上述工作得到的结果表明：