球形Couette流多重周期解的数值的研究.pdfVIP

球形Ceuette流多重周期解的数值研究 293 球形Couette流多重周期解的数值研究。 袁礼 (中国科学院数学与系统科学研究魔计算数学所 科学与工程计算国家重点实验室，北京100080) 摘要本文对内球旋转、外球静止的两同心球之间的不可压缩流进行了数值模拟。 研究了口=006和0．18两个缝隙比情况，雷诺数范围为发生第一次失稳的临界点 至湍流流态之间。通过采用不同的初始条件和施加适当的波形扰动，对口=0．06得 到了多重周期艉，其中一种解的赤道每侧有3对螺旋形Taylor-G6rtler(TG)涡，和 Nakabayashi等人的实验及Dulmas等人的计算结果符合，另外两种新得到的解分别 有2对和1对螺旋形TG涡。还得到了关于赤道不对称的螺旋形rig涡。对口= 0 18得到了有剪切波和Stuart涡的三种周期性流态。研究了这些流态随Re数的演 化。用扰动的基本频率、波数和空间结构刻画了多重周期流态。 关键词球形Couette流非惟一解螺旋形Taylor-C嘉rtler涡。 1引言 两同心球之间的球形Couette流在从层流到湍流的转变过程中蕴涵着丰富多彩的流动结 构和失稳机理。球壳几何上可以看成由极轴附近的平行圆盘和赤道附近的同心圆柱组成，因 此球形Couette流在极轴附近类似于旋转圆盘之间的流动，而在赤道跗近则类似于经典的圆 柱Couette流。球形Couette流的研究既有研究旋转流动的层流一湍流转变机理的理论意义， 又和天体物理、地球物理和工程应用密切相关。1J。 本文考虑内球旋转、外球不动的情况。此时有两个控制参数：雷诺数Re=ntRi／v和缝 隙比p=(R2一R1)／R】，其中v为运动黏性系数，n1为内球旋转角速度，R】和R2分别为内球 和外球的半径。不同的流动结构(如螺旋型TG涡、绕旋波)和它们的特征(波数、基频和旋转 频率)强烈地依赖于缝隙比∥2，“。以轴对称螺圈形TG涡为特征的Taylor不稳定性出现在 不稳定性出现在宽缝情况(口0．3)[4，“。球形Couette流引入注目的特征是具有丰富的涡、 波结构，其中有些是同一＆数下的多重解，包括多重定常解和多重周期解。多重解与封闭域 中Navier—Stokes方程的分叉解有关。实验研究已经发现了多重定常TG涡、多重剪切波和绕 J和连续法计算得 旋波[6～…。数值研究方面，多重定常TG涡已经由时间推进法[”。2 到㈦l“。最近又模拟得到非轴对称流态，如螺旋形TG涡和螺旋波【5’”。“．Sha和 形Couette流多重周期流态的数值模拟尚未见报道。 o国家重点基础规划(G1999032801)和自然科学基金(资助项目。 294近代空气动力学研讨会论文集 本文的目的是研究窄缝和中缝双球的多重周期解[6，…。数值模拟发现缝隙比口=0．06 的窄缝在适当Re数范围至少有三重周期流态，分别对应一对、二对和三对螺旋TG涡(s。： 1，2，3)。Nakabayashi的实验13 将表明多重周期流态很容易从施加波形扰动获得。对缝隙比为J8=018的中缝，得到了剪切 波、Stuart涡和波形外涌结构，和Wimmer的实验基本上符合[…。 数值模拟使用了双时间步推进的人工压缩性方法。差分格式为三阶迎风紧致差分。离散 方程用对角化形式的ADI格式求解。关于方法的细节详见文献[18，12]。本文主要研究第 一次失稳点附近的多重周期流态。 2数学公式和数值方法 考虑两同心球之间的不可压缩流。在一般曲线坐标系下的人工压缩性N—S方程为 等鸲罄+警≯+掣+警≯一o， ㈤ 其中Q=Q／J=(P，“，u，)。／J，Jm=dag(O，l，l，1)7，J为坐标变换的雅可比，P为压 力，“，。和w为直角坐标速度分量，r和t分别为拟时间和物理时间宜，p，e为无黏通量，矗。 p，，O。为黏性通量(如见文献[19])。注意每一物理时间步内需要进行拟时间步的推进，使 得人工压缩性项趋于零以满足零散度条件。 拟时间导数离散采用向后Euler差分，物理时间导数离散采用向后三点二阶精度差分，左 端项的无黏通量采用一阶通量差分裂，离散后的方程用对角化ADI格式求解(例如见文献 『18])