流体的基本计算
1、液体压强计算计算公式；
液体压强；
在液体容器低、内壁、内部中，由液体所产生的液体压强，简称液压。
2、喷嘴射流速度及流量
深度△Z 液体密度ρ出口直径D 流量系数C
出口速度计算公式；
体积流量计算公式；
质量流量计算公式；
3、限孔流场计算
入口直径Di 出口直径Do 压力差△p 流体密度ρ
入口速度计算公式；
出口速度计算公式；
体积流量计算公式；
质量力量计算公式；
4、运动粘度
运动粘度μ密度ρ
运动粘度计算公式；
运动粘度；
运动粘度即流体的运动粘度与同温度下该流体密度ρ之比。
动力粘度；
Μ动力粘度【Pa。s】或【N。S/m2】或【kg/（m。s）】；也被称为动态粘度、绝对粘度或简单粘度，定义为应力与应变速率之比，其数值上等于面积为1m2相距1m的两平板，以1m/s的速度作为相对运动时，因之存在的流体互相作用所产生的内摩擦力。
5、雷诺数
特征速度v 特征长度L 运动粘度V 动力粘度μ密度ρ
雷诺数计算公式；
雷诺数；
一种可用来表征流体流动情况的无量纲数。利用雷诺数可区分为流体的流动是层流或湍流，也可用来确定物体在流体中流动所受的阻力。
6、韦伯数
流体密度ρ特征速度v 特征长度L 秒面张力σ
韦伯数计算公式；
韦伯数
韦伯数是流体力学中的一个无量纲数。当不同的流体之间有交界面时，尤其在多相流中交界面的曲率较大时，它用来分析流体运动。
7、马赫数
流体速度v
马赫数计算公式；
马赫数；
流体力学中表征流体可压缩程度的一个重要的无量纲参数，定义为流场中某点的速度v同该点的当地声速c之比。
8、水力半径和水力直径
流动截面积A
圆周Pw
水力半径计算公式
水力直径计算公式
水力半径；
是水力学中的一个专有名称，指某输水断面的过流面积与输水断水面和接触的边长（圆周）之比，与断面形状有关，常用于计算渠道隧道的输水能力。
水力直径；
是在关内流动中引入的，其目的是为了给非圆管流动取一个合适的特征长度来计算其雷诺数。常用表达式是；2A/P，即二倍的横截面积（A）除以圆周长度（p）。
流体力学复习要点(计算公式)
D D y S x e P gh2 gh1 h2 h1 b L y C C D D y x P hc 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ， 3 m kg 29.1=空气ρ 牛顿内摩擦定律：剪切力： dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度： ρυ μ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱mm 73610/9800012 ===m m N at 2/101325 1m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D +=α sin 1 ) () 2(32121h h h h L e ++= ρ 若01 =h ，则压强为三角形分布，3 2L e y D == ρ 注：①图算法适合于矩形平面；②计算静水压力首先绘制压强分布图， α 且用相对压强绘制。 （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形12 3 bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力z x P P P += 与水平面的夹角 x z P P arct an =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ? ? ??? ??+??+??+??=??+??+??+??=??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ
《流体力学》典型例题20111120解析
《流体力学》典型例题（9大类） 例1~例3——牛顿内摩擦定律（牛顿剪切公式）应用 例4~例5——流体静力学基本方程式的应用——用流体静力学基本方程和等压面计算某点的压强或两点之间的压差。 例6~例8——液体的相对平衡——流体平衡微分方程中的质量力同时考虑重力和惯性力（补充内容） （1）等加速直线运动容器中液体的相对平衡（与坐标系选取有关） （2）等角速度旋转容器中液体的平衡（与坐标系选取有关） 例9——求流线、迹线方程；速度的随体导数（欧拉法中的加速度）；涡量计算及流动有旋、无旋判断 例10~16——速度势函数、流函数、速度场之间的互求 例17——计算流体微团的线变形率、角变形率及旋转角速度 例18~20——动量定理应用（课件中求弯管受力的例子） 例21~22——总流伯努利方程的应用 例23——综合：总流伯努利方程、真空度概念、平均流速概念、流态判断、管路系统沿程与局部损失计算 例题1：如图所示，质量为m ＝5 kg 、底面积为S ＝40 cm ×60 cm 的矩形平板，以U ＝1 m/s 的速度沿着与水平面成倾角θ＝ 30 的斜面作等速下滑运动。已知平板与斜面之间的油层厚度δ ＝1 mm ，假设由平板所带动的油层的运动速度呈线性分布。 求油的动力粘性系数。 U G=mg δ θ 解：由牛顿内摩擦定律，平板所受的剪切应力du U dy τμ μδ =＝ 又因等速运动，惯性力为零。根据牛顿第二定律： 0m ==∑F a ，即： gsin 0m S θτ-?= ()32 4gsin 59.8sin 301100.1021N s m 1406010m U S θδμ--?????==≈????? 粘性是流体在运动状态下，具有的抵抗产生剪切变形速率能力的量度；粘性是流体的一种固有物理属性；流体的粘性具 有传递运动和阻滞运动的双重性。 例题2：如图所示，转轴的直径d ＝0.36 m ，轴承的长度l ＝1 m ，轴与轴承的缝隙宽度δ＝0.23 mm ，缝隙中充满动力粘性系数0.73Pa s μ=?的油，若轴的转速200rpm n =。求克服油的粘性阻力所消耗的功率。 δ d l n 解：由牛顿内摩擦定律，轴与轴承之间的剪切应力 ()60d d n d u y πτμ μδ =＝ 粘性阻力（摩擦力）：F S dl ττπ=?=
第二章计算流体力学的基本知识
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30～40 年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学" 。 从20 世纪60 年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
流体力学例题
第一章 流体及其主要物理性质 例1： 已知油品的相对密度为0.85，求其重度。 解： 例2： 当压强增加5×104Pa 时，某种液体的密度增长0.02%，求该液体的弹性系数。 解： 例3： 已知：A ＝1200cm 2，V ＝0.5m/s μ1＝0.142Pa.s ，h 1＝1.0mm μ2＝0.235Pa.s ，h 2＝1.4mm 求：平板上所受的内摩擦力F 绘制：平板间流体的流速分布图 及应力分布图 解：（前提条件：牛顿流体、层流运 动） 因为 τ1＝τ2 所以 3 /980085.085.0m N ?=?=γδ0=+=?=dV Vd dM V M ρρρρρ d dV V -=Pa dp d dp V dV E p 84105.2105% 02.01111?=??==-==ρρβdy du μ τ=??????? -=-=?2221110 h u h u V μτμτs m h h V h u h u h u V /23.02 112212 2 11 =+= ?=-μμμμμN h u V A F 6.41 1=-==μ τ
第二章 流体静力学 例1： 如图，汽车上有一长方形水箱，高H ＝1.2m ，长L ＝4m ，水箱顶盖中心有一供加水用的通大气压孔，试计算当汽车以加速度为3m/s 2向前行驶时，水箱底面上前后两点A 、B 的静压强（装满水）。 解： 分析：水箱处于顶盖封闭状态，当加速时，液面不变化，但由于惯性力而引起的液体内部压力分布规律不变，等压面仍为一倾斜平面，符合 等压面与x 轴方向之间的夹角 例2： （1）装满液体容器在顶盖中心处开口的相对平衡 分析：容器内液体虽然借离心惯性力向外甩，但由于受容器顶限制，液面并不能形成旋转抛物面，但内部压强分布规律不变： 利用边界条件：r ＝0，z ＝0时，p ＝0 作用于顶盖上的压强： （表压） （2）装满液体容器在顶盖边缘处开口的相对平衡 压强分布规律： =+s gz ax g a tg = θPa L tg H h p A A 177552=??? ?? ?+==θγγPa L tg H h p B B 57602=??? ?? ?-==θγγC z g r p +-?=)2( 2 2ωγg r p 22 2ωγ =C z g r p +-?=)2( 2 2ω γ
计算流体力学常用数值方法简介[1]
计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 一 前 言 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区
流体力学计算公式
1、单位质量力：m F f B B = 2、流体的运动粘度：ρ μ=v （μ[动力]粘度，ρ密度） 3、压缩系数：dp d dp dV V ρρκ?=?-=11（κ的单位是N m 2）体积模量为压缩系数的倒数 4、体积膨胀系数：dT d dT dV V v ρρα?-=?=11（v α的单位是C K ?1,1） 5、牛顿内摩擦定律：为液体厚）为运动速度，以应力表示为y u dy du dy du A T (,μτμ== 6、静止液体某点压强：为该点到液面的距离）h gh p z z g p p ()(000ρρ+=-+= 7、静水总压力： )h (为受压面积，为受压面形心淹没深度为静水总压力，A p ghA A p p c ρ== 8、元流伯努利方程;'2221112w h g p z g u g p z ++=++ρρ('w h 为粘性流体元流单位重量流体由过流断面1-1运动至过流断面2-2的机械能损失，z 为某点的位置高度或位置水头，g p ρ为测压管高度或压强水头，g u ρ2是单位流体具有的动能，u gh g p p g u 22'=-=ρ，u gh C g p p g C u 22'=-=ρC 是修正系数，数值接近于1） 9、总流伯努利方程:w h g v g p z g v g p z +++=++222 221221111αραρ(α为修正系数通常取1） 10、文丘里流量计测管道流量：)21)(41()()(42 122211g d d d k h k g p z g p z k Q -=?=+-+=πμρρμ 11、沿程水头损失一般表达式：g v d l h f 22 λ=（l 为管长，d 为管径，v 为断面平均流速，g
计算流体力学过渡到编程的傻瓜入门教程
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较方便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分方程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子内面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分方程的时候提出的一个问题，用一维无粘可压缩Euler方程就可以描述了。 这里 这个方程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量 ）随空间变化（）的关系。 在CFD中通常把这个方程写成矢量形式 这里 进一步可以写成散度形式
计算流体力学_CFD_的通用软件_翟建华
第26卷第2期河北科技大学学报Vol.26,No.2 2005年6月Journal of Hebei University of Science and T echnology June2005 文章编号:100821542(2005)022******* 计算流体力学(CFD)的通用软件 翟建华 (河北科技大学国际交流与合作处,河北石家庄050018) 摘要:对化学工程领域中的通用CFD(Computational Fluid Dynamics)模拟软件Phoenics,Flu2 ent,CFX等的具体特点和应用情况进行了综述,指出了他们各自的结构特点、特有模块、包含的数学模型和成功应用领域;给出了选用CFD软件平台的7项准则,对今后CFD技术的发展进行了预测,指出,今后CFD研究的主要方向将集中在数学模型开发、工程改造和新设备开发及与工艺软件的匹配连用等方面。 关键词:计算流体力学;模拟软件;CFX;FLUENT;PH OENICS 中图分类号:T Q015.9文献标识码:A Review of commercial CFD software ZH AI Jian2hua (Department of Int ernation Exchange and Cooperation,H ebei University of Science and Technology,Shijiazhuang H ebei 050018,China) Abstr act:The paper summar izes the features and application of the CF D simulation software like Phoenics,F luent and CFX etc in chemical engineering,and discusses their str ucture features,special modules,mathematical models and successful application areas.It also puts forward seven r ules for the good choice of commercial CF D code for the CF D simulation resea rcher s.Based on t he predict ion of the technology development,it points out the possible r esear ch direction for CF D in the future will focus on the development of mathematical model,project transformat ion,new equipment and their matching application with technologi2 cal softwa re. Key words:CF D;simulation software;CF X;FLUENT;P HOENICS CFD(Computational Fluid Dynamics)软件是计算流体力学软件的简称,是用来进行流场分析、计算、预测的专用工具。通过CFD模拟,可以分析并且显示流体流动过程中发生的现象,及时预测流体在模拟区域的流动性能,并通过各种参数改变,得到相应过程的最佳设计参数。CFD的数值模拟,能使我们更加深刻地理解问题产生的机理,为实验提供指导,节省以往实验所需的人力、物力和时间,并对实验结果整理和规律发现起到指导作用。随着计算机软硬件技术的发展和数值计算方法的日趋成熟,出现了基于现有流动理论的商用CFD软件。这使许多不擅长CFD工作的其他专业研究人员能够轻松地进行流体数值计算,从而使研究人员从编制繁杂、重复性的程序中解放出来,以更多的精力投入到研究问题的物理本质、问题提法、边界(初值)条件和计算结果的合理解释等重要方面上,充分发挥商用CFD软件开发人员和其他专业研究人员各自的智力优势,为解决实际工程问题开辟了道路。 CFD研究走过了相当漫长的过程。早期数值模拟阶段,由于缺乏模拟工具,研究者一般根据自身工作性质和研究过程,自行编制模拟程序,其优点是针对性强,对具体问题的解决有一定精度,但是,带来的问题 收稿日期:2004208221;修回日期:2004211221;责任编辑:张军 作者简介:翟建华(19642),男,河北平乡人,教授,主要从事化工CFD、高效传质与分离和精细化工方面的研究。
计算流体力学软件
计算流体力学（CFD）是近代流体力学，数值数学和计算机科学结合的产物，是一门具有强大生命力的边缘科学。它以电子计算机为工具，应用各种离散化的数学方法，对流体力学的各类问题进行数值实验、计算机模拟和分析研究，以解决各种实际问题。 计算流体力学和相关的计算传热学，计算燃烧学的原理是用数值方法求解非线性联立的质量、能量、组分、动量和自定义的标量的微分方程组，求解结果能预报流动、传热、传质、燃烧等过程的细节，并成为过程装置优化和放大定量设计的有力工具。计算流体力学的基本特征是数值模拟和计算机实验，它从基本物理定理出发，在很大程度上替代了耗资巨大的流体动力学实验设备，在科学研究和工程技术中产生巨大的影响。目前比较好的CFD软件有：Fluent、CFX,Phoenics、Star-CD，除了Fluent 是美国公司的软件外，其它三个都是英国公司的产品 ------------------------------------------------------ FLUENT FLUENT是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%。举凡跟流体，热传递及化学反应等有关的工业均可使用。它具有丰富的物理模型、先进的数值方法以及强大的前后处理功能，在航空航天、汽车设计、石油天然气、涡轮机设计等方面都有着广泛的应用。其在石油天然气工业上的应用包括：燃烧、井下分析、喷射控制、环境分析、油气消散/聚积、多相流、管道流动等等。 Fluent的软件设计基于CFD软件群的思想，从用户需求角度出发，针对各种复杂流动的物理现象，FLUENT软件采用不同的离散格式和数值方法，以期在特定的领域内使计算速度、稳定性和精度等方面达到最佳组合，从而高效率地解决各个领域的复杂流动计算问题。基于上述思想，Fluent开发了适用于各个领域的流动模拟软件，这些软件能够模拟流体流动、传热传质、化学反应和其它复杂的物理现象，软件之间采用了统一的网格生成技术及共同的图形界面，而各软件之间的区别仅在于应用的工业背景不同，因此大大方便了用户。其各软件模块包括： GAMBIT——专用的CFD前置处理器，FLUENT系列产品皆采用FLUENT公司自行研发的Gambit 前处理软件来建立几何形状及生成网格，是一具有超强组合建构模型能力之前处理器，然后由Fluent 进行求解。也可以用ICEM CFD进行前处理，由TecPlot进行后处理。 Fluent5.4——基于非结构化网格的通用CFD求解器，针对非结构性网格模型设计，是用有限元法求解不可压缩流及中度可压缩流流场问题的CFD软件。可应用的范围有紊流、热传、化学反应、混合、旋转流（rotating flow）及震波（shocks）等。在涡轮机及推进系统分析都有相当优秀的结果，并且对模型的快速建立及shocks处的格点调适都有相当好的效果。 Fidap——基于有限元方法的通用CFD求解器，为一专门解决科学及工程上有关流体力学传质及传热等问题的分析软件，是全球第一套使用有限元法于CFD领域的软件，其应用的范围有一般流体的流场、自由表面的问题、紊流、非牛顿流流场、热传、化学反应等等。 FIDAP本身含有完整的前后处理系统及流场数值分析系统。对问题整个研究的程序，数据输入与输出的协调及应用均极有效率。 Polyflow——针对粘弹性流动的专用CFD求解器，用有限元法仿真聚合物加工的CFD软件，主要应用于塑料射出成形机，挤型机和吹瓶机的模具设计。 Mixsim——针对搅拌混合问题的专用CFD软件，是一个专业化的前处理器，可建立搅拌槽及混合槽的几何模型，不需要一般计算流力软件的冗长学习过程。它的图形人机接口和组件数据库，让工程师
计算流体力学实例
汽车外部气体流动模拟 振动和噪声控制研究所 1.模型概述 在汽车外部建立一个较大的长方体几何空间，长度约为30m，宽度和高度约为5m，在空间内部挖出汽车形状的空腔，汽车尺寸参照本田CRV为4550mm*1820mm*1685mm。由于汽车向前开进，气体从车头流向车尾，因此将汽车前方空间设为气体入口，后方空间设为气体出口，模拟气体在车外的流动。另外为了节省计算成本将整个模型按1:100的比例缩小，考虑到模型和流体均是对称的，因此仅画出几何模型的一半区域，建立对称面以考虑生成包含理想气体的流体域。在Catia中建立的模型如图1.1所示。 图1.1几何模型 2.利用ICEM CFD进行网格划分 a)导入有Catia生成的stp格式的模型； b)模型修复，删除多余的点、线、面，允许公差设为0.1； c)生成体，由于本模型仅为流体区域，因此将全部区域划分为一个体，选取方法可以 使用整体模型选取； d)为了后面的设置边界方便，因此将具有相同特性的面设为一个part,共设置了in， out，FreeWalls，Symmetry和Body； e)网格划分，设置Max element=2，共划分了1333817个单元，有225390个节点； f)网格输出，设置求解器为ANSYS CFX，输出cfx5文件。 3.利用ANSYS CFX求解 a)生成域，物质选定Air Ideal Gas，参考压强设为1atm，浮力选项为无浮力模型，
域运动选项为静止，网格变形为无；流体模型设定中的热量传输设定为Isothermal，流体温度设定为288k，湍流模型设定为Shear Stress Transport模型，壁面函数 选择Automatic。 b)入口边界设定，类型为Inlet，位置选定在in，质量与栋梁选定Normal Speed，设 定为15m/s，湍流模型设定类型为Intensity and Length Scale=0.05，Eddy Len.Scale=0.1m。 c)出口边界设定，边界类型为Outlet，位置选out。质量与动量选项为Static Pressure，相对压强为0pa。 d)壁面边界设定，边界类型为Wall，位置选在FreeWalls。壁面边界详细信息中指定 WallInfluence On Flow为Free Slip。 e)对称边界设定，边界类型为Symmetry，位置选在Symmetry。 f)汽车外壁面设定，边界类型为Wall，位置设在Body，壁面详细信息选项中指定Wall Influence On Flow为No Slip，即汽车壁面为无滑移壁面。 g)初始条件设定，初始速度分量设为U方向为15m/s，其他两个方向的速度为零。 h)求解设置，残差类型选为RMS，残差目标设定为1e-5，当求解达到此目标时，求解 自动终止。求解之前的模型如图3.1所示。 图3.1求解之前的模型 4.结果后处理 从图4.1中可以看出计算收敛。
流体力学计算软件报告
三维方管内部二次流特征分析 ——基于NUMECA 数值仿真 2120130457 李明月 【摘 要】运用NUMECA 数值仿真的方法，通过在有粘与无粘的工况下三维方管的内部三维流线对比分析，重点在分析粘性工况下方管内部沿流向各截面上的切向速度矢量分布特征和总压系数分布特征对二次流机理进行讨论和分析。 【关键字】数值仿真 二次流 欧拉方程 N-S 方程 压力梯度 0 前言 在边界层内流体质点向着压力梯度相反并与主流运动方向大致垂直的方向流动，称为二次流。几乎所有的过流通到里面都存在着速度和压力分布不均的情况，压力分布不均则产生一个从高压指向低压的作用力，它与惯性力的大小关系是能否形成二次流的关键。而二次流会使叶轮机械叶片的边界层增厚从而导致分离和损失，而二次流在换热器中增强了对流换热，从而强化了传热，故对二次流的成因和特征的研究具有很大的现实意义。而运用NUMECA 软件对一个简单的三维方管在不同工况下进行数值运算，能够直观地观察得到二次流的结果，并对此进行对比和分析，对流体初学者而言，一方面可以熟悉NUMECA 软件的基本操作，一方面可以基于此加深对二次流的理解。 1 几何描述 如图一所示为三维方管的三维图与所需设定的边界条件。在此算例中，最大的特点在于 中部有一个90°的弯道，且出流部分较长。 10m m 30m m 80m m r20m m r10m m 图1 几何模型
2 网格划分与边界条件 在调入IGG data 文件生成几何文件之后，用网格功能中生成网格块的功能用对应网格顶点与几何顶点重合的方式将网格块贴附在几何模型上，再调整网格数量，和Cluster Points 功能调整边界网格大小，使得近壁面的网格较密，使数值计算时能更好地捕捉到近壁面的参数。生成的网格如图2所示。网格生成后一共33×33×129个网格，网格质量为：最小的正交角度为50.68°，最大宽高比为200，最大膨胀比为1.51，多重网格数为3。在边界条件上，管壁设为SOL 类型，另外短管端面设为INL 类型，剩下那一面设为OUT 类型。 3 边界设定及收敛特性 在NUMECA Fine Turbo 里面建立两个工况并命名为一个无粘一个有粘。在无粘的工况下，选择的流动模型为基于Euler 方程的数学模型。在有粘工况下，流动模型选择的是湍流N-S 方程，并且湍流模型为Spalart-Allmaras 模型。两个工况皆为理想气体的定常流动，进口边界设为总量下（total quantities imposed ）马赫数推断（mach number extrapolated ），进口压力为1.3bar ，进口温度为340K 。出口设定为由静压推断（static pressure imposed ），出口压力为1.0bar 。固壁面在欧拉方程下为无粘的欧拉壁，在N-S 方程里为绝热壁。经初始化后选择计算后输出的参数，除了常规的静压静温和速度外，在壁面数据（solid data ）里额外输出一个粘性压力（viscous stress ）。选择500次迭代后，两种工况下的收敛曲线如图3~图6所示。 图2 三维方管网格划分示意图 图3 Euler 方程下残差收敛曲线
流体力学的计算公式
流体力学的计算公式 众所周知，锅炉之类的热力设备可以采用水封管作为限压装置。设备工作时，水封管内的水将设备与大气隔绝，使设备保持一定的压力。水封管内水的多少，即水封管内的实际水位的高低与设备汽压大小有关。汽压大，水位就高；汽压小，水位就低。当设备压力升至额定蒸汽压力时，水封管内的水位应该达到最高水位；而在设备压力稍微大于额定蒸汽压力时，水封管内的水应能立即被冲除掉，使设备内的蒸汽能够迅速得到排放。但是，水封管的内径必须足够大；以保证它的排汽能力大于设备的最大蒸发量，从而防止设备发生超压事故。《蒸汽锅炉安全技术监察规程》第132条规定：“对于额定蒸汽压力小于等于0.1MPa的锅炉可以采用静重式安全阀或水封式安全装置。水封装置的水封管内径不应小于25mm”。《小型和常压热水锅炉安全监察规定》第二十八条要求：“水封管的内径应当根据锅炉的额定容量和压力确定，且内径不得小于25毫米”。 但是，水封管的内径究竟应该根据什么样的公式去计算？这里没有具体加以规定。水封管的内径为何不能小于25mm？这里没有加以说明。 1 确定水封管内径的一种错误观点 有人说，可以按安全阀的排放量计算公式去间接确定水封管的内径。笔者认为这是不对的。理由如下： 其一，按流体力学，安全阀的局部阻碍因为其进出口呈直角型式，边壁是突变的，主流与边壁之间形成大尺度旋涡，蒸汽排放时能量损失很大；而水封管的边壁是渐变的，又不出现减速增压现象的部位，故蒸汽排放时的能量损失很小。其二，安全阀的出口之外肯定还有排汽管，它的排放量计算公式应该考虑这一额外的因素；而水封管却只相当于它的排汽管。其三，安全阀在排汽时，汽流还自始至终受到弹簧或者重锤施加的反方向作用力，能量损失非常之大；而水封管在排汽时却不会受到这样的反作用力，无此项能量损失。显然，安全阀排汽与水封管排汽区别太大，不能按安全阀排放量计算公式去确定水封管的内径。 2 用流体力学理论建立水封管内径计算方法 为了正确确定水封管的应有内径，我们有必要按流体力学的理论来分析一下水封管的流动阻力和能量损失。蒸汽流过水封管时的阻力损失ΔP等于水封管进口压力P1与出口压力P2之差。管道的阻力损失 ΔP＝Pa（1） 式中λ—管道摩擦系数； d —管道内径，m； ρ—流体的密度，k g／m3； u —流体在管内的流速，m／s； L—平直管段的总长度，m； Σξ—管道中各种局部阻力之和。 管道的摩擦系数λ值取决于流体流动的雷诺数Re 和管壁的相对粗糙度。管壁的相对粗糙度等于管壁的绝对粗糙度（即管壁内凸起高度）K与管内径d之比。雷诺数Re =(2) 式中ν－流体运动粘滞系数，m2／s； 其他符号说明同上。 由于水封管进出口压差通常大于40000Pa，排汽时的流速都大于100m／s ，雷诺数Re 都大于80000000，所以蒸汽的流动状态为紊流流动，而且位于阻力平方区。在此区域内，摩擦阻力系数λ值仅取决于管壁的相对粗糙度，可以用尼古拉兹公式进行计算：λ＝（3）
流体力学计算题..
水 水银 题1图 1 2 3 题型一：曲面上静水总压力的计算问题（注：千万注意方向，绘出压力体） 1、AB 曲面为一圆柱形的四分之一，半径R=0.2m ，宽度（垂直纸面）B=0.8m ，水深H=1.2m ，液体密度3 /850m kg =ρ，AB 曲面左侧受到液体压力。求作用在AB 曲面上的水平分力和铅直分力。（10分） 解：（1）水平分力： RB R H g A h P z c x ?-==)2 (ργ…….(3分) N 1.14668.02.0)2 2 .02.1(8.9850=??- ??=，方向向右（2分） 。 （2）铅直分力：绘如图所示的压力体，则 B R R R H g V P z ??? ? ????+-==4)(2πργ……….（3分） 1.1542 8.042.014.32.0)2.02.1(8.98502=???? ? ?????+?-??=，方向向下（2分）。 2．有一圆滚门，长度l=10m ，直径D=4.2m ，上游水深H1=4.2m ，下游水深H2=2.1m ，求作用于圆滚门上的水平和铅直分压力。
解题思路：（1）水平分力： l H H p p p x )(2 1 222121-= -=γ 方向水平向右。 （2）作压力体，如图，则 l D Al V p z 4 432 πγγγ? === 方向垂直向上。 3.如图示，一半球形闸门，已知球门的半径m R 1= ，上下游水位差m H 1= ，试求闸门受到的水平分力和竖直分力的 大小和方向。 解： (1)水平分力： ()2R R H A h P c πγγ?+===左,2R R A h P c πγγ?=' =右 右左P P P x -= kN R H 79.30114.31807.92=???=?=πγ， 方向水平向右。 （２）垂直分力： V P z γ=，由于左、右两侧液体对曲面所形成的压力体均为半球面，且两侧方向相反，因而垂直方向总的压力为0。 4、密闭盛水容器，已知h 1=60cm,h 2=100cm ，水银测压计读值cm h 25=?。试求半径R=0.5m 的半球盖AB 所受总压力的水平分力和铅垂分力。
中科大计算流体力学CFD之大作业一
CFD 实验报告一 姓名： 学号： 一、题目： 利用中心差分格式近似导数22/dx y d ，数值求解常微分方程 x dx y d 2sin 2 2= （10≤≤x ） 00==x y 4 2 s i n 11- ==x y 步长分别取x ?＝0.05, 0.01, 0.001，0.0001。 二、报告要求： 1）列出全部计算公式和步骤； 2）表列出程序中各主要符号和数组意义； 3）绘出数值计算结果的函数曲线，并与精确解比较； 4）比较不同差分格式和不同网格步长计算结果的精度和代价； 5）附源程序。 三、相关差分格式 二阶导数22/dx y d 的三点差分格式有向前差分、向后差分和中心差分，表达 式分别如下： ()()()22122 22122 211 222 222j j j j j j j j j u u u u O x x x u u u u O x x x u u u u O x x x ++--+--+?=+???-+?=+???-+?=+???一阶向前差分：一阶向后差分：二阶中心差分： 代入微分方程可以得到差分方程，表达式分别如下： 212 212 11 2 2=sin 22=sin 22=sin 2j j j j j j j j j j j j u u u x x u u u x x u u u x x ++--+--+?-+?-+?一阶向前差分：一阶向后差分：二阶中心差分： 对于三种差分格式，差分格式可以改写成AY b =的形式，其中A 是相同的，
非齐次项b 不同，如下所示： 2112112112A -????-?? ??=?? -?? ??-?? 系数矩阵 ()()()02112 3221sin 2sin 2sin 2k k y x x b x x x x y ---???? ??? ? ?=????? ?? ?-?? 一阶向前差分 ()()()()2202 322121sin 2sin 2sin 2sin 2k k x x y x x b x x x x y -???-????? ??=????? ???-?? 一阶向后差分 ()()()()2102 232 2211sin 2sin 2sin 2sin 2k k x x y x x b x x x x y --?? ?-????? ??=??????? ?-? ? 二阶中心差分 求解AY b =可以得到各节点y 的值[]T 1 22 1k k Y y y y y --= 。 四、计算公式和步骤； 1.关于精确解的推导： 已知22sin 2d y x dx =，对 x 进行两次积分，得到121 sin 24 y x C x C =-++，再结合 边界条件00 ==x y 和4 2 sin 11-==x y 得到相对应的1C 和2C ，确定最后精确解为： 1 sin 24y x x =-+。 2.关于数值求解方法： 对于方程组AY b =可直接求解，也可以使用追赶法求解，下面介绍简单追赶法求解三对角方程组的过程。
流体力学计算器 源代码
以下为 program.cs的代码 using System; using System.Collections.Generic; using System.Linq; //using System.Threading.Tasks; using System.Windows.Forms; namespace水力计算优化版 { staticclass Program { ///
计算流体力学教学大纲
《计算流体力学》教学大纲 课程编号：00002067 课程中文名称：计算流体力学 课程英文名称：Computational Fluid Mechanics 总学时：48实验学时：0 上机学时：0 学分：3 适用专业：工程力学 一、课程性质、目的和任务（300字内） 计算流体力学是工程力学专业的一门选修课。它的主要任务是通过教学环节，使学生理解和掌握计算流体力学的一些基本理论和基本计算方法。学会运用计算流体力学、计算机解决一些流体数值计算问题。为从事工程技术工作、科学研究及开拓新领域，打下坚实的基础。 通过本课程的学习，使学生掌握计算流体力学有关的基本概念与基本理论，学会将数学、计算机、流体力学知识有机地结合，对科学研究和工程应用中的流动问题进行数值模拟的方法。教学中应贯彻启发式，互动式。着重讲清基本概念、基本理论、分析问题的思路和方法，并配以适当讨论，逐步培养学生具有分析问题和解决问题的能力。指导学生阅读参考书、文献和资料，培养学生自学获取知识的能力。重视实践环节，要求学生上机计算，并分析计算结果。充分利用现代化的教学手段，使学生具有灵活运用知识、进行探索和开拓创新等方面的技能，并具有较高的综合素质。
二、课程教学内容及学时分配 1、理论讲授 3、教学环节的安排及学时分配 教材：1. 吴子牛主编的《计算流体力学基本原理》，科学出版社，2001年第一版。 2. 王福军主编《计算流体动力学》，清华大学出版社，2005年第一版。
参考资料： 1.R. Peyret & T. D. Taylor. Computational Methods for Floud Flow. 2.C.A.J.Fletcher. Computational Techniques for Fluid Dynamics. 3. H K Versteeg的《计算流体动力学导论》，世界图书出版公司，1995年第一版