1、自由对流热转移反应从垂直加热板到地表热通量的振荡马侯赛因和 SK 达斯，孟加拉国达卡和 DAS 里斯，英国巴斯（收稿，1996 年 6 月 3 日1996 年 7 月 12 日修订）总结进行调查的二维不稳定层流沿半无限的竖直板和对流边界的粘性不可压缩流体的层流平均地表热通量关于稳定姿态的小幅度振荡。浮力是有利的，从积极的热通量板表面的流体与往常边界层流的时间周期热通量的相互作用满足审查线性理论。解决方案是利用三种不同的方法，即扩展的一系列扩大低频，高频渐近级数展开法和一般频率的数值有限差分法已经进行了广泛的参数计算，以便找到在波动的幅度和相位角的解决方案。人们已经发现，振幅和相位角，零件表面的
2、剪切应力和表面温度，这三种方法的预测是非常好的在各自有效期的范围。11 介绍在层流边界层理论的领域，莱特希尔1是第一个研究的不稳定平板和圆柱与被迫流动的粘性不可压缩流体的自由流，有小幅度的振荡。自由对流的相似性解决方案从表面均匀热通量与垂直板最初量是研究麻雀格雷格2和获得附近的领先优势有效的解决方案。相应的问题是非定常对流沿垂直板表面温度振荡由纳达和夏尔马3和 Eshghy 等研究。 4。muhuri 和梅蒂5和 Verma6分析了非定常自由对流的振荡表面温度的影响。所有这些调查是基于这样的假设，表面温度执行与时间有关的小振幅振荡意味着温度，和他们进行了通过采用卡门Pohlhausen 近似的
3、积分法。为了获得有效的领先优势，在附近的摄动解下游地区，罗伊7认为是高普朗特数的同类型的问题，威尔克斯8 研究了一个统一的表面热通量垂直板自由对流的问题。在规定的地表热通量的垂直板的情况下，研究由梅尔金和 Mahmood，乔杜里和梅尔金11等人。这些作品被源源不断局限，使解决方案得到了有效的距离，解决方案也给出了中间区域。基于对线性化理论，凯莱赫和杨12研究了层流换热反应逃离对流边界层沿垂直加热板表面温度振荡的问题，最近，侯赛因等调查同类型的问题，详细的表面平均温度，（x），“其中 x 与该板块的领先优势的流向距离成正比。侯赛因等人提出解决方案，在表面热的幅度和相位 ，流向分布频为小型和大型值
4、的传输速率振荡。也有人试图匹配低或高频率的振荡。在相应情况下的自由对流层边界的垂直加热板与非均匀表面热通量还没有被处理之前，应当指出，由于非均匀表面热通量变化可能会出现比物理的表面均匀热通量变化大。重要的是, 要确定具体在多大程度上非均匀表面热通量会影响边界层的响应。目前本文考虑的粘性不可压缩流体的非定常自由对流沿垂直的加热板时，板面的热通量有小幅度的振荡大于平均通量，它本身作为前沿的距离 n 的功率变化。我们调查一般采用（一）延长系列扩展方法低频率范围内，（二）渐近级数展开法在高频率范围，（三）有限差分方法来寻找解决方案的方法。广泛的进行计算参数，以确定表面上的温度和剪切应力的不稳定的影响。
5、2 数学公式在直角坐标系中，一个半无限的垂直板被放置在 y = 0，X_0 的区域，使 x 的测量距离从前沿进入到流体，而且 y 的正常测量从进液板进入到流体。远离表面的环境流体温度为 Too。作为一个积极的表面热通量 q,，结果从板中产生有利的浮力。在 Boussinesq 近似理论下，通常的 Navier-Stokes 方程和能量方程，两三维不可压缩流体，表面热流率是时间相关的情况下，减少到以下边界层方程（Kramer 和排14 ）：其中 u，v 分别为 x 和速度场 y 的组件， 是运动学粘度，T 和 Too 是流体边界层温度和周围流体的温度，g 是重力加速度，fl 是体积膨胀系数， 是
6、热扩散系数。方程（1） - （3）要解决的边界条件传热反应103其中 是振荡频率的表面热通量而 是一个衡量它的幅度。边界条件（4 ）建议式的解决方案。（ 1） - （3 ）可以发现，作为下列表达式的实部（石垣岛）15：其中的组成部分 和 平均流量基本稳定，满足微分方程边界条件u l, vi 和 是非定常流的组成部分，它满足的微分方程受边界条件为了获得相似稳定状态方程（6） - （8），我们将介绍以下组的转换：104 文学硕士侯赛因等。在上面的 满足稳态流的连续性方程，而是 一个平均地表热通量有关的常数。因此，我们获得了方程满足上述方程的边界条件是这里的素数用 表示， 是普朗特数转换（14）带领
7、我们转换以下方程组。（ 10） - （13 ）波动问题的一部分：方程（11）和（12），然后降低到边界条件不稳定的剪应力和表面温度是很难发现的，而这些可以从方程式中得到解决。（15） - （17）和（19） - （21 ）。在这里，我们提出了在振幅和相位的剪应力和表面热流率的解决方案。根据这些定义和 代表的实部和虚部 和 的一部分。传热反应方程（15） - （17）描述了稳定的平均流量和温度场。这些方程的解决方案已经得到乔杜里和梅尔金的11相关的物理参数 pr 和 n 的不同的价值。方程（19） - （21）描述的 波动的解决方案组成部分，这些应通过多种方法解决。在第 3 节我们详细介绍 值，
8、使用了一系列解决方案。由于 和成正比， 见（18），这样的系列解决方案是有效的 x 小的值而 并且为了非常低的频率 与 。第 4 节讨论 的渐近解，这可以作为一个大的距离限制，或解释为高频率的限制。解决方案为中间值 。方程在一般情况下，类似的使用了凯勒盒的方法（见17）。实施这种方法的详情，现在非常标准，并已在13讨论。值得注意的是，当 n =1 的时候。（19）和（20）降低到一对线性常微分方程的解决方案可以通过一个简单的拍摄方法得到的。3 的系列解决方案显然，附近的领先优势，使得在使用结果的基础上，数量有限的影响，忽快忽慢的描述将只在很小的范围内有效。由于 小值也对应非常低的频率，我们期望
9、的流量调整准静态传热边界的波动率。我们扩大 f 和 ，根据代入式（19） - （20），然后等同为像 到零的力量， 我们得到以下对函数 “ ”常微分方程：其中 M = 1，2 ，3，各自的边界条件是.素数再次到 的衍生物。106 文学硕士侯赛因等。可以看出公式（25） - （28 ）是线性的，但加上可独立解决的从一个到另一个。在目前的分析，龙格 - 库塔 - 布彻18的初始值问题求解连同 Nachtsheim Swigert19迭代方案解决式的系统。（ 25） - （28 ）为 。这里的帕德 20也被用来获得更准确地接近了当地的幅度和相位的剪应力和表面热流的波动部分。在第 5 节对详细的数值结
10、果进行了讨论。4 大型的渐近解远离前沿浮力变得越来越重要，直到远离下游的流量将是主要的自由对流，自由流的存在只是略有不安。因此，在本节已经给予了解决方案。（19）和（20） 当大的时候，我们强调的是，这个限制对应的不仅是为 定义的 x 的大值，还有为 x 定义的 的大值。其实，其实，详细检查了凯勒的箱法得到的数值结果显示， 的大值，不稳定的反应，只限于表面附近的薄区域。我们注意到，然而，这个结论没有在更高的电子订单中。因此，我们寻求在高频率范围内的一系列解决方案中，采用零近似的解决方案。出于这个原因，介绍了以下转换：这些换算比例的动机顺序数量级分析（19 ）。然后成为方程（19）和（ 20）和
11、由于这些方程对应的是薄壁层，领先高阶函数的，F 和 ，在这个区域可以表示以下功率具有良好的准确性：传热响应 107其中，根据公式。 （15） - （17 ）基于上述展开式的解决方案。（31 ）和（32 ）可以得到以下形式：当公式（35）代入（31）和（32） 被收集，包括：和质数现在表示关于差异为 Y 相关的边界条件现在解决公式（36）和（ 37），受边界条件（38），我们发现在下面的表达式为 和108 文学硕士侯赛因等。其中 我被评为在第一象限（即 ），并且这里必须指出的是复杂的表达式（39 ）和（40 ）是有效的普朗特数 ，然而，当 是需要的，有时必须采取这些解决方案的限制，为 。5 结果
12、和讨论在目前的分析，一种粘性波动的自由对流流体的解决方案是使不可压缩流体沿垂直受热面小幅度的在地表热通量非均匀稳定的前沿从测量距离的力量变化的热通量振荡。（主导的顺序）稳定问题的解决方案已经由乔杜里和梅尔金11讨论。凯勒在整个频率范围内采用箱法分析了波动问题的部分。该解决方案也已采用渐近法在小频率和大范围的频率使用扰动方法制度。由此获得的结果表示在幅度和相位、剪应力的波动部分以及那些表面温度呈现出不同的表面热通量梯度参数 n 和普朗特数 的影响方面，普朗特数的值被选择代表目前用于核工程冷却剂为液态金属流体（威尔克斯8），例如，锂 0.05 汞 0.01。我们还获得=1.0，0.7，0.1，0.
13、05，0.01 的解决方案。= 1.0 和 n= 1.0 表 1 和表 2 给出的剪应力的波动幅度和相位，通过上述三种方法获得的表面温度的数值。比较结果显示，扰动解和渐近解的差分解决方案很好的吻合。对于 n=0.0，0.25,0.5 和 0.75 =0.7，对应于空气，这些方法获得的数值，这就是上文所述剪应力波动的幅度和相图图形 1 和 2。相应的幅度和相位值表面温度波动图所示 3 和 4。在这些数字的粗曲线。传热反应 109表 1。 =1.0 剪应力的波动幅度和相和 n=1.0a 低频率的系列解决方案b 高频率的渐近解代表凯勒盒的解决方案，圆圈和坚实的盘旋曲线代表的扰动解和渐近解。像以前一样
14、，这些曲线之间的比较确定的扰动解和渐近解与凯勒盒法解决方案在选定的每一个表面热流指数值都非常吻合。图 5 6 和 7 - 8 表示剪应力的波动幅度和相位，以及在 PR=0.7，0.1 ，0.05，0.01 ，N = 0.5 的表面温度。从表 1 和表 2 中，我们可以看到随着频率的增加，无论表面温度梯度规定的普朗特数、剪应力幅值和表面传热，单调减少（从图 1 和 3 中可以看到）。这是由于从表面上看，滞后随频率变化量在相邻流体层的温度滞后。可以看出，在整个频率范围内的表面温度。110 文学硕士侯赛因等。表 2。 =1.0 和 n= 1.0 表面温度波动的幅度和相位a 低频率的系列解决方案b 高频率的渐近解振荡总是导致表面温度波动。相位角， 和 ，在稳态条件下是零，而他们对单调减少渐近值 和 分别为 。我们进一步观察，在低频率范围内的相位角是随着普朗特数的表面温度梯度减小的值（由此可以看出，从图 2 和4）。从图 1 和 3 可以看出，在低频范围，剪应力振幅和波动的表面温度降低时指数表面的热通量增加。同样，从图 2 和 4 可以看出，相角在低频率范围内，随着指数增加的价值，而 值保持常数。传热响应 111
