资源描述：
《流体计算理论基础.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、流体计算理论基础1三大基本方程1.1连续性方程连续性方程也称质量守恒方程，任何流动问题都必须满足质量守恒定律，该定律可表示为：单位吋间若流体不可压缩，密度为常数，于是：uvw0xyz若流体处于稳态，则密度不随吋间变化，可得出：(u)(v)(w)0Xyz1.2动量守恒定律该定律可以表述为：微兀体屮流体的动量对吋间的变化率等于外界作用在该微元体上的各种力之和，该定律实际是牛顿第二定律，按照这一定律，可导出X,y和z三个方向上的动量守恒方程：(u)pXXyxzx・div(uu).FxxxyztPxyyyzy(u)div(uv),Fytyxyz(u)pxz
2、yzzz.div(uw).Fztzxyz式中，p为微元体上的压力，xx,xy和xz等是因分子粘性作用而产牛的作用在微元体表面丄的粘性应力的分量。Fx,Fy和Fz是微元体上的体力，若体力只有重力,且z轴竖直向上，则：FxO,Fy0,Fzgo对于牛顿流体，粘性应力与流体的变形率成比率，有：uuv2div(u);==.()xyyxxxxyxvuwyy2.div(u);xz=zx=.()xzxwvw2div(u);==()yzzyzzxzy其中，为动力粘度，为第二粘度，一般可取~2,将上式代入前式中为:3(u)p.div(uu)div(gradujSutx
3、p(v)div(uv)div(gradv)Svyt(w)p.div(uw)div(gradw)Swzt其中：grad()()/x()/y()/z2为动力粘度(dynamicviscosity),为第二粘度(secondviscosity),一般"J取：(参考3文献：H.Schlichting,BoundaryLayerTheory,8thed,McGrawHill,NewYork,1979)□Su,Sv和Sw为动量守恒方程中的广义源项,SuFxSx,SvFySy,SwFzSz,而其中Sx,Sy和Sz表达式为:UVwS(.)(.)(.)(div(u)
4、)xxxyxxxxUVwS(.)(.)(.)(div(u))yxxyyxyyuVwS(.)(.)(.)(div(u))zXzyzXzz一般来讲,Fx,Fy和Fz是体积力在x,y,，z方向上的分量。Sx,Sy和Sz是小量，对于粘性为常数的不可压缩流体，SxSySz0,动量守恒，简称动量方程，也称N-S方程。关于牛顿体与非牛顿体的定义如下：流体的若为常数，则该类流体为牛顿流体，否则为非牛顿体，空气，水等均为牛顿体；聚合物溶液，含有悬浮粒杂质或纤维的流体为非牛顿体。对于牛顿流体，通常用和［质量］密度的比值代替动力粘度称为运动粘度，单位m2/s。1.3能量
5、守恒方程该方程可以描述为：微元体屮能量的增加率等于进入微元体的净热流量加上体力与面力对微元体所做的功，实际为热力学第一定律。(T)k.div(uT)div(gradT)sTtcpk为流体传热系数，cp为比热容，T为温度，ST为流体内热源及由于粘性作用流体机械能转换为热能的部分，有时简称ST为粘性耗散项。以JL三大基本方程参考：《计算流体动力学分析：CFD软件原理与应用》—王福军2通用控制方程上面的基本方程可以写成下面的通用形式：().div(u)div(gradjSt展开为：()(u)(V)(w.)txyz(.)(.)(,)Sxxyyzz其屮为通用
6、变量，可以代表u,V和w以及T等求解变量，为广义扩散系数，S为广义源项。2几种数值求解方法2」有限差分法主要的思路是用差商代替微商，来近似的表示微分方程.其形式简单，对任意复杂的偏微分方程都可以写成其对应的差分方程,但是微分方程屮各项的物理意义和微分方程所反映的物理定律(如守恒定律)在差分方程屮所表现的特点，在差分方程中没有得到体现.只是微分方程的数学近似，没有反映物理特性，计算结果可能表现出某些不合理现彖.2.2有限元法20世纪60年代出现，离散方程获得的方法主要有:直接刚度法,虚功原理推导，泛函原理推导或加权余量法推导.有限元法的优点是解题能力
7、强,可以较精确的模拟各种复杂的曲线或曲而边界,网格划分比较随意，可以统一•处理多种边界条件，离散方程形式规范，便于编写通用程序，但在应用流体流动和传热原理屮却遇到了一些困难,其原因可归结为按加权余量法推导出的有限元离散方程也只是对原微分方程的数学近似，当处理流动和传热问题的守恒性，强对流，不可压缩条件等方面的要求吋，有限元离散方程屮的各项还无法给出合理的物理解释，对计算屮出现的一些误差也难以进行改进，因此有限元法在流体力学和传热学屮的应用还存在一些问题.2.3有限体积法FVM（或控制体积法CVM）将计算区域划分为网格，并使每个网格点周围有一个互不重
8、复的控制体积，将待解微分方程对每一个控制体积积分，从而得到一组离散方程，其屮的未他数是网格上的因变量，为了求出控制体积的积
