计算流体力学过渡到编程的傻瓜入门教程
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较便。
CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。
回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分程的时候提出的一个问题，用一维无粘可压缩Euler程就可以描述了。
这里
这个程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量）随空间变化（）的关系。
在CFD常把这个程写成矢量形式
这里
进一步可以写成散度形式
一定要熟悉这种矢量形式
以上是控制程，下面说说求解思路。可压缩流动计算中，有限体积（FVM）是最广泛使用的算法，其他算法多多少少都和FVM有些联系或者共通的思路。了解的FVM，学习其他高级点的算法（比如目前比较热门的间断有限元、谱FVM、谱FDM）就好说点了。
针对一个微元控制体，把Euler程在空间积分
用微积分知识可以得到
也就是说控制体气体状态平均值的变化是控制体界面上流通量的结果。因此我们要计算的演化，关键问题是计算控制体界面上的。
FVM就是以这个积分关系式出发，把整个流场划分为多小控制体，每个控制体和围相邻的某个控制体共享一个界面，通过计算每个界面上的通量来得到相邻控制体之间的影响，一旦每个控制体的变化得到，整个流场的变化也就知道了。
所以，再强调一次，关键问题是计算控制体界面上的。
初学者会说，这个不难，把界面上的插值得到，然后就可以计算。有道理！
咱们画个图，有三个小控制体i-1到i+1，中间的“|”表示界面，控制体i
右边的界面用表示，左边的就是。
| i-1 | i | i+1 |
好下个问题：每个小控制体长度都是如插值计算界面上的？最自然的想法就是：取两边的平均值呗，
但是很不幸，这是不行的。
那么换个法？直接平均得到？
还是很不行，这样也不行。
我靠，这是为什么？这明明是符合微积分里面的知识啊？
这个道理有点复杂，说开了去可以讲一本书，可以说从50年代到70年代，CFD科学家就在琢磨这个问题。这里，初学者只需要记住这个结论：对于流动问题，不可以这样简单取平均值来插值或者差分。如果你非要想知道这个究竟，我现在也不想给你讲清楚，因为我眼下的目的是让你快速上手，而且该不刨根问底的时候就不要刨根问底，这也是初学阶段一种重要的学习法。
好了，既然目的只是为了求，我在这里，只告诉你一种计算法，也是非常重要、非常流行的一种法。简单的说，就是假设流动状态在界面
是不连续的，先计算出界面两边的值，和，再由它们用某种法计算出。上述法是非常重要的，是由一个联人Godunov在50年代首创的，后来被发展成为通用Godunov法，著名的ENO/WENO就是其中的一种。
好了，现在的问题是：
1 怎么确定和
2 怎么计算
对于第一个问题，Godunov在他的论文中，是假设每个控制体中是均匀分布的，因此
第二个问题，Godunov采用了精确黎曼解来计算。什么是“精确黎曼解”，就是计算这个激波管问题的精确解。既然有精确解，那还费功夫搞这些FVM算法干什么？因为只有这种简单一维问题有精确解，稍微复杂一点就不行了。精确解也比较麻烦，要分析5种情况，用牛顿法迭代求解（牛顿法是什么？看数值计算的书去，哦，算了，现在暂时可以不必看）。
这是最初Godunov的法，后来在这个思想的基础上，各种变体都出来了。也不过是在这两个问题上做文章，怎么确定，怎么计算。Godunov假设的是每个小控制体是均匀分布，也就是所谓分段常数(piecewise constant)，所以后来有分段线性(picewise linear)或者分段二次分布(picewise parabolic)，到后来ENO/WENO出来，那这个假设的多项式次
数就继续往上走了。都是用多项式近似的，这是数值计算中的一个强大工具，你可以在很多地看到这种近似。
Godunov用的四精确黎曼解，太复杂太慢，也不必要，所以后来就有各种近似黎曼解，最有名的是Roe求解器、HLL求解器和Osher求解器，都是对精确黎曼解做的简化。
这个多项式的阶数是和计算精度密切相关的，阶数越高，误差就越小。不过一般来说，分段线性就能得到不错的结果了，所以工程中都是用这个，Fluent、Fastran以及NASA的CFL3D、OverFlow都是用这个。而黎曼求解器对精度的影响不是那么大，但是对整个算法的物理适用性有影响，也就是说某种近似黎曼求解器可能对某些流动问题不合适，比如单纯的Roe对于钝头体的脱体激波会算得乱七八糟，后来加了熵修正才算搞定。
上次（https://www.360docs.net/doc/4c8136784.html,/node/399）说到了求解可压缩流动的一个重要算法，通用Godunov法。其两个主要步骤就是
1 怎么确定和
2 怎么计算
这里我们给出第一点一个具体的实现法，就是基于原始变量的MUSCL格式（以下简称MUSCL）。它是一种很简单的格式，而且具有足够的精度，NASA著名的CFL3D软件就是使用了这个格式，大家可以去它的主页
（https://www.360docs.net/doc/4c8136784.html,/Cfl3dv6/cfl3dv6.html）上看手册，里面空间离散那一章清楚的写着。
MUSCL假设控制体原始变量（就是）的分布是一次或者二次多项式，如果得到了这个多项式，就可以求出控制体左右两个界面的一侧的值和。
我们以压力p为例来说明怎么构造这个多项式。这里我只针对二次多项式来讲解，你看完之后肯定能自己推导出一次多项式的结果（如果你搞不定，那我对你的智商表示怀疑）。
OK，开始
假设，这个假设不影响最终结论，因为你总可以把一个区间线性的变换到长度为1的区间。
假设压力p在控制体i部的分布是一个二次多项式，控制体i的中心处于处，左右两个界面就是和。
这里先强调一个问题，在FVM中，每个控制体的求解出来的变量实际上是这个控制体的平均值。
所以，
。
我们知道，和，等距网格情况下和处的导数可以近似表示为
那么
由上述三个有关a，b和c的程，我们可以得到
这样就可以得到f(x)的表达式了，由此可以算出和
通常MUSCL格式写成如下形式
对应我们的推导结果（二次多项式假设）。
但是这不是最终形式。如果直接用这个公式，就会导致流场在激波（间断）附近的振荡。因为直接用二次多项式去逼近一个间断，会导致这样的效果。所以科学家们提出要对间断附近的斜率有所限制，因此引入了一个非常重要的修改——斜率限制器。加入斜率限制器后，上述公式就有点变化。
这里是Van Albada限制器
中国石油大学流体力学实验报告
中国石油大学（流体力学）实验报告 实验日期：2012-2-15 成绩： 班级：学号：姓名：教师： 同组者： 实验一、流体静力学实验 一、实验目的 1．掌握用液式测压计测量流体静压强的技能； 2．验证不可压缩流体静力学基本方程，加深对位置水头、压力水头和测压管水头的理解； 3. 观察真空度（负压）的产生过程，进一步加深对真空度的理解； 4．测定油的相对密度； 5．通过对诸多流体静力学现象的实验分析，进一步提高解决实际问题的能力。 二、实验装置 1、在图1-1-1下方的横线上正确填写实验装置各部分的名称 本实验的装置如图所示。 1.测压管； 2.带标尺的测压管； 3.连通管； 4.通气阀； 5.加压打气球； 6.真空测压管； 7 截止阀.；8. U形测压管；9.油柱；
10.水柱；11.减压放气阀 图1-1-1流体静力学实验装置图 2、说明 1．所有测管液面标高均以测压管2标尺零读数为基准； 2．仪器铭牌所注B ?、C ?、D ?系测点B 、C 、D 标高；若同时取标尺零点作为静力学基本方程的基准，则B ?、C ?、D ?亦为B z 、C z 、D z ； 3．本仪器中所有阀门旋柄均以顺管轴线为开。 三、实验原理在横线上正确写出以下公式 1．在重力作用下不可压缩流体静力学基本方程 形式之一： const p =+ γ z （1-1-1a ） 形式之二： h p p γ+=0（1-1b ） 式中z ——被测点在基准面以上的位置高度； p ——被测点的静水压强，用相对压强表示，以下同； 0p ——水箱中液面的表面压强； γ——液体重度； h ——被测点的液体深度。 2. 油密度测量原理 当U 型管中水面与油水界面齐平（图1-1-2），取其顶面为等压面，有 01w 1o p h H γγ==（1-1-2） 另当U 型管中水面和油面齐平（图1-1-3），取其油水界面为等压面，则有 02w o p H H γγ+= 即 02w 2o w p h H H γγγ=-=-（1-1-3）
全日制工程硕士研究生培养方案-北航研究生院-北京航空航天大学
大型飞机高级人才培养班 航空工程全日制工程硕士研究生培养方案 一、适用类别或领域 航空工程（085232） 二、培养目标 材料工程、电子与通信工程、控制工程、航空工程领域全日制工程硕士 (以下简称航空工程等领域全日制工程硕士)是与以上各工程领域任职资格相联系的专业学位，主要为国民经济和国防建设等领域培养应用型、复合型高层次工程技术和工程管理人才。大飞机班旨在探索一条“以国家大型项目人才需求为索引，培养具有献身精神、团结协作精神、开拓创新精神的设计型和复合型人才”的研究生培养新模式，是北航研究生培养体系的一部分。 航空工程等领域全日制工程硕士培养的基本要求是： 1、坚持党的基本路线，热爱祖国、遵纪守法、品行端正、诚实守信、身心健康，具有良好的科研道德和敬业精神。 2、在本领域掌握坚实的基础理论和系统的专门知识，有较宽的知识面和较强的自立能力，具有大飞机设计、制造、运营、管理等领域需求的创造能力和工程实践能力。 3、掌握一门外国语。 三、培养模式及学习年限 1．航空工程等领域全日制工程硕士研究生培养实行导师负责制，或以导师为主的指导小组制，负责制订硕士研究生个人培养计划，选课、组织开题报告、论文中期检查、指导科学研究和学位论文，并与中国商飞、第一飞机设计研究院、西飞公司等航空企业联合培养，实行导师组指导。 2．硕士研究生一般用1学年完成课程学习，课程学习实行学分制，具体学习、考核及管理工作执行《北京航空航天大学研究生院关于研究生课程学习管理规定》。 3．专业实习是全日制工程硕士研究生培养中的重要环节，全日制工程硕士研究生在学期间，应保证不少于0.5年的工程实践。 4．学位论文选题应来源于航空工程等领域工程技术背景。鼓励实行双导师制，其中第一导师为校内导师，校外导师应是与本工程领域相关的专家，也可以根据学生的论文
流体力学实验报告
流体力学 实验指导书与报告 静力学实验 雷诺实验 中国矿业大学能源与动力实验中心
学生实验守则 一、学生进入实验室必须遵守实验室规章制度，遵守课堂纪律，衣着整洁，保持安静，不得迟到早退，严禁喧哗、吸烟、吃零食和随地吐痰。如有违犯，指导教师有权停止基实验。 二、实验课前，要认真阅读教材，作好实验预习，根据不同科目要求写出预习报告，明确实验目的、要求和注意事项。 三、实验课上必须专心听讲，服从指导教师的安排和指导，遵守操作规程，认真操作，正确读数，不得草率敷衍，拼凑数据。 四、预习报告和实验报告必须独自完成，不得互相抄袭。 五、因故缺课的学生，可向指导教师申请一次补做机会，不补做的，该试验以零分计算，作为总成绩的一部分，累计三次者，该课实验以不及格论处，不能参加该门课程的考试。 六、在使用大型精密仪器设备前，必须接受技术培训，经考核合格后方可使用，使用中要严格遵守操作规程，并详细填写使用记录。 七、爱护仪器设备，不准动用与本实验无关的仪器设备。要节约水、电、试剂药品、元器件、材料等。如发生仪器、设备损坏要及时向指导教师报告，属责任事故的，应按有关文件规定赔偿。 八、注意实验安全，遵守安全规定，防止人身和仪器设备事故发生。一旦发生事故，要立即向指导教师报告，采取正确的应急措施，防止事故扩大，保护人身安全和财产安全。重大事故要同时保护好现场，迅速向有关部门报告，事故后尽快写出书面报告交上级有关部门，不得隐瞒事实真相。 九、试验完毕要做好整理工作，将试剂、药品、工具、材料及公用仪器等放回原处。洗刷器皿，清扫试验场地，切断电源、气源、水源，经指导教师检查合格后方可离开。 十、各类实验室可根据自身特点，制定出切实可行的实验守则，报经系(院)主管领导同意后执行，并送实验室管理科备案。 1984年5月制定 2014年4月再修订 中国矿业大学能源与动力实验中心
北京航空航天大学五系流体力学实验染色液流动显示实验报告
研究生《流体力学实验》 ——飞机标模染色液流动显示 实验报告 班级 姓名 实验日期 指导教师 北京航空航天大学流体力学研究所
一、实验目的 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．实验结果和分析
流体力学实践报告
黑龙江科技大学建筑工程二学历实践报告 流体力学实践报告 一、实践概述 在此次实践中,老师给我演示了雷诺试验与伯努利方程试验。下面我就实践的主要内容进行一下总结。 二、雷诺实验 (一)、实验目的 1、观察液体流动时的层流与紊流现象。区分两种不同流态的特征,搞清两种流态产生的条件。分析圆管流态转化的规律,加深对雷诺数的理解。 2、测定颜色水在管中的不同状态下的雷诺数及沿程水头损失。绘制沿程水头损失与断面平均流速的关系曲线,验证不同流态下沿程水头损失的规律就是不同的。进一步掌握层流、紊流两种流态的运动学特性与动力学特性。 3、通过对颜色水在管中的不同状态的分析,加深对管流不同流态的了解。学习古典流体力学中应用无量纲参数进行实验研究的方法,并了解其实用意义。 (二)、实验原理 1、液体在运动时,存在着两种根本不同的流动状态。当液体流速较小时,惯性力较小,粘滞力对质点起控制作用,使各流层的液体质点互不混杂,液流呈层流运动。当液体流速逐渐增大,质点惯性力也逐渐增大,粘滞力对质点的控制逐渐减弱,当流速达到一定程度时,各流层
的液体形成涡体并能脱离原流层,液流质点即互相混杂,液流呈紊流运动。这种从层流到紊流的运动状态,反应了液流内部结构从量变到质变的一个变化过程。 液体运动的层流与紊流两种型态,首先由英国物理学家雷诺进行了定性与定量的证实,并根据研究结果,提出液流型态可用下列无量纲数来判断: Re=Vd/ν Re 称为雷诺数。液流型态开始变化时的雷诺数叫做临界雷诺数。 在雷诺实验装置中,通过有色液体的质点运动,可以将两种流态的根本区别清晰地反映出来。在层流中,有色液体与水互不混惨,呈直线运动状态,在紊流中,有大小不等的涡体振荡于各流层之间,有色液体与水混掺。 2、在如图所示的实验设备图中,取1-1,1-2两断面,由恒定总流的能量方程知: f 2222221111h g 2V a p z g 2V a p z ++γ+=+γ+ 因为管径不变V 1=V 2 ∴=γ +-γ+=)p z ()p z (h 2211f △h 所以,压差计两测压管水面高差△h 即为1-1与1-2两断面间的沿程水头损失,用重量法或体积浊测出流量,并由实测的流量值求得断面平均流速A Q V =,作为lgh f 与lgv 关系曲线,如下图所示,曲线上EC 段与BD 段均可用直线关系式表示,由斜截式方程得: lgh f =lgk+mlgv lgh f =lgkv m h f =kv m m 为直线的斜率
第二章计算流体力学的基本知识
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30～40 年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学" 。 从20 世纪60 年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
北航考研之科研成果及重点实验室汇总
北航考研之科研成果及重点实验室汇总 学校在尖端技术研究领域居于国内高校前列，有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技术研究所、可靠性工程研究所、外国语言研究所、设备工程
流体力学-伯努利方程实验报告
中国石油大学（华东）工程流体力学实验报告 实验日期：2014.12.11成绩： 班级：石工12-09学号：12021409姓名：陈相君教师：李成华 同组者：魏晓彤，刘海飞 实验二、能量方程（伯诺利方程）实验 一、实验目的 1．验证实际流体稳定流的能量方程； 2．通过对诸多动水水力现象的实验分析，理解能量转换特性； 3．掌握流速、流量、压强等水力要素的实验量测技能。 二、实验装置 本实验的装置如图2-1所示。 图2-1 自循环伯诺利方程实验装置 1.自循环供水器； 2.实验台； 3.可控硅无极调速器；4溢流板；5.稳水孔板； 6.恒压水箱； 7.测压机；8滑动测量尺；9.测压管；10.试验管道； 11.测压点；12皮托管；13.试验流量调节阀 说明 本仪器测压管有两种： （1）皮托管测压管（表2-1中标﹡的测压管），用以测读皮托管探头对准点的总水头； （2）普通测压管（表2-1未标﹡者），用以定量量测测压管水头。 实验流量用阀13调节，流量由调节阀13测量。
三、实验原理 在实验管路中沿管内水流方向取n 个过水断面。可以列出进口断面（1）至另一断面（i ）的能量方程式（i =2,3,…,n ） i w i i i i h g v p z g p z -++ + =+ + 1222 2 111 1αγυαγ 取12n 1a a a ==???==，选好基准面，从已设置的各断面的测压管中读出 z+p/r 值，测 出透过管路的流量，即可计算出断面平均流速，从而即可得到各断面测压管水头和总水头。 四、实验要求 1．记录有关常数实验装置编号 No._4____ 均匀段1d = 1.40-210m ?；缩管段2d =1.01-210m ?；扩管段3d =2.00-2 10m ?； 水箱液面高程0?= 47.6-2 10m ?；上管道轴线高程z ?=19 -2 10m ? （基准面选在标尺的零点上） 2．量测（p z γ + ）并记入表2-2。 注：i i i p h z γ =+ 为测压管水头，单位：-2 10m ，i 为测点编号。 3．计算流速水头和总水头。
室外风环境模拟计算报告123
新项目 室外风环境模拟计算报告 计算软件：风模拟分析软件PKPM-CFD 开发单位：中国建筑科学研究院 建研科技股份 合作单位：Software Cradle Co., Ltd. 韵能建筑科技 应用版本：Ver1.00 2015.10.19
室外风环境模拟分析报告 项目名称：新项目 项目地址： 建设单位： 设计单位： 参与单位： 规标准参考依据： 1、《绿色建筑评价标准》（GB/T 50378-2014） 2、《民用建筑设计通则》（GB 50352-2005） 3、《绿色建筑评价技术细则》
一、项目概述 1.1计算模型概况 1.2建筑物概况 图 1 建筑群平面图，红线建筑为目标建筑
二、指标要求 针对室外风环境评价依据为《绿色建筑评价标准》（GB/T 50378-2014）中有关室外风环境的条目要求。 2.1规的评价要求 《绿色建筑评价标准》（GB/T 50378-2014）中有关室外风环境的具体要求如下： 4.2.6 场地风环境有利于室外行走、活动舒适和建筑的自然通风。评分规则如下： 1 冬季典型风速和风向条件下，建筑物周围人行区风速低于5m/s，且室外风速放大系数小于2，得2分；除迎风第一排建筑外，建筑迎风面与背风面表面风压差不超过5Pa，再得1分。 2 过渡季、夏季典型风速和风向条件下，场地人活动区不出现涡旋或无风区，得2分；50%以上可开启外窗室外表面的风压差大于0.5Pa，得1分。 2.2模拟条件设置要求 1、室外风环境模拟的边界条件和基本设置需满足以下规定： 1）计算区域：建筑覆盖区域小于整个计算域面积3%；以目标建筑为中心，半径5H围为水平计算域。建筑上方计算区域要大于3H；H为建筑主体高度； 2）网格划分：建筑的每一边人行高度区1.5m或2m高度应划分10个网格或以上； 3）湍流模型选择：标准k-ε模型。高精度要求时采用Durbin模型或MMK模型。
计算流体力学过渡到编程的傻瓜入门教程
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较方便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分方程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子内面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分方程的时候提出的一个问题，用一维无粘可压缩Euler方程就可以描述了。 这里 这个方程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量 ）随空间变化（）的关系。 在CFD中通常把这个方程写成矢量形式 这里 进一步可以写成散度形式
最新北航计算流体力学第15课
北航计算流体力学第 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 表示内场值。
计算流体力学螺旋管分析报告
重庆大学《计算流体力学与计算传热学基础》上机实验水平螺旋管内的对流换热过程 学生：刘伟文 学号：20123000 指导教师：李隆键 专业：热能与动力工程 重庆大学动力工程学院 二O一五年六月
一、前言 螺旋管在热力、化工、石油及核工业等领域得到了广泛应用，螺旋管换热器也具有结构简单、传热系数高等优点。它的传热系数比直管高，在相同空间里可得到更大的传热面积，布置更长的管道，减少了焊缝，提高了安全性。尽管螺旋管的流体阻力增大，压降增大，但是其传热效率的提高导致能量的节约要高于因阻力增大而消耗的能量。因此，螺旋管在许多行业得到普遍应用而倍受青睐。在工程应用中，由于工艺要求，往往需将流体加热至规定的温度范围，传热是其中的基本单元操作，所以有必要对螺旋管的传热与流动特性进行研究。从理论知识我们知道由于向心力的作用，流体从管中心部分由螺旋管内侧流向外侧壁面，因而造成了螺旋管内侧的低压区。在压差作用下，流体从外侧沿着圆管的上部和下部壁面流回内侧。这种流动是与管的轴向垂直的，也就是与流体的主体流动相垂直，称为二次流。流体的这种二次流与轴向主流复合成螺旋式的前进运动。这样，对于流体的传热传质，不仅可依靠流体的径向扩散，还有径向二次流的作用，相当于边界层进行了破坏，增强了流体传质。 二、GAMBIT建模
1、先建立一个半径为6的圆面。 2、将该圆面向X轴正方向移动120。 3、用圆面sweep形成螺旋柱体。（绕Y轴正方向）
4、重复以上操作，得到如图所示几何体弯管。 5、设置边界层。
并应用至每个截面：
6、设置圆面的网格，选择pave方式，interval size 选择0.6，这样边界层网格与圆面中心网格过渡较平缓。 7、依次建立体网格。 8、检查网格质量。 最差网格为0.41，满足要求。 8、输出网格。
2015北航工程力学考博(航空科学与工程学院)参考书、历年真题、报录比、研究生招生专业目录、复试分数线
2015北京航空航天大学工程力学考博（航空科学与工程学院）参考 书、历年真题、报录比、研究生招生专业目录、复试分数线 一、学院介绍 航空科学与工程学院（以下简称航空学院）是北航最具航空航天特色的学院之一，主要从事大气层内各类航空器（飞机、直升机、飞艇等）、临近空间飞行器、微小型飞行器等的总体、气动、结构、强度、飞行力学与飞行安全、人机环境控制、动力学与控制等方面的基础性、前瞻性、工程性以及新概念、新理论、新方法研究与人才培养工作。 航空学院前身是清华大学航空系，是1952年北航成立时最早的两个系之一，当时称飞机系（设飞机设计和飞机工艺专业），1958年更名为航空工程力学系，1970年更名为五大队，1972年更名为五系，1989年定名为飞行器设计与应用力学系，2003年成立航空科学与工程学院。早期的航空学院荟萃了一批当时国内著名的航空领域的专家，如屠守锷、王德荣、陆士嘉、沈元、王俊奎、吴礼义、张桂联、徐鑫福、徐华舫、何庆芝、伍荣林、史超礼、叶逢培等教授，屠守锷院士（两弹一星元勋）是首任系主任，他们为本院发展奠定了坚实基础。在北航发展史上，航空学院不断输出专业和人才，先后参与组建七系、三系、十四系、宇航学院、飞行学院、无人机所、土木工程系、交通学院等院系。 自建校以来60多年，学院已培养本科毕业生万余人，硕士毕业生两千余人，博士毕业生近千人。毕业生中涌现出王永志、戚发韧、崔尔杰、乐嘉陵、王德臣、张福泽、王浚、钟群鹏、陶宝祺、郭孔辉、赵煦、唐西生、郭孔辉、唐长红等14位两院院士，改革开放后毕业生中也涌现出了“航空报国英模”/原沈飞董事长罗阳、中国商飞董事长金壮龙、第十一届“中国十大杰出青年”/原“神舟”飞船总指挥袁家军、歼15等飞机型号总师孙聪、C919大型客机总师吴光辉以及李玉海、耿汝光、姜志刚、屠恒章、孙聪、方玉峰、王永庆、孙兵、曲景文、李东、余后满、傅惠民、秦福光、陈元先、宋水云、吴宗琼、陈敏、高云峰等一批航空航天院所的年轻总师、总指挥、省市及部门负责人、民营企业家，为我国航空航天、国防事业及国家发展做出突出贡献。 学院作为主力曾先后研制成功我国第一架轻型旅客机“北京一号”、国内第一架高空高速无人侦察机、靶机、蜜蜂系列轻型飞机和第一架共轴式双旋翼直升机等，创造了多项全国第一。学院参与了所有国家重点航空型号以及大部分导弹型号的攻关工作。60多年来，学院取得了上百项国家和省部级教学与科研成果，其中国家级奖20多项。 学院师资力量雄厚，在北航乃至全国同类及相近学院中名列前茅。学院有教授56名（其中博士生导师51名），副教授50名，青年教师中有博士学位的比例为97%。拥有许多国内外著名专家学者，如中国科学院院士高镇同教授、李天教授，中国工程院院士李椿萱教授、王浚教授，“长江学者”特聘教授傅惠民、孙茂、杨嘉陵、高以天、武哲、王晋军、向锦武教授，国家教学名师及“万人计划”王琪教授，杰出青年基金获得者4名，跨/新世纪优秀人才的获得者10名，全国百篇优秀博士学位论文获得者2名；有国家级教学基地2个、国
流体力学报告
流体力学报告 每一门力学学科的建立,都需要建立模型,也就是把实际的问题抽象化,而抽象过程就是把现实中对所研究问题不重要的因素忽略掉,也就是模型假设,从而建立于这个问题相适应的模型进行研究,如果有意义有价值,也就慢慢深入研究,从而形成一门学科,它们都是随社会的发展而发展形成的.比如现如今最前沿的力学学科"纳米力学"就是如此。我们土木工程常说的三大力学有：1.理论力学---分析力学,振动力学,水力学或称为流体力学(这些研究对材料都不太侧重 )2.材料力学---弹性力学,塑性力学(都是又材料特性而分的) 3.结构力学:就是分析复杂的结构的情形。在此我重点叙述我对流体力学这门课学科的学习和认知。 一·流体的基本信息解释： 流体，是与固体相对应的一种物体形态，是液体和气体的总称. 由大量的、不断地作热运动而且无固定平衡位置的分子构成的，它的基本特征是没有一定的形状并且具有流动性。流体都有一定的可压缩性，液体可压缩性很小，而气体的可压缩性较大，在流体的形状改变时，流体各层之间也存在一定的运动阻力(即粘滞性)。当流体的粘滞性和可压缩性很小时，可近似看作是理想流体，它是人们为研究流体的运动和状态而引入的一个理想模型。是液压传动和气压传动的介质。大气和水是最常见的两种流体，大气包围着整个地球，地球表面的70%是水面。大气运动、海水运动(包括波浪、潮汐、中尺度涡旋、环流等)乃至地球深处熔浆的流动都是流体的研究内容。
二·流体力学的阐述： 流体力学是连续介质力学的一门分支，是研究流体(包含气体，液体以及等离子态)现象以及相关力学行为的科学。可以按照研究对象的运动方式分为流体静力学和流体动力学，还可按流动物质的种类分为水力学，空气动力学等等。对流体力学学科的形成作出第一个贡献的是古希腊的阿基米德，他建立了包括物理浮力定律和浮体稳定性在内的液体平衡理论，奠定了流体静力学的基础，特别是从20世纪以来，流体力学已发展成为基础科学体系的一部分，同时又在工业、农业、交通运输、天文学、地学、生物学、医学等方面得到广泛应用。流体力学既包含自然科学的基础理论，又涉及工程技术科学方面的应用。此外，如从流体作用力的角度，则可分为流体静力学、流体运动学和流体动力学;从对不同"力学模型"的研究来分，则有理想流体动力学、粘性流体动力学、不可压缩流体动力学、可压缩流体动力学和非牛顿流体力学等。 三·对流体的研究假设： 连续体假设 物质都由分子构成,尽管分子都是离散分布的,做无规则的热运动.但理论和实验都表明,在很小的范围内,做热运动的流体分子微团的统计平均值是稳定的.因此可以近似的认为流体是由连续物质构成,其中的温度,密度,压力等物理量都是连续分布的标量场。 质量守恒 质量守恒目的是建立描述流体运动的方程组。欧拉法描述为:流进
北航研究生课程实验流体力学重点
实验流体力学 第一章：相似理论和量纲分析 ①流体力学相似？包括几方面内容？有什么意义？ 流体力学相似是指原型和模型流动中，对应相同性质的物理量保持一定的比例关系，且对应矢量相互平行。 内容包括： 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 ρ ρωμ∞∞= ====
北航GPA算法
1、GPA算法 网上流传着各种各样的算法，但是需要强调的是，美国人知道中国的大部分学校不用GPA，因此相当多的学校在网申系统里明确指出，不允许将自己的成绩换算成美国的GPA，比如Caltech，Princeton，Stanford等，这一栏留着不填即可。有的学校则要求按照我们的评定成绩规则填写，北航是用百分制的平均分，那么我们可以填写在保研是用的大学前三年的必修课平均成绩，到时候教务会算好。还有很多的学校没做要求，就按照北航的GPA算法计算。 其实自己填写的GPA只是一个参考，可以写在简历里辅助申请，对方学校会按照他们的标准重新计算。可能具体的教授还会拿出你的某些重要的课程评估你的GPA，因此其实研究哪种算法算的更高没有意义，之所以这么说还有一点原因就是北航的学生在算GPA时必须严格按照自己学校的GPA算法计算。 2、北航的GPA算法 在开成绩单时，学校不给算GPA，但是在成绩单上有GPA的算法，那么在计算GPA是就应该按照这个算法来计算，不可以采用其他的算法。当然也没规定一定严格按照这个算法，但是既然写在了官方的成绩单上，就应该这么算，要不然会有作弊之嫌。 2.1、具体算法如下：85～100/A：4； 70～84/B：3； 60～69/C：2； 不及格/F：0； 按照通过与不通过评分的，算法如下： 通过/P： 3.3； 不通过/F：0； 例如，有三门课，学分分别为1、2、3，得分分别为86、76和通过，那么这三门课的GPA 就是（4×1+3×2+3.3×3）/（1+2+3）=3.32 总体来讲，北航的GPA算法还是很有优势的。 2.2、计算GPA的课程范围： 全部课程，包括所有的必修，任何形式的选修，只要是出现在成绩单上的都要算，大学前三年学过的所有课程都会出现在成绩单上。不要试图去修改成绩，北航也不允许去掉更不允许修改成绩。
流体力学试验
流體力學實驗 老師：Ａ1班→李宗翰老師；Ｂ1班→楊龍杰老師 Ａ2班→蔡欣正老師；Ｂ2班→邵德文老師 時間：Ａ1班→星期五12、13、14節；Ｂ1班→星期三11、12、13節Ａ2班→星期二11、12、13節；Ｂ2班→星期四12、13、14節 上課進度：
成績計算： 1.作業（30﹪）：上課後10分鐘未交報告者扣總分3分！ 當日無故未交者扣總分10分！ 1.課堂（20﹪）：分組合作精神，數據結果及隨堂口試小考。 2.口試（25﹪）：於考前一週公告口試方式。 3.筆試（25﹪）：於考前一週公告考場。 4.上課遲到10分鐘內扣總分3分！ 無故缺課扣總分10分！缺課3次下學期再見！ ※實驗前每組須備有空白數據表格一份，以方便記錄實驗數據※ 規定事項： 一、預習報告：（限用Ａ4大小的紙書寫，不可用打字） 1.封面：包含實驗名稱、組別、班級、姓名、學號、座號。 2.內容：包含實驗目的、實驗原理、實驗步驟及空白數據表格。 3.每人一份，於實驗前由組長收齊交給助教簽章，並於批閱後取回。 二、結論報告：（限用Ａ4大小的紙書寫） 1.個人結報：每人一份，含實驗心得和討論（心得須300字以上）。上課前將 個人結報及前一次實驗領回的預報合訂在一起，交給組長。 2.整組結報：每組一份，含數據、回歸分析結果，回歸分析圖表。 3.回歸分析須有電腦分析報表結果和座標曲線圖，圖可用手畫或電腦處理，若 用手畫請用方格紙，不可用工學院作業紙的背面。 4.未能及時繳交之作業，也一定要儘快繳交，不可缺交。 三、上課期間： 1.在實驗室內不可抽煙、進食及喝水，並注意安全。 2.不可無故離開實驗室，如有需要請先報備，助教會不定時的抽點。 3.組長負責整組的實驗操作、秩序及做完實驗後的清潔。 4.實驗後的數據表格，須在下課前交給助教檢查才算完成，嚴禁抄襲。
计算流体力学结课报告
计算流体力学结课报告200Km/h列车fluent仿真计算 学部：化、环、生学部 学院：化工机械与安全学院 学号：31507095 班级：化1512班 学生姓名：孙金
引言 数值仿真就是对所建立的数值模型进行数值实验和求解的过程。而计算流体力学CFD （Computational Fluid Dynamics）就是在工程仿真实验领域中应用最广泛的一门学科。任何流体运动的规律都是以质量守恒定律、动量守恒定律和能量守恒定律为基础的。这些基本定律可由数学方程组来描述，如欧拉方程、N-S方程。采用数值计算方法，通过计算机求解这些控制流体流动的数学方程，进而研究流体的运动规律这就是CFD研究问题的方法。在实际计算流体力学方面，采用通用的CFD软件来完成工程上的一些流体力学问题，有极为广泛的应用前景。近年来，随着计算机技术以及相关技术的发展，CFD技术已经在工程领域内取得重大的进步，特别是在高速列车的外型设计方面起了很大作用。随着国家经济的发展，国家运输业也有了很大的发展，特别是列车经过几次提速后，高速列车在国家运输行业中所占比例不断提高。高速列车的特点是庞大、细长、在地面轨道上运行，其空气动力学问题非常复杂。空气在列车表面形成空气流场，空气阻力急剧增加，作用在列车的阻力大部分来自压强阻力，而一部分来自表面磨擦阻力，这就使能耗过大，同时列车可能出现较大的空气升力，导致列车产生“飘”的现象，激发列车脱轨事故的发生，因此研究高速列车气动力性能非常重要。用CFD仿真可以详细了解高速列车的空气动力特性，从而设计出阻力小、噪音低等各方面性能完善的高质量列车。本文采用CFD学科中的常用商业软件Fluent仿真一个时速200km/h的二维流线型车头的外流场，对其空气动力性能进行分析，从而得到不同车辆形状其周围流场的不同，进而分析哪种车型更适合。
北航宇航学院航天工程培养方案
宇航学院 航天工程领域（085233） 全日制工程硕士研究生培养方案 一、适用领域 航天工程领域（085233） 二、培养目标 航天工程领域全日制工程硕士是与航天工程领域任职资格相联系的专业性学位，主要为国民经济和国防建设等领域培养应用型、复合型高层次工程技术和工程管理人才。要求掌握航天器总体设计、航天控制技术、航天推进技术的基本概念与理论，能以航天器/空间系统为研究对象，在设计与实现过程运用航天科学的理论与技术，进行系统总体设计、控制系统设计与分析、有效载荷设计与实现、推进系统设计、地球和探测新技术、实验与测试的高层次综合性研究。 三、培养模式及学习年限 航天工程领域全日制工程硕士主要采用校企联合培养，实行校企双方联合导师制，以校内导师指导为主，企业导师参与实践过程、项目研究、课程与论文等多个环节的指导工作。 本领域全日制工程硕士研究生遵循《北京航空航天大学研究生学籍管理规定》，学制一般为2.5年，实行弹性学习年限，一般在1年内完成课程学习，在企业工作时间累积不少于6个月。 全日制专业学位硕士研究生实行学分制，在攻读学位期间，要求在申请硕士学位论文答辩前，依据培养方案，获得知识和能力结构中所规定的各部分学分及总学分；要求全日制专业学位硕士研究生文献综述与开题报告至申请学位论文答辩的时间一般不少于6个月。 四、知识和能力结构 航天工程领域全日制工程硕士研究生培养方案的知识和能力结构由学位理论课程和综合实践环节两部分构成。学位课程的学习是研究生培养环节中的重要内容，学位课程的设置是以全面提高研究生在航天工程领域内的理论及专业知识水平、科学及人文素质、工程能力素质为目标。要秋取得航天工程领域全日制工程硕士学位的研究生必须按培养方案获得表中所规定的各部分学分及总学分，如下表所示。
流体力学计算软件报告
三维方管内部二次流特征分析 ——基于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 方程下残差收敛曲线