计算流体力学作业习题

12014级西安理工大学计算流体力学作业1.写出通用方程，并说明其如何代表各类守恒定律。由守恒型对流-扩散方程：()divUdivTgradSt其中为通用变量；T为广义扩散系数；S为广义原项。若令1;1;0TS时，则得到质量守恒方程（massconservationequation）()()()()0uvwtxyz若令;iu时，则得动量守恒方程（momentumconservationequation）以x方向为例分析，设;uPuSSx，通用方程可化为：()()()()(2)uuuvuwuPudivUtxyzxxxzvuuwFyxyzzx同理可证明y、z方向的动量守恒方程式若令;;TpTTSSC时，则得到能量守恒方程(energyconservationequation)()()hhdivUhdivUdivgradTSt()()TphdivUhdivgradTStC证毕2.用控制体积法离散0)(sdxdTkdxddxdTudtdT，要求对S线性化，据你的理解，谈谈网格如何划分？交界面传热系数何如何计算？边界条件如何处理？根据守恒型对流-扩散方程：()()()uTStxxx，对一维模型进行分析，则有：0)(sdxdTkdxddxdTudtdT2将该一维模型的守恒形式在图A所示的控制容积P在△t时间内做积分。图A()()()()()ettttetttewewwttwTTTTdxuTuTdtKdtSdsdtxx（1）非稳态项选定T随x变化且为阶梯式，既有：()()ettttttPPwTTdxTTx（2）对流项选定T随t的变化规律符合阶梯显示，既有：()()()()ttttewewtuTuTdtuTuTt（3）扩散项()()()()ttttewewtTTTTdttxxxx()()EPeeTTTxx()()PwwTTTwxx（4）原项令S对t和x呈阶梯式变化，既有：ttettwSdsdtSxt综上所述，可以推导出下式：2()()22tttttttttEwEPwPPuuTTTTTKStxx3由图A可知，本次网格划分采用的是外节点法结构化网格划分。对于交界面的传热系数的数值确定，可根据算术平均法（arithmeticmean），在图B中在P、E两点间的λ与x构成线性关系，则可由P,E两点的λ值，确定在e点的传热系数λ值的大小。即：()()()()eeePEeexxxx在计算求解是，若边界为第一边界则可以直接进行迭代计算，若边界为第二、三边界（边界节点的数据为未知数），则采用附加原项计算法进行求解。3.用幂函数格式离散三维通用方程。在直角坐标系下，三维通用方程的离散方程可表述为：PPEEwwaaa4.采用有限体积法离散对流——扩散方程中的对流项时，根据你的理解写出格式的进化过程。由《数值传热学》知，对流-扩散方程表达式：2jjjjuuStuxx其中jjuu为对流项；2jjuxx为扩散项。现以一维对流-扩散方程问题模型方程来阐述对流项格式演变进化过程。()()dddudxdxdxΓ为了分析数值传热问题，人们最早先提出了控制体积中心差分法，即在P点控制容积处做积分，取分段线性型线，最终可演化得：pPEEWWaaa12EeaDFe12WwaDFw()PEWaaaFeFw该类方程的优点在于，连续性方程在数值计算过程中始终得到满足，系数Ea、Wa包括了扩散和对流作用对热传导问题的影响；与流量有关的部分则是界面上分段线型在均匀网格下的表现，很好地体现了对流作用。但是当P2后，中心差分所解得的解将会失去物理意义，因为当P2时，则Ea2，又因为EaWaPa三个系数的值都应当大于零，故在这种情况下使用中心差分格式将会使得计算存在问题。4为了克服由于对流项因为采用中心差分算法引起的问题，进一步提出了对流项的迎风格式算法，在该格式中对流项的一二阶导数均为线性的型线，同时一阶迎风格式离散方程系数Ea、Wa永远大于零，因而无论在何种条件下计算都不会引起解得震荡，其解永远具有物理意义。并且在迎风格式的使用实践，也能为构造更优良的结构网格提供了启示和指导。5.简述压力校正法的基本思想及过程（用详细的方程离散说明）。压力校正法的基本思想：在对于Navier-Stokes方程的离散形式迭代求解的任一上层次，可以给定一个压力场，它可以给是假定的或是上一层次计算所得出的。一个给定的正确的压力场应该使得计算得到速度场满足连续性方程。但是根据这样的给定的压力场计算而得到的速度场，未必能满足连续性方程，因此要对给定的压力场进行修正。图A在时间间隔△t内对主控制体（如上图A所示）做积分，且以0PPt代替t，采用全隐格式，可得：0()()()()0PPewnsxyuuyvvxt将改进后的速度式*//*//()()eeePEnnnPEuudppuudpp代入整理得关于P一阶导数的代数方程：/////PPEEWWNNSSapaPaPaPaPb其中：,,,EeePEWNsaaaaa0****()()()()()PPwesnbuuyvvxt即压力修正算法可以归纳为以下4个基本步骤：5（1）假定一个压力场，记为*P；（2）利用*P，求解动量离散方程，得出相关的速度**uv；（3）利用质量守恒方程来改进压力场，并要求改进后的压力场对应的速度场能满足连续性方程要求；（4）以*/()PP以及*/()uu，*/()vv作为本层次的解并据此开始下层次的计算迭代。6以具体方程式为例详细说明离散方程的迁移特性的概念。我们将中心差分应用于一维非稳态纯对流方程的非守恒形式：0utx有：其中流速u为常数。采用类似的分析方法，对于节点位于1112nnnniiiiutx（i+1）在（n+1）时层有：11122nnnniiiiutx其中：120,nnii所以11(),2niutx而在i-1点处则有：11122nnnniiiiutx因为120,nnii于是得到11(),2niutx。可见在i点的扰动同时沿着相反的两个方向传递，所以对流项的中心差分不具有迁移性。下面对u0的情况来进行分析。对节点i+1，在n时层产生在节点i的扰动对i+1点的影响由下式确定：111111,(0)2nnnnniiiiiutx由此可得11(),niutx而在i-1处则有611122nnnniiiiutx得110ni可见采用一阶迎风格式时，扰动仅仅向着流动的方向传递，故一阶迎风格式具有迁移性。7．以具体方程式为例详细说明离散方程的守恒性的概念。为了便于分析现将一维对流-扩散方程简化为纯对流方程：()0utt再将方程离散为显式格式，然后在一定大小范围内求和。为了讨论书写简便故将对流项中的时间标记删去。11111()2nniiiiiiuutx在如下图所示的均匀网格系统中，任取一段有限区间进行分析，得：2211()lllludxdxtt2211111()()2nnIIiiiiiIiIuutx或2211111()()(){}2IInniiiiiIiIuuxt进一步分析可得：2110()[()()]InniiinutiIxuut上式表明在△t时间内流入与流出某区域中的通量之差等于改时间间隔内该区域中的增量，又由守恒性质可得：2111[()()]()()IiuiinoutiIuuuu78．详细说明差分格式的相容性和收敛性的概念。以一维稳态对流-扩散方程为例，用符号,()inL表示对函数在点（I,n）作某些微分运算的算子。2,,2()()ininLuStxxΓ其中,()0inL是节点（I,n）处的一维模型方程。用符号,()nxtiL表示对ni作某些差分运算的算子，例如：11111,22()2nnnnnnnnniiiiiiixtiiLuStxxΓ于是,()0inL就代表了一维模型方程的显式格式。所谓一个离散方程的截断误差是指其差分方程算子与相应的微分算子的差，记为TE,即：,()()nxtiinTELL在通过Talor级数展开得：2,()()(,)nxtiinTELLOtx由此可见，当时间空间的网格步长趋于零时，如果离散方程的截断误差趋于零，则称该离散方程与微分方程相容。同样在点（I,n）处也存在离散误差：(,)nniiin当时间空间的网格步长趋于零时，如果离散方程的离散误差趋于零，则称该离散方程是收敛的。89．就计算流体力学的内容与应用问题，谈谈自己的一些想法。