计算流体力学作业习题

1、2014级西安理工大学计算流体力学作业1.写出通用方程，并说明其如何代表各类守恒定律。由守恒型对流-扩散方程：其中为通用变量；为广义扩散系数；为广义原项。若令时，则得到质量守恒方程（mass conservation equation）若令时，则得动量守恒方程（momentum conservation equation）以x方向为例分析，设，通用方程可化为：同理可证明y、z方向的动量守恒方程式若令时，则得到能量守恒方程(energy conservation equation)证毕2.用控制体积法离散，要求对S线性化，据你的理解，谈谈网格如何划分？交界面传热系数何如何计算？边界条件如何处理？

2、根据守恒型对流-扩散方程： ，对一维模型进行分析，则有：将该一维模型的守恒形式在图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、界为第二、三边界（边界节点的数据为未知数），则采用附加原项计算法进行求解。3.用幂函数格式离散三维通用方程。在直角坐标系下，三维通用方程的离散方程可表述为：4.采用有限体积法离散对流扩散方程中的对流项时，根据你的理解写出格式的进化过程。由数值传热学知，对流-扩散方程表达式：其中为对流项；为扩散项。现以一维对流-扩散方程问题模型方程来阐述对流项格式演变进化过程。为了分析数值传热问题，人们最早先提出了控制体积中心差分法，即在P点控制容积处做积分，取分段线性型线，最终可演化得： 该类方程的优点在于，连续性方程在数值计算过程中始终得到满足，系数、包括了扩散和对流作用对热传导问题的影响；与流量有关的部分

4、则是界面上分段线型在均匀网格下的表现，很好地体现了对流作用。但是当2后，中心差分所解得的解将会失去物理意义，因为当2时，则0的情况来进行分析。对节点i+1，在n时层产生在节点i的扰动对i+1点的影响由下式确定：由此可得 而在i-1处则有得 可见采用一阶迎风格式时，扰动仅仅向着流动的方向传递，故一阶迎风格式具有迁移性。7以具体方程式为例详细说明离散方程的守恒性的概念。为了便于分析现将一维对流-扩散方程简化为纯对流方程：再将方程离散为显式格式，然后在一定大小范围内求和。为了讨论书写简便故将对流项中的时间标记删去。在如下图所示的均匀网格系统中，任取一段有限区间进行分析，得：或 进一步分析可得：上式表明在t时间内流入与流出某区域中的通量之差等于改时间间隔内该区域中的增量，又由守恒性质可得：8详细说明差分格式的相容性和收敛性的概念。以一维稳态对流-扩散方程为例，用符号表示对函数在点（I,n）作某些微分运算的算子。其中是节点（I,n）处的一维模型方程。用符号表示对作某些差分运算的算子，例如：于是就代表了一维模型方程的显式格式。所谓一个离散方程的截断误差是指其差分方程算子与相应的微分算子的差，记为TE,即：在通过Talor级数展开得：由此可见，当时间空间的网格步长趋于零时，如果离散方程的截断误差趋于零，则称该离散方程与微分方程相容。同样在点（I,n）处