用变分法与差分法计算热对流稳定性的研究中的本征值问题.pdfVIP

用变分法和差分法计算热对流稳 定性研究中的本征值问题 酆庆增 (北京大学，力学与工程科学系，湍流研究国家重点实验室) 摘要 本文用变分法和差分法二种方法计算了旋转流体热对流稳定性研究中遇到 的本征值问题。对两者的计算结果和特点进行了比较。 jJ_，1 1．引论 在Boussinesq近似F研究轴对称各器内旋转流体的热对流稳定性，最 后归结为解如下的偏微分方程(1j： 式中”+为无量纲切向小扰动振幅，Re为R掣n。lds数，Pr为Prandfl数，Ra 兢 为Rayleigh数，o有负实部，则表明扰动衰减，不会出现热对流。L为偏 微分算子： L：乓+；昙一乓一丢 (2)～ 。E2 {延 芘2 E2 式中E为无量纲径向坐标I=r／h，I为无量纲轴向坐标t=Wh，h为 容器特征尺度。 作为例子我们考虑圆台形容器，其高为h，顶和底的半径分别为h、2h。 客器绕对称轴以角速度。旋转。顶面温度为T0，底面温度为T1，删面温度 在T．与To间呈线性变化。流动参数为： (3) Re=oh、，v．R丑=Ⅸgh4盯i-丁拶(vK)，Pr=v／K， 式中V为运动学粘性系数．x为导温系数。 我们把v·按算子L的本征函数系展开。下面用变分法和差分法两种方 法求解本征值问题。变分法的特点是用较少项数就能得到大致可靠的结果。 一般第一本征值的误差较小，以后各次本征值的误差会逐渐增加。如增加 项数．则推导、计算工作会急剧增加a差分法直接从微分方程出发进行计 算。只要网格间距足够小，就可以得到相当精确的结果。但每次都要解很 大阶数的代数本征值问题。 2周变分法求解算子L的本征值问题 算子L的本征函数中和本征值x由以下方程决定 L巾=一^巾， (4) 边界条件为 中=O， 当‘=0，‘：l和l=O，‘=2一‘． (5) 方程(4)可以改写为如下的形式 奏l簧]+壶I墓]一；+砭+=。 cs， 按照Courant和Hilbert‘。，以上本征值问题相当于如下泛函的极值问题 川减柏势t㈤91：唰k㈤ 泛函(7)的变分为零，即6 J=O，的必要条件是(7)中的被积函数F(I， ‘，士，中E，由二)满足Ealer方程： 差～；十蒡～；一■。o (8) 容易证明·从(8)可推出(6)。(6)有无穷多个本征值，最大孤征值 趋向无穷大‘二1e在计算实践中最重要的是前面几个本征值和本i正函数。 我们把巾按满足边条件(5)的线性独立函数v。(滓2，2，3……)的完备 集展开 65 十=∑bi■ (9) i 将(9)代入(7)，由泛函变分为零的条件得到 detlAij·^B