计算流体力学分析.docVIP

计算流体力学作业 第一部分 用微元法推导出连续性方程，N-S方程，以及牛顿流体的化简形式 解：由质量守恒定律，单位时间内流体微元体中质量的增加，等于同一时间间隔内流入微元体的净质量，得质量守恒方程： 0 即为连续性方程： 0 动量守恒定律：微元体中流体的动量对时间的变化率等于外界作用在该微元体上各种力之和。 以Fx为例，Fx分为体积力（Body forces），与表面力（Surface forces），表面力又分为压力（Pressure）以及粘性应力（Viscous），其中粘性应力分为普通的（如τxx）与剪切力（如τxy，τxy表示在y方向上作用于与x所垂直的平面的力） x方向的体积力为： 在abcd平面上，x方向的表面力为 与abcd相平行，在efgh平面上，x方向的表面力为 注意二者方向相反，综合为，为x轴正方向 同理可得另外两个平行平面的的合力。 另： x方向的压力为： x负方向的压力为 压力合力为x方向： 综上， 微元体的质量： 微元体加速度： 则 将式子左侧变形为全微分形式 其中， 且 则 如此可得x方向的N-S方程： 同理得y方向： Z方向： 对于不可压缩流体： 对于牛顿流体： 将上述式子带入N-S方程，联立可得： 其中，; 2 用有限体积法离散二维对流扩散方程 , 将得出代数方程： 并确定式中系数。 解：对扩散方程在上图控制容积内积分，有：s 由图可知，控制容积的边界面积（长度） 则上述方程式可以写为 计算通过控制容积截面的场变量值或者导数值： 穿过东侧边界的对流量： 穿过东侧边界的扩散量： 穿过西侧边界的对流量： 穿过西侧边界的扩散量： 穿过北侧边界的对流量： 穿过北侧边界的扩散量： 穿过南侧边界的对流量： 穿过南侧边界的扩散量： 源项： 此时，另， 则： 带入通过各界面的扩散量及对流量，然后积分方程变为： 按照节点变量整理，最终可得： 即 式中，,,,, , 2设某场变量经过对流扩散过程从一维区域的x 0点输运到x L点，输运方程为 其中：流体密度 1.0kg/m3，L 2.0m,扩散系数 0.2kg/ m.s 边界条件, 求（1）当流速u 0.1m/s，对对流项采用中心差分格式将区域离散成十个节点时，用有限体积法求在区域内的分布并作图。 （2）当流速u 2.5m/s，对对流项采用迎风格式将区域离散成二十个节点时，用有限体积法求在区域内的分布并作图。 解： 求解区域的离散网格系统如图所示： 每一个控制容积的长度 0.2m,此时,,,,对所有控制容积都成立，中间节点2~9的离散方程为，边界点1与边界点10需要特殊处理。 对节点1所在的控制容积积分，有 （此时点1为P） 则 在西侧界面时，有 则 按照节点场变量整理，有 即 式中 ,, , 当 常数时，表达式中（。 同理，对节点5所在的控制容积积分，得： 按节点场变量整理，有 即 式中,, , 综上，我们得到下表所列的离散方程系数和等效源项计算公式 节点 1 0 D-F/2 - 2D+F 2D+F φA 2到9 D+F/2 D-F/2 0 0 10 D+F/2 0 - 2D-F 2D-F φB 第一种计算情况时，,带入上表，可得离散方程系数 节点 1 0 0.95 1.1 -2.1 3.05 2到9 1.05 0.95 0 0 2.0 10 1.05 0 0 -2.1 3.15 将方程写成增广矩阵形式（未知数为各节点处的值） 3.05 -1.05 0 0 0 0 0 0 0 0 2.1 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2 -1.05 0 0 0 0 0 0 0 0 0 -0.95 2.95 0 用追赶法解此方程组，得： 0.97094 0.820348 0.684099 0.560826 0.449293 0.348382 0.257082 0.174477 0.099739 0.032119 %%%%图有问题 （2）此时 当流动为正方向（从左往右）时，采用迎风格式即为： 此时离散方程为 按照节点场变量整理为： 在边界点1处，对流项采用迎风格式： 同理，在边界点20处： 按照节点场变量整理为的形式，系数值见下表 节点 S_P 1 0 D -（2D+F） 2~1