计算流体力学参考书介绍
https://www.doczj.com/doc/9416020295.html, & https://www.doczj.com/doc/9416020295.html, yjs808
一、中文参考书：
1.中国人民解放军总装备部军事训练教材编辑工作委员会，计算流体力学及应用，国防工业出版社，2003年。（*****重点推荐）
该书比较系统，有限基本解法（亚超声速流的有限基本解方法），有限差分法和有限体积法均有介绍，对于应用的网格生成技术（网格生成基本方法，网格分区与重叠网格技术，非结构网格生成技术），非线性速势方程解法，不可压NS方程计算（压力泊松方程方法，压力修正方法（SIMPLE方法），虚拟压缩性方法），可压缩的Euler和NS方程计算（MacCormack 显式差分法，Beam-Warming因式分解格式，Jameson有限体积格式，TVD格式，NND格式，ENO及ENN格式，几种常用的隐式（时间方向离散）算法（近似因式分解方法，对角化算法，LU-ADI算法，LU-SGS算法），湍流模型）以及稀薄气体的蒙特卡罗的数值模拟都有专题，知识介绍结构清晰，有助于了解计算流体力学的整体知识框架。
2.苏铭德，黄素逸，计算流体力学基础，清华大学出版社，1997年。（*****重点推荐）
该书虽说取名为基础，但如果能全部看明白，倒也很不容易。该书分为若干个专题：数值模拟专题，数值计算方法，流场的数值计算。里面从发展型方程的有限差分法讲起，讲到偏微分方程的几种常用数值方法（特征线法，有限元法，泊松方程的直接方法，对流扩散方程的有限解析法，发展方程的谱方法），再转到流场的数值计算（无粘性流体流动的数值计算（等熵流动的数值计算，Burgers方程及其求解，激波的捕捉，Riemann问题的解和Godunov格式，多维气体流动的数值计算），粘性流体流动的数值计算（可压粘性流动数值计算的MacCormack格式、Beam-Warming格式、反扩散和NND格式、通量分裂法和推进迭代法；用流函数-旋度方程、有限解析法、有限差分法、推进迭代法、谱方法等求解不可压粘性流动方程）上应用；该书总体来讲还算是浅显易懂；通过看这本书我们能学习到很多东西，比如各种数值方法（有限差分法，有限解析法，有限元法，谱方法，边界元法）的区别和联系小知识等等。
3.张涵信，沈孟育，计算流体力学-差分方法的原理和应用，国防工业出版社，2003年。（*****重点推荐）
该书首先从流体力学的各级近似方程开始，从分析方程数学性质，到定解条件的提法；再用模型方程及其差分计算的分析理论，引出一系列的经典的差分格式（MacCormack格式，基于Runge-Kutta法的显式耗散格式，Beam-Warming隐式耗散格式，Jameson隐式耗散格式，混合反扩散格式）以及无波动的不含自由参数的差分格式；由低阶格式的模拟能力限制提出了建立高阶格式的必要、建立高阶精度格式的原则，介绍了一些典型的高精度差分格式（三阶ENN格式，ENO格式，广义紧致格式，时空守恒方法）；由低维流动问题的差分格式到多维流动问题的差分格式（半离散化的多维差分格式，Runge-Kutta方法，交替方向隐式(ADI)方法，时间分裂格式（分数步格式），隐式近似因式分解方法）；最后讲述了边界条件的数学处理，网格生成技术以及一些实用的数值计算规律。这本书有点像工具书，内容全面而实用。
4.刘儒勋，舒其望，计算流体力学的若干新方法，科学出版社，2003年。（****一般重点推荐）
该书深入地介绍了当代偏微分方程数值方法的新发展，特别是其中最具有理论意义和实用价值的间断解或弱解的数值方法，以及我们面临的新课题和新挑战。对于研究比较前沿课题和需要在CFD方面进阶的同学，这本书算是一本很不错的参考书。该书主要内容包括：预备篇（有关数学模型的简单介绍（线性试验模型，非线性发展方程试验模型，流体动力学的某些数值模拟实验模型），关于有限差分方法和其他数值方法的某些重要概念（差分格式的相容性、收敛性和稳定性，有限差分法和其他数值方法的守恒性、单调性问题，差分格式和其他数值方法的迎风性设计，间断解或者弱解的一些重要概念，Riemann问题的数值方法），有限差分方法的数值耗散、数值频散和数值群速度效应（有限差分格式的修正方程和格式余项效应分析方法，数值模拟中所发生的奇怪现象的解释，差分格式的改造、改进和优化）；间断解问题篇（间断分解和Riemann问题的特征表示，Riemann间断分解问题的Roe解法（基本思想，Roe参向量和近似Riemann解算子，Roe方法提高求解精度的关键和存在的问题），积分平均型间断解问题方法的设计和讨论（积分平均型守恒、单调格式的一般构思，Van Leer 的MUSCL方法的构造，Collela和Woodward的PPM，积分平均格式构造的进一步讨论）；高分辨率方法篇（TVD方法的构造和控制函数（从Lax-Wendroff格式谈起，TVD格式的构造与讨论，限制函数或限制器的简单讨论），ENO和WENO格式（一维、二维标量方程的ENO和WENO 格式，时间离散的Runge-Kutta方法，方程组问题）；有限体积法篇（非结构网格的生成与构造，非结构网格有限体积法）；非标准有限方法篇（混合有限元简介，运动有限元方法，间断有限元方法，时空有限元方法）；运动界面追踪问题的数值方法篇（VOF方法和运动界面的重构方法，等值面（Level Set）函数方法）。
5.李万平，计算流体力学基础，华中科技大学出版社，2004年。（****一般重点推荐）
短小精悍是该书的一大特点，适合于初学者结合其他书籍作辅助性阅读。该书偏重于流体力学的有限体积法介绍，通过阅读An Introduction to Computational Fluid Dynamics: The Finite Volume Method (1995 H. K. VERSTEEG）不难发现，这本中文书很大部分是对这本英文书的翻译，书中的其他部分则是《计算流体力学基础》（苏铭德，黄素逸）中有关发展型方程差分方法的章节。该书内容主要包括：流体流动的数学模型，湍流与湍流模型，发展方程的有限差分法，有限体积法（扩散问题的有限体积法，对流扩散问题的有限体积法，定常流动的压强速度耦合算法，非定常流动的有限体积法），网格生成技术以及数值模拟二维不可压缩粘性流体绕圆柱的流动。
6.吴江航，韩庆书，计算流体力学的理论，方法及应用，科学出版社，1988年。（****一般重点推荐）
该书比较老，不过里面的有些内容还是比较好的，能做很好的辅助性阅读，比如，离散近似的伪物理效应及其修正等等，数值方法主要讲述有限差分法，有限元法，有限分析法。该书主要内容包括：流体力学基本方程组，流体力学的模型方程及其数学物理性质，有限差分近似及其数学性质，数值解算的数学理论基础，离散近似的伪物理效应（“逆风”效应与“迎风”修正，数值耗散与数值频散，能量增长与反常能量谱转移效应，伪湍流现象），流
体力学模型方程的几种差分格式（过程的稳定性和定解条件的恰当性（小扰动方程的由来），对流方程的几种差分格式，扩散方程和对流扩散方程的差分格式，KdV方程的差分格式，双曲型方程组特征型和守恒型的差分格式），非线性方程的差分格式（无粘性方程的差分格式，粘性方程的差分格式，隐式格式，线性化方法，精解差分格式），求解多维初值问题的分步方法和交替方向的隐式格式，计算气体与浅水波运动物理解的差分格式（单调差分格式，一维、二维Godunov格式，活动网格的Godunov格式），不可压粘性流的差分解（网格步长与差分格式的选取，定常Navier-Stokes方程的差分格式，边界格式，压力的泊松方程格式，人工压缩性方法的交替方向隐式（ADI）格式，不定常的Navier-Stokes方程的差分格式，MAC方法，罚函数方法），流体力学中的有限元方法，二维不可压粘性流与浅水流动的有限元解，对流扩散问题与不可压粘性流的有限分析方法，对流扩散问题与不可压粘性流的分布杂交方法，流体力学中的Green函数方法。
7.傅德薰，马延文，计算流体力学，高等教育出版社，2002年。（***一般推荐）
该书给我最大的印象就是严重偏向于可压缩流体的计算，虽说目录章节框架看起来比较系统，但内容读起来有点乱有点儿晦涩难懂，不易阅读。内容主要包括：流体力学方程及模型方程，偏微分方程的数值解法（有限差分法，偏微分方程的全离散（非定常流动时间方向的离散方法（时间导数简单离散法，Taylor级数展开法，Runge-Kutta方法，多步法）），有限体积法，有限元方法，谱方法），高精度有限差分法及数值解的行为分析（模型方程及半离散化方程，高精度差分逼近式，数值解的精度及分辨率分析，数值解的耗散效应与色散效应，数值解的群速度，时间离散的色散与耗散效应），代数方程的求解（Gauss消去法，追赶法，交替方向追赶法，非线性方程的求解，时间关系法及局部时间步长法，多重网格技术），可压缩流体力学方程组的离散（一维流体力学方程及Jacobian系数矩阵的分裂，一维Euler 方程的离散，Godunov间断分解法，Roe格式与Roe分解，多维问题的差分逼近，粘性项的差分逼近），激波高分辨率差分格式（数值解中的非物理振荡，一阶、二阶TVD格式，MUSCL 格式），不可压Navier-Stokes方程的差分逼近，网格技术。
8.吴子牛，计算流体力学基本原理，科学出版社，2001年。（**一般推荐）
拜读作者的前言及后语，终于明白作者的原意是尽量避开较为烦琐的各种计算格式介绍，重点介绍基本原理；可仔细阅读书的内容之后，却给人一种头重脚轻虎头蛇尾的感觉。流体力学中的各种形式的方程倒是介绍了一大堆，这对于写论文需要这个的同学倒是一个不错的参考；等讲到构造计算格式的基本原理的时候，书也快结束了，你的感觉才刚刚开始，却发现书已经结束了，但里面的基本原理也许真的像作者的原意那样具有重要的启发性意义吧。
……
二、英文参考书
1.John D. Anderson，Computational Fluid Dynamics: The Basics with Applications，McGraw-Hill Press, 1995。（*****重点推荐）
该书在国内有清华大学的影印版。很多大学在讲授CFD课程时
都喜欢推荐这本书给学生作开始读物，该书在基础方面的确很不
错，前面讲述了相关的流体力学基础知识，应该比国内很多中文书
都讲的好，概念清楚易懂，后面介绍了计算流体力学的一些基础及
基本应用。该书中有很多知识点作者都是用框架图画出的，使读者
对该领域有个总体框架的印象。该书对于已经有一定基础的人，参
考价值不是很大。
2.John C. Tannehill, Dale A. Anderson, Richard H. Pletcher, Computational Fluid Mechanics and Heat Transfer, Taylor & Franics Publisher, 1997。（*****重点推荐）
该书比较系统，从偏微分方程入手，讲到基本的离散方法，再讲到模型方程的数值方法；然后列出流体力学和传热学中控制方程，紧接着重点按照几种物理近似讲述各种近似方程的数值解法：无粘流的数值方法，边界层类型方程的数值方法，抛物化NS方程的数值方法以及NS方程的数值方法，最后介绍了网格生成方面的知识。层次清晰，根据不同的流体方程讲述不同的数值方法。
3.Joel H. Ferziger, Milovan Peric, Computational Methods for Fluid Dyanmics, Springer-Verlag, 2002。（*****重点推荐）
该书层次相当地清晰，涵盖了CFD流程中的大部分内容，从数值方法（主要为有限差分法和有限体积法）入手，讲述不同的物理方程的数值解法，较为经典。该书附有网址，提供了一些可以下载的代码使用。
4.J.Blazek, Computational Fluid Dynamics: Pinciples and Applications, Elsevier, 2001。（*****重点推荐）
该书首先讲述流体力学控制方程，然后讲述求解控制方程的原理，接着讲述空间离散（结构网格和非结构网格的有限体积方法），时间离散，湍流模型，边界条件，加速技术，计算格式相容性精度以及稳定性要求与分析方法，最后讲述了网格生成原理。随书光盘提供了一些网格生成，以及一维二维Euler方程求解器的代码。
5.T. J. Chung, Computational Fluid Dyanmics, Cambridge University Press, 2002。（*****重点推荐）
该书也是重点讲述不同的数值方法在不同物理方程数值计算
中的应用，重点介绍有限差分法，有限体积法，有限元法。内容
很全，除了传统的计算流体力学，还包括多相流计算，化学反应
流与燃烧计算，计算声学，混合模型的辐射传热学计算，磁流体
力学计算，相对大气流动计算等。有点像工具书。
6.C. A. J. Fletcher, Computational Techniques for Fluid Dynamics 1 Fundamental and General Techniques & 2 Specific Techniques for Different Flow Categories, Springer-Verlag, 1987。（*****重点推荐）
该书分为上下两个分册，第一册重点讲述偏微分方程数值方法基础，讲述典型的模型方程的数值计算方法，稳态问题的数值方法，发展问题的数值方法；第二册重点讲述那些数值方法在不同流动下的应用，无粘流动，边界层流动，物理近似的NS方程计算，不可压粘性流动，可压缩粘性流动等等。该书知识框架明确，不像很多中文书说不清楚CFD到底包括那些过程，书中包括大量的代码是该书的另一大特色。
7.H. K. Versteeg, An Introduction to Computational Fluid Dynamics: The Finite Volume Method, Longman Scientific & Technical,1995。（****一般重点推荐）
该书比较简短，主要介绍有限体积法在计算流体力学中的应
用，从模型方程开始，到流动方程的数值计算，一步步地深入，
限于它的篇幅，内容不太全面。
8.O. C. Zienkiewicz, The Finite Element Method Volume 3: Fluid Dynamics, Butterworth Heinemann, 2000。（****一般重点推荐）
该书是有限元法牛人O.C.Zienkiewicz的经典之作The
Finite Element Method三卷本中的第三本，数值方法显然会
主要讲解有限元法在计算流体力学中的应用。没有仔细读过，
不太好评论，估计大牛的书应该错不了。
9.K A. Hoffmann, Computational Fluid Dynamics Volume 1, 2, 3, A Publication of Engineering Education System, 1989, 1993, 1998, 2000。（****一般重点推荐）
这套书是工具书，阵容庞大，有点儿偏向于偏微分方程的数值解法，当然包括了CFD 中的绝大部分内容，第三分册主要讲述湍流。
10.Dimitris Drikakis, William Rider, High-Resolution Methods for Incompressible and Low-Speed Flows, Springer, 2005。（****一般重点推荐）
该书主要讲述不可压低速流动的高分辨率计算方法，书中的内容比较新，对于做这方面研究的人是很好的参考。
……
三、CFD相关的书籍
1.Joe F Thompson, Handbook of Grid Generation, CRC Press, 1998。（*****重点推荐）
网格生成手册，显然有工具书的味道，书的内容相当庞大，有这一本在手，网格生成方面的大部分问题应该能在上面找到。
2.M. Farraskhalvat, Basic Structured Grid Generation with An Introduction to Unstructured Grid Generation, Butterworth Heinemann, 2003。（*****重点推荐）
该书写的相当精致工整，看起来很舒服，里面附有一些代码。对于不太复杂的网格生成该书是比较好的参考书。
3. Alfio Quarteroni, Numerical Mathematics, Springer, 2000。（*****重点推荐）
该书主要讲述数值计算方法，对于CFD中的代数方程组求解方法有一定的全面的介绍。
4.Richard Barrett etc. Templates for the Solution of Linear Systems Building Blocks for Iterative Methods, CFD Online。（*****重点推荐）
该文章简短地讲述几乎所有可能用到的线性方程组的数值解法。
... ...
D D y S x e P gh2 gh1 h2 h1 b L y C C D D y x P hc 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ， 3 m kg 29.1=空气ρ 牛顿内摩擦定律：剪切力： dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度： ρυ μ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱mm 73610/9800012 ===m m N at 2/101325 1m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D +=α sin 1 ) () 2(32121h h h h L e ++= ρ 若01 =h ，则压强为三角形分布，3 2L e y D == ρ 注：①图算法适合于矩形平面；②计算静水压力首先绘制压强分布图， α 且用相对压强绘制。 （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形12 3 bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力z x P P P += 与水平面的夹角 x z P P arct an =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ? ? ??? ??+??+??+??=??+??+??+??=??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ
流体力学的发展和现状 作为物理的一部分，流体力学在很早以前就得到发展。在19世纪，流体力学沿着两个方面发展，一方面，将流体视为无粘性的，有一大批有名的力学数学家从事理论研究，对数学物理方法和复变函数的发展，起了相当重要的作用; 另一方面，由于灌溉、给排水、造船，及各种工业中管道流体输运的需要，使得工程流体力学，特别是水力学得到高度发展。将二者统一起来的关键是本世纪初边界层理论的提出，其中心思想是在大部分区域，因流体粘性起的作用很小，流体确实可以看成是无粘的。这样，很多理想流体力学理论就有了应用的地方。但在邻近物体表面附近的一薄层中，粘性起着重要的作用而不能忽略。边界层理论则提供了一个将这两个区域结合起来的理论框架。边界层这样一个现在看来是显而易见的现象，是德国的普朗特在水槽中直接观察到的。这虽也是很多人可以观察到的，却未引起重视，普朗特的重大贡献就在于他提出了处理这种把两个物理机制不同的区域结合起来的理论方法。这一理论提出后，在经过约10年的时间，奠定了近代流体力学的基础。 流体力学又是很多工业的基础。最突出的例子是航空航天工业。可以毫不夸大地说，没有流体力学的发展，就没有今天的航空航天技术。当然，航空航天工业的需要，也是流体力学，特别是空气动力学发展的最重要的推动力。就以亚音速的民航机为例，如果坐在一架波音747飞机上，想一下这种有400多人坐在其中，总重量超过300吨，总的长宽有大半个足球场大的飞机，竟是由比鸿毛还轻的空气支托着，这是任何人都不能不惊叹流体力学的成就。更不用说今后会将出现更大、飞行速度更快的飞机。 同样，也不可能想象，没有流体力学的发展，能设计制造排水量超过50万吨的船舶，能建造长江三峡水利工程这种超大规模工程，能设计90万kW汽轮机组，能建造每台价值超过10亿美元的海上采油平台，能进行气候的中长期预报，等等。甚至天文上观测到的一些宇宙现象，如星系螺旋结构形成的机理，也通过流体力学中形成的理论得到了解释。近年来从流体力学的角度对鱼类游动原理的研究，发现了采用只是摆动尾部（指身体大部不动）来产生推进力的鱼类，最好的尾型应该是细长的月牙型。这正是经过几亿年进化而形成的鲨鱼和鲸鱼的尾型，而这些鱼类的游动能力在鱼类中是最好的。这就为生物学进化方面提供了说明，引起了生物学家的很大兴趣。 所以很明显，流体力学研究，既对整个科学的发展起了重要的作用，又对很多与国计民生有关的工业和工程，起着不可缺少的作用。它既有基础学科的性质，又有很强的应用性，是工程科学或技术科学的重要组成部分。今后流体力学的发展仍应二者并重。 本世纪的流体力学取得多方面的重大进展，特别是在本世纪下半叶，由于实验测试技术、数值计算手段和分析方法上的进步，在多种非线性流动以及力学和其他物理、化学效应相耦合的流动等方面呈现了丰富多采的发展态势。 在实验方面，已经建立了适合于研究不同马赫数、雷诺数范围典型流动的风洞、激波管、弹道靶以及水槽、水洞、转盘等实验设备，发展了热线技术、激光技术、超声技术和速度、温度、浓度及涡度的测量技术，流动显示和数字化技术的迅猛发展使得大量数据采集、处理和分析成为可能，为提供新现象和验证新理论创造了条件。 流体力学是在人类同自然界作斗争，在长期的生产实践中，逐步发展起来的。早在几千年前，劳动人民为了生存，修水利，除水害，在治河防洪，农田灌溉，河道航运，水能利用等方面总结了丰富的经验。我国秦代李冰父子根据“深淘滩，低作堰”的工程经验，修建设计的四川都江堰工程具有相当高的科学水平，反映出当时人们对明渠流和堰流的认识已经达
〈〈产品数据管理（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 液体和气体统称为流体，流体的基本特性是具有流动性，只要剪应力存在流动就持续进行，流体在静止时不能承受剪应力。 2 流体连续介质假设：把流体当做是由密集质点构成的，内部无空隙的连续体来研究。 3 流体力学的研究方法：理论、数值、实验。 4 作用于流体上面的力 （1）表面力：通过直接接触，作用于所取流体表面的力。 作用于A 上的平均压应力 作用于A 上的平均剪应力 应力 法向应力 切向应力 （2）质量力：作用在所取流体体积内每个质点上的力，力的大小与流体的质量成比例。（常见的质量力：重力、惯性力、非惯性力、离心力） ΔF ΔP ΔT A ΔA V τ 法向应力 周围流体作用 的表面力 切向应力 A P p ??=A T ??=τA 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 =2m s
单位为 5 流体的主要物理性质 （1） 惯性：物体保持原有运动状态的性质。质量越大，惯性越大，运动状态越难改变。 常见的密度（在一个标准大气压下）： 4℃时的水 20℃时的空气 （2） 粘性 牛顿内摩擦定律： 流体运动时，相邻流层间所产生的切应力与剪切变形的速率成正比。即 以应力表示 τ—粘性切应力，是单位面积上的内摩擦力。由图可知 —— 速度梯度，剪切应变率（剪切变形速度） 粘度 μ是比例系数，称为动力黏度，单位“pa ·s ”。动力黏度是流体黏性大小的度量，μ值越大，流体越粘，流动性越差。 运动粘度 单位：m2/s 同加速度的单位 说明： 1）气体的粘度不受压强影响，液体的粘度受压强影响也很小。 2）液体 T ↑ μ↓ 气体 T ↑ μ↑ 3 /1000m kg =ρ3 /2 .1m kg =ρdu T A dy μ=? h u u+du U y dy x dt dr dy du ?=?=μμτdu u dy h =ρμ ν=
流体力学 实验指导书与报告 静力学实验 雷诺实验 中国矿业大学能源与动力实验中心
学生实验守则 一、学生进入实验室必须遵守实验室规章制度，遵守课堂纪律，衣着整洁，保持安静，不得迟到早退，严禁喧哗、吸烟、吃零食和随地吐痰。如有违犯，指导教师有权停止基实验。 二、实验课前，要认真阅读教材，作好实验预习，根据不同科目要求写出预习报告，明确实验目的、要求和注意事项。 三、实验课上必须专心听讲，服从指导教师的安排和指导，遵守操作规程，认真操作，正确读数，不得草率敷衍，拼凑数据。 四、预习报告和实验报告必须独自完成，不得互相抄袭。 五、因故缺课的学生，可向指导教师申请一次补做机会，不补做的，该试验以零分计算，作为总成绩的一部分，累计三次者，该课实验以不及格论处，不能参加该门课程的考试。 六、在使用大型精密仪器设备前，必须接受技术培训，经考核合格后方可使用，使用中要严格遵守操作规程，并详细填写使用记录。 七、爱护仪器设备，不准动用与本实验无关的仪器设备。要节约水、电、试剂药品、元器件、材料等。如发生仪器、设备损坏要及时向指导教师报告，属责任事故的，应按有关文件规定赔偿。 八、注意实验安全，遵守安全规定，防止人身和仪器设备事故发生。一旦发生事故，要立即向指导教师报告，采取正确的应急措施，防止事故扩大，保护人身安全和财产安全。重大事故要同时保护好现场，迅速向有关部门报告，事故后尽快写出书面报告交上级有关部门，不得隐瞒事实真相。 九、试验完毕要做好整理工作，将试剂、药品、工具、材料及公用仪器等放回原处。洗刷器皿，清扫试验场地，切断电源、气源、水源，经指导教师检查合格后方可离开。 十、各类实验室可根据自身特点，制定出切实可行的实验守则，报经系(院)主管领导同意后执行，并送实验室管理科备案。 1984年5月制定 2014年4月再修订 中国矿业大学能源与动力实验中心
湍流模型发展综述 摘要：在概述了湍流问题的基础上，本文简要介绍了湍流的四种模型，对湍流模型在不同情况下的模拟能力进行了对比，最后简述了湍流模型的发展方向。 关键词：湍流模型；Navier-Stokes方程组；J-K模型 Abstract：On the basis of introducing the problems of turbulence, this paper briefly analyzed four kinds of turbulence models and compared their ability of simulation in different situations. At last, the paper expounded the development direction of the turbulence model. Key words：Turbulence model; Navier-Stokes equations; J-K model 一、引言 湍流又称紊流,是自然界中常见的一种很不规则的流动现象。当粘性阻尼无法消除惯性的影响时,自然界中的绝大部分流动都是湍流。 湍流运动的实验研究表明，虽然湍流结构十分复杂，但它仍然遵循连续介质的一般动力学规律，湍流流动的各物理量的瞬时值也应该服从一般的N-S方程。对粘性流体服从的N-S方程进行时均化，就可以得到雷诺平均方程。与定常的N-S方程相比，不同之处是在该式右边多了九项与脉动量有关的项，这脉动量的乘积的平均值与密度的乘积是湍流流动中的一种应力，称为湍流应力或雷诺应力。其中，法向雷诺应力和切向雷诺应力各有三个。 湍流问题就是在给定的边界条件下解雷诺方程。由于雷诺平均方程中未知数个数远多于方程个数而出现了方程不封闭的问题，这就需要依据各种半经验理论提出相应的补充方程式，即各种湍流模型。一般按照所用湍流量偏微分方程的物理含义或者数量进行区分，分别称为梅罗尔—赫林方法和雷诺方法。而后者又将湍流模型分成四类。（1）零方程模型；（2）一方程模型；（3）二方程模型；（4）应力方程模型。下面就对这些模型进行简单的描述。 二、湍流模型简介 1、零方程模型 最初的湍流模型只考虑了一阶湍流计算统计量的动力学微分方程，即平均方程，没有引进高阶统计量的微分方程，因而称之为一阶封闭模式或零方程模型。零方程模型又称为代数模型，代数模型又可以分成以下几种模型：（1）Cebeci —Smith 模型，（2）Baldwin—Lomax 模型，（3）Johnson—King 模型。 其中，B-L与C-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方法，也是最成熟、最常用的方法。它将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的 导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求出差分万程组 的解，就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题 的近似数值解法。
流体力学-笔记参考书籍： 《全美经典-流体动力学》 《流体力学》张兆顺、崔桂香 《流体力学》吴望一 《一维不定常流》 《流体力学》课件清华大学王亮主讲 目录： 第一章绪论 第二章流体静力学 第三章流体运动的数学模型 第四章量纲分析和相似性 第五章粘性流体和边界层流动 第六章不可压缩势流 第七章一维可压缩流动 第八章二维可压缩流动气体动力学 第九章不可压缩湍流流动 第十章高超声速边界层流动 第十一章磁流体动力学 第十二章非牛顿流体 第十三章波动和稳定性
第一章 绪论 1、牛顿流体： 剪应力和速度梯度之间的关系式称为牛顿关系式，遵守牛顿关系式的流体是牛顿流体。 2、理想流体：无粘流体，流体切应力为零，并且没有湍流？。此时，流体内部没有内摩擦，也就没有内耗散和损失。 层流：纯粘性流体，流体分层，流速比较小； 湍流：随着流速增加，流线摆动，称过渡流，流速再增加，出现漩涡，混合。因 为流速增加导致层流出现不稳定性。 定常流：在空间的任何点，流动中的速度分量和热力学参量都不随时间改变， 3、欧拉描述：空间点的坐标； 拉格朗日：质点的坐标； 4、流体的粘性引起剪切力，进而导致耗散。 5、无黏流体—无摩擦—流动不分离—无尾迹。 6、流体的特性：连续性、易流动性、压缩性 不可压缩流体：0D Dt ρ= const ρ=是针对流体中的同一质点在不同时刻保持不变，即不可压缩流体的密度在任何时刻都保持不变。是一个过程方程。 7、流体的几种线 流线：是速度场的向量线，是指在欧拉速度场的描述； 同一时刻、不同质点连接起来的速度场向量线； (),0dr U x t dr U ??=
迹线：流体质点的运动轨迹，是流体质点运动的几何描述； 同一质点在不同时刻的位移曲线； 涡线：涡量场的向量线，(),,0U dr x t dr ωωω=????= 涡线的切线和当地的涡量或准刚体角速度重合，所以，涡线是流体微团 准刚体转动方向的连线，形象的说：涡线像一根柔性轴把微团穿在一起。 第二章流体静力学 1、压强：0lim A F dF p A dA ?→?==? 静止流场中一点的应力状态只有压力。 2、流体的平衡状态： 1）、流体的每个质点都处于静止状态，==整个系统无加速度； 2）、质点相互之间都没有相对运动，==整个系统都可以有加速度； 由于流体质点之间都没有相对运动，导致剪应力处处为零，故只有： 体积力(重力、磁场力)和表面力(压强和剪切力)存在。 3、表面张力：两种不可混合的流体之间的分界面是曲面，则在曲面两边存在一 个压强差。 4、正压流场：流体中的密度只是压力(压强)的单值函数。() dp p ρ? 5、涡量不生不灭定理 拉格朗日定理：理想正压流体在势力场中运动时，如某一时刻连续流场无旋，则 流场始终无旋。0,,ndA U ωω?==??? 有斯托克斯公式得：00,A l U x ndA δωΓ=?=?=??
水 水银 题1图 1 2 3 题型一：曲面上静水总压力的计算问题（注：千万注意方向，绘出压力体） 1、AB 曲面为一圆柱形的四分之一，半径R=0.2m ，宽度（垂直纸面）B=0.8m ，水深H=1.2m ，液体密度3/850m kg =ρ，AB 曲面左侧受到液体压力。求作用在AB 曲面上的水平分力和铅直分力。（10分） 解：（1）水平分力：RB R H g A h P z c x ?-==)2 (ργ…….(3分) N 1.14668.02.0)22 .02.1(8.9850=??- ??=，方向向右（2分）。 （2）铅直分力：绘如图所示的压力体，则 B R R R H g V P z ??? ? ????+-==4)(2πργ……….（3分） 1.15428.04 2.014.32.0)2.02.1(8.98502=???? ? ?????+?-??=，方向向下（2分） 。 2．有一圆滚门，长度l=10m ，直径D=4.2m ，上游水深H1=4.2m ，下游水深H2=2.1m ，求作用于圆滚门上的水平和铅直分压力。
解题思路：（1）水平分力： l H H p p p x )(2 1222121-=-=γ 方向水平向右。 （2）作压力体，如图，则 l D Al V p z 4 432 πγγγ? === 方向垂直向上。 3.如图示，一半球形闸门，已知球门的半径m R 1= ，上下游水位差m H 1= ，试求闸门受到的水平分力和竖直分力的大小和方向。 解： (1)水平分力： ()2R R H A h P c πγγ?+===左,2R R A h P c πγγ?='=右 右左P P P x -= kN R H 79.30114.31807.92=???=?=πγ， 方向水平向右。 （２）垂直分力： V P z γ=，由于左、右两侧液体对曲面所形成的压力体均为半球面，且两侧方向相反，因而垂直方向总的 压力为0。 4、密闭盛水容器，已知h 1=60cm,h 2=100cm ，水银测压计读值cm h 25=?。试求半径R=0.5m 的半球盖AB 所受总压力的水平分力和铅垂分力。
课后习题 第一章 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基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
3347流体力学全国自考 第一章绪论 1、液体和气体统称流体，流体的基本特性是具有流动性。流动性是区别固体和流体的力学特性。 2、连续介质假设：把流体当作是由密集质点构成的、内部无空隙的连续踢来研究。 3、流体力学的研究方法：理论、数值和实验。 4、表面力：通过直接接触，作用在所取流体表面上的力。 5、质量力：作用在所取流体体积内每个质点上的力，因力的大小与流体的质量成比例，故称质量力。重力是最常见的质量力。 6、与流体运动有关的主要物理性质：惯性、粘性和压缩性。 7、惯性：物体保持原有运动状态的性质；改变物体的运功状态，都必须客服惯性的作用。 8、粘性：流体在运动过程中出现阻力，产生机械能损失的根源。粘性是流体的内摩擦特性。粘性又可定义为阻抗剪切变形速度的特性。 9、动力粘度：是流体粘性大小的度量，其值越大，流体越粘，流动性越差。 10、液体的粘度随温度的升高而减小，气体的粘度随温度的升高而增大。 11、压缩性：流体受压，分子间距离减小，体积缩小的性质。 12、膨胀性：流体受热，分子间距离增大，体积膨胀的性质。 13、不可压缩流体：流体的每个质点在运动过程中，密度不变化的流体。 14、气体的粘度不受压强影响，液体的粘度受压强影响也很小。 第二章流体静力学 1、精致流体中的应力具有一下两个特性： 应力的方向沿作用面的内法线方向。 静压强的大小与作用面方位无关。 2、等压面：流体中压强相等的空间点构成的面；等压面与质量力正交。 3、绝对压强是以没有气体分子存在的完全真空为基准起算的压强、 4、相对压强是以当地大气压强为基准起算的压强。 5、真空度：若绝对压强小于当地大气压，相对压强便是负值，有才呢个·又称负压，这种状态用真空度来度量。 6、工业用的各种压力表，因测量元件处于大气压作用之下，测得的压强是改点的绝对压强超过当地大气压的值，乃是相对压强。因此，先跪压强又称为表压强或计示压强。 7、z+p/ρg=C: z为某点在基准面以上的高度,可以直接测量,称为位置高度或位置水头.。 p/ρg=h p，称为测压管高度或压强水头,其物理意义是单位重量的液体具有的压强势能,简称压能。 z+p/ρg称为测压管水头，是单位重量液体具有的总势能，其物理意义是静止液体中各点单位重量液体具有的总势能相等。 第三章流体动力学基础 1、描述流体运动的两种方法：拉格朗日法和欧拉法。 2、拉格朗日法：从整个流体运动是无数个质点运动的综合出发，以个别质点为观察对象来描述，再讲每个质点的运动情况汇总起来，就描述了流体的整个流动。 3、欧拉法：以流体运动的空间点作为观察对象，观察不同时刻各空间点上流体质点的运动，再将每个时刻的情况汇总起来，就描述了整个运动。
流体力学综述 奇妙的流体力学：流体力学是研究流体在受到一系列力和边界条件作用时流体的运动和内部应力的科学。生活中我们与流体力学息息相关，很多貌似很奇怪、很难理解的现象，如河流中沙丘的形成、浴帘效应、杯中水的涡旋等现象，都可以用流体力学的只是来解释 流体力学的特点： 1、研究对象范围广：大到宇宙中的天体星云，小到人体内的毛细血管大气运动、沙漠迁移、河流泥沙运动、管道中煤粉输送、化工中气体催化剂的运动海水运动(包括波浪、潮汐、中尺度涡旋、环流等)乃至地球深处熔浆的流动，都属于流体力学的研究范畴。 2、研究历史悠久：流体力学是在人类同自然界作斗争和在生产实践中逐步发展起来的。公元前2000余年中国有大禹治水疏通江河的传说；秦朝李冰父子带领劳动人民修建的都江堰，至今还在发挥着作用；大约与此同时，古罗马人建成了大规模的供水管道系统等等，这都是人们对于自然的探索与改造。箭弩的发明反映了原始人对箭头的流线型降低摩阻及尾翅的稳定性问题的探索。 3、对整个自然科学贡献大：流体力学为自然科学的研究提供里一个完整的体系，并且对整个自然科学的一些根本性的东西产生了重要的影响。作为与量子力学、相对论相齐名的一个重大科学理论，混沌理论自产生以来产生的巨大影响同时也被广泛应用于各领域。 流体力学的研究内容： 1、流体静力学：主要研究流体处于静止状态时的力和平衡关系。 ? 浮力规律的探讨—阿基米德（十七世纪以前） 对流体力学学科的形成作出第一个贡献的是古希腊的阿基米德，静力学和流体静力学的奠基人，他建立了包括物理浮力定律和浮体稳定性在内的液体平衡理论，奠定了流体静力学的基础。他的著作《论浮体》相当详细地讨论了正回旋抛物体在流体中的稳定性，研究了不同的高与底的比、具有不同的比重及在流体中处于不同位置时这种立体的性态。 ? 帕斯卡原理（静压传递原理）（1651-1654） 其基本内容是加在密闭液体任何一部分上 const 2g 2 =+g v p ρ 的压强，必然按照其原来的大小由液体向各个方向传递。帕斯卡的著作《论液体的平衡和空气的重力》代表了十七世纪力学发展的里程碑。静压传递原理是液压与气压传动的基础，在液压千斤顶、液压机等机械上有非常广泛的应用。 F1／F2（F 为施加的力）＝S2／S2（S 指大小活塞的面积） 2、 理想流体的运动学和动力学 ：流体动力学中主要研究无粘性不可压缩流体在绕过物体时的流动和管内流动规律的一个分支，又称经典流体动力学。这一学科分支的任务是求解流场中的速度、压力分布和物体受力。它忽略了真实流体的粘性和压缩性，也不考虑表面张力，从而大大简化了复杂的流体动力学问题，故常作为近似处理许多工程问题的依据。 牛顿、伯努利、欧拉、拉格朗日、拉普拉斯、纳维、柯西、泊松、圣维南、斯托克斯、雷诺等科学家对流体力学的发展都做出了巨大贡献。 ? 牛顿-站在巨人的肩上：他的著作《自然哲学的数学原理》发表于1687年，牛顿内摩擦定律-第一个系统研究流体力学的人。他所研究的流体我们定义为“牛顿流体”，是指在受力
试卷 1. 简述计算流体力学的特点及其应用领域。 CFD 是以计算机作为模拟手段，运用一定的计算技术寻求流体力学各种复杂问题的离散化数值解。它的主要特征：(1)数值解而不是解析解；(2)计算技术起关键作用；(3)与计算机的发展紧密相关。（成本较低，适用范围宽，可靠性差，表达困难）应用领域：航空、航天、气象、船舶、武器装备、 水利、化工、建筑、机械、汽车、海洋、体育、环境、卫 生等 2. 等步长网格分布情况下u x ??的一阶向前差分、22u x ??的二阶中心差分表达式。（P89） 一阶向前差分：1,,,()i j i j i j u u u x x x +-?=+O ???（） 二阶中心差分：21,,1,2,22 2()()i j i j i j i j u u u u x x x +--+?=+O ???（） 3. 简答题 1) 什么是差分方程的相容性? 差分方程与微分方程的差别是截断误差R 。必要时通过缩小空间步长（网格尺寸）h 和时间步长t ，这一误差应可缩小至尽可能小。当h->0和t->0时，若R->0，则差分方程趋于微分方程，表示这两个方程是一致的。这时称该差分方程与微分方程是相容的。 2) 什么是差分解的收敛性? 当微分方程在离散为差分方程来求解，当步长h 0→时，存在着差分方程的解 n y 能够收敛到微分方程的准确解y()n x ，这就是差分方法的收敛性。 收敛性定义：对于任意节点的0n x x nh =+，如果数值解n y 当h 0→（同时n →∞）时趋向于准确解y()n x ，则称该方法是收敛的。 3) 什么是差分解的稳定性? 数值计算时，除计算机舍入误差（字长有限）外，初始条件或方程中某些常数项 也有可能给的不尽精确。舍入误差和这些误差在计算过程中可能一步步积累与传 递，误差的传递，有时可能变大，有时可能变小。某一步舍入误差放大或缩小的
计算流体力学的发展及应用 计算流体力学的发展： 20世纪30年代，由于飞机工业的需要、要求用流体力学理论来了解和指导飞机设计，当时由于飞行速度很低，可以忽略粘性和旋涡，因此流动的模型为拉普拉斯方程，研究工作的重点是椭圆型方程的数值解。利用复变函数理论和解的迭加方法来求解析解。随着飞机外形设计越来越复杂，出现了求解奇异边界积分方程的方法。以后为了考虑粘性效应，有了边界层方程的数值计算方法，并发展成以位势方程为外流方程，与内流边界层方程相结合，通过迭代求解粘性干扰流场的计算方法。同一时期许多数学家研究了偏微分方程的数学理论，Courant，Fredric等人研究了偏微分方程的基本特性、数学提法的适定性、物理波的传播特性等问题，发展了双曲型偏微分方程理论。以后，Courant，Fredric，Lowy等人发表了经典论文，证明了连续的椭圆型、抛物型和双曲型方程组解的存在性和唯一性定理，并针对线性方程的初值问题，首先将偏微分方程离散化，然后证明了离散系统收敛到连续系统，最后利用代数方法确定了差分解的存在性；他们还给出了著名的稳定性判别条件：CFL条件。这些工作是差分方法的数学理论基础。20世纪40年代，V onNeumann，Richmyer，Hopf，Lax和其他一些学者建立了非线性双曲型方程守恒定律的数值方法理论，为含有激波的气体流动数值模拟打下了理论基础。
在20世纪50年代，仅采用当时流体力学的方法，研究比较复杂的非线性流动现象是不够的，特别是不能满足高速发展起来的宇航飞行器绕流流场特性研究的需要。针对这种情况，一些学者开始将基于双曲型方程数学理论基础的时问相关方法用于求解宇航飞行器的气体的定常绕流场问题，这种方法虽然要求花费更多的计算机时，但因数学提法适定，又有较好的理论基础，且能模拟流体运动的非定常过程，所以在60年代这是应用范围较广的一般方法。以后由Lax、Kais和其他著者给出的非定常偏微分方程差分逼近的稳定性理论，进一步促进了时间相关方法。当时还出现了一些针对具体问题发展起来的特殊算法。 进人2O世纪80年代以后，计算机硬件技术有了突飞猛进的发展，计算机逐渐进人人们的实践活动范围。随着计算方法的不断改进和数值分析理论的发展高精度模拟已不再是天方夜谭。同时随着人类生产实践活动的不断发展，科学技术的日新月异，一大批高新技术产业对计算流体力学提出了新的要求，同时也为计算流体力学的发展提供了新的机遇。实践与理论的不断互动，形成计算流体力学的新热点、新动力，从而推动计算流体力学不断向前发展。首先，在计算模型方面，又提出了一些新的模型，如新的大涡模拟模型、考虑壁面曲率等效应的新的湍流模式、新的多相流模式、新的飞行器气动分析与热结构的一体化模型等这就使得计算流体力学的计算模型由最初的Euler和Ⅳ—s方程，扩展到包括湍流、两相流、化学非平衡、太阳风等问题研究模型
借宝地写几个小短文，介绍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】的经典结论进行对比，评述本
(计算流体力学概述) CFD仿真 3月20日309 计算流体力学概述 流体力学，是研究流体（液体和气体）的力学运动规律及其应用的学科。主要研究在各种力的作用下，流体本身的状态，以及流体和固体壁面、流体和流体间、流体与其他运动形态之间的相互作用的力学分支。流体力学是力学的一个重要分支，它主要研究流体本身的静止状态和运动状态，以及流体和固体界壁间有相对运动时的相互作用和流动的规律。在生活、环保、科学技术及工程中具有重要的应用价值。 计算流体力学的发展 计算流体动力学（Computational Fluid Dynamics）简写为CFD，是20世纪60年代起伴随计算科学与工程（Computational Science and Engineering, 简称CSE）迅速崛起的一门学科分支，经过半个世纪的迅猛发展，这门学科已经是相当的成熟了，一个重要的标志就是近几十年来，各种CFD通用软件的陆续出现，成为商品化软件，服务于传统的流体力学和流体工程领域，如航空、航天、船舶、水利等。随着CFD通用软件的性能日益完善，应用的范围也不断的扩大，在化工、冶金、建筑、环境等相关领域中也被广泛应用。 现代流体力学研究方法包括理论分析，数值计算和实验研究三个方面。这些方法针对不同的角度进行研究，相互补充。理论分析研究能够表述参数影响形式，为数值计算和实验研究提供了有效的指导；试验是认识客观现实的有效手段，验证理论分析和数值计算的正确性；计算流体力学通过提供模拟真实流动的经济手段补充理论及试验的空缺。 更重要的是，计算流体力学提供了廉价的模拟、设计和优化的工具，以及提供了分析三维复杂流动的工具。在复杂的情况下，测量往往是很困难的，甚至是不可能的，而计算流体力学则能方便的提供全部流场范围的详细信息。与试验相比，计算流体力学具有对于参数没有什么限制，费用少，流场无干扰的特点。出于计算流体力学如此的优点，我们选择它来进行模
《流体力学考》考点重点知识归纳 1.流体元：就有线尺度的流体单元，称为流体“质元”，简称流体元。流体元可看做大量流体质点构成的微小单元。 2.流体质点：（流体力学研究流体在外力作用下的宏观运动规律） （1）流体质点无线尺度，只做平移运动 （2）流体质点不做随即热运动，只有在外力的作用下作宏观运动； （3）将以流体质点为中心的周围临街体积的范围内的流体相关特性统计的平均值作为流体质点的物理属性； 3.连续性介质模型的内容：根据流体指点概念和连续介质模型，每个流体质点具有确定的宏观物理量，当流体质点位于某空间点时，若将流体质点的物理量，可以建立物理的空间连续分布函数，根据物理学基本定律，可以建立物理量满足的微分方程，用数学连续函数理论求解这些方程，可获得该物理量随空间位置和时间的连续变化规律。 4.连续介质假设：假设流体是有连续分布的流体质点组成的介质。 5.牛顿的粘性定律表明：牛顿流体的粘性切应力与流体的切变率成正比，还表明对一定的流体，作用于流体上的粘性切应力由相邻两层流体之间的速度梯度决定的，而不是由速度决定的： 6.牛顿流体：动力粘度为常数的流体称为牛顿流体。 7.分子的内聚力：当两层液体做相对运动时，两层液体的分子的平均距离加大，分子间的作用力变现为吸引力，这就是分子的内聚力。 液体快速流层通过分子内聚力带动慢流层，漫流层通过分子的内聚力阻滞快流层的运动，表现为内摩擦力。、 流体在固体表面的不滑移条件：分子之间的内聚力将流体粘附在固体表面，随固体一起运动或静止。 8.温度对粘度的影响：温度对流体的粘度影响很大。液体的粘度随温度升高而减小，气体的粘度则相反，随温度的升高而增大。 压强对粘性的影响：压强的变化对粘度几乎没有什么影响，只有发生几百个大气压的变化时，粘度才有明显改变，高压时气体和液体的粘度增大。 9.描述流体运动的两种方法 拉格朗日法：拉格朗日法又称为随体法。它着眼于流体质点，跟随流体质点一起运动，记录流体质点在运动过程中会各种物理量随所到位置和时间的变化规律，跟中所有质点便可了解整个流体运动的全貌。 欧拉法：欧拉法又称当地法。它着眼于空间点，把流体的物理量表示为空间位置和时间的函数。空间点的物理量是指，某个时刻占据空间点的。 流体质点的物理量，不同时刻占据该空间点的流体质点不同。 10.速度场：速度场是由流体空间各个坐标点的速度矢量构成的场。速度场不仅描述速度矢量的空间分布，还可描述这种分布随时间的变化。 11.毛细现象：玻璃管内的液体在表面张力的作用下液面升高或降低的现象称为毛细现象； 12.迹线：流体质点运动的轨迹。在流场中对某一质点作标记，将其在不同时刻的所在位置点连成线就是该流体质点的迹线。 13.定常流动：流动参数不随时间变化的流动。反之，流体参数随时间变化的流动称为不定长流动。 14.流线：流线是指示某一时刻流场中各点速度矢量方向的假象曲线。