流体计算理论基础.docVIP

流体计算理论基础 1三大基本方程 1.1连续性方程 连续性方程也称质量守恒方程，任何流动问题都必须满足质量守恒定律，该定 律可表示为：单位吋间若流体不可压缩，密度为常数，于是： u v w 0 x y z 若流体处于稳态，则密度不随吋间变化，可得出： (u) ( v) ( w) 0 X y z 1.2动量守恒定律 该定律可以表述为：微兀体屮流体的动量对吋间的变化率等于外界作用在该微 元体上的各种力之和，该定律实际是牛顿第二定律，按照这一定律，可导出X, y和z三个方向上的动量守恒方程： (u) p XX yx zx ? div( uu) . Fx x x y z t P xy yy zy (u)div( uv) ,Fy t y x y z (u) p xz yz zz . div( uw) . Fz t z x y z 式中，p为微元体上的压力，xx, xy和xz等是因分子粘性作用而产牛的 作用在微元体表 面丄的粘性应力 的分量。Fx, Fy和Fz是微元体上的体力，若体力只有重力, 且z轴竖直向上，则：Fx O,Fy 0, Fz go 对于牛顿流体，粘性应力与流体的变形率成比率，有： u u v 2 div(u); = = .()xyyx xx x y x v u w yy 2 . div(u); xz= zx=.() x z x w v w 2 div(u); = = ()yzzy zz x z y 其中， 为动力粘度， 为第二粘度，一般可取 ~2,将上式代入前式中为: 3 (u) p . div( uu) div( gradujSu t x p ( v)div( uv) div( gradv)Sv y t ( w) p . div( uw) div( gradw)Sw z t 其中： grad() ()/ x ()/ y ()/ z 2 为动力粘度(dynamic viscosity), 为第二粘度(second viscosity),一般J取： (参考3 文献：H.Schlichting,Boundary Layer Theory,8th ed,McGraw Hill, New York, 1979)□ Su,Sv和Sw为动量守恒方程中的广义源项,Su FxSx, Sv FySy, Sw FzSz, 而其中Sx, Sy和Sz表达式为: U V w S ( .)( .)( .)( div(u)) xxx y x x x x U V w S ( .)( .)( .)( div(u)) y x x y y x y y u V w S (.)( .)(.)(div(u)) z X z y z X z z 一般来讲, Fx, Fy和Fz是体积力在 x, y, ，z方向上的分量。 Sx, Sy 和 Sz 是 小量，对于粘性为常数的不可压缩流体，Sx Sy Sz 0,动量守恒，简称动量 方程，也称N-S方程。关于牛顿体与非牛顿体的定义如下： 流体的若 为常数，则该类流体为牛顿流体，否则为非牛顿体，空气，水等 均为牛顿体；聚合物溶液，含有悬浮粒杂质或纤维的流体为非牛顿体。 对于牛顿流体，通常用和［质量］密度的比值代替动力粘度 称为运动粘度，单位m2/s。 1.3能量守恒方程 该方程可以描述为：微元体屮能量的增加率等于进入微元体的净热流量加上体 力与面力对微元体所做的功，实际为热力学第一定律。 (T)k . div( uT) div(gradT)sT tcp k为流体传热系数，cp为比热容，T为温度，ST为流体内热源及由于粘性作用 流体机械能转换为热能的部分，有时简称ST为粘性耗散项。 以JL三大基本方程参考： 《计算流体动力学分析：CFD软件原理与应用》—王福军 2通用控制方程 上面的基本方程可以写成下面的通用形式： TOC \o 1-5 \h \z ().div( u ) div( grad jS t 展开为： ()(u ) ( V ) ( w .) t x y z (.)(.)(,)S x x y y z z 其屮 为通用变量，可以代表u, V和w以及T等求解变量， 为广义扩散系 数，S为广义源项。 2几种数值求解方法 2」有限差分法 主要的思路是用差商代替微商，来近似的表示微分方程.其形式简单，对任意复杂 的偏微分方程都可以写成其对应的差分方程,但是微分方程屮各项的物理意义和 微分方程所反映的物理定律(如守恒定律)在差分方程屮所表现的特点，在差分方 程中没有得到体现.只是微分方程的数学近似，没有反映物理特性，计算结果可能 表现出某些不合理现彖. 2.2有限元法 20世纪60年代出现，离散方程获得的方法主要有:直接刚度法,虚功原理推导，泛 函原理推导或加权余量法推导. 有限元法的优点是解题能力强,可以较精确的模拟各种复杂的曲线或