一、判断题 1、 根据牛顿内摩擦定律，当流体流动时，流体内部内摩擦力大小与该处的流速大小成正比。 2、 一个接触液体的平面壁上形心处的水静压强正好等于整个受压壁面上所有各点水静压强的平均 值。 3、 流体流动时，只有当流速大小发生改变的情况下才有动量的变化。 4、 在相同条件下，管嘴出流流量系数大于孔口出流流量系数。 5、 稳定（定常）流一定是缓变流动。 6、 水击产生的根本原因是液体具有粘性。 7、 长管是指运算过程中流速水头不能略去的流动管路。 8、 所谓水力光滑管是指内壁面粗糙度很小的管道。 9、 外径为D ，内径为d 的环形过流有效断面，其水力半径为4 d D -。 10、 凡是满管流流动，任何断面上的压强均大于大气的压强。 二、填空题 1、某输水安装的文丘利管流量计，当其汞-水压差计上读数cm h 4=?，通过的流量为s L /2，分析 当汞水压差计读数cm h 9=?，通过流量为 L/s 。 2、运动粘度与动力粘度的关系是 ，其国际单位是 。 3、因次分析的基本原理是： ；具体计算方法分为两种 。 4、断面平均流速V 与实际流速u 的区别是 。 5、实际流体总流的伯诺利方程表达式为 ， 其适用条件是 。 6、泵的扬程H 是指 。 7、稳定流的动量方程表达式为 。 8、计算水头损失的公式为 与 。 9、牛顿内摩擦定律的表达式 ，其适用范围是 。 10、压力中心是指 。 一、判断题 ×√×√× ×××√× 二、填空题 1、 3 L/s 2、 ρμν=，斯(s m /2 ) 3、 因次和谐的原理，п定理 4、 过流断面上各点的实际流速是不相同的，而平均流速在过流断面上是相等的 5、 22222212111 122z g v a p h g v a p z +++=++-γγ，稳定流，不可压缩流体，作用于流体上的质量力只有重力，所取断面为缓变流动 6、 单位重量液体所增加的机械能 7、 ∑?=F dA uu cs n ρ
第二章 2-1.已知某种物质的密度ρ=2.94g/cm3，试求它的相对密度d。 解：d=ρ/ρw=2.94（g/cm3）/1（g/cm3）=2.94 2-2.已知某厂1号炉水平烟道中烟气组分的百分数为α（co2）=13.5%，a(SO2)=0.3%，a(O2)=5.2%，a(N2)=76%，a(H2O)=5%。试求烟气的密度。 2-3.上题中烟气的实测温度t=170℃，实测静计压强Pe=1432Pa，当地大气压强 Pa=10058Pa。试求工作状态下烟气的密度和运动粘度。
2-4.当压强增量为50000Pa时，某种液体的密度增长0.02%，试求该液体的体积模量。 2-5.绝对压强为3.923×10^5Pa的空气的等温体积模量和等熵体积模量各等于多少？ 2-6. 充满石油的油槽内的压强为4.9033×10^5Pa,今由槽中排出石油40kg，使槽内压强降到9.8067×10^4Pa，设石油的体积模量K＝1.32×10^9 Pa。试求油槽的体积。 2-7. 流量为50m3/h，温度为70℃的水流入热水锅炉，经加热后水温升到90℃，而水的体胀系数αV=0.000641/℃，问从锅炉中每小时流出多少立方米的水？
2-8. 压缩机压缩空气，绝对压强从9.8067×104Pa升高到5.8840×105Pa，温度从20℃升高到78℃，问空气体积减少了多少？ 2-9. 动力粘度为2.9×10^-4Pa·S，密度为678kg/m3的油，其运动粘度等于多少？解：V=u/ρ=2.9×10^-4/678=4.28×10^-7m2/s 2-10. 设空气在0℃时的运动粘度ν0=13.2×10-6m2/s,密度ρ0=1.29kg/m3。试求在150℃时空气的动力粘度。 2-11. 借恩氏粘度计测得石油的粘度为8.5oE,如石油的密度为ρ=850kg/m3,试求石油的动力粘度。
1．流体的连续介质模型：研究流体的宏观运动，在远远大于分子运动尺度的范围里考察流体运动，而不考虑个别分子的行为，因此我们可以把流体视为连续介质。 它有如下性质： （1）流体是连续分布的物质，它可以无限分割为具有均布质量的宏观微元体。 （2）不发生化学反应和离解等非平衡热力学过程的运动流体中，微元体内流体状态服 从热力学关系 （3）除了特殊面外，流体的力学和热力学状态参数在时空中是连续分布的，并且通常 认为是无限可微的 2．应力：有限体的微元面积上单位面积的表面力称为表面力的局部强度，又称为应力，定义如下：=n T A F A δδδlim 0→ 3．流体的界面性质：微元界面两侧的流体的速度和温度相等，应力向量的大小相等.方向相反或应力分量相等。 4．流体具有易流行和压缩性。 5．应力张量具有对称性。 6．欧拉描述法：在任意指定的时间逐点描绘当地的运动特征量（如速度、加速度）及其它的物理量的分布（如压力、密度等）。 7．拉格朗日描述法：从某个时刻开始跟踪质点的位置、速度、加速度和物理参数的变化，这种方法是离散质点的运动描述法称为拉格朗日描述法。 8.流线：速度场的向量线，该曲线上的任意一点的切向量与当地的的速度向量重合。 迹线：流体质点点的运动迹象。 差别：迹线是同一质点在不同时刻的位移曲线。 流线是同一时刻、不同质点连接起来的速度场向量线。 流线微分方程：ω dz v dy u dx == 迹线微分方程：t x U i i ??= 9．质点加速度：质点速度向量随时间的变化率。 U U t U a )(??+??= 质点加速度=速度的局部导数+速度的迁移导数。 物理量的质点导数=物理量的局部导数+物理量的对流导数。
〈〈产品数据管理（PDM技术》课程简介 课程代码:AM011 课程简介： 本门课程将讲授PDM技术的基本概念、理论方法、系统结构和PDM^r业实施案例以及典型PDM^统介绍等相关专题，以满足我国企业信息化工程对大量复合型人才的需求 本课程的主要任务是： 1、掌握PDMJ术的发展与应用； 2、掌握PDMJ术的基本理论和方法； 3、掌握PDMK统体系结构和主要功能； 4、掌握PD"对象的建模方法和对象模型； 5、了解PDMK统实施方法； 6、了接国内外著名PDMS用系统。 本课程是一门实用性和系统性很强的课程，包含了机械工程和工业工程等领域高级技术 人员必须掌握的基本知识和内容。课程学习的目的是使学生掌握 PD M 基本理论知识和方法，为今后从事企业信息化工作，特别是从事产品数字化设计、制造与管理工作打下坚实的理论基础。 This course is the basic course on product development, it covers the following topics: Development and applications of PDM technology, Supporting technologies of PDM, Product data management technology, Product development lifecycle management technology, PDM implementation methodology, Introduction to SIPM/PDM.
院（系）公章: 撰写人:
一．填空题（共30分，每小题2分） 1.均质不可压缩流体的定义为。 2. 在常压下,液体的动力粘度随温度的升高而 。 3.在渐变流过流断面上，动压强分布规律的表达式为。 5.只要比较总流中两个渐变流断面上单位重量流体的大小，就能判别出流动方向。 6.产生紊流附加切应力的原因是。 7.在静止流体中，表面力的方向是沿作用面的方向。 8.圆管紊流粗糙区的沿程阻力系数λ与有关。 9.渐变流流线的特征是。 10.任意空间点上的运动参数都不随时间变化的流动称为。 11.局部水头损失产生的主要原因是。 12.直径为d的半满管流的水力半径R=。 13.平面不可压缩流体的流动存在流函数的条件是流速 x u和 y u满足方程。 14.弗劳德数Fr表征惯性力与之比。 15.在相同的作用水头下,同样口径管嘴的出流量比孔口的出流量。 二．（14分）如图所示，一箱形容器，高 1.5 h m =，宽（垂直于纸面）2 b m =，箱内充满水，压力表的读数为2 20/ kN m，用一半径1 r m =的园柱封住箱的一角，求作用在园柱面上的静水总压力的大小与方向。 三．（14分）如图所示，一水平放置的管道在某混凝土建筑物中分叉。已知主管直径3 D m =，主管流量3 35/ Q m s =，分叉管直径2 d m =，两分叉管流量均为2 Q，分叉管转角0 60 θ=，1-1断面中点的压强2 294/ p kN m =，不计水头损失，求水流对支座的作用力。 题二图题三图
四.（14分）如图所示，长 50L m =、直径0.21D m =的自流管，将水自水池引至吸水井中，然后用水泵送至水塔。已知泵的吸水管直径0.2d mm =，管长6l m =， 泵的抽水量30.064/Q m s =，滤 水网的局部阻力系数12 6.0ξξ==，弯头的局部阻力系数 30.3ξ=，自流管和吸水管的沿程阻力系数0.02λ=。假定自流 管中通过的流量等于泵的抽水量。试求：（1）水池水面与吸水 井的水面高差h ；（2）水泵的安装高度2s H m =时，水泵进口断 面的真空度。 五．（14分）油在管中以1/v m s =的速度向下流动，油 的密度3920/kg m ρ=，管长3l m =，管径25d mm =，水银 压差计测得9h cm =.试求：(1)油在管中的流态；（2）油的运动 粘度ν；（3）若保持相同的平均速度反方向运动时，压差计的读数有何变化（水银密度313600/kg m ρ'=） 六．（14分）用文丘里流量计量测空气（绝热指数 1.4k =，气体常数287/R J kg K =?）的质量流量。已知进口断面的直径为1400d mm =，绝对压强为21140/p kN m =，温度为0118T C =，喉部断面2150d mm =，22116/p kN m =，假定流动为一元恒定等熵气流，求通过流量计的质量流量。 参考答案： 一．填空题（共30分，每小题2分） 1.均质不可压缩流体的定义为 ρ=c 。 2.在常压下,液体的动力粘度随温度的升高而 降低 。 3.在渐变流过流断面上，动压强分布规律的表达式为P/ρg +z=0 。 5.只要比较总流中两个渐变流断面上单位重量流体的 总机械能 大小，就能判别出流动方向。 题四图 题五图
计算流体力学课程总结 计算流体动力学(computational Fluid Dynamics,简称CFD)是通过计算机数值 计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。是用电子计算机和离散化的数值方法对流体力学问题进行数值模拟和分析的一个分支。 流体力学和其他学科一样，是通过理论分析和实验研究两种手段发展起来的。很早就已有理论流体力学和实验流体力学两大分支。理论分析是用数学方法求出问题的定量结果。但能用这种方法求出结果的问题毕竟是少数，计算流体力学正是为弥补分析方法的不足而发展起来的。计算流体力学是目前国际上一个强有力的研究领域,是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术，广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计、HAVC&R 等诸多工程领域。 计算流体力学的任务是流体力学的数值模拟。数值模拟是“在计算机上实现的一 个特定的计算，通过数值计算和图像显示履行一个虚拟的物理实验——数值实验“。 数值模拟包括以下几个部分。首先，要建立反映问题(工程问题、物理问题等)本质数 学模型。其次，数学模型建立以后需要解决的问题是寻求高效率、高准确度的计算方法。再次，在确定了计算方法和坐标系统后，编制程序和进行计算式整个工作的主体。最后，当计算工作完成后，流畅的图像显示是不可缺少的部分。 还有一个就是CFD的基本思想问题，它就是把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通 过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求 解代数方程组获得场变量的近似值。 经过四十多年的发展，CFD出现了多种数值解法。这些方法之间的主要区别在于 对控制方程的离散方式。根据离散的原理不同，CFD大体上可分为三个分支： ?有限差分法(Finite Different Method，FDM) ?有限元法(Finite EIement Method，FEM) ?有限体积法(Finite Volume Method，FVM) 有限差分法是应用最早、最经典的CFD方法，也是最成熟、最常用的方法。它将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的 导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求出差分万程组 的解，就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题 的近似数值解法。
《工程流体力学》复习题及参考答案 整理人:郭冠中内蒙古科技大学能源与环境学院热能与动力工程09级1班 使用专业:热能与动力工程 一、名词解释。 1、雷诺数２、流线 3、压力体 4、牛顿流体５、欧拉法 6、拉格朗日法 ７、湿周 8、恒定流动 9、附面层10、卡门涡街11、自由紊流射流 12、流场 13、无旋流动1４、贴附现象１5、有旋流动1６、自由射流 17、浓差或温差射流 18、音速1９、稳定流动２0、不可压缩流体21、驻点2２、 自动模型区 二、是非题。 1.流体静止或相对静止状态的等压面一定是水平面。（ ) 2.平面无旋流动既存在流函数又存在势函数。( ) 3.附面层分离只能发生在增压减速区。（ ) 4.等温管流摩阻随管长增加而增加，速度和压力都减少。 ( ） 5.相对静止状态的等压面一定也是水平面。 ( ） 6.平面流只存在流函数,无旋流动存在势函数。 （） 7.流体的静压是指流体的点静压。 （) 8.流线和等势线一定正交。（） 9.附面层内的流体流动是粘性有旋流动。 （） 10.亚音速绝热管流摩阻随管长增加而增加，速度增加，压力减小。( ） 11.相对静止状态的等压面可以是斜面或曲面。（ ) 12.超音速绝热管流摩阻随管长增加而增加,速度减小，压力增加。( ） 13.壁面静压力的压力中心总是低于受压壁面的形心。（ ) 14.相邻两流线的函数值之差，是此两流线间的单宽流量。 ( ） 15.附面层外的流体流动时理想无旋流动。（） 16.处于静止或相对平衡液体的水平面是等压面。（） 17.流体的粘滞性随温度变化而变化，温度升高粘滞性减少;温度降低粘滞性增大。 （） 18.流体流动时切应力与流体的粘性有关，与其他无关。（） 三、填空题。
计算流体力学教案 Teaching plan of computational fluid mechanics
计算流体力学教案 前言：本文档根据题材书写内容要求展开，具有实践指导意义，适用于组织或个人。便于学习和使用，本文档下载后内容可按需编辑修改及打印。 一、流体地基本特征 1.物质地三态 在地球上，物质存在地主要形式有：固体、液体和气体。 流体和固体地区别:从力学分析地意义上看，在于它们对外力抵抗地能力不同。 固体：既能承受压力，也能承受拉力与抵抗拉伸变形。 流体：只能承受压力，一般不能承受拉力与抵抗拉伸变形。 液体和气体地区别：气体易于压缩；而液体难于压缩； 液体有一定地体积，存在一个自由液面；气体能充满任意形状地容器，无一定地体积，不存在自由液面。 液体和气体地共同点：两者均具有易流动性，即在任何 微小切应力作用下都会发生变形或流动，故二者统称为流体。 2.流体地连续介质模型
微观：流体是由大量做无规则运动地分子组成地，分子之间存在空隙，但在标准状况下，1cm3液体中含有3.3×1022个左右地分子，相邻分子间地距离约为3.1×10-8cm。1cm3气体中含有2.7×1019个左右地分子，相邻分子间地距离约为3.2×10-7cm。 宏观：考虑宏观特性，在流动空间和时间上所采用地一切特征尺度和特征时间都比分子距离和分子碰撞时间大得多。 （1）概念 连续介质（continuum/continuous medium）：质点连续充满所占空间地流体或固体。 连续介质模型（continuum continuous medium model）：把流体视为没有间隙地充满它所占据地整个空间地一种连续介质，且其所有地物理量都是空间坐标和时间地连续函数地一种假设模型：u =u（t,x,y,z）。 （2）优点 排除了分子运动地复杂性。物理量作为时空连续函数，则可以利用连续函数这一数学工具来研究问题。 3.流体地分类
课后习题 第一章 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基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
工程流体力学考试试卷 一． 解答下列概念或问题 （15分） 1． 恒定流动 2． 水力粗糙管 3． 压强的表示方法 4． 两流动力学相似条件 5． 减弱水击强度的措施 二． 填空 （10分） 1．流体粘度的表示方法有（ ）粘度、（ ）粘度和（ ）粘度。 2．断面平均流速表达式V =（ ）；时均流速表达式υ=（ ）。 3．一两维流动y 方向的速度为),,(y x t f y =υ，在欧拉法中y 方向的加速度为y a =（ ）。 4．动量修正因数（系数）的定义式0α=（ ）。 5．雷诺数e R =（ ），其物理意义为（ ）。 三． 试推求直角坐标系下流体的连续性微分方程。 （15分） 四． 已知平面不可压缩流体流动的流速为y x x x 422-+=υ， y xy y 22--=υ （20分） 1． 检查流动是否连续； 2． 检查流动是否有旋；
3．求流场驻点位置； 4．求流函数。 五．水射流以20s m/的速度从直径mm d100 =的喷口射出，冲击一对称叶片，叶片角度 θ，求：（20分） 45 = 1．当叶片不动时射流对叶片的冲击力； 2．当叶片以12s m/的速度后退而喷口固定不动时，射流对叶片的冲击力。 第（五）题图
六． 求如图所示管路系统中的输水流量V q ，已知H =24， m l l l l 1004321====， mm d d d 100421===， mm d 2003=， 025.0421===λλλ，02.03=λ，30=阀ξ。（20分） 第（六）题图 参考答案 一．1．流动参数不随时间变化的流动； 2．粘性底层小于壁面的绝对粗糙度（?<δ）； 3．绝对压强、计示压强（相对压强、表压强）、真空度； 4．几何相似、运动相似、动力相似； 5．a)在水击发生处安放蓄能器；b)原管中速度0V 设计的尽量小些；c)缓慢关闭；d)采用弹性管。 二．1．动力粘度，运动粘度，相对粘度； 第2 页 共2 页
第二章 2-1.已知某种物质的密度ρ=cm3，试求它的相对密度d。 解：d=ρ/ρw=（g/cm3）/1（g/cm3）= 2-2.已知某厂1号炉水平烟道中烟气组分的百分数为α（co2）=%，a(SO2)=%，a(O2)=%，a(N2)=76%，a(H2O)=5%。试求烟气的密度。 2-3.上题中烟气的实测温度t=170℃，实测静计压强Pe=1432Pa，当地大气压强 Pa=10058Pa。试求工作状态下烟气的密度和运动粘度。
2-4.当压强增量为50000Pa时，某种液体的密度增长%，试求该液体的体积模量。 2-5.绝对压强为×10^5Pa的空气的等温体积模量和等熵体积模量各等于多少 2-6. 充满石油的油槽内的压强为×10^5Pa,今由槽中排出石油40kg，使槽内压强降到×10^4Pa，设石油的体积模量K＝×10^9 Pa。试求油槽的体积。 2-7. 流量为50m3/h，温度为70℃的水流入热水锅炉，经加热后水温升到90℃，而水的体胀系数αV=℃，问从锅炉中每小时流出多少立方米的水
2-8. 压缩机压缩空气，绝对压强从×104Pa升高到×105Pa，温度从20℃升高到78℃，问空气体积减少了多少 2-9. 动力粘度为×10^-4Pa·S，密度为678kg/m3的油，其运动粘度等于多少 解：V=u/ρ=×10^-4/678=×10^-7m2/s 2-10. 设空气在0℃时的运动粘度ν0=×10-6m2/s,密度ρ0=m3。试求在150℃时空气的动力粘度。 2-11. 借恩氏粘度计测得石油的粘度为,如石油的密度为ρ=850kg/m3,试求石油的动力 粘度。
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分程的时候提出的一个问题，用一维无粘可压缩Euler程就可以描述了。 这里 这个程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量）随空间变化（）的关系。 在CFD常把这个程写成矢量形式 这里 进一步可以写成散度形式
《计算流体力学》课程大作业 ——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟 张伊哲 航博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.写出静止流体中的应力张量,解释其中非0项的意义. 无粘流体或静止流场中，由于不存在切向应力，即p ij =0（i ≠j ），此时有 P =00000 0xx yy zz p p p ??????????=000000p p p -????-????-??=-p 00000011????1?????? = -p I 式中I 为单位张量，p 为流体静压力。 流体力学中，常将应力张量表示为 p =-+P I T （2-9） 式中p 为静压力或平均压力，由于其作用方向与应力定义的方向相反，所以取负值；T 称为偏应力张量，即 T =xx xy xz yx yy yz zx zy zz τττττττττ?????????? （2-10） 偏应力张量的分量与应力张量各分量的关系为：i ＝j 时，p ij 为法向应力，τii = p ij - p ；当i ≠j 时p ij 为粘性剪切应力，τij =p ij 。τii =0的流体称为非弹性流体或纯粘流体，τii ≠0的流体称为粘弹性流体。 3.分析可压缩(不可压缩)流体和可压缩(不可压缩)流动的关系. 当气体速度流动较小（马赫数小于0.3)时，其密度变化不大，或者说对气流速度的变化不十分敏感，气体的压缩性没有表现出来。因此，在处理工程实际问题时，可以把低速气流看成是不可压缩流动，把气体可以看作是不可压缩流体。而当气体以较大的速度流动时，其密度要发生明显的变化，则此时气体的流动必须看成是可压缩流动。 流场任一点处的流速v 与该点（当地)气体的声速c 的比值，叫做该点处气流的马赫数，用符号Ma 表示： Ma /v c v == （4－20） 当气流速度小于当地声速时，即Ma<1时，这种气流叫做亚声速气流；当气流速度大于当地声速时，即Ma>l 时，这种气流称为超声速气流；当气流速度等于当地声速时，即Ma=l 时，这种气流称为声速气流。以后将会看到，超声速气流和亚声速气流所遵循的规律有着本质的不同。 马赫数与气流的压缩性有着直接的联系。由式(4－11)可得 所以有 222Ma d ρv dv dv ρc v v =-=-。 （4－21） 当Ma≤0.3时，dρ/ρ≤0.09dv /v 。由此可见，当速度变化一倍时，气体的密度仅仅改变9%以下，一般可以不考虑密度的变化，即认为气流是不可压缩的。反之，当Ma>0.3时，气流必须看成是可压缩的。 4.试解释为什么有时候飞机飞过我们头顶之后才能听见飞机的声音. 5.试分析绝能等熵条件下截面积变化对气流参数(v ,p ,ρ,T )的影响.
流体力学试卷及答案
1．绝对压强p abs与相对压强p 、真空度p v、当地大气压p a之间的关系是： A. p abs =p+p v； B. p=p abs－p a C. p v= p a－p abs D. p=p abs+p a 2．如图所示 A. p0=p a； B. p0>p a； C. p0
f水银；D、不一定。 5．流动有势的充分必要条件是( )。 A. 流动是无旋的； B. 必须是平面流动； C. 必须是无旋的平面流动； D. 流线是直线的流动。 6．雷诺数 Re 反映了( )的对比关系 A.粘滞力与重力 B.重力与惯性力 C. 惯性力与粘滞力 D. 粘滞力与动水压力7．一密闭容器内下部为水，上部为空气，液面下4.2m处测压管高度为2.2m，设当地大气压为1个工程大气压，则容器内气体部分的相对压强为___ 水柱（）。 A. 2m B. 1m C. 8m D. －2m 8．如图所示，下述静力学方程哪个正确？B
9．下列压强分布图中哪个是错误的？B 10．粘性流体总水头线沿程的变化是( ) 。 A. 沿程下降 B. 沿程上升 C. 保持水平 D. 前三种情况都有可能 一．名词解释（共10小题，每题2分，共20分） 1．粘滞性——流体在受到外部剪切力作用时发生变形(流 动)，其内部相应要产生对变形的抵抗，并以内摩擦力的形 式表现出来，这种流体的固有物理属性称为流体的粘滞性 或粘性 2．迹线——流体质点的运动轨迹曲线 流线——同一瞬时，流场中的一条线，线上每一点切线方向与流体在该点的速度矢量方向一致 3．层流——流体运动规则、稳定，流体层之间没有宏观的横向掺混 4．量纲和谐——只有量纲相同的物理量才能相加减，所以正确的物理关系式中各加和项的量纲必须是相同的，等式两边的量纲也必然是相同的 5．偶极流——由相距2a的点源与点汇叠加后，令a趋近于零得到的流动 6．排挤厚度——粘性作用造成边界层速度降低，相比理想流体有流量损失，相当于中心区理想流体的流通面积减少，计算时将平板表面上移一个厚度，此为排挤厚度7．顺压力梯度——沿流动方向压力逐渐降低，边界层的流动受压力推动不会产生分离8．时均速度——湍流的瞬时速度随时间变化，瞬时速度的时间平均值称为时均速度 9．输运公式——将系统尺度量转换成与控制体相关的表达式
工程流体力学_第四版_孔珑_作业答案_详 解
第二章 2-1.已知某种物质的密度ρ=2.94g/cm3，试求它的相对密度d。 解：d=ρ/ρw=2.94（g/cm3）/1（g/cm3）=2.94 2-2.已知某厂1号炉水平烟道中烟气组分的百分数为α（co2）=13.5%，a(SO2)=0.3%，a(O2)=5.2%，a(N2)=76%，a(H2O)=5%。试求烟气的密度。 2-3.上题中烟气的实测温度t=170℃，实测静计压强Pe=1432Pa，当地大气压强 Pa=10058Pa。试求工作状态下烟气的密度和运动粘度。
2-4.当压强增量为50000Pa时，某种液体的密度增长0.02%，试求该液体的体积模量。
2-5.绝对压强为3.923×10^5Pa的空气的等温体积模量和等熵体积模量各等于多少？ 2-6. 充满石油的油槽内的压强为4.9033×10^5Pa,今由槽中排出石油40kg，使槽内压强降到9.8067×10^4Pa，设石油的体积模量K＝1.32×10^9 Pa。试求油槽的体积。 2-7. 流量为50m3/h，温度为70℃的水流入热水锅炉，经加热后水温升到90℃，而水的体胀系数αV=0.000641/℃，问从锅炉中每小时流出多少立方米的水？
2-8. 压缩机压缩空气，绝对压强从9.8067×104Pa升高到5.8840×105Pa，温度从20℃升高到78℃，问空气体积减少了多少？ 2-9. 动力粘度为2.9×10^-4Pa·S，密度为678kg/m3的油，其运动粘度等于多少？ 解：V=u/ρ=2.9×10^-4/678=4.28×10^-7m2/s 2-10. 设空气在0℃时的运动粘度ν0=13.2×10-6m2/s,密度ρ0=1.29kg/m3。试求在150℃时空气的动力粘度。 2-11. 借恩氏粘度计测得石油的粘度为8.5oE,如石油的密度为ρ=850kg/m3,试求石油的动力粘度。
高等流体力学 第一章 流体力学的基本概念 连续介质：流体是由一个紧挨着一个的连续的质点所组成的，没有任何空隙的连续体，即所 谓的连续介质。 流体质点：是指微小体积内所有流体分子的总和。 欧拉法质点加速度：时变加速度与位变加速度和 z u u y u u x u u t u dt du a x z x y x x x x x ??+??+??+??== 质点的随体导数：质点携带的物理量随时间的变化率称为质点的随体导数，用dt d 表示。在欧拉法描述中的任意物理量Q 的质点随体导数表述如下： x k k Q u t Q dt dQ ??+??= 式中Q 可以是标量、矢量、张量。质点的随体导数公式对任意物理量都成立，故将质点的 随体导数的运算符号表示如下： x k k u t dt d ??+??= 其中 t ?? 称为局部随体导数，x k k u ??称为对流随体导数，即在欧拉法描述的流动中，物理 量的质点随体导数等于局部随体导数与对流随体导数之和。 体积分的随体导数：质点携带的物理量随时间的变化率称为质点的随体导数。则在由流体质点组成的流动体积V 中标量函数Φ（x, t ）随时间的变化率就是体积分的随导函数。 由两部分组成①函数Φ 对时间的偏导数沿体积V 的积分，是由标量场的非恒定性引起的。②函数Φ通过表面S 的通量。由体积V 的改变引起的。 ()dV divv dt d dV v div t dS u dV t dV dt d v v n s v v ?? ? ???Φ+Φ=??????Φ+?Φ?=Φ+?Φ?=Φ??????????????()dV adivv dt da dV av div t a dS au dV t a adV dt d v v n s v v ?? ????+=??????+??=+??=?????????????? 变形率张量： 11ε 12ε13ε D ij = 21ε 22ε 23ε 31ε 32ε 33ε
《高等流体力学》考试大纲 一、考试性质 《高等流体力学》是我校相关专业博士入学专业基础课考试科目。 二、考试形式与试卷结构 1、答卷方式：闭卷，笔试 2、答题时间；180分钟 3、题型比例 概念20% 计算与应用80% 4、参考书目 《高等流体力学》高学平，天津大学出版社，2005. 《高等工程流体力学》张鸣远等，西安交通大学出版社，2006. 三、考试要点 1、流体力学的基本概念 连续介质、欧拉法质点加速度、质点随体导数、体积分的随体导数、变形率张量、旋转角速度、判断有旋流与无旋流、涡量与速度环量的关系、应力张量的概念（包括切应力的特性、压应力的特性）、牛顿流体的本构方程（本构方程的概念、切应力和法向应力与变形的关系）。 2、流体运动的基本方程 微分形式的连续方程的表达形式、不可压缩流体的确切定义、理解其含义。N-S方程的各种表示形式、流体的能量包括哪几种形式，
并对各种形式进行解释，写出单位质量流体能量的表达式、流体运动微分形式的基本方程组有哪些方程组成，通常有几个未知量，方程组是否封闭、对于不可压缩流体，如何求解速度场、压强场以及温度场，说明其求解步骤。 3、势流运动 势流运动控制方程及求解步骤；势流求解常用的方法有哪些。速度势函数与流函数；复势与复速度；恒定平面势流的解析方法有哪几种途径；保角变换法的思路。 4、粘性流体运动 基本方程及求解途径；黏性流体运动的基本性质；黏性流体运动的解析解（如两平行板间的层流、普阿塞流的流速分布的推导）、小雷诺数流动近似解的思路；边界层的概念；边界层厚度（名义厚度、位移厚度）；边界层方程的相似性解的概念；边界层的分离现象。5、紊流运动 紊流的特征及分类；壁面剪切紊流的发生过程及紊流结构；时间平均法和系综平均法的概念。紊流运动方程—雷诺方程的推导思路，雷诺方程的形式及与N-S方程的区别，雷诺应力项的意义。紊流模型的用途，紊流模型通常有哪几类（零方程模型、一方程模型、二方程模型、其他模型）；紊流动能k、能量耗散率ε。 6、涡旋运动 涡旋的运动学性质、涡旋运动的基本方程；涡旋的形成。
高等工程流体力学 粘性流动 康顺 华北电力大学能源与动力工程系学院 Kangs@https://www.doczj.com/doc/8b13140238.html,
内容提纲 ?边界层及其方程 ?层流边界层流动转捩 ?湍流边界层结构 ?流动分离、二次流动与旋涡 能源动力领域流动问题的主要特征 ?全三维 ?非定常 ?粘性 ?高雷诺数，边界层 ?边界层：层流、转捩、湍流（紊流），分离流动，旋涡运动 叶轮机械（透平和压气机等）大多由单个或多个级组成。每个级含有一 排静子叶片列和一排转子叶片列。在级内的气流场中，一般至少有以下 几种流动现象发生：1、前缘马蹄涡；2、通道涡；3、顶部间隙涡；4、 边界层转捩；5、叶片尾迹；6、旋涡、尾迹等与叶片列周期性非定常相 互作用。 ?激波、激波与边界层相互作用
边界层流动 边界层 边界层概念：粘性很小的流体以大雷诺数运动时，在大部分流场上可以略去粘性的作用；但在物面附近的很薄的一层流体内必须考虑粘性作用。这一薄层流体称为边界层。 平板边界层示意图有边界的流动图谱 如右上图所示：流动分为三个区：边界层，尾迹区，位流区（外部势流区） 二维平板的边界层微分方程 设直匀流以零迎角平行流过一块长度为的平板,如左下图所示，人为规定，当某个y处的速度达到层外自由流的99％时，这一点到物体表面的距离（即y）称为边界层在改点的厚度，记为。显然，边界层的厚度是与X有关的，所以可以写成。 平板边界层 边界层的厚度很小，满足此关系式： 在忽略质量力的前提下，粘性平面不可压流的运动方程加上连续方程是: 用边界层条件式上式，y的数值限制在边界层之内，即 υ ∞l δδ（x） δ（x）l δ（x）<< 22 22 22 22 1 () 1 () u u u p u u u t x y x x y p u t x y y x y u x y υν ρ υυυυυ υν ρ υ ? ?????? ++=-++? ??????? ? ??????? ++=-++? ??????? ? ?? +=? ???? l δ（x）<<0yδ ≤≤
《计算流体力学》教学大纲 课程编号：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年第一版