计算流体动力学在化工中的应用
计算流体力学在化工中的应用
摘要：计算流体力学（CFD）用于求解固定几何形状设备内的流体的动量、热量和质量方程以及相关的其它方程，是化学工程师用于分析问题和解决问题强有力的和用途广泛的工具。本文综述了CFD 在化学工程领域的应用进展及发展趋势。
关键词：计算流体力学；流体流动；化学工程；数值模拟
计算流体力学（Computational Fluid Dynamics, CFD）是流体力学的一个分支，用于求解固定几何形状空间内的流体的动量、热量和质量方程以及相关的其他方程，并通过计算机模拟获得某种流体在特定条件下的有关数据[1]。CFD最早运用于汽车制造业、航天事业及核工业，解决空气动力学中的流体力学问题。CFD计算相对于实验研究，具有成本低、速度快、资料完备、可以模拟真实及理想条件等优点，从而使CFD成为研究各种流体现象，设计、操作和研究各种流动系统和流动过程的有利工具。20世纪60年代末，CFD技术已经在流体力学各相关行业得到了广泛的应用，化学工程的模拟计算始于20世纪90年代后期，如今CFD已经成为研究化工领域中流体流动和传质重要工具。CFD可以用于各种化工装置的模拟、分析及预测，如模拟搅拌槽混合设备的设计、放大；可以预测流体流动过程中的传质、传热，如模拟加热器中的传热效果，蒸馏塔中的两相传质流动状态；可以描述化学反应及反应速率，进行反应器模拟，如模拟出燃烧反应器、生化反应器中的反应速率；还可有效模拟分离、过滤及干燥等设备及装置内流体的流动。
一、CFD在化学工程中的基本原理
CFD是通过数值计算方法来求解化工中几何形状空间内的动量、热量、质量方程等流动主控方程，从而发现化工领域中各种流体的流动现象和规律，其主要以化学方程式中的动量守恒定律、能量守恒定律及质量守恒方程为基础。一般情况下，CFD的数值计算方法主要包括有限差分法、有限元法及有限体积法[2]。CFD是多领域交叉的学科，涉及计算机科学、流体力学、偏微分方程的数学理论、计算几何、数值分析等学科。这些学科的交叉融合,相互促进和支持,推动着这些学科的深入发展。CFD的真实模拟，其主要目的是对流体流动进行预测，以获得流体流动的信息，从而有效控制化工领域中的流体流动。随着信息技术的发展，市场上也出现了CFD软件，其具有对流场进行分析、计算、预测的功能，计算流体力学软件操作简单，界面直观形象，有利于化学工程师对流体进行准确的计算。常见的CFD软件有：FLUENT、PHOENICS、CFX、STAR-CD、FIDIP等。其中FLUENT由美国FLUENT公司于1983年推出的，是目前功能最全面、适用性最广、国内使用最广泛的CFD软件之一。
二、CFD在化学工程中的应用
1.在搅拌槽中的应用分析
搅拌槽由于其内部流动的复杂性，搅拌混合目前尚未形成完善的理论体系，对搅拌槽等混合设备的放大设计，经验成分往往多于理论计算。在工业实际中，特别是快速反应体系或高黏度非牛顿物系，工业规模的反应器存在不同程度的非均匀性，随着规模的增大，这种不均匀性更加严重，经验放大设计方法的可靠性受到前所未有的挑战，因此对搅拌槽内部流场有必要进行更深入的研究。
自从Harvey等[3]用计算机对搅拌槽内的流场进行二维模拟以来，近年来利用CFD的方法研究搅拌槽内的流场发展很快，使用这种方法不仅可以节省大量的研究经费，而且还可以获得通过实验手段所不能得到的数据。Sun等[4]利用CFD中的相关湍流模型计算了搅拌槽内的气液两相流动，并且对其进行了三维模拟，通过实验研究表明，计算流体力学的数值模拟能有效的计算搅拌器上部的气体部分，但其模拟数值也存在一定的缺陷，即不能有效模拟搅拌器底部区域。Wang等[5]以欧拉-欧拉方法为基础，采用相关湍流模型对搅拌槽中液-液-固三相体系各相的流场分布进行了CFD数值模拟，结果表明固体颗粒对液液两相分布有很大的影响，液相分布与实验结果吻合较好，固相分布结果与实验数据还存在一些差异，但是随在叶轮转速的增加也趋向一致。
此外，CFD与多普勒激光测速仪（Laser Doppler Velocimetry, LDV）有效结合可以更深入的研究搅拌装置，因为LDV测量数据可以验证CFD计算结果，并且使用LDV 测定特定点的速度也可作为CFD计算的边界条件。
2.在换热器中的应用分析
换热设备在化学工程中被广泛使用，详细、准确地预测壳程的流动、传热特性对设计经济和可靠的换热器以及评价现有管壳式换热器的性能对工业应用十分重要。针对管壳式换热器几何结构复杂，流动和传热的影响因素很多等特点，运用CFD对管壳式换热器的壳侧流场进行计算机模拟，可以对其他方法难以掌握的壳侧瞬态的温度场和速度场有所了解，利于换热器的机理分析和结构优化。
熊智强等[6]利用CFD技术对管壳式换热器弓形折流板附近流场进行了数值模拟，发现在弓形折流板背面，有部分区域的流速较低，一定程度上存在着流动死区，采用在弓形折流板上开孔的方法后，CFD计算结果显示其传热效率提高，壳侧压降减。
3.在反应工程中的应用研究
CFD在化学反应工程领域也得到了广泛应用。Meier等[7]使用CFD设计和优化新型高温太阳能化学反应器。CFD模拟提供了计算速度、温度和压力场,以及粒子运动轨迹,而这些数据在高辐射和高温条件下无法测量得到。Holgren和Anderson[8]通过CFD软件采用有限元法模拟了用整块催化剂进行的催化燃烧、部分氧化和液相加氢等化学反应，并与实验结果进行了比较。
4.在精馏塔中的应用研究
CFD是研究精馏塔内气液两相流动和传质的重要工具，通过CFD模拟可获得塔
内气液两相微观的流动状况。在板式塔板上的气液传质方面，Mudde等[9]应用低雷诺数的模型对鼓泡塔反应器的持液量和速度分布进行了模拟，在塔气相负荷、塔径、塔高和气液系统的参数大范围变化的情况下，模拟结果和现实的数据能够较好的吻合。Vivek等[10]以欧拉-欧拉方法为基础，充分考虑了塔壁对塔内流体的影响，用CFD软件模拟计算了矩形鼓泡塔内气液相的分散性能，以及气泡数量、大小和气相速度之间的关系，取得了很好的效果。
三、结束语
由于化工过程中经常会出现流体，所以CFD在化学工程领域得到了广泛的应用。所有涉及流体流动、热交换、分子输运、燃烧等现象的问题，几乎都可以通过CFD的方法进行分析和模拟。但CFD还不是一种很成熟的技术，在处理复杂的物理现象、湍流和反应等现象，难以找到合适的模型，对计算机配置要求也高，对于许多问题所应用的数学模型也还不够精确。即使是所谓的通用CFD软件，也不是适合于所有流体力学问题，需要使用者根据研究的对象做认真的选择。即使如此，CFD 已经成为化工过程研究中不可缺少的工具，随着现代计算机硬件和软件技术的发展，CFD将会在化工领域得到更广泛的应用。
参考文献：
[1]姚朝晖，周强. 计算流体力学入门[M]. 北京：清华大学出版社，2010.
[2]王福军. 计算流体动力学分析—CFD软件原理与应用[M]. 北京：清华大学出版社，2004.
[3] Harvey P S, Greaves M. Turbulent flow in an agitated vessel. PartⅡ:Numerical solution and model
prediction [J]. Transaction of the Institute of Chemical Engineers, 1982, 60: 201-210.
[4] Sun H Y, Mao Z S, Yu G Z. Experimental and numerical study of gas hold-up in surface aerated
slimed tanks [J]. Chemical Engineering Science, 2006, 61:4098-4110.
[5] Wang F, Mao Z S, Wang Y F, et al. Measurement of phase holdups in liquid-liquid-solid three-
phase stirred tanks and CFD simulation [J]. Chemical Engineering Science, 2006, 61:7535-7550.
[6]熊智强，喻九阳，曾春.折流板开孔改进管壳式换热器性能的CF'D分析[J]. 武汉化工学院学
报，2006, 38(4): 67-70.
[7] Meier A, et al. Modeling of a novel high-temperature solar chemical reactor [J].Chemical
Engineering Science, 1996, 51 (11): 3181-3186.
[8] Holmgern A, Anderson B. Mass transfer in monolith catalysts–CO oxidation experiments and
simulations [J]. Chemical Engineering Science, 1998, 53 (13): 2285-2298.
[9] Mudde R F. 2D and 3D simulations of an internal airlift loop reactor on the basis of a two-fluid
model[J]. Chemical Engineering Science, 2001, 56 (21): 6351-6358.
[10]Vivek V. B.Dynamics of gas–liquid flow in a rectangular bubble column: experiments and
single/multi-group CFD simulations [J]. Chemical Engineering Science, 2002, 57(22): 4715-4736.
计算流体动力学分析-CFD软件原理与应用_王福军--阅读笔记
计算流体动力学(简称CFD)是建立在经典流体动力学与数值计算方法基础之上的一门新型独立学科，通过计算机数值计算和图像显示的方法，在时间和空间上定量描述流场的数值解，从而达到对物理问题研究的目的。它兼有理论性和实践性的双重特点。 第一章节 流体流动现象大量存在于自然界及多种工程领域中，所有这些过程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。本章向读者介绍这些守恒定律的数学表达式，在此基础上提出数值求解这些基本方程的思想，阐述计算流体力学的任务及相关基础知识，最后简要介绍目前常用的计算流体动力学商用软件。 计算流体动力学((Computational Fluid Dynamics简称CFD)是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。CFD的基本思想可以归结为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看做是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒方程)控制卜对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理量，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 1.1.2计算流体动力学的工作步骤 采用CFD的方法对流体流动进行数值模拟，通常包括如下步骤: (1)建立反映工程问题或物理问题本质的数学模型。具体地说就是要建立反映问题各个量之间关系的微分方程及相应的定解条件，这是数值模拟的出发点。没有正确完善的数 学模型，数值模拟就毫无意义。流体的基本控制方程通常包括质量守恒方程、动量守恒方程、能量守恒方程，以及这些方程相应的定解条件。 (2}}寻求高效率、高准确度的计算方法，即建立针对控制方程的数值离散化方法，如有限差分法、有限元法、有限体积法等。这里的计算方法不仅包括微分方程的离散化方法及求解方法，还包括贴体坐标的建立，边界条件的处理等。这些内容，可以说是c}}的核心。 (3})编制程序和进行计算。这部分工作包括计算网格划分、初始条件和边界条件的输入、控制参数的设定等。这是整个工作中花时间最多的部分。由于求解的问题比较复杂，比如Na}ier-Stakes方程就是一个讨，分复杂的非线性方程，数值求解方法在理论上不是绝对完善的，所以需要通过实验加以验证。正是从这个意义上讲.数值模拟又叫数值试验。应该指出，这部分工作不是轻而易举就可以完成的。 4})显示计算结果。计算结果一般通过图表等方式显示，这对检查和判断分析质量和结果有重要参考意义。 以上这些步骤构成了CFD数值模拟的全过程。其中数学模型的建立是理论
流体力学复习要点(计算公式)
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 +++=++ρρ
《计算流体动力学分析》学习报告
《计算流体动力学分析》学习报告 计算流体力学基础： 本章主要讲解流体动力学的核心思想以及流体动力学的控制方程。 1、计算流体动力学（Computational Fluid Dynamic ）基本思想：把原来在时间和空间上的连续的物理量，用一系列离散点上的变量值来代替，通过一定的原则和方式建立变量之间的代数方程式，求解之后获得变量的近似值。 2、CFD 控制方程： 质量守恒方程 0)·=?+??u t ρρ（ 动量守恒方程（Navier-Stokes 方程） Fz z y x z u w div t w F z y x y u v div t v F z y x x u u div t u zz zx zx y zy yy xy x zx yx xx +??+??+??+??-=+??+??+??+??+??-=+??+??+??+??+??-=+??τττρρρτττρρρτττρρρ)()()()()()( 能量守恒方程 T p S gradT c k div T u div t +=+??)()(T ( ρρ） S T 为粘性耗散项。 方程含有u ，v ，w ，p ，T 和ρ六个未知量，所以还需要一个方程组，才能使其封闭，而这个方程组就是联系P 和ρ的状态方程组：P=（ρ，T ）。 组分质量守恒方程（在一个系统中，可能存在质的交换，或者存在化学组分时使用。） ()s s s s S c grad D div c u div t +=+??)()(c (s ρρρ ） 为便于对控制方程进行计算和分析，对CFD 控制方程写成通用格式： ()S z z y y x x z w y v x u t S grad div u div t +??Γ??+??Γ??+??Γ??=??+??+??+??+Γ=+??)()()()()()())()(φφφφρφρφρρφφφρρφ 依次为瞬态项，对流项，扩散项和源项。 3、湍流控制方程 三维的N-S 方程无论对于层流还是湍流都是是使用的，但由于直接求解三维瞬态的控制方程，对计算机的内存和速度要求很高，因此在工程上广为采用的方法是对瞬态的N-S 方程进行实践平均处理，同时补充反应湍流特性的其他方程，例如湍动能方程以及湍流耗散率方程
流体力学计算公式
C3.6.2 达西摩擦因子 为了确定λ与Re 的关系,人们作了大量实验和理论研究，下面介绍有代表性的结果。 1．尼古拉兹实验 尼古拉兹（J.Nikuradse，1932）分析了达西的圆管沿程阻力实验数据后，发现壁面粗糙度对λ的影响很大，决定用人工粗糙度方法实现对粗糙度的控制。他用当地黄砂砂粒经筛选后分类均匀粘贴在管内壁上，相对粗糙度ε/d 从1/30—1/1014分6种，测得λ与Re 的关系，得到尼古拉兹图（图C3.6.1）。 2． 常用计算公式 从尼古拉兹图中看到在不同Re 数和ε/d 值的区域，λ有不同的变化规律。 图C3.6.1
（1）层流区 由泊肃叶定律推导的沿程水头损失（C3.4.10）式可得 代入达西公式（C3.6.3）式，可得层流区λ的解析式 上式表明层流区λ与管壁粗糙度无关，写成常用对数形式为 上式在双对数坐标系中是一条直线，与尼古拉兹图吻合。 （2）过渡区 该区是层流向湍流的转捩区（2000 CFD 流体动力学软件介绍 CFD—计算流体动力学，因历史原因，国一直称之为计算流体力学。其结构为： 提出问题—流动性质（流、外流；层流、湍流；单相流、多相流；可压、不可压等等），流体属性（牛顿流体：液体、单组分气体、多组分气体、化学反应气体；非牛顿流体） 分析问题—建模—N-S方程（连续性假设），Boltzmann方程（稀薄气体流动），各类本构方程与封闭模型。 解决问题—差分格式的构造/选择，程序的具体编写/软件的选用，后处理的完成。 成果说明—形成文字，提交报告，赚取应得的回报。 CFD实现过程： 1.建模——物理空间到计算空间的映射。 主要软件： 二维： AutoCAD： 大家不要小看它，非常有用。一般的网格生成软件建模都是它这个思路，很少有参数化建模的。相比之下AutoCAD的优点在于精度高，草图处理灵活。可以这样说，任何一个网格生成软件自带的建模工具都是非参数化的，而对于非参数化建模来说，AutoCAD应该说是最好的，毕竟它发展了很多很多年！ 三维： CATIA：航空航天界CAD的老大，法国人的东西，NB，实体建模厉害，曲面建模独步武林。本身可以生成有限元网格，前几天又发布了支持ICEM-CFD的插件ICEM-CFD CAA V5。有了它和ICEM-CFD，可以做任何建模与网格划分！ UG：总觉得EDS脑袋进水了，收了I-deas这么久了，也才发布个几百M的UG NX 2.0，还被大家争论来争论去说它如何的不好用！其实，软件本身不错，大公司用得也多，可是就这么打市场，早晚是走下坡路。按CAD建模的功能来说它排不上第一，也不能屈居第二，尤其是加上了I－DEAS更是如虎添翼。现 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 工程流体力学公式总结 第二章流体得主要物理性质 ?流体得可压缩性计算、牛顿内摩擦定律得计算、粘度得三种表示方法。１.密度ρ= m／V 2.重度γ= G /V 3．流体得密度与重度有以下得关系:γ= ρg或ρ= γ/ g ４.密度得倒数称为比体积,以υ表示υ= １/ ρ= V/m 5.流体得相对密度:d = γ流/γ水= ρ流/ρ水 6．热膨胀性 ７.压缩性、体积压缩率κ 8．体积模量 ９.流体层接触面上得内摩擦力 10．单位面积上得内摩擦力(切应力)(牛顿内摩擦定律) 1１.、动力粘度μ： 12.运动粘度ν:ν=μ/ρ 1３．恩氏粘度°E:°E = t 1 /t 2 第三章流体静力学 ?重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体得压强计算、流体静压力得计算（压力体)。 1.常见得质量力: 重力ΔW = Δmg、 直线运动惯性力ΔFI ＝Δm·ａ 离心惯性力ΔＦR ＝Δm·rω２、 ２.质量力为F。:F= ｍ·aｍ= m(ｆｘi+f yj＋ｆzk) am ＝F/m = f xi＋f yj＋fｚk为单位质量力,在数值上就等于加速度 实例:重力场中得流体只受到地球引力得作用，取z轴铅垂向上,xoy为水平面，则单位质量力在x、y、z轴上得分量为 fｘ= ０,ｆy＝0 , fz=-mg/ｍ= －ｇ式中负号表示重力加速度g与坐标轴z方向相反 ３流体静压强不就是矢量,而就是标量,仅就是坐标得连续函数。即:p=p(x,y,z),由此得静压强得全微分为： 4.欧拉平衡微分方程式 单位质量流体得力平衡方程为: 课后习题 第一章 1.计算流体动力学的基本任务是什么 计算流体动力学是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。 2.什么叫控制方程？常用的控制方程有哪几个？各用在什么场合？ 流体流动要受物理守恒定律的支配，基本的守恒定律包括：质量守恒定律、动量守恒定律、能量守恒定律。如果流动包含有不同组分的混合或相互作用，系统还要遵守组分守恒定律。如果流动处于湍流状态，系统还要遵守附加的湍流输运方程。控制方程是这些守恒定律的数学描述。 常用的控制方程有质量守恒方程、动量守恒方程、能量守恒方程、组分质量守恒方程。质量守恒方程和动量守恒方程任何流动问题都必须满足，能量守恒定律是包含有热交换的流动系统必须满足的基本定律。组分质量守恒方程，在一个特定的系统中，可能存在质的交换，或者存在多种化学组分，每种组分都需要遵守组分质量守恒定律。 4.研究控制方程通用形式的意义何在？请分析控制方程通用形式中各项的意义。 建立控制方程通用形式是为了便于对各控制方程进行分析，并用同一程序对各控制方程进行求解。 各项依次为瞬态项、对流项、扩散项、源项。 6.CFD商用软件与用户自行设计的CFD程序相比，各有何优势？常用的商用CFD软件有哪些？特点如何？ 由于CFD的复杂性及计算机软硬件条件的多样性，用户各自的应用程序往往缺乏通用性。 CFD商用软件的特点是 功能比较全面、适用性强。 具有比较易用的前后处理系统和其他CAD及CFD软件的接口能力，便于用户快速完成造型、网格划分等工作。 具有比较完备的容错机制和操作界面，稳定性高。 可在多种计算机、多种操作系统，包括并行环境下运行。 常用的商用CFD软件有PHOENICS、CFX、SRAR-CD、FIDAP、FLUENT。PHOENICS除了通用CFD软件应该拥有的功能外，PHOENICS软件有自己独特的功能：开放性、CAD接口、运动物体功能、多种模型选择、双重算法选择、多模块选择。 CFX除了可以使用有限体积法外，还采用基于有限元的有限体积法。用于模拟流体流动、传热、多相流、化学反应、燃烧问题。其优势在于处理流动物理现象简单而几何形状复杂的问题。 SRAR-CD基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。 非流体、热动专业CFD新手入门 首先掌握流体力学基本原理，丁祖荣主编的流体力学这本教材，仔细看两天，这样就会知道gambit中为什么会有边界层设置，边界层厚度如何设置；雷诺系数如何确定来判断层流与湍流；马赫数如何确定来判断流体是可压还是不可压，这样就能解决Fluent，是基于压力还是基于密度求解。能够对实际中一些看似简单的流体现象有深刻的认识，能够准确判断是定常流还是非定常流。 CFD网格划分 网格划分对于初学者所接触案例，其实非常简单。但实际工程中，大项目，特别涉及到整套工程，如环保，飞机，网格质量与数量都要求非常高，往往服务器类的PC才能解决问题，所谓的内存128G，CPU四核主频3.0以上。初学者，简单的管道，一般的机器还是没问题。有机械三维软件基础的，对于gambit建模就非常容易了。往往大项目，复杂的结构gambit 建模显得力不从心，所以对于流体工作者来说，学习三维软件对于建模有莫大的帮助，如Proe。 1.1Gambit介绍 网格的划分使用Gambit软件，首先要启动Gambit，在Dos下输入Gambit 计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 一 前 言 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区 借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。 言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分程的时候提出的一个问题，用一维无粘可压缩Euler程就可以描述了。 这里 这个程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量）随空间变化（）的关系。 在CFD常把这个程写成矢量形式 这里 进一步可以写成散度形式 流体力学的计算公式 众所周知，锅炉之类的热力设备可以采用水封管作为限压装置。设备工作时，水封管内的水将设备与大气隔绝，使设备保持一定的压力。水封管内水的多少，即水封管内的实际水位的高低与设备汽压大小有关。汽压大，水位就高；汽压小，水位就低。当设备压力升至额定蒸汽压力时，水封管内的水位应该达到最高水位；而在设备压力稍微大于额定蒸汽压力时，水封管内的水应能立即被冲除掉，使设备内的蒸汽能够迅速得到排放。但是，水封管的内径必须足够大；以保证它的排汽能力大于设备的最大蒸发量，从而防止设备发生超压事故。《蒸汽锅炉安全技术监察规程》第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 所受总压力的水平分力和铅垂分力。 第一章 绪论 单位质量力： 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 -=-= 压强单位换算：水银柱水柱m m 73610/9800012===m m N at 2/1013251m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D += 1 ) () 2(32121h h h h L e ++= 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 arctan =θ 潜体和浮体的总压力： 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 +++=++ρρ 1738年瑞士数学家：伯努利在名著《流体动力学》中提出了伯努利方程。 1755年欧拉在名著《流体运动的一般原理》中提出理想流体概念，并建立了理想流体基本方程和连续方程，从而提出了流体运动的解析方法，同时提出了速度势的概念。 1781年拉格朗日首先引进了流函数的概念。 1826年法国工程师纳维，1845年英国数学家、物理学家斯托克思提出了著名的N-S方程。 1876年雷诺发现了流体流动的两种流态：层流和紊流。 1858年亥姆霍兹指出了理想流体中旋涡的许多基本性质及旋涡运动理论，并于1887年提出了脱体绕流理论。 19世纪末，相似理论提出，实验和理论分析相结合。 1904年普朗特提出了边界层理论。 20世纪60年代以后，计算流体力学得到了迅速的发展。流体力学内涵不断地得到了充实与提高。 理想势流伯努利方程 （3-14） 或（3-15） 物理意义:在同一恒定不可压缩流体重力势流中，理想流体各点的总比能相等即在整个势流场中，伯努利常数C 均相等。 （应用条件：“”所示） 符号说明 二、沿流线的积分 1.只有重力作用的不可压缩恒定流，有 2.恒定流中流线与迹线重合： 沿流线（或元流）的能量方程： （3-16） 注意：积分常数C，在非粘性、不可压缩恒定流流动中，沿同一流线保持不变。一般不同流线各不相同（有旋流）。 （应用条件：“”所示，可以是有旋流） 流速势函数（势函数）观看录像>> ?存在条件：不可压缩无旋流，即或 必要条件存在全微分d 直角坐标 （3-19） 式中：——无旋运动的流速势函数，简称势函数。 ?势函数的拉普拉斯方程形式 对于不可压缩的平面流体流动中，将（3-19）式代入连续性微分方程（3-18），有： 或（3-20） 适用条件：不可压缩流体的有势流动。 点击这里练习一下 极坐标 （3-21） 流函数 以下为 program.cs的代码 using System; using System.Collections.Generic; using System.Linq; //using System.Threading.Tasks; using System.Windows.Forms; namespace水力计算优化版 { staticclass Program { /// 计算流体力学入门第九章库特流代码 fortan90版 ! -------------------------------------------------- ! Silverfrost FTN95 for Microsoft Visual Studio ! Free Format FTN95 Source File ! -------------------------------------------------- program piple implicitnone real,dimension(21)::u real,dimension(21)::uu real,dimension(21,3)::cf integer::i real::s=0.0 real::err=1 ! judgement of wheather stop or not do i=1,21 u(i)=0 enddo dowhile(err>1e-8) u(1)=0.0 u(21)=1.0 uu(1)=0.0 uu(21)=1.0 cf(:,1)=-0.5 cf(:,2)=2.0 do i=2,20 cf(i,3)=0.5*(u(i+1)+u(i-1)) enddo cf(20,3)=cf(20,3)+0.5 do i=3,20 cf(i,2)=cf(i,2)-(cf(i,1)*cf(i-1,1))/cf(i-1,2) cf(i,3)=cf(i,3)-(cf(i-1,3)*cf(i,1))/cf(i-1,2) enddo uu(20)=cf(20,3)/cf(20,2) do i=19,1,-1 uu(i)=(cf(i,3)+0.5*uu(i+1))/cf(i,2) enddo uu(1)=0 do i=1,21 s=s+abs(uu(i)-u(i)) enddo u=uu err=s s=0.0 print*,err enddo print*,uu read*,i endprogram piple ! -------------------------------------------------- ! Silverfrost FTN95 for Microsoft Visual Studio ! Free Format FTN95 Source File ! -------------------------------------------------- program piple 第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1 计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。 从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解析解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解析解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力 第1章CFD 基础 计算流体动力学(computational fluid dynamics，CFD)是流体力学的一个分支，它通过计算机模拟获得某种流体在特定条件下的有关信息，实现了用计算机代替试验装置完成“计算试验”，为工程技术人员提供了实际工况模拟仿真的操作平台，已广泛应用于航空航天、热能动力、土木水利、汽车工程、铁道、船舶工业、化学工程、流体机械、环境工程等领域。 本章介绍CFD一些重要的基础知识，帮助读者熟悉CFD的基本理论和基本概念，为计算时设置边界条件、对计算结果进行分析与整理提供参考。 1.1 流体力学的基本概念 1.1.1 流体的连续介质模型 流体质点(fluid particle)：几何尺寸同流动空间相比是极小量，又含有大量分子的微元体。 连续介质(continuum/continuous medium)：质点连续地充满所占空间的流体或固体。 连续介质模型(continuum/continuous medium model)：把流体视为没有间隙地充满它所占据的整个空间的一种连续介质，且其所有的物理量都是空间坐标和时间的连续函数的一种假设模型：u =u(t,x,y,z)。 1.1.2 流体的性质 1. 惯性 惯性(fluid inertia)指流体不受外力作用时，保持其原有运动状态的属性。惯性与质量有关，质量越大，惯性就越大。单位体积流体的质量称为密度(density)，以r表示，单位为kg/m3。对于均质流体，设其体积为V，质量为m，则其密度为 m ρ=(1-1) V 对于非均质流体，密度随点而异。若取包含某点在内的体积V?，其中质量m ?，则该 点密度需要用极限方式表示，即 0lim V m V ρ?→?=? (1-2) 2. 压缩性 作用在流体上的压力变化可引起流体的体积变化或密度变化，这一现象称为流体的可压缩性。压缩性(compressibility)可用体积压缩率k 来量度 d /d /d d V V k p p ρρ =-= (1-3) 式中：p 为外部压强。 在研究流体流动过程中，若考虑到流体的压缩性，则称为可压缩流动，相应地称流体为可压缩流体，例如高速流动的气体。若不考虑流体的压缩性，则称为不可压缩流动，相应地称流体为不可压缩流体，如水、油等。 3. 粘性 粘性(viscosity)指在运动的状态下，流体所产生的抵抗剪切变形的性质。粘性大小由粘度来量度。流体的粘度是由流动流体的内聚力和分子的动量交换所引起的。粘度有动力粘度μ和运动粘度ν之分。动力粘度由牛顿内摩擦定律导出： d d u y τμ= (1-4) 式中：τ为切应力，Pa ；μ为动力粘度，Pa ?s ；d /d u y 为流体的剪切变形速率。 运动粘度与动力粘度的关系为 μ νρ = (1-5) 式中：ν为运动粘度，m 2/s 。 在研究流体流动过程中，考虑流体的粘性时，称为粘性流动，相应的流体称为粘性流体；当不考虑流体的粘性时，称为理想流体的流动，相应的流体称为理想流体。 根据流体是否满足牛顿内摩擦定律，将流体分为牛顿流体和非牛顿流体。牛顿流体严格满足牛顿内摩擦定律且μ保持为常数。非牛顿流体的切应力与速度梯度不成正比，一般又分为塑性流体、假塑性流体、胀塑性流体3种。 塑性流体，如牙膏等，它们有一个保持不产生剪切变形的初始应力0τ，只有克服了这个初始应力后，其切应力才与速度梯度成正比，即 0d d u y ττμ=+ (1-6) 假塑性流体，如泥浆等，其切应力与速度梯度的关系是CFD—计算流体动力学软件介绍
流体力学计算公式
流体力学公式总结
计算流体力学课后题作业
CFD计算流体动力学入门教程选择
计算流体力学常用数值方法简介[1]
计算流体力学过渡到编程的傻瓜入门教程
流体力学的计算公式
流体力学计算题..
流体力学复习要点(计算公式)
流体主要计算公式
流体力学计算器 源代码
计算流体力学入门 第九章 库特流代码 fortan90版
第二章计算流体力学的基本知识
CFD基础(流体力学)