流体力学基本概念
精品文档
**流函数：由连续性方程导出的、其值沿流线保持不变的标量函数。
**粘性：在运动状态下，流体内部质点间或流层间因相对运动而产生内摩擦力以抵抗剪切变形，这种性质叫做粘性。粘性的大小用黏度表示，是用来表征液体性质相关的阻力因子。粘度又分为动力黏度.运动黏度和条件粘度。
**内摩擦力：流体内部不同流速层之间的黏性力。
**牛顿流体：剪切变形率与切应力成线性关系的流体(水，空气)。
**非牛顿流体：黏度系数在剪切速率变化时不能保持为常数的流体(油漆，高分子溶液)。
**表面张力：1.表面张力作用于液体的自由表面上。2.气体不存在表面张力。3.表面张力是液体分子间吸引力的宏观表现。4.表面张力沿表面切向并与界线垂直。5.液体表面上单位长度所受的张力。6.用σ表示，单位为N/m。
**流线：表示某瞬时流动方向的曲线，曲线上各质点的流速矢量皆与该曲线相切。性质：a、同一时刻的不同流线，不能相交。b、流线不能是折线，而是一条光滑的曲线。c、流线簇的疏密反映了速度的大小。
**过流断面：与元流或总流的流向相垂直的横断面称为过流断面。(元流：在微小流管内所有流体质点所形成的流动称为元流。总流：若流管的壁面是流动区域的周界，将流管内所有流体质点所形成的流动称为总流。)
**流量：单位时间内通过某一过流断面的流体体积称为该过流断面的体积流量，简称流量。
**控制体：被流体所流过的，相对于某个坐标系来说，固定不变的任何体积称之为控制体。控制体的边界面，称之为控制面。控制面总是封闭表面。占据控制体的诸流体质点随着时间而改变。
**边界层：水和空气等黏度很小的流体，在大雷诺数下绕物体流动时，黏性对流动的影响仅限于紧贴物体壁面的薄层中，而在这一薄层外黏性影响很小，完全可以忽略不计，这一薄层称为边界层。
**边界层厚度：边界层内、外区域并没有明显的分界面，一般将壁面流速为零与流速达到来流速度的99％处之间的距离定义为边界层厚度。
**边界层的基本特征：(1) 与物体的特征长度相比，边界层的厚度很小。(2) 边界层内沿厚度方向，存在很大的速度梯度。(3) 边界层厚度沿流体流动方向是增加的，由于边界层内流体质点受到黏性力的作用，流动速度降低，所以要达到外部势流速度，边界层厚度必然逐渐增加。(4) 由于边界层很薄，可以近似认为边界层中各截面上的压强等于同一截面上边界层外边界上的压强值。 (5) 在边界层内，黏性力与惯性力同一数量级。 (6) 边界层内的流态，也有层流和紊流两种流态。
**滞止参数：设想某断面的流速以等熵过程减小到零，此断面的参数称为滞止参数。
**滞止参数性质：（1）在等熵流动中，滞止参数值不变；（2）在等熵流动中，速度增大，参数值降低；（3）气流中最大音速是滞止音速；（4）在有摩擦的绝热过程中，机械能转化为内能，总能量不变。
**层流：是指流体质点不互相混杂，流体质点作有条不紊的有序的直线运动。特点：1.有序性。2.水头损失与流速的一次方成正比。3.在流速较小且雷诺数Re较小时发生。4.层流遵循牛顿内摩擦定律，粘性抑制或约束质点作横向运动。
**湍流：黏性流体质点互相掺混，局部压强、速度等随时间和空间有随机脉动的流动。
**雷诺数：临界流速v与过流断面的特性几何尺寸(管径)d、流体的动力粘度μ和密度ρ有关，这四个量可以组成一个特征数(量纲一的量或无量纲数)称雷诺数 Re 。
**雷诺应力：紊流时均流动中由于流速脉动引起质点间的动量交换而产生的附加应力。
**马赫数：流场中某点的速度与该点处的声速之比。
流体与气体：两者均具有易流动性，即在任何微小切应力作用下都会发生变形或流动，故二者统称为流体。区别：气体易于压缩；而液体难于压缩。液体有一定的体积，存在自由液面；气体能充满任意形状的容器，无一定体积，不存在自由液面。
牛顿内摩擦定律：流体内摩擦力的大小与流体的性质有关，与流体的速度梯度和接触面积成正比。（切应力与剪切变形速度成正比）
实际流体：自然界中存在的具有粘性的流体。理想流体：假想的完全没有粘性的流体。利用理想流体的概念可以在研究上大简化问题，找出规律后再考虑粘性的影响进行修正，这种修正多数借助实验。
表面力：作用在隔离表面上的力，其大小和受力作用的表面面积成正比，包括垂直于作用面的压力和平行于作用面的切力。应力：单位面积上的表面力。质量力：作用在隔离体内每个流体质点上的力，其大小是和流体的质量成正比的，因为在均质流体中必然和体积相关，因此又称体积力，主要包括重力和惯性力。
连续介质：质点连续地充满所占空间的流体或固体。连续介质模型：把流体看作是全部充满、内部没有任何间隙的质点所组成的一种连续介质，且其所有的物理量都是空间坐标和时间的连续函数的一种假设模型。
恒定流：若流场中各空间点上的任何流动要素均不随时间变化，则称流动为恒定流，也称为定常流。
非恒定流：若流场中各空间点上的其中任何一个流动要素随时间变化，则称流动为非恒定流，也称为非定常流。
收集于网络，如有侵权请联系管理员删除
一、中国科学院数学与系统科学研究院简介 中国科学院数学与系统科学研究院由中科院数学研究所、应用数学研究所、系统科学研究所及计算数学与科学工程计算研究所四个研究所整合而成，此外还拥有科学与工程计算国家重点实验室、中科院管理决策与信息系统重点实验室、中科院系统控制重点实验室、中科院数学机械化重点实验室、华罗庚数学重点实验室、随机复杂结构与数据科学重点实验室，以及中科院晨兴数学中心和中科院预测科学研究中心等。2010年11月成立国家数学与交叉科学中心，旨在从国家层面搭建一个数学与其它学科交叉合作的高水平研究平台。数学与系统科学研究院拥有完整的学科布局，研究领域涵盖了数学与系统科学的主要研究方向。共有16个硕士点和13个博士点（二级学科），分布在经济学、数学、系统科学、统计学、计算机科学与技术、管理科学与工程六个一级学科中，可以在此范围内招收和培养硕士与博士研究生。在2006年全国学科评估中，我院数学学科的整体评估得分为本学科的最高分数。数学与系统科学研究院硕士招生类别为硕士研究生、硕博连读生和专业学位硕士研究生。2019年共计划招收122名。 二、中国科学院大学计算数学专业招生情况、考试科目
三、中国科学院大学计算数学专业分数线 2018年硕士研究生招生复试分数线 2017年硕士研究生招生复试分数线 四、中国科学院大学计算数学专业考研参考书目 616数学分析 现行（公开发行）综合性大学（师范大学）数学系用数学分析教程。 801高等代数 [1] 北京大学编《高等代数》，高等教育出版社，1978年3月第1版，2003年7月第3版，2003年9月第2次印刷. [2] 复旦大学蒋尔雄等编《线性代数》，人民教育出版社，1988. [3] 张禾瑞，郝鈵新，《高等代数》，高等教育出版社， 1997. 五、中国科学院大学计算数学专业复试原则 在中国科学院数学与系统科学研究院招生工作小组领导下，按研究所成立招收硕士研究生复试小组，设组长1人、秘书1人。 复试总成绩按百分制计算，其中专业知识成绩占60％，英语听力及口语测试成绩占20％，综合素质成绩占20%。 在面试环节，每位考生有5分钟自述，考查内容主要包括专业知识、外语（口语）水平
〈〈产品数据管理（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.
院（系）公章: 撰写人:
《计算流体动力学分析》学习报告 计算流体力学基础： 本章主要讲解流体动力学的核心思想以及流体动力学的控制方程。 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 方程进行实践平均处理，同时补充反应湍流特性的其他方程，例如湍动能方程以及湍流耗散率方程
计算流体力学课程总结 计算流体动力学(computational Fluid Dynamics,简称CFD)是通过计算机数值 计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。是用电子计算机和离散化的数值方法对流体力学问题进行数值模拟和分析的一个分支。 流体力学和其他学科一样，是通过理论分析和实验研究两种手段发展起来的。很早就已有理论流体力学和实验流体力学两大分支。理论分析是用数学方法求出问题的定量结果。但能用这种方法求出结果的问题毕竟是少数，计算流体力学正是为弥补分析方法的不足而发展起来的。计算流体力学是目前国际上一个强有力的研究领域,是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术，广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计、HAVC&R 等诸多工程领域。 计算流体力学的任务是流体力学的数值模拟。数值模拟是“在计算机上实现的一 个特定的计算，通过数值计算和图像显示履行一个虚拟的物理实验——数值实验“。 数值模拟包括以下几个部分。首先，要建立反映问题(工程问题、物理问题等)本质数 学模型。其次，数学模型建立以后需要解决的问题是寻求高效率、高准确度的计算方法。再次，在确定了计算方法和坐标系统后，编制程序和进行计算式整个工作的主体。最后，当计算工作完成后，流畅的图像显示是不可缺少的部分。 还有一个就是CFD的基本思想问题，它就是把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通 过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求 解代数方程组获得场变量的近似值。 经过四十多年的发展，CFD出现了多种数值解法。这些方法之间的主要区别在于 对控制方程的离散方式。根据离散的原理不同，CFD大体上可分为三个分支： ?有限差分法(Finite Different Method，FDM) ?有限元法(Finite EIement Method，FEM) ?有限体积法(Finite Volume Method，FVM) 有限差分法是应用最早、最经典的CFD方法，也是最成熟、最常用的方法。它将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的 导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求出差分万程组 的解，就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题 的近似数值解法。
【2012年】《液压与气压传动》继海宋锦春高常识-第1-7章课后答案【最新经典版】 1.1 液体传动有哪两种形式？它们的主要区别是什么？ 答：用液体作为工作介质来进行能量传递的传动方式被称之为液体传动。按照其工作 原理的不同，液体传动又可分为液压传动和液力传动，其中液压传动是利用在密封容器 液体的压力能来传递动力的；而液力传动则的利用液体的动能来传递动力的。 1.2 液压传动系统由哪几部分组成？各组成部分的作用是什么？ 答：（1）动力装置：动力装置是指能将原动机的机械能转换成为液压能的装置，它是 液压系统的动力源。 （2）控制调节装置：其作用是用来控制和调节工作介质的流动方向、压力和流量，以 保证执行元件和工作机构的工作要求。 （3）执行装置：是将液压能转换为机械能的装置，其作用是在工作介质的推动下输出 力和速度（或转矩和转速），输出一定的功率以驱动工作机构做功。 （4）辅助装置：除以上装置外的其它元器件都被称为辅助装置，如油箱、过滤器、蓄 能器、冷却器、管件、管接头以及各种信号转换器等。它们是一些对完成主运动起辅助作
用的元件，在系统中是必不可少的，对保证系统正常工作有着重要的作用。（5）工作介质：工作介质指传动液体，在液压系统常使用液压油液作为工作介质。 1.3 液压传动的主要优缺点是什么？ 答：优点：（1）与电动机相比，在同等体积下，液压装置能产生出更大的动力，也就 是说，在同等功率下，液压装置的体积小、重量轻、结构紧凑，即：它具有大的功率密度 或力密度，力密度在这里指工作压力。 （2）液压传动容易做到对速度的无级调节，而且调速围大，并且对速度的调节还可 以在工作过程中进行。 （3）液压传动工作平稳，换向冲击小，便于实现频繁换向。 （4）液压传动易于实现过载保护，能实现自润滑，使用寿命长。 （5）液压传动易于实现自动化，可以很方便地对液体的流动方向、压力和流量进行调 节和控制，并能很容易地和电气、电子控制或气压传动控制结合起来，实现复杂的运动和 操作。 （6）液压元件易于实现系列化、标准化和通用化，便于设计、制造和推广使用。答：缺点：（1）由于液压传动中的泄漏和液体的可压缩性使这种传动无法保证严格
课后习题 第一章 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基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
自然环境和工程装置中的流动常常是湍流流动，模拟任何实际过程首先遇到的就是湍流问题，而湍流问题本身又是流体力学理论上的难题。 对湍流最根本的模拟方法是在湍流尺度的网格尺寸内求解瞬态的三维N-S 方程的全模拟方法，此时无需引进任何模型。然而由于计算方法及计算机运算水平的限制，该种方法不易实现。另一种要求稍低的方法是亚网格尺寸度模拟即大涡模拟（LES ），也是由N-S 方程出发，其网格尺寸比湍流尺度大，可以模拟湍流发展过程的一些细节，但由于计算量仍然很大，只能模拟一些简单的情况，直接应用于实际的工程问题也存在很多问值题[1]。目前数模拟主要有三种方法：1.平均N-S 方程的求解，2.大涡模拟（LES ），3.直接数值模拟（DNS ），而模拟的前提是建立合适的湍流模型。 2、基本湍流模型 常用的湍流模型有： 零方程模型：C-S 模型，由Cebeci-Smith 给出；B-L 模型，由Baldwin-Lomax 给出。一方程模型：来源由两种，一种从经验和量纲分析出发，针对简单流动逐步发展起来，如Spalart-Allmaras(S-A)模型；另一种由二方程模型简化而来，如Baldwin-Barth(B-B)模型。二方程模型：应用比较广泛的两方程模型有Jones 与Launder 提出的标准k-e 模型，以及k-omega 模型。 2.1 零方程模型 上世纪30年代发展的一系列湍流的半经验理论，如Prandtl 的混合长度理论、Taylor 的涡量输运理论、von Karman 的相似性理论等，本质上即是零方程湍流模型。零方程模型直接建立雷诺应力与平均速度之间的代数关系，由于不涉及代数关系故称为另方程模型： ''m u u v y ρρε?-=? 其中m ε称为涡粘系数，他与分子的运动粘性系数ν有相同的量级。对于一般的三维的情况，上式可写为： '' 223 i j m ij ij u v S K ρεδ-=- K 为单位质量的湍流脉动动能。为了发展上述方法，需要建立m ε与平均速度之间的关系。1925年，普朗特沿这一方向做了重要工作，提出可混合长度理论，混合长度理论认为，存在这样的长度l ,在此长度内流体质点运动是自由的（不与
第一章绪论1-1． 20℃的水 2.5m 3，当温度升至80℃时，其体积增加多少？[ 解 ] 温度变化前后质量守恒，即1V12V2 又20℃时，水的密度80℃时，水的密度1998.23kg / m3 2971.83kg / m3 V2 1V 1 2.5679m3 2 则增加的体积为V V2 V1 0.0679 m3 1-2．当空气温度从0℃增加至 20℃时，运动粘度增加15%，重度减少 10% ，问此时动力粘度增加多少（百分数）？ [ 解 ] (1 0.15) 原 (1 0.1) 原 1.035 原原 1.035 原 原 1.035 原原 0.035 原原 此时动力粘度增加了 3.5% 1-3．有一矩形断面的宽渠道，其水流速度分布为u 0.002 g( hy 0.5y2 ) /，式中、分别为水的密度和动力粘度，h 为水深。试求h 0.5m 时渠底（y=0）处的切应力。 [ 解 ] du 0.002 g (h y) / dy du 0.002 g(h y) dy 当h =0.5m，y=0时 0.002 1000 9.807(0.50) 9.807Pa 1-4．一底面积为 45× 50cm2，高为 1cm 的木块，质量为 5kg，沿涂有润滑油的斜面向下作等速运动，木块运动速度 u=1m/s，油层厚 1cm，斜坡角 22.620（见图示），求油的粘度。 u
[ 解 ] 木块重量沿斜坡分力 F 与切力 T 平衡时，等速下滑 mg sin T A du dy mg sin 5 9.8 sin 22.62 A u 0. 4 0.45 1 0.001 0.1047 Pa s 1-5．已知液体中流速沿 y 方向分布如图示三种情况，试根据牛顿内摩擦定律 du ，定性绘出切应力 dy 沿 y 方向的分布图。 y y y u u u u u u [ 解 ] y y y = 0 = 1-6．为导线表面红绝缘，将导线从充满绝缘涂料的模具中拉过。已知导线直径 0.9mm ，长度 20mm ，涂料 的粘度 =0.02Pa ． s 。若导线以速率 50m/s 拉过模具，试求所需牵拉力。 （1.O1N ） [ 解 ] A dl 3.14 0.8 10 3 20 10 3 5.024 10 5 m 2
科学与工程计算国家重点实验室 简介 中国科学院科学与工程计算国家重点实验室（简称LSEC）是在已故著名数学家、中国计算数学的奠基人和开拓者冯康院士的倡导、并亲自筹备和组织下，由原中科院计算中心从事计算数学研究的部分课题组成的。实验室筹建于1990年，1993年10月经中科院验收后正式投入运行，1994年向国内外开放，1995年9月和 2005年3月两次通过国家验收。 实验室主要开展科学与工程计算中具有重要意义的基础理论研究，解决科学与工程领域中的重大计算问题，着重研究计算方法的构造、理论分析及实现。研究内容包括：动力系统与数值方法，研究各类保结构算法的理论、算法的构造和数值试验；有限元边界元方法，针对具有应用背景的椭圆边值问题及其它相关问题，提出适合于这些问题的有限元边界元新型高性能计算方法；非线性最优化，主要研究求解非线性规划的新算法以及算法的收敛性；计算流体力学,研究非定常不可压N-S方程和可压缩流的计算方法；并行计算方法和科学计算可视化；非均匀多孔介质中渗流问题的多尺度计算方法。 实验室主任是陈志明研究员。实验室学术委员会主任是中国工程院院士崔俊芝。 实验室建设以来在动力系统几何算法，非线性优化，有限元边界元，数理方程反问题，计算流体力学，并行算法，科学计算可视化等方面取得了大量的研究成果，十分突出的是关于哈密尔顿系统的辛几何算法的研究。其成果荣获“国家自然科学一等奖”。实验室在设备研制方面也取得了显著的成绩。 实验室现有科研人员19人，中科院院士2人（石钟慈、林群）,中国工程院院士1人（崔俊芝），其中研究员16人，此外，实验室还获得多项其它重要奖项，其中石钟慈院士在 2000年获“何梁何利科学与技术进奖”,林群院士获2001年获捷克科学院“数学科学成就荣誉奖”、2004年获“何梁何利科学与技术进奖”。实验室十分重视队伍建设和人才培养工作，尤其注重青年学术骨干的培养和引进。目前通过中科院“百人计划”已引进３位年轻的学科带头人，其中实验室主任陈志明研究员被国家科技部任命为９７３计划项目“高性能科学计算研究”首席科学家，一批优秀青年学术骨干脱颖而出，他们在各自的研究领域取得了可喜的成果，并因此获得了荣誉。例如，袁亚湘研究员曾获1995年首届“冯康科学计算奖”、1996年度“中国青年科学家奖”、“国家杰出青年科学基金”、1998年度“全国十大杰出青年”称号；2005年度“北京市科学技术一等奖”；张林波研究员曾获1995年度“中科院青年科学家二等奖”、1997年度“中科院优秀青年”奖、2000年度“国家科技进步奖二等奖”；白中治研究员获得1998年度“中科院自然科学三等奖”、1999年度“中科院青年科学家二等奖”、“中科院优秀青年”称号、2005年度“国家杰出青年科学基金”；许学军研究员获2000年度“钟家庆数学奖”；陈志明研究员获2000年度“国家杰出青年科学基金”、2001年度“第四届冯康科学计算奖”、2003年度“第七届中科院杰出青年”称号、2004年度“新世纪百千万人才工程国家级人选”、2005年度“海外青年学者合作研究基金”；周爱辉研究员获2004年度“国家杰出青年科学基金”。
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分程的时候提出的一个问题，用一维无粘可压缩Euler程就可以描述了。 这里 这个程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量）随空间变化（）的关系。 在CFD常把这个程写成矢量形式 这里 进一步可以写成散度形式
重庆大学《计算流体力学与计算传热学基础》上机实验水平螺旋管内的对流换热过程 学生：刘伟文 学号：20123000 指导教师：李隆键 专业：热能与动力工程 重庆大学动力工程学院 二O一五年六月
一、前言 螺旋管在热力、化工、石油及核工业等领域得到了广泛应用，螺旋管换热器也具有结构简单、传热系数高等优点。它的传热系数比直管高，在相同空间里可得到更大的传热面积，布置更长的管道，减少了焊缝，提高了安全性。尽管螺旋管的流体阻力增大，压降增大，但是其传热效率的提高导致能量的节约要高于因阻力增大而消耗的能量。因此，螺旋管在许多行业得到普遍应用而倍受青睐。在工程应用中，由于工艺要求，往往需将流体加热至规定的温度范围，传热是其中的基本单元操作，所以有必要对螺旋管的传热与流动特性进行研究。从理论知识我们知道由于向心力的作用，流体从管中心部分由螺旋管内侧流向外侧壁面，因而造成了螺旋管内侧的低压区。在压差作用下，流体从外侧沿着圆管的上部和下部壁面流回内侧。这种流动是与管的轴向垂直的，也就是与流体的主体流动相垂直，称为二次流。流体的这种二次流与轴向主流复合成螺旋式的前进运动。这样，对于流体的传热传质，不仅可依靠流体的径向扩散，还有径向二次流的作用，相当于边界层进行了破坏，增强了流体传质。 二、GAMBIT建模
1、先建立一个半径为6的圆面。 2、将该圆面向X轴正方向移动120。 3、用圆面sweep形成螺旋柱体。（绕Y轴正方向）
4、重复以上操作，得到如图所示几何体弯管。 5、设置边界层。
并应用至每个截面：
6、设置圆面的网格，选择pave方式，interval size 选择0.6，这样边界层网格与圆面中心网格过渡较平缓。 7、依次建立体网格。 8、检查网格质量。 最差网格为0.41，满足要求。 8、输出网格。
《计算流体力学》课程大作业 ——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟 张伊哲 航博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】的经典结论进行对比，评述本
计算流体力学结课报告200Km/h列车fluent仿真计算 学部：化、环、生学部 学院：化工机械与安全学院 学号：31507095 班级：化1512班 学生姓名：孙金
引言 数值仿真就是对所建立的数值模型进行数值实验和求解的过程。而计算流体力学CFD （Computational Fluid Dynamics）就是在工程仿真实验领域中应用最广泛的一门学科。任何流体运动的规律都是以质量守恒定律、动量守恒定律和能量守恒定律为基础的。这些基本定律可由数学方程组来描述，如欧拉方程、N-S方程。采用数值计算方法，通过计算机求解这些控制流体流动的数学方程，进而研究流体的运动规律这就是CFD研究问题的方法。在实际计算流体力学方面，采用通用的CFD软件来完成工程上的一些流体力学问题，有极为广泛的应用前景。近年来，随着计算机技术以及相关技术的发展，CFD技术已经在工程领域内取得重大的进步，特别是在高速列车的外型设计方面起了很大作用。随着国家经济的发展，国家运输业也有了很大的发展，特别是列车经过几次提速后，高速列车在国家运输行业中所占比例不断提高。高速列车的特点是庞大、细长、在地面轨道上运行，其空气动力学问题非常复杂。空气在列车表面形成空气流场，空气阻力急剧增加，作用在列车的阻力大部分来自压强阻力，而一部分来自表面磨擦阻力，这就使能耗过大，同时列车可能出现较大的空气升力，导致列车产生“飘”的现象，激发列车脱轨事故的发生，因此研究高速列车气动力性能非常重要。用CFD仿真可以详细了解高速列车的空气动力特性，从而设计出阻力小、噪音低等各方面性能完善的高质量列车。本文采用CFD学科中的常用商业软件Fluent仿真一个时速200km/h的二维流线型车头的外流场，对其空气动力性能进行分析，从而得到不同车辆形状其周围流场的不同，进而分析哪种车型更适合。
第一章绪论 3 1-1. 20C的水2.5m，当温度升至80C时，其体积增加多少？ ［解］温度变化前后质量守恒，即V 2V 3 又20C时，水的密度i 998.23kg /m 3 80C 时，水的密度 2 971.83kg/m3 V2— 2.5679m3 2 3 则增加的体积为V V V i 0.0679m 1-2.当空气温度从0C增加至20C时，运动粘度增加15%，重度减少10%，问此时动力粘度增加多少(百分数)？ ［解］(1 0.15)原(1 0.1)原 1.035原原1.035原 原 1.035原原 0.035 原原 此时动力粘度增加了 3.5% 1-3?有一矩形断面的宽渠道，其水流速度分布为u 0.002 g(hy 0.5y2)/ ，式中、分别为水的密度和动力粘度，h为水深。试求h 0.5m时渠底(y=0)处的切应力。 ［解］——0.002 g(h y)/ dy 0.002 g(h y) dy 当h =0.5m , y=0 时 0.002 1000 9.807(0.5 0) 9.807Pa 1-4.一底面积为45 x 50cm2,高为1cm的木块，质量为5kg，沿涂有润滑油的斜面向下作等速运动，木块运动速度u=1m/s，油层厚1cm,斜坡角22.620(见图示)，求油的粘度。
［解］木块重量沿斜坡分力F与切力T平衡时，等速下滑 mg sin du T A dy mg sin A U 5 9.8 sin 22.62 1 0.4 0.45 - 0.001 0.1047 Pa s 1-5.已知液体中流速沿y方向分布如图示三种情况，试根据牛顿内摩擦定律 沿y方向的分布图。 3 3 5 2 ［解］ A dl 3.14 0.8 10 20 10 5.024 10 m 石，定性绘出切应力 1-6 ?为导线表面红绝缘，将导线从充满绝缘涂料的模具中拉过。已知导线直径 的粘度=0.02Pa. s。若导线以速率50m/s拉过模具，试求所需牵拉力。 0.9mm，长度20mm，涂料 (1.O1N) y
中科院海洋研究所硕士研究生入学考试 《流体力学》考试大纲 本流体力学考试大纲适用于中国科学院研究生院力学专业的硕士研究生入学考试。流体力学是现代力学的重要分支，是许多学科专业的基础理论课程，本科目的考试内容主要包括流体的物理性质、流体运动学、动力学和静力学，无粘不可压缩、可压缩流动，粘性不可压缩流动及湍流、流体波动和漩涡理论等方面。要求考生对其基本概念有较深入的了解，能够熟练地掌握基本方程的推导，并具有综合运用所学知识分析问题和解决问题的能力。 一、考试内容： （一）流体的物理性质 固液气体的宏观性质与微观结构，连续介质假设及其适用条件，流体的物理性质(粘性、可压缩性与热膨胀性、输运性质、表面张力与毛细现象) ，质量力与表面力。 （二）流体运动学 流体运动的描述(拉格朗日描述与欧拉描述及其间的联系、物质导数与随体导数、迹线、流线及脉线)，流场中的速度分解，涡量，涡量场，涡线、涡管、涡通量，涡管强度及守恒定理。 （三）流体动力学 连续性方程(雷诺输运定理)，动量方程(流体的受力、应力张量)，能量方程(热力学定律)，本构关系，状态方程，流体力学方程组及定解条件，正交曲线坐标系，量纲分析与流动相似理论，流体力学中的无量纲量及其物理意义、相似原理的应用。 （四）流体静力学 控制方程，液体静力学规律，自由面的形状，非惯性坐标系中的静止液体。 （五）无粘流动的一般理论 无粘流动的控制方程，Bernoulli方程，Bernoulli方程和动量定理的应用。 （六）无粘不可压缩流体的无旋流动 控制方程及定解条件，势函数及无旋流动的性质，平面定常无旋流动(流函数、源汇、点涡、偶极子、镜像法、保角变换)，无旋轴对称流动，非定常无旋流动。 （七）液体表面波 控制方程(小振幅水波) 及定解条件，平面单色波，水波的色散和群速度，水波的能量及其传输，速度与压力场特性，表面张力波及分层流体的重力内波，非线性水波理论。 （八）旋涡运动 涡量动力学方程和涡量的产生，涡量场(空间特性、时间特性)，典型的涡模型。 （九）粘性不可压缩流动 控制方程及定解条件，定常的平行剪切流动(Couette流动、Poiseuille流动等)，非定常的平行剪切流动(Stokes第一和第二问题、管道流动的起动问题)，圆对称的平面粘性流动(圆柱Couette流及其起动过程)，小雷诺数粘性流动。 （十）层流边界层和湍流 边界层的概念，层流边界层方程(Blasius平板边界层)，边界层的分离，湍流的发生，层流到湍流的转捩，雷诺方程和雷诺应力。 （十一）无粘可压缩流动 声速和马赫数，膨胀波、弱压缩波的形成及其特点，一维等熵流(定常和非定常)，激波(正激波和斜激波)，拉瓦尔喷管流动的特征。 二、考试要求：
计算流体力学实验报告——热传导方程求解 姓名：梁庆 学号：0808320126 指导老师：江坤 日期：2010/12/30
基于FTCS格式热传导方程求解程序设计 摘要 计算流体力学是通过数值方法求解流体力学控制方程，得到流场的定量描述，并以预测流体运动规律的学科。在CFD中，我们将流体控制方程中积分微分项，近似的表示为离散的代数形式，使得积分或微分形式的控制方程转化为离散的代数方程组；然后通过计算机求解这些代数方程，从而得到流场在空间和时间点上的数值解。 基于以上思路，我们利用FTCS格式差分，工程上常用的热传导方程，并编制计算机求解程序，解出其数值解。并通过Matlab绘制，求解结果，分别以二维，三维的形式，给出求解结果，本实验通过求解的数值解，制作了1秒内长度为1的距离内，热传导情况动画，以备分析所用。 关键词FTCS 有限差分热传导方程
一、 问题重述 编制一个可以有限差分程序，实现求解热传导方程。 非定常热传导方程： 22(0) u u t t γγ??=>?? 初边值问题的有限差分求解。初始条件和边界条件为： (,0)() (0,)()0(1,)()0 u x f x u t a t u t b t =?? ==??==? 其中1γ=，初值条件为：000.3()1 0.30.71010 0.7 1.0 3 3x f x x x ? <
计算数学研究方向 网上摘抄：计算数学研究方向及网上资料 计算数学目的为物理学和工程学作计算。主要研究方向包括： 数值泛函分析；连续计算复杂性理论；数值偏微与有限元；非线性数值代数及复动力系统； 非线性方程组的数值解法；数值逼近论；计算机模拟与信息处理等；工程问题数学建模与计算等等。 目前发展最好的方向已经与应用数学的CAGD 方向合二为一。现在最热的方向应该是微分方程的数值求解、数值代数和流形学习，数值计算名校：西安交通大学、北京大学、大连理工大学 从计算数学的字面来看，应该与计算机有密切的联系，也强调了实践对于计算数学的重要性。 也许Parlett 教授的一段话能最好地说明这个问题： How could someone as brilliant as von Neumann think hard about a subject as mundane as triangular factoriz-ation of an invertible matrix and not perceive that, with suitable pivoting, the results are impressively
good Partial answers can be suggested-lack of hands-on experience, concentration on the inverse rather than on the solution of Ax = b -but I do not find them adequate. Why did Wilkinson keep the QR algorithm as a backup to a Laguerre-based method for the unsymmetric eigenproblem for at least two years after the appearance of QR Why did more than 20 years pass before the properties of the Lanczos algorithm were understood I believe that the explanation must involve the impediments to comprehension of the effects of finite-precision arithmetic. ( 引自既然是计算数学专业的学生，就不能对自己领域内的专家不有所了解。早些年华人在计算数学领域里面占有一席之地是因为冯康院士独立于西方，创立了有限元方法，而后又提出辛算法。这里只是列出几位比较年轻的华人计算数学专家，因为他们代表了当前计算数学的研究热点，也反映华人对计算数学的发展的贡献。 侯一钊（加州理工） 研究方向：计算流体力学、多尺度计算与模拟、多相流 鄂维南（Princeton 大学） 北京大学长江学者，研究方向：多尺度计算与模拟 包刚（Michigan 州立大学） 吉林大学长江学者，研究方向：光学与电磁场中的计算等 金石(Wisconsin 大学)
第一章习题答案 选择题（单选题） 1.1 按连续介质的概念，流体质点是指：（d ） （a ）流体的分子；（b ）流体内的固体颗粒；（c ）几何的点；（d ）几何尺寸同流动空间相比是极小量，又含有大量分子的微元体。 1.2 作用于流体的质量力包括：（c ） （a ）压力；（b ）摩擦阻力；（c ）重力；（d ）表面张力。 1.3 单位质量力的国际单位是：（d ） （a ）N ；（b ）Pa ；（c ）kg N /；（d ）2/s m 。 1.4 与牛顿内摩擦定律直接有关的因素是：（b ） （a ）剪应力和压强；（b ）剪应力和剪应变率；（c ）剪应力和剪应变；（d ）剪应力和流速。 1.5 水的动力黏度μ随温度的升高：（b ） （a ）增大；（b ）减小；（c ）不变；（d ）不定。 1.6 流体运动黏度ν的国际单位是：（a ） （a ）2/s m ；（b ）2/m N ；（c ）m kg /；（d ）2/m s N ?。 1.7 无黏性流体的特征是：（c ） （a ）黏度是常数；（b ）不可压缩；（c ）无黏性；（d ）符合 RT p =ρ 。 1.8 当水的压强增加1个大气压时，水的密度增大约为：（a ） （a ）1/20000；（b ）1/10000；（c ）1/4000；（d ）1/2000。 1.9 水的密度为10003 kg/m ，2L 水的质量和重量是多少？ 解：10000.0022m V ρ==?=（kg ） 29.80719.614G mg ==?=（N ） 答：2L 水的质量是2kg ，重量是19.614N 。 1.10 体积为0.53m 的油料，重量为4410N ，试求该油料的密度是多少？ 解：44109.807 899.3580.5 m G g V V ρ= ===（kg/m 3） 答：该油料的密度是899.358kg/m 3 。 1.11 某液体的动力黏度为0.005Pa s ?，其密度为8503 /kg m ，试求其运动黏度。
课程综合作业课程名称：计算流体力学 专业班级：研究方向： 学生姓名：学号： 完成日期：
计算流体力学课程综合报告 1.简介 计算流体动力学（Computational Fluid Dynamics，简称CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。其基本思想为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看作是在流动基本方程（质量守恒方程、动量守恒方程、能量守恒方程）控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量（速度、压力、温度、浓度等）的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理星，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 2.计算流体动学的特点: ①流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和边界条件复杂，很难求得解析解，而用CFD方法则有可能找出满足工程需要的数值解。 ②可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较。 ③它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件。 ④数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差。 ⑤它不像物理模型实验一开始就能给出流动现象并定性地描述，往往需要由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模型进行验证。