计算流体力学课程总结 计算流体动力学(computational Fluid Dynamics,简称CFD)是通过计算机数值 计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。是用电子计算机和离散化的数值方法对流体力学问题进行数值模拟和分析的一个分支。 流体力学和其他学科一样，是通过理论分析和实验研究两种手段发展起来的。很早就已有理论流体力学和实验流体力学两大分支。理论分析是用数学方法求出问题的定量结果。但能用这种方法求出结果的问题毕竟是少数，计算流体力学正是为弥补分析方法的不足而发展起来的。计算流体力学是目前国际上一个强有力的研究领域,是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术，广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计、HAVC&R 等诸多工程领域。 计算流体力学的任务是流体力学的数值模拟。数值模拟是“在计算机上实现的一 个特定的计算，通过数值计算和图像显示履行一个虚拟的物理实验——数值实验“。 数值模拟包括以下几个部分。首先，要建立反映问题(工程问题、物理问题等)本质数 学模型。其次，数学模型建立以后需要解决的问题是寻求高效率、高准确度的计算方法。再次，在确定了计算方法和坐标系统后，编制程序和进行计算式整个工作的主体。最后，当计算工作完成后，流畅的图像显示是不可缺少的部分。 还有一个就是CFD的基本思想问题，它就是把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通 过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求 解代数方程组获得场变量的近似值。 经过四十多年的发展，CFD出现了多种数值解法。这些方法之间的主要区别在于 对控制方程的离散方式。根据离散的原理不同，CFD大体上可分为三个分支： ?有限差分法(Finite Different Method，FDM) ?有限元法(Finite EIement Method，FEM) ?有限体积法(Finite Volume Method，FVM) 有限差分法是应用最早、最经典的CFD方法，也是最成熟、最常用的方法。它将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的 导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求出差分万程组 的解，就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题 的近似数值解法。
大型飞机高级人才培养班 航空工程全日制工程硕士研究生培养方案 一、适用类别或领域 航空工程（085232） 二、培养目标 材料工程、电子与通信工程、控制工程、航空工程领域全日制工程硕士 (以下简称航空工程等领域全日制工程硕士)是与以上各工程领域任职资格相联系的专业学位，主要为国民经济和国防建设等领域培养应用型、复合型高层次工程技术和工程管理人才。大飞机班旨在探索一条“以国家大型项目人才需求为索引，培养具有献身精神、团结协作精神、开拓创新精神的设计型和复合型人才”的研究生培养新模式，是北航研究生培养体系的一部分。 航空工程等领域全日制工程硕士培养的基本要求是： 1、坚持党的基本路线，热爱祖国、遵纪守法、品行端正、诚实守信、身心健康，具有良好的科研道德和敬业精神。 2、在本领域掌握坚实的基础理论和系统的专门知识，有较宽的知识面和较强的自立能力，具有大飞机设计、制造、运营、管理等领域需求的创造能力和工程实践能力。 3、掌握一门外国语。 三、培养模式及学习年限 1．航空工程等领域全日制工程硕士研究生培养实行导师负责制，或以导师为主的指导小组制，负责制订硕士研究生个人培养计划，选课、组织开题报告、论文中期检查、指导科学研究和学位论文，并与中国商飞、第一飞机设计研究院、西飞公司等航空企业联合培养，实行导师组指导。 2．硕士研究生一般用1学年完成课程学习，课程学习实行学分制，具体学习、考核及管理工作执行《北京航空航天大学研究生院关于研究生课程学习管理规定》。 3．专业实习是全日制工程硕士研究生培养中的重要环节，全日制工程硕士研究生在学期间，应保证不少于0.5年的工程实践。 4．学位论文选题应来源于航空工程等领域工程技术背景。鼓励实行双导师制，其中第一导师为校内导师，校外导师应是与本工程领域相关的专家，也可以根据学生的论文
第一章绪论 1.工程流体力学的研究对象：工程流体力学以流体（包括液体和气体）为研究对象，研究流体宏观 的平衡和运动的规律，流体与固体壁面之间的相互作用规律，以及这些规律在工程实际中的应用。 第二章流体的主要物理性质 1.★流体的概念：凡是没有固定的形状，易于流动的物质就叫流体。 2.★流体质点：包含有大量流体分子，并能保持其宏观力学性能的微小单元体。 3.★连续介质的概念：在流体力学中，把流体质点作为最小的研究对象，从而把流体看成是： 1)由无数连续分布、彼此无间隙地； 2)占有整个流体空间的流体质点所组成的介质。 4.密度：单位体积的流体所具有的质量称为密度，以ρ表示。 5.重度：单位体积的流体所受的重力称为重度，以γ表示。 6.比体积：密度的倒数称为比体积，以υ表示。它表示单位质量流体所占有的体积。 7.流体的相对密度：是指流体的重度与标准大气压下4℃纯水的重度的比值，用d表示。 8.★流体的热膨胀性：在一定压强下，流体体积随温度升高而增大的性质称为流体的热膨胀性。 9.★流体的压缩性：在一定温度下，流体体积随压强升高而减少的性质称为流体的压缩性。 10.可压缩流体：ρ随T 和p变化量很大，不可视为常量。 11.不可压缩流体：ρ随T 和p变化量很小，可视为常量。 12.★流体的粘性：流体流动时，在流体内部产生阻碍运动的摩擦力的性质叫流体的粘性。 13.牛顿内摩擦定律：牛顿经实验研究发现，流体运动产生的内摩擦力与沿接触面法线方向的速度变 化（即速度梯度）成正比，与接触面的面积成正比，与流体的物理性质有关，而与接触面上的压强无关。这个关系式称为牛顿内摩擦定律。 14.非牛顿流体：通常把满足牛顿内摩擦定律的流体称为牛顿流体，此时不随dυ/d n而变化，否则称 为非牛顿流体。 15.动力粘度μ：动力粘度表示单位速度梯度下流体内摩擦应力的大小，它直接反映了流体粘性的 大小。 16.运动粘度ν：在流体力学中，动力粘度与流体密度的比值称为运动粘度，以ν表示。 17.实际流体：具有粘性的流体叫实际流体（也叫粘性流体）， 18.理想流体：就是假想的没有粘性（μ= 0）的流体 第三章流体静力学 1.★流体的平衡：（或者说静止）是指流体宏观质点之间没有相对运动，达到了相对的平衡。 2.★绝对静止：流体对地球无相对运动，也称为重力场中的流体平衡。 3.★相对平衡：流体整体对地球有相对运动，但流体对运动容器无相对运动，流体质点之间也无相 对运动，这种静止或叫流体的相对静止★：体积力：作用于流体的每一个流体质点上，其大小与流体所具有的质量成正比的力。在均质流体中，质量力与受作用流体的体积成正比，因此又叫。 4.★表面力：表面力是作用于被研究流体的外表面上，其大小与表面积成正比的力。 5.★压强：在静止或相对静止的流体中，单位面积上的内法向表面力称为压强。 6.等压面：在静止流体中，由压强相等的点所组成的面。 7.★位置水头（位置高度）：流体质点距某一水平基准面的高度。 8.压强水头（压强高度）：由流体静力学基本方程中的p/（ρg）得到的液柱高度。 9.★静力水头：位置水头z和压强水头p/（ρg）之和。 10.压强势能：流体静力学基本方程中的p/ρ项为单位质量流体的压强势能。
研究生《流体力学实验》 ——飞机标模染色液流动显示 实验报告 班级 姓名 实验日期 指导教师 北京航空航天大学流体力学研究所
一、实验目的 1. 掌握染色流动显示技术的基本原理、应用方法和实验过程中应注意的技术问题。 2. 了解战斗机典型绕流现象和特性，包括机翼前缘涡（边条涡）、机头涡的形态、特征、涡 系间相互作用，以及攻角影响等，并分析这些流动现象对飞机气动性能的影响。 二、基本原理 流动显示技术是显示技术包括方法、设备、记录手段、图像处理和数据分析等方面，逐渐形成专门的实验技术。 水洞中常用的流动显示技术有氢气泡方法和染色方法等（属于示踪粒子方法），配以激光片光源等辅助手段可以得到很多有意义的细节结果。染色线流动显示是在在被观测的流场中设置若干个点，在这些点上不断释放某种颜色的液体，它随流过该点的流体微团一起往下游流去，流过该点的所有流体微团组成了可视的染色线。染料选取应注意：1.所选取的染料应使染色线扩散慢、稳定性好；2.染色液应与水流具有尽可能相同的密度（与酒精混合）； 3. 染料颜色与流场背景形成强的反差（荧光染料）注入方式；4.在绕流物体表面开孔；5.直接注入流场中所需要观测的位置。 本实验选用飞机标模，利用染色液方法观察其绕流的典型流动现象，重点关注机头涡、边条涡及其对基本翼（主翼也称后翼）流动的影响。 三、实验装置及模型 1.实验模型 飞机标模由机身、机翼、尾翼构成，见图2。机身为尖拱型头部加圆柱形后体，机翼为大后掠边条加中度后略三角翼主翼，尾翼包括水平尾翼和垂直尾翼（单立尾）。各部分表面都布有染色液出孔。
2.实验风洞 北航1.2米多用途低速串联水平回流式水洞。该水洞实验段尺寸大、流场品质高，与同类设备比较，不但在国内领先，而且达到国际先进水平。设备主实验段1.2米×1米×16米（高×宽×长），流速范围0.1～1.0米/秒。主实验段主要流场品质：湍流度0.27％～0.45％，截面速度不均匀度：0.46％。 四、实验步骤 1.实验准备，将染色液注入系统； 2.开启水洞，水流速度稳定到10cm/s； 3.调整攻角； 4.待流场稳定后，调节染色液流量，得到清晰的流动结构显示形态； 5.待流动稳定后，观察稳定的流态，拍摄照片； 6. 将攻角分别调整到0 o，5o，10o，15o，20o，25o，30o，35o，40o，45o，50o，55o，60o，重复步骤5，直到所要求的攻角状态实验全部完成。 五、实验结果报告 1．实验条件： ①水温t=20o C； 水的运动粘性系数υ=0.878×10-6米2秒； 附：水的运动粘性系数随温度的变化： ②水流速度 U = 0.1 米/秒； ③特征长度C=0.115m （C为模型机翼平均弦长） 计算：雷诺数 Re = UC /υ= 1.310×104； 2．实验结果和分析
2012～2013学年第1学期 12级研究生《计算流体力学》结课作业 适用专业：供热供燃气通风及空调工程 一、结合某一具体学科，阐述纯理论方法、实验方法及数值方法在科学研究中的各自优缺点，在此基础上论述数值模拟方法的发展前景。（不少于4千字）。 流体力学是力学的一个重要分支, 是研究流体(液体和气体)的力学运动规律及其应用的学科, 主要研究在各种力的作用下，流体本身的静止状态和运动状态特征，以及流体和相邻固体界面有相对运动时的相互作用和流动规律。在人们的生活和生产活动中随时随地都可遇到流体，流体力学与人类的日常生活和生产事业密切相关。按其研究内容的侧重点不同，分为理论流体力学和工程流体力学。其中理论流体力学主要采用严密的数学推理方法，力求准确性和严密性，工程流体力学侧重于解决工程实际中出现的问题，而不追求数学上的严密性。当然由于流体力学研究的复杂性，在一定程度上，两种方法都必须借助于实验研究，得出经验或半经验的公式。 在实际工程的诸多领域流体力学都起着十分重要的作用。如气象、水利的研究，船舶、飞行器、叶轮机械和核电站的设计及其运行，可燃气体或炸药的爆炸，都广泛地用到流体力学知识。许多现代科学技术所关心的问题既受流体力学的指导，同时也促进了流体力学自身的不断发展。1950年后，计算机的发展给予流体力学以极大的推动作用。 目前，解决流体力学问题的方法主要有实验方法、理论分析方法和数值方法三种。 实验方法 同物理学、化学等学科一样，流体力学的研究离不开实验，尤其是对新的流体运动现象的研究。实验能显示运动特点及其主要趋势，有助于形成概念，检验理论的正确性。二百年来流体力学发展史中每一项重大进展都离不开实验。流体力学实验研究方法有实物实验、比拟研究和模型研究三类：实物实验是用仪器实测原型系统的流动参数，适用于较小的原型；比拟实验是利用电场和磁场来模拟流场，实施起来限制条件较多；模型研究是实验流体力学最常用的研究方法。 实验研究的一般过程是：在相似理论的指导下建立实验模型，用流体测量技术测量流动参数，处理和分析实验数据。建立实验模型要求模型与原型满足相似理论，即满足两个流场
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30～40 年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学" 。 从20 世纪60 年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
流体力学总结 第一章 流体及其物理性质 1. 流体：流体是一种受任何微小剪切力作用都能连续变形的物质，只要这种力继续作用， 流体就将继续变形，直到外力停止作用为止。流体一般不能承受拉力，在静止状态下也不能承受切向力，在任何微小切向力的作用下，流体就会变形，产生流动 2. 流体特性：易流动(易变形)性、可压缩性、粘性 3. 流体质点：宏观无穷小、微观无穷大的微量流体。 4. 流体连续性假设：流体可视为由无数连续分布的流体质点组成的连续介质。稀薄空气和 激波情况下不适合。 5. 密度0lim V m m V V δδρδ→== 重度0lim V G G g V V δδγρδ→=== 比体积1v ρ= 6. 相对密度：是指某流体的密度与标准大气压下4C 时纯水的密度（1000）之比 w w S ρ ρρ= 为4C 时纯水的密度 13.6Hg S = 7. 混合气体密度1 n i i i ρρα == ∑ 8. 体积压缩系数：温度不变，单位压强增量引起的流体体积变化率。体积压缩系数的倒数 为体积模量1 P P K β= 1p V p V δβδ=- 1 1 0 1.4p p T Q p p βγβγ→= === 9. 温度膨胀系数：压强不变，单位温升引起的流体体积变化率。 1T V T V δβδ= 1 T p T β→= 10. 不可压缩流体：流体受压体积不减少，受热体积不膨胀，密度保持为常数，液体视为不 可压缩流体。气体流速不高，压强变化小视为不可压缩流体
11. 牛顿内摩擦定律： du dy τμ = 黏度du dy τ μ= 流体静止粘性无法表示出来，压强对黏 度影响较小，温度升高，液体黏度降低，气体黏度增加 μ υρ = 。满足牛顿内摩擦定律的流体为牛顿流体。 12. 理想流体：黏度为0，即0μ=。完全气体：热力学中的理想气体 第二章 流体静力学 1. 表面力：流体压强p 为法向表面应力，内摩擦τ是切向表面应力（静止时为0）。 2. 质量力（体积力）：某种力场对流体的作用力，不需要接触。重力、电磁力、电场力、 虚加的惯性力 3. 单位质量力：x y z F f f i f j f k m ==++ ，单位与加速度相同2m s 4. 流体静压强： 1）流体静压强的方向总是和作用面相垂直且指向该作用面，即沿着作用面的内法线方向 2）在静止流体内部任意点处的流体静压强在各个方向都是相等的。 x y z n p p p p === 5. 流体平衡微分方程式（欧拉平衡方程） 101010 x y z p f x p f y p f z ρρρ?- =??-=??-=? 10 p p p f p p i j k x y z ρ???-?=?= ++??? 6. 压差方程 ()x y z dp f dx f dy f dz ρ=++ 7. 势函数 ()()() ,,x y z f f f x y z πππ?-?-?-= ==??? ()dp d ρπ=-
北航考研之科研成果及重点实验室汇总 学校在尖端技术研究领域居于国内高校前列，有40多项科研成果具有开辟意义;该校研制发射成功的多种型号飞行器填补了国内多项空白，如中国第一架轻型旅客机“北京一号”、亚洲第一枚探空火箭“北京二号”、中国第一架无人驾驶飞机“北京五号”、“蜜蜂”系列飞机、共轴式双旋翼无人驾驶直升机等。自2001年至2013年，北航共获国家三大科技奖励49项;特别在2005年至2013年，该校连续7年获得7项国家级科技奖励一等奖。 2013年学校科研经费到款23.23亿元，6项成果获得国家奖励，3位教授及其团体获国家自然科学二等奖，获批5项973项目、12项自然科学基金重点项目，以及237项面上项目，重大工程项目进展顺利并获嘉奖。以“3D打印”为代表的技术创新取得重大进展，“月宫一号”实验装置取得实质进展，物理科学与核能工程学院参与的发现四夸克物质Zc(3900)被评为2013美国《物理》杂志年度研究热点、生物与医学工程学院两篇文章分别名列ESI 近两年热点论文和材料领域近十年高引用论文。 标志性成果 昆虫飞行的空气动力学和飞行力学 实时三维图形平台BH-GRAPH 航空航天、先进制造等复杂工程系统 过渡金属及其化合物纳米材料的可控制备、微结构及特性研究 六方铁磁体的稳定磁结构耦合及其可控磁功能特性 科研机构 截止2013年，北京航空航天大学拥有1个国家实验室、9个国家级重点实验室、4个国家级工程研究中心、42个省部级重点实验室、3个省部级工程中心和3个中关村开放实验室。 国家实验室 航空科学与技术国家实验室 国家重点实验室 虚拟现实新技术国家重点实验室、软件开发环境国家重点实验室、国家计算流体力学实验室、航空发动机气动热力国家科技重点实验室、“863”高技术CIMS设计自动化工程实验室、惯性技术国家级重点实验室、可靠性与环境工程实验室、飞行器控制一体化技术实验室、国家空管新航行系统技术重点实验室。 省部级重点实验室 航空可靠性综合航空科技重点实验室、数字化设计与制造技术北京市重点实验室、网络技术北京市重点实验室、计算机新技术实验室、材料力学实验室、流体力学教育部重点实验室、先进仿真技术航空科技重点实验室、航空电子航空科技重点实验室、特种功能材料与薄膜技术北京市重点实验室、聚合物基复合材料北京市高技术实验室、“复杂系统分析与管理决策”教育部重点实验室、“城市运行应急保障模拟技术”北京市重点实验室等。 研究所 航空探测研究所、A TE技术研究所、可靠性工程研究所、外国语言研究所、设备工程
流体力学 实验指导书与报告 静力学实验 雷诺实验 中国矿业大学能源与动力实验中心
学生实验守则 一、学生进入实验室必须遵守实验室规章制度，遵守课堂纪律，衣着整洁，保持安静，不得迟到早退，严禁喧哗、吸烟、吃零食和随地吐痰。如有违犯，指导教师有权停止基实验。 二、实验课前，要认真阅读教材，作好实验预习，根据不同科目要求写出预习报告，明确实验目的、要求和注意事项。 三、实验课上必须专心听讲，服从指导教师的安排和指导，遵守操作规程，认真操作，正确读数，不得草率敷衍，拼凑数据。 四、预习报告和实验报告必须独自完成，不得互相抄袭。 五、因故缺课的学生，可向指导教师申请一次补做机会，不补做的，该试验以零分计算，作为总成绩的一部分，累计三次者，该课实验以不及格论处，不能参加该门课程的考试。 六、在使用大型精密仪器设备前，必须接受技术培训，经考核合格后方可使用，使用中要严格遵守操作规程，并详细填写使用记录。 七、爱护仪器设备，不准动用与本实验无关的仪器设备。要节约水、电、试剂药品、元器件、材料等。如发生仪器、设备损坏要及时向指导教师报告，属责任事故的，应按有关文件规定赔偿。 八、注意实验安全，遵守安全规定，防止人身和仪器设备事故发生。一旦发生事故，要立即向指导教师报告，采取正确的应急措施，防止事故扩大，保护人身安全和财产安全。重大事故要同时保护好现场，迅速向有关部门报告，事故后尽快写出书面报告交上级有关部门，不得隐瞒事实真相。 九、试验完毕要做好整理工作，将试剂、药品、工具、材料及公用仪器等放回原处。洗刷器皿，清扫试验场地，切断电源、气源、水源，经指导教师检查合格后方可离开。 十、各类实验室可根据自身特点，制定出切实可行的实验守则，报经系(院)主管领导同意后执行，并送实验室管理科备案。 1984年5月制定 2014年4月再修订 中国矿业大学能源与动力实验中心
计算流体力学结课报告200Km/h列车fluent仿真计算 学部：化、环、生学部 学院：化工机械与安全学院 学号：31507095 班级：化1512班 学生姓名：孙金
引言 数值仿真就是对所建立的数值模型进行数值实验和求解的过程。而计算流体力学CFD （Computational Fluid Dynamics）就是在工程仿真实验领域中应用最广泛的一门学科。任何流体运动的规律都是以质量守恒定律、动量守恒定律和能量守恒定律为基础的。这些基本定律可由数学方程组来描述，如欧拉方程、N-S方程。采用数值计算方法，通过计算机求解这些控制流体流动的数学方程，进而研究流体的运动规律这就是CFD研究问题的方法。在实际计算流体力学方面，采用通用的CFD软件来完成工程上的一些流体力学问题，有极为广泛的应用前景。近年来，随着计算机技术以及相关技术的发展，CFD技术已经在工程领域内取得重大的进步，特别是在高速列车的外型设计方面起了很大作用。随着国家经济的发展，国家运输业也有了很大的发展，特别是列车经过几次提速后，高速列车在国家运输行业中所占比例不断提高。高速列车的特点是庞大、细长、在地面轨道上运行，其空气动力学问题非常复杂。空气在列车表面形成空气流场，空气阻力急剧增加，作用在列车的阻力大部分来自压强阻力，而一部分来自表面磨擦阻力，这就使能耗过大，同时列车可能出现较大的空气升力，导致列车产生“飘”的现象，激发列车脱轨事故的发生，因此研究高速列车气动力性能非常重要。用CFD仿真可以详细了解高速列车的空气动力特性，从而设计出阻力小、噪音低等各方面性能完善的高质量列车。本文采用CFD学科中的常用商业软件Fluent仿真一个时速200km/h的二维流线型车头的外流场，对其空气动力性能进行分析，从而得到不同车辆形状其周围流场的不同，进而分析哪种车型更适合。
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较方便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分方程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子内面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分方程的时候提出的一个问题，用一维无粘可压缩Euler方程就可以描述了。 这里 这个方程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量 ）随空间变化（）的关系。 在CFD中通常把这个方程写成矢量形式 这里 进一步可以写成散度形式
流体力学知识点总结 第一章 绪论 1 液体和气体统称为流体，流体的基本特性是具有流动性，只要剪应力存在流动就持续进行，流体在静止时不能承受剪应力。 2 流体连续介质假设：把流体当做是由密集质点构成的，内部无空隙的连续体来研究。 3 流体力学的研究方法：理论、数值、实验。 4 作用于流体上面的力 （1）表面力：通过直接接触，作用于所取流体表面的力。 作用于A 上的平均压应力 作用于A 上的平均剪应力 应力 法向应力 切向应力 （2）质量力：作用在所取流体体积内每个质点上的力，力的大小与流体的质量成比例。（常见的质量力： 重力、惯性力、非惯性力、离心力） 单位为 5 流体的主要物理性质 （1） 惯性：物体保持原有运动状态的性质。质量越大，惯性越大，运动状态越难改变。 常见的密度（在一个标准大气压下）： 4℃时的水 20℃时的空气 （2） 粘性 ΔF ΔP ΔT A ΔA V τ 法向应力周围流体作用 的表面力 切向应力 A P p ??=A T ??=τA F A ??=→?lim 0δA P p A A ??=→?lim 0为A 点压应力，即A 点的压强 A T A ??=→?lim 0τ 为A 点的剪应力 应力的单位是帕斯卡（pa ） ，1pa=1N/㎡，表面力具有传递性。 B F f m =u u v v 2m s 3 /1000m kg =ρ3 /2.1m kg =ρ
牛顿内摩擦定律： 流体运动时，相邻流层间所产生的切应力与剪切变形的速率成正比。即 以应力表示 τ—粘性切应力，是单位面积上的内摩擦力。由图可知 —— 速度梯度，剪切应变率（剪切变形速度） 粘度 μ是比例系数，称为动力黏度，单位“pa ·s ”。动力黏度是流体黏性大小的度量，μ值越大，流体越粘，流动性越差。 运动粘度 单位：m2/s 同加速度的单位 说明： 1）气体的粘度不受压强影响，液体的粘度受压强影响也很小。 2）液体 T ↑ μ↓ 气体 T ↑ μ↑ 无黏性流体 无粘性流体，是指无粘性即μ=0的液体。无粘性液体实际上是不存在的，它只是一种对物性简化的力学模型。 （3） 压缩性和膨胀性 压缩性：流体受压，体积缩小，密度增大，除去外力后能恢复原状的性质。 T 一定，dp 增大，dv 减小 膨胀性：流体受热，体积膨胀，密度减小，温度下降后能恢复原状的性质。 P 一定，dT 增大，dV 增大 A 液体的压缩性和膨胀性 液体的压缩性用压缩系数表示 压缩系数：在一定的温度下，压强增加单位P ，液体体积的相对减小值。 由于液体受压体积减小，dP 与dV 异号，加负号，以使к为正值；其值愈大，愈容易压缩。к的单位是“1/Pa ”。（平方米每牛） 体积弹性模量K 是压缩系数的倒数，用K 表示，单位是“Pa ” 液体的热膨胀系数：它表示在一定的压强下，温度增加1度，体积的相对增加率。 du T A dy μ =? dt dr dy du ? =?=μ μτdu u dy h =ρ μν= dP dV V dP V dV ? -=-=1/κρ ρ κ d dP dV dP V K =-==1
北航计算流体力学第 15课
进口 出口
n n n 外边界 l l 外流边界形状 n n n 周期边 进口边界 出口边界 （b ）叶栅流 n n n n 进口边界 出口边界 （a ）通道流 固体壁 内流边界形状
二．几个重要概念 边界条件的定义： 边界条件表示求解域外的信息（扰动）对求解域边界的影响。 确定边界条件的原则： 1．若一信息由边界传入求解域，就应指定该信息的边界条件（第一原则）； 2．若一信息由求解域内传出边界，则不应指定该信息的边界条件（第二原则）。 由第一原则确定的边界条件称为解析边界条件； 由第二原则确定不给边界条件，但在数值求解中必须补充的边界条件称为数值边界条件。 由于信息传播的方式由方程的类型所决定，所以边界条件如何确定是由方程的类型所决定的。 又由于信息（扰动）是沿特征线传播的，所以边界条件的确定与特征线与边界交汇的方式有关。
进口 出口 三．进口与出口条件 （一） 一维Euler 方程 0t x U F += 式中： U u e ρρ?? ??=?? ???? ()2u F u p e p u ρρ???? =+????+?? 补充状态方程 21 12 p e u ργ= +- 1．进口边界（用下标 “in ”表示） 1）超音流（u a ∞>） 3个解析边界条件均由来流条件决定，即 in u u ∞= ，in ρρ∞= ，in p p ∞= 2）亚音流（u a ∞<） 2个解析边界条件，1个数值边界条件 in u u ∞= ，in ρρ∞= ，in inner p p = 下标inner 表示内场值。
流体力学总结 [题型]:简答题 流体静压强实验的操作步骤, 答案:(1)搞清仪器组成及其用法;(2)检查仪器是否密封，将阀门关闭，加压后检查测管液面高程是否恒定，若下降，表明漏气，应查明原因并加以处理;(3)量测点静压强(各点压强用厘米水柱高表示);(4)打开排气阀，记录水箱液面标高?0和各测压管液面标高?H (?H =0);(5) 关闭排气阀，用加压器缓慢加压，U形管出现压差?h。在加压的同时，观察左侧? A1、A2、B1、B2管的液柱上升情况，由于水箱内部的压强向各个方向传递，在左侧的测压管中，可以看到由于A、B两点在水箱内的淹没深度h不同，在压强向各点传递时，先到A点后到B点。在测压管中反应出的是A1管的液柱先上升，而B1管的液柱滞后一点也在上升，当停止加压时，A1、B1两点在同一水平面上， A2、B2两点与水箱内液面在同一水平面上，测记?0及各?H(此过程反复进行3 次;(6)打开排气阀，使液面恢复到同一水平面上。关闭排气阀，打开密闭容器底部的水门，放出一部分水，造成容器内压力下降，观察U形管中液柱的变化情况，测记?0及各?H(此过程反复进行3次)。 难度:1 分数:100 所属知识点: 知识体系/建筑系/土木专业/《流体力学》/流体动力学 [题型]:简答题 雷诺方程演示实验的操作步骤, 答案:(1)测记本实验的有关常数;(2)观察两种流态;(3)打开开关3使水箱充水至溢流水位，经稳定后，微微开启调节阀9，并注入颜色水于实验管内，使颜色水
流成一直线。通过颜色水质点的运动观察管内水流的层流流态，然后逐步开大调节阀，通过颜色水直线的变化观察层流转变到紊流的水力特征，待管中出现完全紊流后，再逐步关小调节阀，观察由紊流转变为层流的水力特征;(4)测定下临界雷诺数;(5) 将调节阀打开，使管中呈完全紊流，再逐步关小调节阀使流量减小。当流量调节到使颜色水在全管刚呈现出一稳定直线时，即为下临界状态;(6)待管中出现临界状态时，用体积法或重量法测定流量;(7)根据所测流量计算下临界雷诺数，并与公认值(2320)比较，偏离过大，需重测;(8)重新打开调节阀，使其形成完全紊流，按照上述步骤重复测量不少于三次;(9)同时用水箱中的温度计测记水温，从而求得水的运动粘度。难度:1 分数:100 所属知识点: 知识体系/建筑系/土木专业/《流体力学》/流体动力学 [题型]:简答题 流体力学综合实验的操作步骤, 答案:(1)测记本实验的有关常数;(2)打开电子调速器开关，使恒压水箱充水，排除实验管道中的滞留气体，待水箱溢流后，检查泄水阀全关时，各测压管液面是否齐平，若不平，则需排气调平;(3)打开泄水阀至最大开度，待流量稳定后，测记测压管读数，同时用体积法测记流量;(4)改变泄水阀开度3,4次，分别测记测压管读数及流量;(5) 实验完成后关闭泄水阀，检查测压管液面齐平后再关闭进水阀。 难度:2 分数:100 所属知识点: 知识体系/建筑系/土木专业/《流体力学》/流体动力学 [题型]:简答题 能量方程演示实验的操作步骤,
2015北京航空航天大学工程力学考博（航空科学与工程学院）参考 书、历年真题、报录比、研究生招生专业目录、复试分数线 一、学院介绍 航空科学与工程学院（以下简称航空学院）是北航最具航空航天特色的学院之一，主要从事大气层内各类航空器（飞机、直升机、飞艇等）、临近空间飞行器、微小型飞行器等的总体、气动、结构、强度、飞行力学与飞行安全、人机环境控制、动力学与控制等方面的基础性、前瞻性、工程性以及新概念、新理论、新方法研究与人才培养工作。 航空学院前身是清华大学航空系，是1952年北航成立时最早的两个系之一，当时称飞机系（设飞机设计和飞机工艺专业），1958年更名为航空工程力学系，1970年更名为五大队，1972年更名为五系，1989年定名为飞行器设计与应用力学系，2003年成立航空科学与工程学院。早期的航空学院荟萃了一批当时国内著名的航空领域的专家，如屠守锷、王德荣、陆士嘉、沈元、王俊奎、吴礼义、张桂联、徐鑫福、徐华舫、何庆芝、伍荣林、史超礼、叶逢培等教授，屠守锷院士（两弹一星元勋）是首任系主任，他们为本院发展奠定了坚实基础。在北航发展史上，航空学院不断输出专业和人才，先后参与组建七系、三系、十四系、宇航学院、飞行学院、无人机所、土木工程系、交通学院等院系。 自建校以来60多年，学院已培养本科毕业生万余人，硕士毕业生两千余人，博士毕业生近千人。毕业生中涌现出王永志、戚发韧、崔尔杰、乐嘉陵、王德臣、张福泽、王浚、钟群鹏、陶宝祺、郭孔辉、赵煦、唐西生、郭孔辉、唐长红等14位两院院士，改革开放后毕业生中也涌现出了“航空报国英模”/原沈飞董事长罗阳、中国商飞董事长金壮龙、第十一届“中国十大杰出青年”/原“神舟”飞船总指挥袁家军、歼15等飞机型号总师孙聪、C919大型客机总师吴光辉以及李玉海、耿汝光、姜志刚、屠恒章、孙聪、方玉峰、王永庆、孙兵、曲景文、李东、余后满、傅惠民、秦福光、陈元先、宋水云、吴宗琼、陈敏、高云峰等一批航空航天院所的年轻总师、总指挥、省市及部门负责人、民营企业家，为我国航空航天、国防事业及国家发展做出突出贡献。 学院作为主力曾先后研制成功我国第一架轻型旅客机“北京一号”、国内第一架高空高速无人侦察机、靶机、蜜蜂系列轻型飞机和第一架共轴式双旋翼直升机等，创造了多项全国第一。学院参与了所有国家重点航空型号以及大部分导弹型号的攻关工作。60多年来，学院取得了上百项国家和省部级教学与科研成果，其中国家级奖20多项。 学院师资力量雄厚，在北航乃至全国同类及相近学院中名列前茅。学院有教授56名（其中博士生导师51名），副教授50名，青年教师中有博士学位的比例为97%。拥有许多国内外著名专家学者，如中国科学院院士高镇同教授、李天教授，中国工程院院士李椿萱教授、王浚教授，“长江学者”特聘教授傅惠民、孙茂、杨嘉陵、高以天、武哲、王晋军、向锦武教授，国家教学名师及“万人计划”王琪教授，杰出青年基金获得者4名，跨/新世纪优秀人才的获得者10名，全国百篇优秀博士学位论文获得者2名；有国家级教学基地2个、国
《计算流体力学》课程大作业 ——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟 张伊哲 航博101 1、 引言和综述 2、 问题的提出，怎样使用涡量-流函数方法建立差分格式 3、 程序说明 4、 计算结果和讨论 5、 结论 1引言 虽然不可压缩流动的控制方程从形式上看更为简单，但实际上，目前不可压缩流动的数值方法远远不如可压缩流动的数值方法成熟。 考虑不可压缩流动的N-S 方程： 01()P t νρ??=? ? ??+??=-?+???? U U UU f U (1.1) 其中ν是运动粘性系数，认为是常数。将方程组写成无量纲的形式： 01()Re P t ??=?? ??+??=-?+????U U UU f U (1.2) 其中Re 是雷诺数。 从数学角度看，不可压缩流动的控制方程中不含有密度对时间的偏导数项，方程表现出椭圆-抛物组合型的特点；从物理意义上看，在不可压缩流动中，压力这一物理量的波动具有无穷大的传播速度，它瞬间传遍全场，以使不可压缩条件在任何时间、任何位置满足，这就是椭圆型方程的物理意义。这就造成不可压缩的N-S 方程不能使用比较成熟的发展型．．．偏微分方程的数值求解理论和方法。 如果将动量方程和连续性方程完全耦合求解，即使使用显示的离散格式，也将会得到一个刚性很强的、庞大的稀疏线性方程组，计算量巨大，更重要的问题是不易收敛。因此，实际应用中，通常都必须将连续方程和动量方程在一定程度上解耦。 目前，求解不可压缩流动的方法主要有涡量-流函数法，SIMPLE 法及其衍生的改进方法，有限元法，谱方法等，这些方法各有优缺点。其中涡量-流函数法是解决二维不可压缩流动的有效方法。作者本学期学习了研究生计算流体课程，为了熟悉计算流体的基本方法，选择使用涡量-流函数法计算不可压缩方腔驱动流问题，并且对于不同雷诺数下的解进行比较和分析，得出一些结论。 本文接下来的内容安排为：第2节提出不可压缩方腔驱动流问题，并分析该问题怎样使用涡量-流函数方法建立差分格式、选择边界条件。第3节介绍程序的结构。第4节对于不同雷诺数下的计算结果进行分析，并且与U.GHIA 等人【1】的经典结论进行对比，评述本
实验流体力学 第一章：相似理论和量纲分析 ①流体力学相似？包括几方面内容？有什么意义？ 流体力学相似是指原型和模型流动中，对应相同性质的物理量保持一定的比例关系，且对应矢量相互平行。 内容包括： 1.几何相似—物体几何形状相似，对应长度成比例； 2.动力相似—对应点力多边形相似，同一性质的力对应成比例并相互平行 （加惯性力后，力多边形封闭）； 3.运动相似—流场相似，对应流线相似，对应点速度、加速度成比例。 ②什么是相似参数？举两个例子并说明其物理意义 必须掌握的相似参数：Ma ，Re ，St 。知道在什么流动条件下必须要考虑这些相似参数。 相似参数又称相似准则，是表征流动相似的无量纲特征参数 。 1.两物理过程或系统相似则所有对应的相似参数相等。例如：假定飞机缩比模型风洞试验可以真正模拟真实飞行，则原型和模型之间所有对应的相似参数都相等，其中包括C L , C D , C M ： S V L C L 22 1 ρ= S V D C D 22 1 ρ= Sb V M C M 22 1 ρ= 风洞试验可以测得CL, CD, CM 值，在此基础上，将真实飞行条件带入CL, CD, CM 表达式，可以求得真实飞行的升力、阻力和力矩等气动性能参数。 2.所有对应的相似参数相等且单值条件相似则两个物理过程或系统相似。例如：对于战斗机超音速风洞试验，Ma 和Re 是要求模拟的相似参数，但通常在常规风动中很难做到。 由于对于此问题，Ma 影响更重要，一般的方案是保证Ma 相等，对Re 数影响进行修正。 ; Re V p Ma a RT a V L l St V ρ ρωμ∞∞= ====
《流体力学考》考点重点知识归纳 1.流体元：就有线尺度的流体单元，称为流体“质元”，简称流体元。流体元可看做大量流体质点构成的微小单元。 2.流体质点：（流体力学研究流体在外力作用下的宏观运动规律） （1）流体质点无线尺度，只做平移运动 （2）流体质点不做随即热运动，只有在外力的作用下作宏观运动； （3）将以流体质点为中心的周围临街体积的范围内的流体相关特性统计的平均值作为流体质点的物理属性； 3.连续性介质模型的内容：根据流体指点概念和连续介质模型，每个流体质点具有确定的宏观物理量，当流体质点位于某空间点时，若将流体质点的物理量，可以建立物理的空间连续分布函数，根据物理学基本定律，可以建立物理量满足的微分方程，用数学连续函数理论求解这些方程，可获得该物理量随空间位置和时间的连续变化规律。 4.连续介质假设：假设流体是有连续分布的流体质点组成的介质。 5.牛顿的粘性定律表明：牛顿流体的粘性切应力与流体的切变率成正比，还表明对一定的流体，作用于流体上的粘性切应力由相邻两层流体之间的速度梯度决定的，而不是由速度决定的： 6.牛顿流体：动力粘度为常数的流体称为牛顿流体。 7.分子的内聚力：当两层液体做相对运动时，两层液体的分子的平均距离加大，分子间的作用力变现为吸引力，这就是分子的内聚力。 液体快速流层通过分子内聚力带动慢流层，漫流层通过分子的内聚力阻滞快流层的运动，表现为内摩擦力。、 流体在固体表面的不滑移条件：分子之间的内聚力将流体粘附在固体表面，随固体一起运动或静止。 8.温度对粘度的影响：温度对流体的粘度影响很大。液体的粘度随温度升高而减小，气体的粘度则相反，随温度的升高而增大。 压强对粘性的影响：压强的变化对粘度几乎没有什么影响，只有发生几百个大气压的变化时，粘度才有明显改变，高压时气体和液体的粘度增大。 9.描述流体运动的两种方法 拉格朗日法：拉格朗日法又称为随体法。它着眼于流体质点，跟随流体质点一起运动，记录流体质点在运动过程中会各种物理量随所到位置和时间的变化规律，跟中所有质点便可了解整个流体运动的全貌。 欧拉法：欧拉法又称当地法。它着眼于空间点，把流体的物理量表示为空间位置和时间的函数。空间点的物理量是指，某个时刻占据空间点的。 流体质点的物理量，不同时刻占据该空间点的流体质点不同。 10.速度场：速度场是由流体空间各个坐标点的速度矢量构成的场。速度场不仅描述速度矢量的空间分布，还可描述这种分布随时间的变化。 11.毛细现象：玻璃管内的液体在表面张力的作用下液面升高或降低的现象称为毛细现象； 12.迹线：流体质点运动的轨迹。在流场中对某一质点作标记，将其在不同时刻的所在位置点连成线就是该流体质点的迹线。 13.定常流动：流动参数不随时间变化的流动。反之，流体参数随时间变化的流动称为不定长流动。 14.流线：流线是指示某一时刻流场中各点速度矢量方向的假象曲线。