第13章-计算流体力学CFD2.pptVIP

第13章 计算流体力学CFD(2) 粘性激波层 如果沿流向没有局部的逆向分离流动，包含流向导数的粘性项（如 ）小到可以忽略，而且流动是定常的，则导出的方程称为抛物化N-S (PNS)方程。 “抛物化”粘性流动(抛物型方程) 粘性激波层 PNS方程的优点在于： （1）比N-S方程更简单，包含的项更少。 （2）可以采用向下游推进的方法求解。 “抛物化”粘性流动(抛物型方程) 粘性激波层 PNS方程的缺点在于：不能计算沿流向存在分离区的粘性流动。因为涉及流向导数的粘性项被忽略了，而这些导数项反映了流动中粘性效应向上游反馈信息的物理机制。 “抛物化”粘性流动(抛物型方程) 3.4.2.3 非定常热传导(抛物型方程) 非定常热传导(抛物型方程) 热传导方程： 是抛物型方程。其中，?为热扩散率，其物理意义是流体吸收热流量的能力。 两个半无限长壁面之间有传热流体，初始时，壁面和流体的温度都为T1 非定常热传导(抛物型方程) 假设t=0时刻，右边壁面温度突然增加到T2，流体温度发生非定常变化，其瞬时温度分布由热传导方程确定。 非定常热传导(抛物型方程) 右图显示了不同时刻温度在x方向上的分布。 非定常热传导(抛物型方程) 3.4.3 椭圆型方程 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 对于椭圆型方程，特征线是虚数，所以与特征线有关的解法对于椭圆型方程来说是不适用的。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 对于椭圆型方程，没有有限的影响区域和依赖区域，信息可以向任何方向传播到任何地方。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 如右图所示，平面上一点P，假设问题的定义域是abcd所围阴影区域，P点位于abcd内。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 P点能影响阴影区域内的每一个点。反过来，P点的解要受整个封闭边界abcd的影响。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 P点的解需要与流场中所有的点同时求解。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 椭圆型方程的求解不能采用双曲型或抛物型方程所能采用的“推进”解法。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 因为区域内的解依赖于所有边界，必须在整个边界abcd上给定边界条件。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 这些边界条件可以有以下几种形式： 1）在边界上指定未知函数u和v，这种边界条件称为Dirichlet条件。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 这些边界条件可以有以下几种形式： 2）在边界上指定未知函数的导数，例如?u/?x，这种边界条件称为Neumann条件。 椭圆型方程 与二维椭圆型方程的解有关的区域和边界 这些边界条件可以有以下几种形式： 3）Dirichlet条件和Neumann条件的混合条件。 3.4.3.1 定常亚声速无粘流动(椭圆型方程) 定常亚声速无粘流动(椭圆型方程) 绕翼型的低速亚声速无粘流动 在亚声速流动中，扰动以声速传播，理论上可以向上游无限制地传播。 绕翼型的低速亚声速无粘流动 在无粘的亚声速流动中，有限的扰动可以沿着任何方向传播到无穷远。 定常亚声速无粘流动(椭圆型方程) 双曲型方程 由双曲型方程决定的流场可以“推进”求解。计算从给定的初始条件（如y轴）开始，沿着x轴方向一步步推进，逐步求解流场。 与双曲型方程的解有关的区域和边界（二维定常流） 3.4.1.1 定常无粘超声速流动(双曲型方程) 定常无粘超声速流动(双曲型方程) 对于欧拉方程，如果当地马赫数是超声速时，不论方程写成守恒形式还是非守恒形式，定常流动的方程都属于双曲型。 特征线法的初值线（二维定常流） 主流方向为x方向，在物体上游的ab线作为初值线，从初值线开始沿着x方向逐步向下游推进，数值求解控制方程。 特征线法的初值线（二维定常流） 定常无粘超声速流动(双曲型方程) 对于三维定常超声速无粘流动，特征曲线将变成特征曲面。 与双曲型方程的解有关的区域和边界（三维定常流） 定常无粘超声速流动(双曲型方程) 取流动的方向作为X方向。 与双曲型方程的解有关的区域和边界（三维定常流） 定常无粘超声速流动(双曲型方程) 过P点向前的特征面所包含的阴影区域为P点的影响区域。 与双曲型方程的解有关的区域和边界（三维定常流） 定常无粘超声速流动(双曲型方程) 过P点向后的特征面在yz平面截取的区域（用交叉线表示）内的初值才能对P点有影响。 与双曲型方程的解有关的区域和边界（三维定常流） 定常无粘超声速流动(双曲型方程) 可以从YZ平面内给定的初值开