流体的基本计算 The final edition was revised on December 14th, 2020.
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，管道内压力P，能否求管道中流体的流速和流量？怎么求 已知管道直径D，管道内压力P，还不能求管道中流体的流速和流量。你设想管道末端有一阀门，并关闭的管内有压力P，可管内流量为零。管内流量不是由管内压力决定，而是由管内沿途压力下降坡度决定的。所以一定要说明管道的长度和管道两端的压力差是多少才能求管道的流速和流量。 对于有压管流，计算步骤如下： 1、计算管道的比阻S，如果是旧铸铁管或旧钢管，可用舍维列夫公式计算管道比阻s=0.001736/d^5.3 或用s=10.3n2/d^5.33计算，或查有关表格； 2、确定管道两端的作用水头差H=P/(ρg)，)，H 以m为单位；P为管道两端的压强差(不是某一断面的压强），P以Pa为单位； 3、计算流量Q：Q = (H/sL)^(1/2) 4、流速V=4Q/(3.1416d^2) 式中：Q―― 流量，以m^3/s为单位；H――管道起端与末端的水头差，以m^为单位；L――管道起端至末端的长度，以m为单位。 管道中流量与压力的关系 管道中流速、流量与压力的关系 流速:V=C√(RJ)=C√[PR/(ρgL)] 流量：Q=CA√(RJ)=√[P/(ρgSL)] 式中：C――管道的谢才系数；L――管道长度；P――管道两端的压力差；R――管道的水力半径；ρ――液体密度；g――重力加速度；S――管道的摩阻。 管道的内径和压力流量的关系 似呼题目表达的意思是：压力损失与管道内径、流量之间的关系，如果是这个问题，则正确的答案应该是：压力损失与流量的平方成正比，与内径5.33方成反比，即流量越大压力损失越大，管径越大压力损失越小，其定量关系可用下式表示： 压力损失（水头损失）公式（阻力平方区） h=10.3*n^2 * L* Q^2/d^5.33 上式严格说是水头损失公式，水头损失乘以流体重度后才是压力损失。式中n――管内壁粗糙度；L――管长；Q――流量；d――管内径 在已知水管:管道压力0.3Mp、管道长度330、管道口径200、怎么算出流速与每小时流量？ 管道压力0.3Mp、如把阀门关了，水流速与流量均为零。（应提允许压力降） 管道长度330、管道口径200、缺小单位，管道长度330米？管道内径200为毫米？其中有无阀门与弯头，包括其形状与形式。 水管道是钢是铸铁等其他材料，其内壁光滑程度不一样。 所以无法计算。 如果是工程上大概数，则工程中水平均流速大约在0.5－－1米/秒左右，则每小时的流量为：0.2×0.2×0.785×1（米/秒，设定值）×3600＝113（立方/小时） 管道每米的压力降可按下式计算：
本节概要 本节讨论喷管内流量、流速的计算。工程上通常依据已知工质初态参数和背压，即喷管出口截面处的工作压力，并在给定的流量等条件下进行喷管设计计算，以选择喷管的外形及确定其几何尺寸；有时也需就已有的喷管进行校核计算，此时喷管的外形和尺寸已定，须计算在不同条件下喷管的出口流速及流量。在喷管的计算中要注意到背压对确定喷管出口截面上压力的作用。 本节内容 4.8.1 流速计算及其分析 4.8.2 临界压力比 4.8.3 流量计算及分析 4.8.4 例题 本节习题 4-24、4-25、4-26、4-27、4-29 下一节 流速计算及其分析 1.喷管出口截面的流速计算 2.压力比对流速的影响 …喷管出口截面的流速计算 据能量方程,气体在喷管中绝热流动时任一截面上的流速可由下式计算： (4-28) 因此，出口截面上流速： (4-28a) 或(4-28b)
在入口速度较小时，上式中可忽略不计，于是： (4-28c) (4-28)各式表明，气流的出口流速取决于气流在喷管中的绝热焓降。值得注意的是，上述各式中焓的单位是J/kg。 如果理想气体可逆绝热流经喷管，可据初态参数(p1,T1)及速度求取滞止参数， 然后结合出口截面参数如p2按可逆绝热过程方程式求出T2从而计算h2再求得；对水蒸汽 可逆绝热流经喷管，可以利用h-s图，根据进口蒸汽的状态查得初态点1，通过点1作垂线与喷管出口截面上压力p2相交，得出状态点2，从点1和2可查出h1和h2，代入式(4-28)即可求出出口流速。 ☆ 式子对理想气体和实际气体均适用；与过程是否可逆无关，但不可逆绝 热流动，若用可逆的关系求出h2在求得的需修正，若h2是不可逆过程终态的焓，则求出的不需修正。 式的适用范围是什么？是否与过程的可逆与否有关？与工质的性质有关？ 返回
计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 一 前 言 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区
课后习题 第一章 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基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
第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软件，是一个专业化的前处理器，可建立搅拌槽及混合槽的几何模型，不需要一般计算流力软件的冗长学习过程。它的图形人机接口和组件数据库，让工程师
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分程的时候提出的一个问题，用一维无粘可压缩Euler程就可以描述了。 这里 这个程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量）随空间变化（）的关系。 在CFD常把这个程写成矢量形式 这里 进一步可以写成散度形式
三维方管内部二次流特征分析 ——基于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 方程下残差收敛曲线
压力与流速的计算公式 没有“压力与流速的计算公式”。流体力学里倒是有一些类似的计算公式，那是附加了很多苛刻的条件的，而且适用的范围也很小。 1，压力与流速并不成比例关系，随着压力差、管径、断面形状、有无拐弯、管壁的粗糙度、是否等径/流体的粘度属性……，无法确定压力与流速的关系。 2，如果你要确保流速，建议你安装流量计和调节阀。也可以考虑定容输送。 要使流体流动，必须要有压力差（注意：不是压力！），但并不是压力差越大流速就一定越大。当你把调节阀关小后，你会发现阀前后的压力差更大，但流量却更小。 管道的水力计算包括长管水力计算和短管水力计算。区别是后者在计算时忽略了局部水头损失，只考虑沿程水头损失。（水头损失可以理解为固体相对运动的摩擦力） 以常用的长管自由出流为例，则计算公式为 H=(v^2*L)/(C^2*R), 其中H为水头，可以由压力换算， L是管的长度， v是管道出流的流速， R是水力半径R=管道断面面积/内壁周长=r/2， C是谢才系数C=R^(1/6)/n， n是糙率，其大小视管壁光洁程度，光滑管至污秽管在0.011至0.014之间取 列举五种判别明渠水流三种流态的方法 [ 标签：明渠,水流,方法 ] （1）明渠水流的分类 明渠恒定均匀流 明渠恒定非均匀流 明渠非恒定非均匀流 明渠非恒定均匀流在自然界是不可能出现的。 明渠非均匀流根据其流线不平行和弯曲的程度，又可以分为渐变流和急变流。 （2）明渠梯形断面水力要素的计算公式： 水面宽度 B = b+2 mh （5—1） 过水断面面积 A =（b+ mh）h （5—2） 湿周（5—3） 水力半径（5—4）
式中：b为梯形断面底宽，m为梯形断面边坡系数，h为梯形断面水深。 （3）当渠道的断面形状和尺寸沿流程不变的长直渠道我们称为棱柱体渠道。 （4）掌握明渠底坡的定义，明渠有三种底坡：正坡（i＞0）平坡（i=0）和逆坡（i＜0。 明渠均匀流特性和计算公式 （1）明渠均匀流的特征： a）均匀流过水断面的形状、尺寸沿流程不变，特别是水深h沿程不变，这个水深也称为正常水深。 b）过水断面上的流速分布和断面平均流速沿流程不变。 c）总水头线坡度、水面坡度、渠底坡度三者相等，J = Js = I。 即水流的总水头线、水面线和渠底线三条线平行。 从力学意义上来说：均匀流在水流方向上的重力分量必须与渠道边界的摩擦阻力相等才能形成均匀流。因此只有在正坡渠道上才可能形成均匀流。 （2）明渠均匀流公式 明渠均匀流计算公式是由连续性方程和舍齐公式组成的，即 Q = A v （5—5）（5—5） 也可表示为：（5—7） 曼宁公式为（5—8） 式中K是流量模数，它表示当底坡为i = 1的时候，渠道中通过均匀流的流量。 水在管道内的流速与水所受的压力有关系吗？ [ 标签：管道流速,流速,关系 ] 水在一根管道内的流速与他所受的压力有什么关系？加上管道对水的阻力之后呢？ 管道的水力计算包括长管水力计算和短管水力计算。区别是后者在计算时忽略了局部水头损失，只考虑沿程水头损失。（水头损失可以理解为固体相对运动的摩擦力） 以常用的长管自由出流为例，则计算公式为 H=(v^2*L)/(C^2*R), 其中H为水头，可以由压力换算，
以下为 program.cs的代码 using System; using System.Collections.Generic; using System.Linq; //using System.Threading.Tasks; using System.Windows.Forms; namespace水力计算优化版 { staticclass Program { ///
流体力学基本概念和基础知识（部分） 1．什么是粘滞性？什么是牛顿内摩擦定律？不满足牛顿内摩擦定律的流体是牛顿流体还是非牛顿流体？ 流体内部质点间或流层间因相对运动而产生内摩擦力以反抗相对运动的性质 dy du A T μ= 满足牛顿内摩擦定律的流体是牛顿流体 请阐述液体、气体的动力粘滞系数随着温度、压强的变化规律。 水的黏滞性随温度升高而减小；空气的黏滞性随温度的升高而增大。（动力粘度μ体现黏滞性）通常的压强对流体的黏滞性影响不大，但在高压作用下，气液的动力黏度随压强的升高而增大。 2．在流体力学当中，三个主要的力学模型是指哪三个？并对其进行说明。 连续介质（对流体物质结构的简化）、无黏性流体（对流体物理性质的简化）、不可压流体（对流体物理性质的简化） 3．什么是理想流体？ 不考虑黏性作用的流体，称为无黏性流体（或理想流体） 4．什么是实际流体？ 考虑黏性流体作用的实际流体 5．什么是不可压缩流体？ 流体在流动过程中，其密度变化可以忽略的流动，称为不可压缩流动。 6．为什么流体静压强的方向必垂直作用面的内法线？ 流体在静止时不能承受拉力和切力，所以流体静压强的方向必然是沿着作用面的内法线方向 7．为什么水平面必是等压面？
由于深度相等的点，压强也相同，这些深度相同的点所组成的平面是一个水平面，可见水平面是压强处处相等的面，即水平面必是等压面。 8．什么是等压面？满足等压面的三个条件是什么？ 在同一种液体中,如果各处的压强均相等由各压强相等的点组成的面称为等压面。满足等压面的三个条件是同种液体连续液体静止液体。 9．什么是阿基米德原理？ 无论是潜体或浮体的压力体均为物体浸入液体的体积，也就是物体排开液体的体积。 10．潜体或浮体在重力G和浮力P的作用，会出现哪三种情况？ 重力大于浮力，物体下沉至底。重力等于浮力，物体在任一水深维持平衡。重力小于浮力，物体浮出液体表面，直至液体下部分所排开的液体重量等于物体重量为止。 11．等角速旋转运动液体的特征有那些？ （1）等压面是绕铅直轴旋转的抛物面簇；（2）在同一水平面上的轴心压强最低，边缘压强最高。 12．什么是绝对压强和相对压强？两者之间有何关系？通常提到的压强是指绝对压强还是相对压强？1个标准大气压值以帕（Pa）、米水柱（mH2O）、毫米水银柱（mmHg）表示，其值各为多少？ 绝对压强：以毫无一点气体存在的绝对真空为零点起算的压强。相对压强：当地同高程的大气压强ap为零点起算的压强。压力表的度数是相对压强，通常说的也是相对压强。1atm=101325pa=10.33mH2O=760mmHg. 13．什么叫自由表面？和大气相通的表面叫自由表面。 14．什么是流线？什么是迹线？流线与迹线的区别是什么？ 流线是某一瞬时在流场中画出的一条空间曲线，此瞬时在曲线上任一点的切线方向与该点的速度方向重合，这条曲线叫流线。区别：迹线是流场中流体质点在一段时间过程中所走过的轨迹线。流线是由无究多个质点组成的，它是表示这无究多个流
计算流体力学入门第九章库特流代码 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)? 计算流体力学(Computational Fluid Dynamics)是流体力学的一个新兴的分支，是一个采用数值方法利用计算机来求解流体流动的控制偏微分方程组，并通过得到的流场和其它物理场来研究流体流动现象以及相关的物理或化学过程的学科。事实上，研究流动现象就是研究流动参数如速度、压力、温度等的空间分布和时间变化，而流动现象是由一些基本的守恒方程（质量、动量、能量等）控制的，因此，通过求解这些流动控制方程，我们就可以得到流动参数在流场中的分布以及随时间的变化，这听起来似乎十分简单。但遗憾的是，常见的流动控制方程如纳维一斯托克斯(Navier-Stokes)方程或欧拉(Euler)方程都是复杂的非线性的偏微分方程组，以解析方法求解在大多数情况下是不可能的。实际上，对于绝大多数有实际意义的流动，其控制方程的求解通常都只能采用数值方法的求解。因此，采用CFD方法在计算机上模拟流体流动现象本质上是流动控制方程（多数情况下是纳维一斯托克斯方程或欧拉方程）的数值求解，而CFD软件本质上就是一些求解流动控制方程的计算机程序。 二、计算流体力学的控制方程 计算流体力学的控剖方程就是流体流动的质量、动量和能量守恒方程。守恒方程的常见的推导方法是基于流体微元的质量、动量和能量衡算。通过质量衡算可以得到连续性方程，通过动量守恒可以得到动量方程，通过能量衡算可以得到能量方程。式(1)一(3)是未经任何简化的流动守恒微分方程，即纳维一斯托克斯方程( N-S方程)。 N-S方程可以表示成许多不同形式，上面的N-S方程是所谓的守恒形式，
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） 流函数
我所知道的计算流体力学(CFD)大牛们 (1) Jameson的故事 Jameson是当今CFD届的超级大牛。偶的超级偶像哦。Jameson是个英国人，出生在军人世家。从小随老爹驻守印度。于是长大了也抗起枪到海外保卫日不落帝国，军衔是Second Lieutenant。无奈“日不落”已落，皇家陆军已经不需要他了。大概有什么立功表现把，退役后就直接进了剑桥大学。在那里拿到博士学位。辗转间从英国来到了美国，从工厂又到了学校。成了Princeton的教授。在那里提出了著名的中心差分格式和有限体积法。就是在这里，发表了他那篇著名的中心差分离散的有限体积法。中心差分格式，大家都知道，是二阶，但是稳定范围特别小，Pe不能超过2，于是就得加人工粘性（一听这名字，数学家就倔嘴巴，不科学嘛），这是大学生都知道的事，怎么加就是学问了。Jameson用二阶项做背景粘性，用四阶项抑制激波振荡（也亏他想得出来），配合他提出的有限体积法，获得了极大的成功，很快风靡世界，工程界几乎无一例外在使用他的方法，原因很简单，他的方法乐百氏，而且又有相当精度。从此大行于市，座上了P大的航空系系主任，也确立了CFD界第一大牛人的地位。 Jameson发文章有个特点，喜欢发在小会议上或者烂杂志上，反正是SCI检索不到地方。包括后来关于非结构网格，多重网格等等经典的开创性文章，都是这样。（如果按照清华的唯SCI论的评判标准，我估计在清华最多只能给他评一个副教授当当。）牛牛的人总是遭人忌妒，哪里都这样。看着Jameson的有限体积方法这么受欢迎，有些人就红眼了。于是说，有限体积方法不错，可惜只适合于定常问题计算，非定常计算就不怎么样嘛。Jameson那里能容忍别人对他的得意之做胡说。于是，灵机一动，想出了一个双时间尺度的方法，引进一个非物理时间，把非定常问题变成了一个定常问题计算，还真好使，又风靡世界，从此天下太平。97年，Jameson年龄到了，就从P大退休了，结果又被聘请到Standford大学当Thomas V. Jones Professor搞起了湍流来。前不久偶导师见他回来，对欧们边摇头边说，“几年不见，老得快不行了”，言下之意，我们如果想多活几年，不要去搞什么湍流。 (2) Steven A. Orszag的故事 Steven A. Orszag是一个天才级别的人物啦。在直接数值模拟，谱方法，湍流模型等等许多方面都有开创性的贡献。天才嘛，总是有缺陷的，不是生活不能自理，就是不懂得处理人际关系。前者还好办，只是lp不舒服，后者嘛，让同事和同行不舒服，可麻烦就大了。 不幸的是，Orszag属于后者。对于他的恃才傲物，有人早就恨得牙根痒痒，报复的机会终于来了。 三十年前，湍流模型的先驱们，是通过数值试验，再连懵带猜的确定下了双方程湍流模型的参数。20年前，Orszag突发奇想，能否用RNG（重整化群理论）从理论上推导这些参数呢？RNG理论在相变上取得了很大的成功，发明者也在81年获得了Nobel奖。牛人就是牛人很快居然真从理论上推出了这些参数。这下湍流模型界可炸开了锅，这岂不是要砸掉很多人的饭碗？这不等于说那些老家伙几十年前的工作一钱不值么？这帮大学霸可不是省油的灯。环顾地球之大，Orszag居然找不到一本杂志愿意接受他这篇文章。Orszag这个郁闷呀，这个生