计算流体力学课程总结 计算流体动力学(computational Fluid Dynamics,简称CFD)是通过计算机数值 计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。是用电子计算机和离散化的数值方法对流体力学问题进行数值模拟和分析的一个分支。 流体力学和其他学科一样，是通过理论分析和实验研究两种手段发展起来的。很早就已有理论流体力学和实验流体力学两大分支。理论分析是用数学方法求出问题的定量结果。但能用这种方法求出结果的问题毕竟是少数，计算流体力学正是为弥补分析方法的不足而发展起来的。计算流体力学是目前国际上一个强有力的研究领域,是进行传热、传质、动量传递及燃烧、多相流和化学反应研究的核心和重要技术，广泛应用于航天设计、汽车设计、生物医学工业、化工处理工业、涡轮机设计、半导体设计、HAVC&R 等诸多工程领域。 计算流体力学的任务是流体力学的数值模拟。数值模拟是“在计算机上实现的一 个特定的计算，通过数值计算和图像显示履行一个虚拟的物理实验——数值实验“。 数值模拟包括以下几个部分。首先，要建立反映问题(工程问题、物理问题等)本质数 学模型。其次，数学模型建立以后需要解决的问题是寻求高效率、高准确度的计算方法。再次，在确定了计算方法和坐标系统后，编制程序和进行计算式整个工作的主体。最后，当计算工作完成后，流畅的图像显示是不可缺少的部分。 还有一个就是CFD的基本思想问题，它就是把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通 过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求 解代数方程组获得场变量的近似值。 经过四十多年的发展，CFD出现了多种数值解法。这些方法之间的主要区别在于 对控制方程的离散方式。根据离散的原理不同，CFD大体上可分为三个分支： ?有限差分法(Finite Different Method，FDM) ?有限元法(Finite EIement Method，FEM) ?有限体积法(Finite Volume Method，FVM) 有限差分法是应用最早、最经典的CFD方法，也是最成熟、最常用的方法。它将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的 导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。求出差分万程组 的解，就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题 的近似数值解法。
自然环境和工程装置中的流动常常是湍流流动，模拟任何实际过程首先遇到的就是湍流问题，而湍流问题本身又是流体力学理论上的难题。 对湍流最根本的模拟方法是在湍流尺度的网格尺寸内求解瞬态的三维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 ,在此长度内流体质点运动是自由的（不与
计算流体力学软件CFD在燃烧器设计中的应用探讨[摘要]本文通过对目前燃烧器的现状与技术发展的研究，探讨计算流体力学 软件CFD在燃烧器设计中应用的必要性和可行性，以CFD(计算流体力学)软件为工具，以普通大气式燃烧器为研究对象，采用实验和理论相结合的方法，充分利用现代计算机技术，达到降低燃烧器设计成本和研制费用的目的。 [关键词]燃烧器数值模拟计算流体力学 一、燃烧器的发展现状 1.部分预混式燃烧器的产生及其原理 燃烧的方法被分为扩散式燃烧、部分预混式燃烧和完全预混式燃烧。扩散式燃烧易产生不完全燃烧产物，燃烧温度很低，并未充分利用燃气的能量；而一旦预先混入一部分空气后火焰就会变的清洁，燃烧温度也可以提高，燃烧较充分。完全预混燃烧（无焰燃烧）要求事先按照化学当量比将燃气和空气均匀混合（实际应用中空气系数要大于1），燃烧充分，火焰温度很高，但稳定性较差，易回火。所以民用燃具多采用部分预混式燃烧。 1855年工程师本生发明了一种燃烧器，能从周围大气中吸入一些空气和燃气预混，在燃烧时形成不发光的蓝色火焰，这就是实验室常用的本生灯（单火孔燃烧器）。这种燃烧技术就被称作部分预混式燃烧。 本生灯燃烧所产生的火焰为部分预混层流火焰（俗称本生火焰）。它由内焰，外焰及燃烧区域外围肉眼看不见的高温区组成。火焰一般呈锥体状。燃气—空气的混合气体先在内锥燃烧，中间产物及未燃尽的部分便从锥内向外流出，且混合气体出流的速度与内锥表面火焰向内传播速度相互平衡，此外便形成一个稳定的焰面，呈蓝色。而未燃烧尽的混合气体残余物继续与大气中的空气进行二次混合燃烧，形成火焰外锥。如图1所示，完成燃烧后产生高温co2和水进而在外焰的外侧形成外焰膜（肉眼看不见的高温层）: 图1. 本生灯示意图 如果混合气流是处于层流状态，则外焰面呈较光滑的锥形；如果处于紊流状态，则外焰面产生褶皱，直至产生强烈扰动，气团不断飞散、燃尽。
黑龙江科技大学建筑工程二学历实践报告 流体力学实践报告 一、实践概述 在此次实践中,老师给我演示了雷诺试验与伯努利方程试验。下面我就实践的主要内容进行一下总结。 二、雷诺实验 (一)、实验目的 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 为直线的斜率
计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 一 前 言 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30～40 年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学" 。 从20 世纪60 年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
计算流体力学软件Fluent在烟气脱硫中的应用 ０引言 污染最为有效的方法之一，而石灰石—石膏湿烟气脱硫是目前能大规模控制燃煤造成ＳＯ ２ 法脱硫技术以其脱硫效率高、吸收剂来源丰富、成本低廉、技术成熟和运行可靠等优点获得广泛应用．从气液两相流体力学和化学反应动力学的观点看，脱硫吸收塔内流体流动的目的是强化气液两相的混合和质量传递、延长气液两相在塔内的接触时间、增大气液两相的接触面积并尽量减小吸收塔的阻力．合理的塔内流场分布对提高脱硫效率、降低脱硫投资和运行成本都具有重要意义． 目前，国内外对烟气脱硫吸收塔进行大量研究，主要采用实验方法，如研究塔的阻力特性、液滴运动速度沿塔高变化和ＴＣＡ塔内温度场分布等，这些研究对指导工业应用具有重要意义，但其结果往往只针对特定的设备或结构，具有较大的局限性．随着计算机技术的迅速发展，计算流体力学（ＣｏｍｐｕｔａｔｉｏｎａｌＦｌｕｉｄＤｙｎａｍｉｃ，ＣＦＤ）已成为研究三维流动的重要方法：周山明等［４］利用ＦＬＵＥＮＴ计算空塔和喷淋状态下的塔热态流场，结果表明脱硫吸收塔入口处流场变化最剧烈、压降损失最大，并根据计算结果改造来流烟道；孙克勤等采用混合网格和随机颗粒生成模型对烟气脱硫吸收塔的热态流场进行数值模拟；郭瑞堂等采用ＦＬＵＥＮＴ结合非稳态反应传质－反应理论对湿法脱硫液柱冲的吸收进行数值模拟． 击塔内的流场和ＳＯ ２ 本文尝试应用ＦＬＵＥＮＴ对某脱硫吸收塔内烟气脱硫过程进行初步数值模拟，通过对内部流场进行分析验证本文模拟的合理性，进而对脱硫过程中脱硫吸收塔内是否存在湿壁现象进行深入分析研究． １基于ＲＡＮＳ求解器的ＣＦＤ数值模拟 方法 １．１控制方程 时均的不可压缩连续性方程和Ｎ Ｓ方程 （ＲＡＮＳ方程）如下： １．２湍流模型和多相流模型
谈流体力学的研究内容及发展简史 流体力学是力学的一个独立分支，是一门研究流体的平衡和流体机 械运动规律及其实际应用的技术科学，在许多工业部门中都有着广泛应 用，航空工业中飞机的制造离不开空气动力学；造船工业部门要用到水 动力学，与土建类各专业有着更加密切的关系，了解流体动力学的研究 内容及发展简史对学习流体力学知识具有的一定的引导作用，为以后的 学习铺设台阶，引起学习的兴趣。 流体力学的研究内容 流体是气体和液体的总称。在人们的生活和生产活动中随时随地都 可遇到流体，所以流体力学是与人类日常生活和生产事业密切相关的。 大气和水是最常见的两种流体，大气包围着整个地球，地球表面的70% 是水面。大气运动、海水运动(包括波浪、潮汐、中尺度涡旋、环流等) 乃至地球深处熔浆的流动都是流体力学的研究内容。 流体力学既包含自然科学的基础理论，又涉及工程技术科学方面的 应用。此外，如从流体作用力的角度，则可分为流体静力学、流体运动 学和流体动力学；从对不同“力学模型”的研究来分，则有理想流体动力 学、粘性流体动力学、不可压缩流体动力学、可压缩流体动力学和非牛 顿流体力学等。 在流体力学中为简化计算，对流体模型做出了假设：质量守恒;动量 守恒;能量守恒。 在流体力学中常会假设流体是不可压缩流体，也就是流体的密 度为一定值。液体可以算是不可压缩流体，气体则不是。有时也会 假设流体的黏度为零，此时流体即为非粘性流体。气体常常可视为 非粘性流体。若流体黏度不为零，而且流体被容器包围（如管子）， 则在边界处流体的速度为零。 流体的主要物理性质： 1、流体：只能承受压力，一般不能承受拉力与抵抗拉伸变形。液体 有一定的体积，存在一个自由液面；气体能充满任意形状的容器，无一 定的体积，不存在自由液面。 2、流体的连续介质模型 微观：流体是由大量做无规则运动的分子组成的，分子之间存在空隙，但在标准状况下，1cm3液体中含有3.3×1022个左右的分子，相邻分子间的距离约为3.1×10-8cm。1cm3气体中含有2.7×1019个左右的分子，相邻分子间的距离约为3.2×10-7cm。 宏观：考虑宏观特性，在流动空间和时间上所采用的一切特征尺度和特征时间都
借宝地写几个小短文，介绍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、输出网格。
流体力学报告 每一门力学学科的建立,都需要建立模型,也就是把实际的问题抽象化,而抽象过程就是把现实中对所研究问题不重要的因素忽略掉,也就是模型假设,从而建立于这个问题相适应的模型进行研究,如果有意义有价值,也就慢慢深入研究,从而形成一门学科,它们都是随社会的发展而发展形成的.比如现如今最前沿的力学学科"纳米力学"就是如此。我们土木工程常说的三大力学有：1.理论力学---分析力学,振动力学,水力学或称为流体力学(这些研究对材料都不太侧重 )2.材料力学---弹性力学,塑性力学(都是又材料特性而分的) 3.结构力学:就是分析复杂的结构的情形。在此我重点叙述我对流体力学这门课学科的学习和认知。 一·流体的基本信息解释： 流体，是与固体相对应的一种物体形态，是液体和气体的总称. 由大量的、不断地作热运动而且无固定平衡位置的分子构成的，它的基本特征是没有一定的形状并且具有流动性。流体都有一定的可压缩性，液体可压缩性很小，而气体的可压缩性较大，在流体的形状改变时，流体各层之间也存在一定的运动阻力(即粘滞性)。当流体的粘滞性和可压缩性很小时，可近似看作是理想流体，它是人们为研究流体的运动和状态而引入的一个理想模型。是液压传动和气压传动的介质。大气和水是最常见的两种流体，大气包围着整个地球，地球表面的70%是水面。大气运动、海水运动(包括波浪、潮汐、中尺度涡旋、环流等)乃至地球深处熔浆的流动都是流体的研究内容。
二·流体力学的阐述： 流体力学是连续介质力学的一门分支，是研究流体(包含气体，液体以及等离子态)现象以及相关力学行为的科学。可以按照研究对象的运动方式分为流体静力学和流体动力学，还可按流动物质的种类分为水力学，空气动力学等等。对流体力学学科的形成作出第一个贡献的是古希腊的阿基米德，他建立了包括物理浮力定律和浮体稳定性在内的液体平衡理论，奠定了流体静力学的基础，特别是从20世纪以来，流体力学已发展成为基础科学体系的一部分，同时又在工业、农业、交通运输、天文学、地学、生物学、医学等方面得到广泛应用。流体力学既包含自然科学的基础理论，又涉及工程技术科学方面的应用。此外，如从流体作用力的角度，则可分为流体静力学、流体运动学和流体动力学;从对不同"力学模型"的研究来分，则有理想流体动力学、粘性流体动力学、不可压缩流体动力学、可压缩流体动力学和非牛顿流体力学等。 三·对流体的研究假设： 连续体假设 物质都由分子构成,尽管分子都是离散分布的,做无规则的热运动.但理论和实验都表明,在很小的范围内,做热运动的流体分子微团的统计平均值是稳定的.因此可以近似的认为流体是由连续物质构成,其中的温度,密度,压力等物理量都是连续分布的标量场。 质量守恒 质量守恒目的是建立描述流体运动的方程组。欧拉法描述为:流进
第26卷第2期河北科技大学学报Vol.26,No.2 2005年6月Journal of Hebei University of Science and T echnology June2005 文章编号:100821542(2005)022******* 计算流体力学(CFD)的通用软件 翟建华 (河北科技大学国际交流与合作处,河北石家庄050018) 摘要:对化学工程领域中的通用CFD(Computational Fluid Dynamics)模拟软件Phoenics,Flu2 ent,CFX等的具体特点和应用情况进行了综述,指出了他们各自的结构特点、特有模块、包含的数学模型和成功应用领域;给出了选用CFD软件平台的7项准则,对今后CFD技术的发展进行了预测,指出,今后CFD研究的主要方向将集中在数学模型开发、工程改造和新设备开发及与工艺软件的匹配连用等方面。 关键词:计算流体力学;模拟软件;CFX;FLUENT;PH OENICS 中图分类号:T Q015.9文献标识码:A Review of commercial CFD software ZH AI Jian2hua (Department of Int ernation Exchange and Cooperation,H ebei University of Science and Technology,Shijiazhuang H ebei 050018,China) Abstr act:The paper summar izes the features and application of the CF D simulation software like Phoenics,F luent and CFX etc in chemical engineering,and discusses their str ucture features,special modules,mathematical models and successful application areas.It also puts forward seven r ules for the good choice of commercial CF D code for the CF D simulation resea rcher s.Based on t he predict ion of the technology development,it points out the possible r esear ch direction for CF D in the future will focus on the development of mathematical model,project transformat ion,new equipment and their matching application with technologi2 cal softwa re. Key words:CF D;simulation software;CF X;FLUENT;P HOENICS CFD(Computational Fluid Dynamics)软件是计算流体力学软件的简称,是用来进行流场分析、计算、预测的专用工具。通过CFD模拟,可以分析并且显示流体流动过程中发生的现象,及时预测流体在模拟区域的流动性能,并通过各种参数改变,得到相应过程的最佳设计参数。CFD的数值模拟,能使我们更加深刻地理解问题产生的机理,为实验提供指导,节省以往实验所需的人力、物力和时间,并对实验结果整理和规律发现起到指导作用。随着计算机软硬件技术的发展和数值计算方法的日趋成熟,出现了基于现有流动理论的商用CFD软件。这使许多不擅长CFD工作的其他专业研究人员能够轻松地进行流体数值计算,从而使研究人员从编制繁杂、重复性的程序中解放出来,以更多的精力投入到研究问题的物理本质、问题提法、边界(初值)条件和计算结果的合理解释等重要方面上,充分发挥商用CFD软件开发人员和其他专业研究人员各自的智力优势,为解决实际工程问题开辟了道路。 CFD研究走过了相当漫长的过程。早期数值模拟阶段,由于缺乏模拟工具,研究者一般根据自身工作性质和研究过程,自行编制模拟程序,其优点是针对性强,对具体问题的解决有一定精度,但是,带来的问题 收稿日期:2004208221;修回日期:2004211221;责任编辑:张军 作者简介:翟建华(19642),男,河北平乡人,教授,主要从事化工CFD、高效传质与分离和精细化工方面的研究。
计算流体力学结课报告200Km/h列车fluent仿真计算 学部：化、环、生学部 学院：化工机械与安全学院 学号：31507095 班级：化1512班 学生姓名：孙金
引言 数值仿真就是对所建立的数值模型进行数值实验和求解的过程。而计算流体力学CFD （Computational Fluid Dynamics）就是在工程仿真实验领域中应用最广泛的一门学科。任何流体运动的规律都是以质量守恒定律、动量守恒定律和能量守恒定律为基础的。这些基本定律可由数学方程组来描述，如欧拉方程、N-S方程。采用数值计算方法，通过计算机求解这些控制流体流动的数学方程，进而研究流体的运动规律这就是CFD研究问题的方法。在实际计算流体力学方面，采用通用的CFD软件来完成工程上的一些流体力学问题，有极为广泛的应用前景。近年来，随着计算机技术以及相关技术的发展，CFD技术已经在工程领域内取得重大的进步，特别是在高速列车的外型设计方面起了很大作用。随着国家经济的发展，国家运输业也有了很大的发展，特别是列车经过几次提速后，高速列车在国家运输行业中所占比例不断提高。高速列车的特点是庞大、细长、在地面轨道上运行，其空气动力学问题非常复杂。空气在列车表面形成空气流场，空气阻力急剧增加，作用在列车的阻力大部分来自压强阻力，而一部分来自表面磨擦阻力，这就使能耗过大，同时列车可能出现较大的空气升力，导致列车产生“飘”的现象，激发列车脱轨事故的发生，因此研究高速列车气动力性能非常重要。用CFD仿真可以详细了解高速列车的空气动力特性，从而设计出阻力小、噪音低等各方面性能完善的高质量列车。本文采用CFD学科中的常用商业软件Fluent仿真一个时速200km/h的二维流线型车头的外流场，对其空气动力性能进行分析，从而得到不同车辆形状其周围流场的不同，进而分析哪种车型更适合。
计算流体力学（CFD）是近代流体力学，数值数学和计算机科学结合的产物，是一门具有强大生命力的边缘科学。它以电子计算机为工具，应用各种离散化的数学方法，对流体力学的各类问题进行数值实验、计算机模拟和分析研究，以解决各种实际问题。 计算流体力学和相关的计算传热学，计算燃烧学的原理是用数值方法求解非线性联立的质量、能量、组分、动量和自定义的标量的微分方程组，求解结果能预报流动、传热、传质、燃烧等过程的细节，并成为过程装置优化和放大定量设计的有力工具。计算流体力学的基本特征是数值模拟和计算机实验，它从基本物理定理出发，在很大程度上替代了耗资巨大的流体动力学实验设备，在科学研究和工程技术中产生巨大的影响。目前比较好的CFD软件有：Fluent、CFX,Phoenics、Star-CD，除了Fluent 是美国公司的软件外，其它三个都是英国公司的产品 ------------------------------------------------------ FLUENT FLUENT是目前国际上比较流行的商用CFD软件包，在美国的市场占有率为60%。举凡跟流体，热传递及化学反应等有关的工业均可使用。它具有丰富的物理模型、先进的数值方法以及强大的前后处理功能，在航空航天、汽车设计、石油天然气、涡轮机设计等方面都有着广泛的应用。其在石油天然气工业上的应用包括：燃烧、井下分析、喷射控制、环境分析、油气消散/聚积、多相流、管道流动等等。 Fluent的软件设计基于CFD软件群的思想，从用户需求角度出发，针对各种复杂流动的物理现象，FLUENT软件采用不同的离散格式和数值方法，以期在特定的领域内使计算速度、稳定性和精度等方面达到最佳组合，从而高效率地解决各个领域的复杂流动计算问题。基于上述思想，Fluent开发了适用于各个领域的流动模拟软件，这些软件能够模拟流体流动、传热传质、化学反应和其它复杂的物理现象，软件之间采用了统一的网格生成技术及共同的图形界面，而各软件之间的区别仅在于应用的工业背景不同，因此大大方便了用户。其各软件模块包括： GAMBIT——专用的CFD前置处理器，FLUENT系列产品皆采用FLUENT公司自行研发的Gambit 前处理软件来建立几何形状及生成网格，是一具有超强组合建构模型能力之前处理器，然后由Fluent 进行求解。也可以用ICEM CFD进行前处理，由TecPlot进行后处理。 Fluent5.4——基于非结构化网格的通用CFD求解器，针对非结构性网格模型设计，是用有限元法求解不可压缩流及中度可压缩流流场问题的CFD软件。可应用的范围有紊流、热传、化学反应、混合、旋转流（rotating flow）及震波（shocks）等。在涡轮机及推进系统分析都有相当优秀的结果，并且对模型的快速建立及shocks处的格点调适都有相当好的效果。 Fidap——基于有限元方法的通用CFD求解器，为一专门解决科学及工程上有关流体力学传质及传热等问题的分析软件，是全球第一套使用有限元法于CFD领域的软件，其应用的范围有一般流体的流场、自由表面的问题、紊流、非牛顿流流场、热传、化学反应等等。 FIDAP本身含有完整的前后处理系统及流场数值分析系统。对问题整个研究的程序，数据输入与输出的协调及应用均极有效率。 Polyflow——针对粘弹性流动的专用CFD求解器，用有限元法仿真聚合物加工的CFD软件，主要应用于塑料射出成形机，挤型机和吹瓶机的模具设计。 Mixsim——针对搅拌混合问题的专用CFD软件，是一个专业化的前处理器，可建立搅拌槽及混合槽的几何模型，不需要一般计算流力软件的冗长学习过程。它的图形人机接口和组件数据库，让工程师
流体力学与传热学 流体静力学:研究静止流体中压强分布规律及对固体接触面的作用问题 流体动力学:研究运动流体中各运动参数变化规律，流体与固体作用面的相互作用力的问题 传热学研究内容：研究热传导和热平衡规律的科学上篇：流体力学基础 第一章流体及其主要力学性质 第一节流体的概念 一流体的概述 ⒈流体的概念：流体是液体和气体的统称 ⒉流体的特点:易流动性—在微小剪切力的作用下，都将连续不断的产生变形（区 别于固体的特点） ⑴液体：具有固定的体积；在容器中能够形成一定的自由表面；不可压缩性 ⑵气体;没有固定容积;总是充满所占容器的空间;可压缩性
二连续介质的模型 ⒈连续介质的概念 所谓连续介质即是将实际流体看成是一种假想的,由无限多流体质点所组成的稠密而无间隙的连续介质.而且这种连续介质仍然具有流体的一切基本力学性质. ⒉连续介质模型意义 所谓流体介质的连续性,不仅是指物质的连续不间断,也指一些物理性质的连续不间断性.即反映宏观流体的密度,流速,压力等物理量也必定是空间坐标的连续函数(可用连续函数解决流体力学问题）
第二节流体的性质 一密度—--表征流体质量性质 ⒈密度定义:单位体积内所具有的流体质量 ⑴对于均质流体:ρ=m/v 式中ρ-流体的密度（㎏/m 3） m-流体的质量（㎏） v —流体的体积（m 3) ⑵对于非均质流体:ρ=⒉比体积（比容)：单位质量流体所具有的体积（热力学和 气体动力学概念） ⑴对于均质流体：v=V/m=1/ρ（m 3/㎏) 3．液体的密度在一般情况下，可视为不随温度或压强而变化；但气体的密度则随温度和压强可发生很大的变化。 二流体的压缩性和膨胀性 dv dm v m v =??→?0lim
三维方管内部二次流特征分析 ——基于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 方程下残差收敛曲线
计算流体力学实验报告——热传导方程求解 姓名：梁庆 学号：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 ? <
########学院 计算流体力学与传热学 学号： 专业： 学生姓名： 任课教师：教授 2013年12月
目录 第一章验证显式格式的稳定性 (4) 1.1 概述 (4) 1.2 数学推导 (4) 1.3 问题描述 (4) 1.4 数值模拟 (4) 1.5 结果及分析 (5) 第二章判断肋片可以按一维问题处理的主要依据 (6) 2.1 概述 (6) 2.2 问题描述及算法 (6) 2.3 数值模拟 (7) 2.4 结果及分析 (8) 第三章三层墙导热 (9) 3.1 概述 (9) 3.2 问题描述 (9) 3.3 TDMA算法 (9) 3.4 结果 (10) 第四章一维无源稳态对流扩散问题 (11) 4.1 公式及初值 (11) 4.2 情况一 (11) 4.3 情况二 (12) 4.4 情况三 (13)
第五章用ADI算法计算长方肋内的温度分布 (14) 5.1 问题描述 (14) 5.2 初始参数 (14) 5.3 情况一，一列列扫 (14) 5.4 情况二，一行行扫 (14) 5.5 情况三，采用ADI算法 (15) 5.6 结果分析 (15) 参考文献 (16)
第一章 验证显式格式的稳定性 1.1 概述 将一维非稳态热传导方程用显式格式差分化为代数方程，在求解的迭代过程中必须满足一定的条件，才能使方程收敛且结果正确。此处即验证β≤?。 1.2 数学推导 方程： 22T t T x α??=?? （1） 显式离散格式： 此处时间向前差分，空间中心差分 111 22n n n n n i i i i i T T T T T t x α+-+--+=?? 1112(2)n n n n n i i i i i t T T T T T x α +-+?-=-+? 令β=2 t x α ??则： 111(2)n n n n n i i i i i T T T T T β+-+-=-+ （2） 误差也应该满足上式，故： ()()1()()()2()()i i i i i Ikx Ikx Ik x x Ikx Ik x x n n n n n T e T e T e T e T e ψψβψψψ----?--+?+??-=-+?? ()()()1()12()()()i i i i Ikx Ikx Ik x x Ik x x n n n n T e T e T e T e ψβψβψψ----?-+?+??=-++?? ()()1()12()()i i i Ikx Ikx Ikx n n Ik x Ik x n T e T e e e T e ψβψβψ---+-??=-++ ()()1() 121() n Ik x Ik x n T e e T ψββψ+-??=-++≤ 因此 β≤?。即当β≤? 时方程（2）才会有收敛的解。 1.3 问题描述 在验证过程中同时可模拟一个实际问题，即冬季里墙壁中的温度分布。此时室内壁温设为Tl=30.0℃，室外壁温Tr=-25.0℃，墙壁以11号楼为例，L=1m ，热扩散系数ɑ=alfa=1.33e-6m 2/s 然后分别取β=0.4，n=10和β=0.6，n=10两种情况，看最后的结果是否收敛和正确。 1.4 数值模拟