Orr-Sommerfeld方程临界Reynolds的数值求解

摘要

本文采用两点中心差分格式的BOX方法,以分离特征函数实部虚部方式对O-S方程做了数值离散.采用块三角法求解O-S方程,使用牛顿法迭代得到系数矩阵的特征关系。以平板Poiseuille流动和平板Blasius流动为例,应用BOX方法,求解了对应的中性稳定性曲线.根据方程迭代计算过程中的“零”迭代现象,给出了一个简单判断临界Reynolds的方法。