计算流体力学作业习题浅析.docxVIP

PAGE \* MERGEFORMAT9 2014级西安理工大学计算流体力学作业 1.写出通用方程，并说明其如何代表各类守恒定律。 由守恒型对流-扩散方程： 其中为通用变量；为广义扩散系数；为广义原项。 若令时，则得到质量守恒方程（mass conservation equation） 若令时，则得动量守恒方程（momentum conservation equation） 以x方向为例分析，设，通用方程可化为： 同理可证明y、z方向的动量守恒方程式 若令时，则得到能量守恒方程(energy conservation equation) 证毕 2.用控制体积法离散，要求对S线性化，据你的理解，谈谈网格如何划分？交界面传热系数何如何计算？边界条件如何处理？ 根据守恒型对流-扩散方程： ，对一维模型进行分析，则有： 将该一维模型的守恒形式在图A所示的控制容积P在△t时间内做积分。 图A （1）非稳态项 选定T随x变化且为阶梯式，既有： （2）对流项 选定T随t的变化规律符合阶梯显示，既有： （3）扩散项 （4）原项 令S对t和x呈阶梯式变化，既有： 综上所述，可以推导出下式： 由图A可知，本次网格划分采用的是外节点法结构化网格划分。对于交界面的传热系数的数值确定，可根据算术平均法（arithmetic mean），在图B中在P、E两点间的λ与x构成线性关系，则可由P,E两点的λ值，确定在e点的传热系数λ值的大小。即： 在计算求解是，若边界为第一边界则可以直接进行迭代计算，若边界为第二、三边界（边界节点的数据为未知数），则采用附加原项计算法进行求解。 3.用幂函数格式离散三维通用方程。 在直角坐标系下，三维通用方程的离散方程可表述为： 4.采用有限体积法离散对流——扩散方程中的对流项时，根据你的理解写出格式的进化过程。 由《数值传热学》知，对流-扩散方程表达式： 其中为对流项；为扩散项。 现以一维对流-扩散方程问题模型方程来阐述对流项格式演变进化过程。 为了分析数值传热问题，人们最早先提出了控制体积中心差分法，即在P点控制容积处做积分，取分段线性型线，最终可演化得： 该类方程的优点在于，连续性方程在数值计算过程中始终得到满足，系数、包括了扩散和对流作用对热传导问题的影响；与流量有关的部分则是界面上分段线型在均匀网格下的表现，很好地体现了对流作用。但是当2后，中心差分所解得的解将会失去物理意义，因为当2时，则2，又因为三个系数的值都应当大于零，故在这种情况下使用中心差分格式将会使得计算存在问题。 为了克服由于对流项因为采用中心差分算法引起的问题，进一步提出了对流项的迎风格式算法，在该格式中对流项的一二阶导数均为线性的型线，同时一阶迎风格式离散方程系数、永远大于零，因而无论在何种条件下计算都不会引起解得震荡，其解永远具有物理 意义。并且在迎风格式的使用实践，也能为构造更优良的结构网格提供了启示和指导。 5.简述压力校正法的基本思想及过程（用详细的方程离散说明）。 压力校正法的基本思想： 在对于Navier-Stokes方程的离散形式迭代求解的任一上层次，可以给定一个压力场，它可以给是假定的或是上一层次计算所得出的。一个给定的正确的压力场应该使得计算得到速度场满足连续性方程。但是根据这样的给定的压力场计算而得到的速度场，未必能满足连续性方程，因此要对给定的压力场进行修正。 图A 在时间间隔△t内对主控制体（如上图A所示）做积分，且以代替，采用全隐格式，可得： 将改进后的速度式 代入整理得关于P一阶导数的代数方程： 其中： 即压力修正算法可以归纳为以下4个基本步骤： （1）假定一个压力场，记为； （2）利用，求解动量离散方程，得出相关的速度； （3）利用质量守恒方程来改进压力场，并要求改进后的压力场对应的速度场能满足连续性方程要求； （4）以以及，作为本层次的解并据此开始下层次的计算迭代。 6以具体方程式为例详细说明离散方程的迁移特性的概念。 我们将中心差分应用于一维非稳态纯对流方程的非守恒形式： 有： 其中流速u为常数。 采用类似的分析方法，对于节点位于（i+1）在（n+1）时层有： 其中： 所以 而在i-1点处则有： 因为于是得到。 可见在i点的扰动同时沿着相反的两个方向传递，所以对流项的中心差分不具有迁移性。 下面对u0的情况来进行分析。对节点i+1，在n时层产生在节点i的扰动对i+1点的影响由下式确定： 由此可得 而在i-1处则有 得 可见采用一阶迎风格式时，扰动仅仅向着流动的方向传递，故一阶迎风格式具有迁移性。 7．以具体方程式为例详细说明离散方程的守恒性的概念。 为了便于分析现将一维对流-?