第11讲_相似解法.pptVIP

SIMILARITY SOLUTION相似解法 PARALLEL FLOW OVER A FLAT PLATE平板流问题 何谓解的相似性？ 相似解法的方法及思路 针对二维边界层层流流动问题，速度或温度函数f(x,y)若能转化为f(η)，则η称为相似变量 相似变量η是关于x,y的一个组合 寻找相似变量的方法有许多，可以采用一种直观的方法 下面所采用的这种直观方法是建立在一种对流线或等温线的思考基础之上 边界层流动可看成从板端起头的若干流线或等温线组成 尽管处于板的不同位置(x)，但同一流线或等温线具有相同的无量纲速度或无量纲温度 等流线或等温线可用y/g(x)描述，g(x)从物理意义上可看成与速度边界层或温度边界层厚度相关的函数 故：η=y/g(x) 相似解存在条件及相似变量的求解 温度场相似解存在及相似变量求解举例 Pr1时Nu数表达式 本讲总结 相似解法是一种基于物理背景的数学解法 相似解法是具有坐标变换思想的一种方法 相似解法的坐标与流线和等温线的坐标有一定关系 相似解法将偏微分方程的求解转化为常微分方程的求解 相似解法的结果可通过数值积分最终数字化 Nu数表达式是精确解的近似化结果 * * 测量所得的边界层内速度、温度以及浓度分布如图所示，使人直觉判断所有这些分布本质上具有相似性，如果坐标重新选择，这种相似性就会明显的显现出来。 当边界层方程具有相似性解时，其流速u(x,y)、温度T(x,y)和浓度ω(x,y)的分布具有以下性质：如果把任意x断面的u(x,y)、 T(x,y) 或ω(x,y)采用无量纲坐标表示时，其分布图形均相同。 从数学上说，这就是坐标变换。边界层数学描述，通过这样的坐标变换，会从偏微分方程形式转换为常微分方程形式。当然，方程的非线性性质仍然没有改变。 对于二维定常问题，引入流函数ψ(x,y)，并定义： 如果假设流函数可写成如下形式： 则流速分量为： U为纵掠平板的来流速度 代入层流边界层动量微分方程式中可得： 边界条件为： 可见，只有 为常数时，方程才是常微分方程，即f只是η的函数 若令： 则可解得： 相似解存在的条件 由此可得相似变量为 若壁面温度在板的起始位置开始维持常数，引入无量纲温度： 则能量方程及相应边界条件为： 按照η=y/g(x)的假设，可以有： 代入偏微分方程可得： 很清楚，只有： 此时f 才仅为η的函数，常微分方程才可得解。 由此积分可得 则g(x)有较简单的形式， 从而，相似变量 从所得结果来看，g(x)确实与边界层厚度相关 也可写成 TEMPERATURE DISTRIBUTION SOLUTION 按照这样的相似变量， 能量方程及边界条件变为： 积分后得： VELOCITY DISTRIBUTION SOLUTION 其中， 为无量纲速度 其解称之为著名的布拉修斯解 Blasius 是通过级数展开进行求解的，而后人则是用迭代法或用计算机数值积分法求解 TEMPERATURE DISTRIBUTION 代入速度分布，温度分布可通过数值积分方式得出，如图所示。 换热过程微分方程式改写为 于是 Nusselt Number solution 换热问题的求解最终有赖于对流换热系数的求解，一般表达为无量纲的形式 无量纲的Nu数的求解最终转化为无量纲温度梯度的求解 按照温度分布表达式，可以得出： 此无量纲温度梯度只与Pr数有关 Nusselt Number solution 将不同Pr数下无量纲温度梯度值进行计算，然后拟合，可得无量纲温度与Pr数的对应关系 Nusselt Number solution 拟合可得： 可看出二者的接近程度 两者对比如下表所示： Nusselt Number solution Local Nusselt number From then Pr1 积分 f(0)=0 适用于液态金属，Pr0.05 * * * * *