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 )(??+??= 质点加速度=速度的局部导数+速度的迁移导数。 物理量的质点导数=物理量的局部导数+物理量的对流导数。
《计算流体力学》课程大作业 作业内容：3-4人为小组完成数值模拟，在第8次课上每组进行成果展示，并在课程结束后每组上交一份纸质版报告。 数值模拟实现形式：自编程或者使用任意的开源、商业模型。 成果展示要求：口头讲述和幻灯片结合的方式，每组限时10分钟（8分钟讲述，2分钟提问和讨论）。 报告要求：按照期刊论文的思路和格式进行撰写（包括但不限于如下内容：摘要、绪论\引言、数值模型简介、数值结果分析\讨论、结论、参考文献）。 （以下题目二选一） 题目一：固定单方柱扰流问题 根据文章《Interactions of tandem square cylinders at low Reynolds numbers》中的实验进行数值模拟，完成但不局限于如下工作： （1）根据Fig. 2 中的雷诺数和方柱排列形式，进行相同雷诺数不同间距比情况下的方柱绕流数值模拟，并做出流线图和Fig.2中的结果对比。 （2）根据Fig. 3 中的雷诺数和方柱排列形式，进行相同雷诺数后柱不同转角情况下的方柱绕流数值模拟，并做出流线图和Fig.3中的结果对比。 （3）根据Fig. 12, 13 中的雷诺数和方柱间距比的设置进行数值模拟，作出频率、斯特劳哈尔数、阻力系数随雷诺数变化的折线并与图中对应的折线画在同一坐标系下比较。 （中共有4条折线，对应4种不同的方柱排列形式下的物理参数随雷诺数变化的规律，仅需选取单柱模型和其中一种双柱模型进行数值模拟，共计16个工况）。 题目二：溃坝问题 根据文章《Experimental investigation of dynamic pressure loads during dam break》中的实验进行数值模拟，完成但不局限于如下工作： （1）分别完成二维、三维的溃坝的数值建模，讨论二维、三维模型的区别。 （2）分别将二维、三维溃坝的数值模拟结果和Fig. 7,10中各时刻的自由面形态进行对比，并分别观测溃坝前端水舌的位置随时间的变化，其结果和Fig. 12 种的各试验结果放在同一坐标系下进行对比。 （3）根据实验设置数值观测点，分别观测与实验测点相对应的数值观测点上的水体高度、压力随时间的变化曲线，并和Fig.16, 18,21,30,31,32,33,35中的实验结果进行对比。
高等计算流体力学讲义（2） 第二章 可压缩流动的数值方法 §1. Euler 方程的基本理论 0 概述 在计算流体力学中，传统上，针对可压缩Navier －Stokes 方程的无粘部分和粘性部分分别构造数值方法。其中最为困难和复杂的是无粘部分的离散方法；而粘性项的离散相对简单，一般采用中心差分离散。所以，本章主要研究无粘的Euler 方程的解法。在推广到Navier －Stokes 方程时，只需在Euler 方程的基础上，加上粘性项的离散即可。Euler 方程是一种典型的非线性守恒系统。下面我们将讨论一般的非线性守恒系统以及Euler 方程的一些数学理论，作为研究数值方法的基础。 1非线性守恒系统和Euler 方程 一维一阶非线性守恒系统(守恒律)可写为下列一般形式 =??+??x F t U ，0,>∈t R x （1） 其中U 称为守恒变量，是有m 个分量的列向量，即T m u u u U ),...,(21=。T m f f f F ),...,(21=称为通量函数，是U 的充分光滑的函数，且满足归零条件，即： 0)(lim =→U F U 即通量是对守恒变量的输运，守恒变量为零时，通量也为零。 守恒律的物理意义 设U 的初始值为：0(,0)(),U x U x x =∈R 。如果0()U x 在x ∈R 中有紧支集（即0U 在有限区域以外恒为零），则0(,)()U x t dx U x dx =??R R 。即此时虽然(,)U x t 的分布可以随时 间变化，但其总量保持守恒。 多维守恒律可以写为 )(=++??+??k H j G i F t U （2） 守恒律的空间导数项可以写为散度形式。 守恒系统（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基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
计算流体力学软件CFD在燃烧器设计中的应用探讨[摘要]本文通过对目前燃烧器的现状与技术发展的研究，探讨计算流体力学 软件CFD在燃烧器设计中应用的必要性和可行性，以CFD(计算流体力学)软件为工具，以普通大气式燃烧器为研究对象，采用实验和理论相结合的方法，充分利用现代计算机技术，达到降低燃烧器设计成本和研制费用的目的。 [关键词]燃烧器数值模拟计算流体力学 一、燃烧器的发展现状 1.部分预混式燃烧器的产生及其原理 燃烧的方法被分为扩散式燃烧、部分预混式燃烧和完全预混式燃烧。扩散式燃烧易产生不完全燃烧产物，燃烧温度很低，并未充分利用燃气的能量；而一旦预先混入一部分空气后火焰就会变的清洁，燃烧温度也可以提高，燃烧较充分。完全预混燃烧（无焰燃烧）要求事先按照化学当量比将燃气和空气均匀混合（实际应用中空气系数要大于1），燃烧充分，火焰温度很高，但稳定性较差，易回火。所以民用燃具多采用部分预混式燃烧。 1855年工程师本生发明了一种燃烧器，能从周围大气中吸入一些空气和燃气预混，在燃烧时形成不发光的蓝色火焰，这就是实验室常用的本生灯（单火孔燃烧器）。这种燃烧技术就被称作部分预混式燃烧。 本生灯燃烧所产生的火焰为部分预混层流火焰（俗称本生火焰）。它由内焰，外焰及燃烧区域外围肉眼看不见的高温区组成。火焰一般呈锥体状。燃气—空气的混合气体先在内锥燃烧，中间产物及未燃尽的部分便从锥内向外流出，且混合气体出流的速度与内锥表面火焰向内传播速度相互平衡，此外便形成一个稳定的焰面，呈蓝色。而未燃烧尽的混合气体残余物继续与大气中的空气进行二次混合燃烧，形成火焰外锥。如图1所示，完成燃烧后产生高温co2和水进而在外焰的外侧形成外焰膜（肉眼看不见的高温层）: 图1. 本生灯示意图 如果混合气流是处于层流状态，则外焰面呈较光滑的锥形；如果处于紊流状态，则外焰面产生褶皱，直至产生强烈扰动，气团不断飞散、燃尽。
一.选择题 1.牛顿内摩擦定律适用于（）。 A．任何流体B．牛顿流体C．非牛顿流体 2.液体不具有的性质是（）。 A．易流动性B．压缩性C．抗拉性D．粘滞性 3连续介质假定认为流体（）连续。 A．在宏观上B．在微观上C．分子间D．原子间 4.在国际单位制中流体力学基本量纲不包括（）。 A．时间B．质量C．长度D．力． 5.在静水中取一六面体，作用在该六面体上的力有（） A．切向力、正压力B．正压力C．正压力、重力D．正压力、切向力、重力 6．下述哪些力属于质量力( ) A．惯性力B．粘性力C．弹性力D．表面张力E．重力 7．某点存在真空时，（）（） A．该点的绝对压强为正值B．该点的相对压强为正值c.该点的绝对压强为负值D．该点的相对压强为负值 8.流体静压强的（）。 A．方向与受压面有关B．大小与受压面积有关B．大小与受压面方位无关 9.流体静压强的全微分式为（）。 A．B．C． 10.压强单位为时，采用了哪种表示法（）。 A．应力单位B．大气压倍数C．液柱高度 11.密封容器内液面压强小于大气压强，其任一点的测压管液面（）。A．高于容器内液面B．低于容器内液面C．等于容器内液面 12.流体运动的连续性方程是根据（）原理导出的。 A.动量守恒 B. 质量守恒 C.能量守恒 D. 力的平衡 13. 流线和迹线重合的条件为（）。
A.恒定流 B.非恒定流 C.非恒定均匀流 14.总流伯努利方程适用于（）。 A.恒定流 B.非恒定流 C.可压缩流体 15. 总水头线与测压管水头线的基本规律是：（）、（） A.总水头线总是沿程下降的。 B.总水头线总是在测压管水头线的上方。 C.测压管水头线沿程可升可降。 D.测压管水头线总是沿程下降的。 16 管道中液体的雷诺数与（）无关。 A. 温度 B. 管径 C. 流速 D. 管长 17.. 某圆管直径d=30mm，其中液体平均流速为20cm/s。液体粘滞系数为0.0114cm3/s，则此管中液体流态为（）。 A. 层流 B. 层流向紊流过渡 C.紊流 18.等直径圆管中紊流的过流断面流速分布是（）A呈抛物线分布B. 呈对数线分布 C.呈椭圆曲线分布 D. 呈双曲线分布 19.等直径圆管中的层流，其过流断面平均流速是圆管中最大流速的（） A 1.0倍B.1/3倍C. 1/4倍D. 1/2倍 20.圆管中的层流的沿程损失与管中平均流速的（）成正比. A. 一次方 B. 二次方 C. 三次方 D. 四次方 21..圆管的水力半径是( ) A. d/2 B. d/3 C. d/4 D. d/5. 22谢才公式中谢才系数的单位是（）A. 无量纲B. C. D. . 23. 判断层流和紊流的临界雷诺数是（） A.上临界雷诺数 B.下临界雷诺数 C.上下临界雷诺数代数平均 D.上下临界雷诺数几何平均 24.. 对于管道无压流，当充满度分别为（）时，其流量和速度分别达到最大。 A. 0.5, 0.5 B. 0.95, 0.81 C. 0.81, 081 D. 1.0, 1.0 25.对于a, b, c三种水面线，下列哪些说法是错误（）（） A．所有a、c型曲线都是壅水曲线，即，水深沿程增大。B．所有b型曲线都是壅水曲线，即，水深沿程增大。C．所有a、c型曲线都是降水曲线，即，水深沿程减小。C．所有b型曲线都是降水曲线，即，水深沿程减
2014级西安理工大学计算流体力学作业 1.写出通用方程，并说明其如何代表各类守恒定律。 由守恒型对流-扩散方程： ()()() div U div T grad S t φφρφρφφ?+=+? 其中φ为通用变量；T φ为广义扩散系数；S φ为广义原项。 若令1;1;0T S φφφ===时，则得到质量守恒方程（mass conservation equation ） ()()()() 0u v w t x y z ρρρρ????+++=???? 若令;i u φ=时，则得动量守恒方程（momentum conservation equation ） 以x 方向为例分析，设;u P u S S x φφ?==- ?，通用方程可化为： ()()()()(2)u uu vu wu P u divU t x y z x x x ρρρρλη???????+++=-++??????? z v u u w F y x y z z x ηηρ???????????? ??+++++?? ? ????????????????? 同理可证明y 、z 方向的动量守恒方程式 若令;;T p T T S S C φφλ φ===时，则得到能量守恒方程(energy conservation equation) ()()() ()h h div Uh div U div gradT S t ρρρλφ?+=-+++? ()()()T p h div Uh div gradT S t C ρλ ρ?+=+? 证毕 2.用控制体积法离散 0)(=+++s dx dT k dx d dx dT u dt dT ，要求对S 线性化，据你的理解，谈谈网格如何划分？交界面传热系数何如何计算？边界条件如何处理？ 根据守恒型对流-扩散方程： ()()()u T S t x x x ρφρ?φ ????' +=+????，对一维模型 进行分析，则有： 0)(=+++s dx dT k dx d dx dT u dt dT
计算流体力学教案 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 提出问题 [问题描述] Sod 激波管问题是典型的一类Riemann 问题。如图所示，一管道左侧为高温高压气体，右侧为低温低压气体，中间用薄膜隔开。t=0 时刻，突然撤去薄膜，试分析其他的运动。 Sod 模型问题：在一维激波管的左侧初始分布为：0 ,1 ,1111===u p ρ，右侧分布为：0 ,1.0 ,125.0222===u p ρ，两种状态之间有一隔膜位于5.0=x 处。隔膜突然去掉，试给出在14.0=t 时刻Euler 方程的准确解，并给出在区间10≤≤x 这一时刻u p , ,ρ的分布图。 2 一维Euler 方程组 分析可知，一维激波管流体流动符合一维Euler 方程，具体方程如下： 矢量方程： 0U f t x ??+=?? (0.1) 分量方程： 连续性方程、动量方程和能量方程分别是： 2 22,,p u ρ
() ()()()2 000u t x u u p t x x u E p E t x ρρρρ???+ =?????????++=? ??????+?????+ =????? (0.2) 其中 22v u E c T ρ?? =+ ?? ? 对于完全气体，在量纲为一的形式下，状态方程为： ()2 p T Ma ργ∞ = (0.3) 在量纲为一的定义下，定容热容v c 为： () 21 1v c Ma γγ∞= - (0.4) 联立(1.2)，(1.3)，(1.4)消去温度T 和定容比热v c ，得到气体压力公式为： ()2112p E u γρ??=-- ??? (0.5) 上式中γ为气体常数，对于理想气体4.1=γ。 3 Euler 方程组的离散 3.1 Jacibian 矩阵特征值的分裂 Jacibian 矩阵A 的三个特征值分别是123;;u u c u c λλλ==+=-，依据如下算法将其分裂成正负特征值： () 12 222 k k k λλελ±±+= (0.6) 3.2 流通矢量的分裂 这里对流通矢量的分裂选用Steger-Warming 分裂法，分裂后的流通矢量为 ()()()()()()()12312322232121212122f u u c u c u u c u c w γλλλργλλλγλλγλ?? ? -++ ?=-+-++ ? ? ? -+-+++ ??? ＋＋＋＋＋＋＋ ＋＋ ＋＋ (0.7)
计算流体力学软件Fluent在烟气脱硫中的应用 ０引言 污染最为有效的方法之一，而石灰石—石膏湿烟气脱硫是目前能大规模控制燃煤造成ＳＯ ２ 法脱硫技术以其脱硫效率高、吸收剂来源丰富、成本低廉、技术成熟和运行可靠等优点获得广泛应用．从气液两相流体力学和化学反应动力学的观点看，脱硫吸收塔内流体流动的目的是强化气液两相的混合和质量传递、延长气液两相在塔内的接触时间、增大气液两相的接触面积并尽量减小吸收塔的阻力．合理的塔内流场分布对提高脱硫效率、降低脱硫投资和运行成本都具有重要意义． 目前，国内外对烟气脱硫吸收塔进行大量研究，主要采用实验方法，如研究塔的阻力特性、液滴运动速度沿塔高变化和ＴＣＡ塔内温度场分布等，这些研究对指导工业应用具有重要意义，但其结果往往只针对特定的设备或结构，具有较大的局限性．随着计算机技术的迅速发展，计算流体力学（ＣｏｍｐｕｔａｔｉｏｎａｌＦｌｕｉｄＤｙｎａｍｉｃ，ＣＦＤ）已成为研究三维流动的重要方法：周山明等［４］利用ＦＬＵＥＮＴ计算空塔和喷淋状态下的塔热态流场，结果表明脱硫吸收塔入口处流场变化最剧烈、压降损失最大，并根据计算结果改造来流烟道；孙克勤等采用混合网格和随机颗粒生成模型对烟气脱硫吸收塔的热态流场进行数值模拟；郭瑞堂等采用ＦＬＵＥＮＴ结合非稳态反应传质－反应理论对湿法脱硫液柱冲的吸收进行数值模拟． 击塔内的流场和ＳＯ ２ 本文尝试应用ＦＬＵＥＮＴ对某脱硫吸收塔内烟气脱硫过程进行初步数值模拟，通过对内部流场进行分析验证本文模拟的合理性，进而对脱硫过程中脱硫吸收塔内是否存在湿壁现象进行深入分析研究． １基于ＲＡＮＳ求解器的ＣＦＤ数值模拟 方法 １．１控制方程 时均的不可压缩连续性方程和Ｎ Ｓ方程 （ＲＡＮＳ方程）如下： １．２湍流模型和多相流模型
《计算流体力学》课程大作业 ——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟 张伊哲 航博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.模型概述 在汽车外部建立一个较大的长方体几何空间，长度约为30m，宽度和高度约为5m，在空间内部挖出汽车形状的空腔，汽车尺寸参照本田CRV为4550mm*1820mm*1685mm。由于汽车向前开进，气体从车头流向车尾，因此将汽车前方空间设为气体入口，后方空间设为气体出口，模拟气体在车外的流动。另外为了节省计算成本将整个模型按1:100的比例缩小，考虑到模型和流体均是对称的，因此仅画出几何模型的一半区域，建立对称面以考虑生成包含理想气体的流体域。在Catia中建立的模型如图1.1所示。 图1.1几何模型 2.利用ICEM CFD进行网格划分 a)导入有Catia生成的stp格式的模型； b)模型修复，删除多余的点、线、面，允许公差设为0.1； c)生成体，由于本模型仅为流体区域，因此将全部区域划分为一个体，选取方法可以 使用整体模型选取； d)为了后面的设置边界方便，因此将具有相同特性的面设为一个part,共设置了in， out，FreeWalls，Symmetry和Body； e)网格划分，设置Max element=2，共划分了1333817个单元，有225390个节点； f)网格输出，设置求解器为ANSYS CFX，输出cfx5文件。 3.利用ANSYS CFX求解 a)生成域，物质选定Air Ideal Gas，参考压强设为1atm，浮力选项为无浮力模型，
域运动选项为静止，网格变形为无；流体模型设定中的热量传输设定为Isothermal，流体温度设定为288k，湍流模型设定为Shear Stress Transport模型，壁面函数 选择Automatic。 b)入口边界设定，类型为Inlet，位置选定在in，质量与栋梁选定Normal Speed，设 定为15m/s，湍流模型设定类型为Intensity and Length Scale=0.05，Eddy Len.Scale=0.1m。 c)出口边界设定，边界类型为Outlet，位置选out。质量与动量选项为Static Pressure，相对压强为0pa。 d)壁面边界设定，边界类型为Wall，位置选在FreeWalls。壁面边界详细信息中指定 WallInfluence On Flow为Free Slip。 e)对称边界设定，边界类型为Symmetry，位置选在Symmetry。 f)汽车外壁面设定，边界类型为Wall，位置设在Body，壁面详细信息选项中指定Wall Influence On Flow为No Slip，即汽车壁面为无滑移壁面。 g)初始条件设定，初始速度分量设为U方向为15m/s，其他两个方向的速度为零。 h)求解设置，残差类型选为RMS，残差目标设定为1e-5，当求解达到此目标时，求解 自动终止。求解之前的模型如图3.1所示。 图3.1求解之前的模型 4.结果后处理 从图4.1中可以看出计算收敛。
流体力学-大作业
————————————————————————————————作者: ————————————————————————————————日期：
一．选择题 1．牛顿内摩擦定律适用于()。 A.任何流体B．牛顿流体 C.非牛顿流体 2．液体不具有的性质是()。 A．易流动性Ｂ.压缩性Ｃ.抗拉性 D.粘滞性 3连续介质假定认为流体（)连续。 A．在宏观上 B.在微观上 C.分子间Ｄ．原子间 4.在国际单位制中流体力学基本量纲不包括（）。 A．时间 B.质量 C.长度D．力． 5.在静水中取一六面体，作用在该六面体上的力有（) Ａ．切向力、正压力Ｂ．正压力C．正压力、重力D．正压力、切向力、重力 ６. 下述哪些力属于质量力（) A．惯性力Ｂ.粘性力 C.弹性力Ｄ.表面张力Ｅ．重力 7.某点存在真空时,( ）(） A.该点的绝对压强为正值 B.该点的相对压强为正值c.该点的绝对压强为负值 D.该点的相对压强为负值 ８．流体静压强的( )。 A.方向与受压面有关Ｂ.大小与受压面积有关B．大小与受压面方位无关 9．流体静压强的全微分式为（)。 A. B. C． １0．压强单位为时，采用了哪种表示法（)。 A．应力单位Ｂ.大气压倍数C．液柱高度 11.密封容器内液面压强小于大气压强，其任一点的测压管液面( ）。 A.高于容器内液面B．低于容器内液面 C.等于容器内液面 12．流体运动的连续性方程是根据( ）原理导出的。 A.动量守恒 B. 质量守恒 C.能量守恒 D. 力的平衡 13．流线和迹线重合的条件为（)。
A．恒定流B.非恒定流C.非恒定均匀流 １４．总流伯努利方程适用于（）。 A.恒定流 B.非恒定流Ｃ.可压缩流体 15. 总水头线与测压管水头线的基本规律是：( ）、( ） A.总水头线总是沿程下降的。B．总水头线总是在测压管水头线的上方。 C.测压管水头线沿程可升可降。 D.测压管水头线总是沿程下降的。 １6 管道中液体的雷诺数与（)无关。 A.温度B．管径Ｃ. 流速D．管长 １7．. 某圆管直径d=3０mｍ，其中液体平均流速为20ｃｍ/s。液体粘滞系数为0.0１14cm3/s，则此管中液体流态为( )。 A. 层流 B. 层流向紊流过渡C．紊流 1８．等直径圆管中紊流的过流断面流速分布是（)A呈抛物线分布B．呈对数线分布 C.呈椭圆曲线分布D．呈双曲线分布19．等直径圆管中的层流,其过流断面平均流速是圆管中最大流速的（） A １．０倍 B.1/３倍Ｃ.1/４倍D. 1/2倍 20.圆管中的层流的沿程损失与管中平均流速的（）成正比. A. 一次方 B.二次方 C. 三次方Ｄ. 四次方 ２1..圆管的水力半径是（) A. d/2Ｂ.d/3 C. d/４D. d／5. 2２谢才公式中谢才系数的单位是（）A.无量纲B．Ｃ.D．. 23．判断层流和紊流的临界雷诺数是（） Ａ.上临界雷诺数 B.下临界雷诺数 C.上下临界雷诺数代数平均 D.上下临界雷诺数几何平均 24..对于管道无压流,当充满度分别为( )时，其流量和速度分别达到最大。Ａ．0.５,0.５B．0.95,0．８1 C.0.８1, ０8１ D. １.0，1.０ 25.对于ａ, ｂ，c三种水面线，下列哪些说法是错误( ）（) A．所有a、ｃ型曲线都是壅水曲线，即，水深沿程增大。B．所有
2012～2013学年第1学期 12级研究生《计算流体力学》结课作业 适用专业：供热供燃气通风及空调工程 一、结合某一具体学科，阐述纯理论方法、实验方法及数值方法在科学研究中的各自优缺点，在此基础上论述数值模拟方法的发展前景。（不少于4千字）。 流体力学是力学的一个重要分支, 是研究流体(液体和气体)的力学运动规律及其应用的学科, 主要研究在各种力的作用下，流体本身的静止状态和运动状态特征，以及流体和相邻固体界面有相对运动时的相互作用和流动规律。在人们的生活和生产活动中随时随地都可遇到流体，流体力学与人类的日常生活和生产事业密切相关。按其研究内容的侧重点不同，分为理论流体力学和工程流体力学。其中理论流体力学主要采用严密的数学推理方法，力求准确性和严密性，工程流体力学侧重于解决工程实际中出现的问题，而不追求数学上的严密性。当然由于流体力学研究的复杂性，在一定程度上，两种方法都必须借助于实验研究，得出经验或半经验的公式。 在实际工程的诸多领域流体力学都起着十分重要的作用。如气象、水利的研究，船舶、飞行器、叶轮机械和核电站的设计及其运行，可燃气体或炸药的爆炸，都广泛地用到流体力学知识。许多现代科学技术所关心的问题既受流体力学的指导，同时也促进了流体力学自身的不断发展。1950年后，计算机的发展给予流体力学以极大的推动作用。 目前，解决流体力学问题的方法主要有实验方法、理论分析方法和数值方法三种。 实验方法 同物理学、化学等学科一样，流体力学的研究离不开实验，尤其是对新的流体运动现象的研究。实验能显示运动特点及其主要趋势，有助于形成概念，检验理论的正确性。二百年来流体力学发展史中每一项重大进展都离不开实验。流体力学实验研究方法有实物实验、比拟研究和模型研究三类：实物实验是用仪器实测原型系统的流动参数，适用于较小的原型；比拟实验是利用电场和磁场来模拟流场，实施起来限制条件较多；模型研究是实验流体力学最常用的研究方法。 实验研究的一般过程是：在相似理论的指导下建立实验模型，用流体测量技术测量流动参数，处理和分析实验数据。建立实验模型要求模型与原型满足相似理论，即满足两个流场
第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、高效传质与分离和精细化工方面的研究。
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.简介 计算流体动力学（Computational Fluid Dynamics，简称CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。其基本思想为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看作是在流动基本方程（质量守恒方程、动量守恒方程、能量守恒方程）控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量（速度、压力、温度、浓度等）的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理星，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 2.计算流体动学的特点: ①流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和边界条件复杂，很难求得解析解，而用CFD方法则有可能找出满足工程需要的数值解。 ②可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较。 ③它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件。 ④数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差。 ⑤它不像物理模型实验一开始就能给出流动现象并定性地描述，往往需要由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模型进行验证。