9.粘性流体力学基本方程组

流体力学大作业

《计算流体力学》课程大作业 作业内容：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 02 00Re L V L V L V μρμρ= = ) (/)(00003 000020T T C L V L V T T C V Ec w p w p - =-= ρρ 热传递中流体压缩性的影响，也就是推进功与对流热之比。00 0Pr K C p μ= 表示流体的物性的影响，表征温度场和速度场的相似程度。边界层特征厚度dy u u h e e ?- =0 * )1(ρρδ 边界层的存在而使自由流流线向外推移的距离。 θ δ* =H 能够反映速度剖面的形状，H 值越小， 剖面越饱满。动量积分方程：不可压流二维 f e w e e C u dx du u H dt d ==++2)2(ρτθθ /2 普朗特方程的导出，相似解的概念，布拉休斯解的主要结论 ?????????????+??+??-=??+??+????+??+??-=??+??+??=??+ ??)(1)(1022222222y v x v y p y v v x v u t v y u x u x p y u v x u u t u y v x u νρνρ 将方程无量纲化： ./,/,/,/*2***L tU t u p p U u u L x x ====ρ ν/Re UL =，Re /1*≈δ ,/,/,,**L L y U u v L y u v δδ=?==?= 分析：当Re 趋于很大时，**y p ??是大量，则**y p ??＝0，根据量纲分析，去掉小量化为有量纲形式则可得到普朗特边界层方程： ???? ?? ??? =????+??-=??+??+??=??+??01022y p y u x p y u v x u u t u y v x u υρ 相似解的概念：对不同x 截面上的速度剖面u(x,y)都可以通过调整速度u 和坐标y 的尺度因子,使他们重合在一起。外部势流速度Ue(x)作为u 的尺度因子，g(x)作为坐标y 的尺度因子。则无量纲坐标)(x g y ，无量纲速度)(x u u e ，则 对所有不同的x 截面其速度剖面的形状将会相 同。即= )(])(,[111x u x g y x u e ) (] ) (,[222x u x g y x u e 布拉修斯解（零攻角沿平板流动的解）的主要结论： x x Re 721.1* =δx x Re 664.0=θ 591.2/*==θδH 壁面切应力为： x y w U y u Re 1332.0)(2 0∞ ==??=ρμτ 壁面摩擦系数为：x w f u C Re 1664.022 ==∞ρτ 平均为：l l f Df dx C l C Re 1328.110? == 湍流的基本概念及主要特征，湍流脉动与分子随机运动之间的差别湍流是随机的，非定常的，三维的有旋流动，随机背后还存在拟序结构。特征：随机脉动耗散性，有涡性（大涡套小涡）。 湍流脉动：不断成长、分裂和消失的湍流微团；漩涡的裂变造成能量的传递；漩涡运动与边界条件有密切关系，漩涡的最小尺度必大于分子的自由程。分子随机运动：是稳定的个体；碰撞时发生能量交换；平均自由程λ与平均速度 和边界条件无关。层流稳定性的基本思想：在临界雷诺数以下时，流动本身使得流体质点在外力的作用下具有一定的稳定性，能抵抗微弱的扰动并使之消失，因而能保持层流；当雷诺数超过临界值后，流动无法保持稳定，只要存在微弱的扰动便会迅速发展，并逐渐过渡到湍流。平板边界层稳定性研究得到的主要结果：1.雷诺数达到临界雷诺数时流动开始不稳定，成为不稳定点，而转捩点则对应与更高的雷诺数。2.导致不稳定扰动最小波长 δ δλ65.17min ≈=*，可见不稳定波是一种 波长很长的扰动波，约为边界层厚度的6倍。3. 不稳定扰动波传播速度远小于边界层外部势流速度，其最大的扰动波传播速度 4.0/=∞U c r 。当雷诺数相当大时，中性稳定线的上下两股趋于水平轴。判别转捩的试验方法： 升华法（主要依据：湍流的剪切应力大小）热膜法（主要依据：层流和湍流边界层内 气流脉动和换热能力的差别）液晶法（主要依 据：湍流传热和层流传热能力之间的差异）湍流的两种统计理论：1. 湍流平均量的半经验分 析（做法：主要研究各个参数的平均量以及它们之间的相互关系，如平均速度，压力，附面层厚度等。2. 湍流相关函数的统计理论分析（做法；将流体视为连续介质，将各物理量如：流速，压力，温度等脉动值视为连续的随机函数， 并通过各脉动值的相关函数和谱函数来描述湍流结构。）耗散涡、含能涡的尺度耗散涡为小尺 度涡，它的尺度受粘性限制，但必大于分子自由行程。控制小尺度运动的参数包括单位质量的能量消耗量ε和运动粘性系数ν。因此，由 量纲分析，小涡各项尺度为：长度尺度 4/13)(ενη=时间尺度2/1)(εντ=速度尺度4/1)(νε=v 耗散雷诺数 1Re →=νη v d 可知：小尺度涡体的湍流 脉动是粘性主宰的耗散流动，因此这一尺度的 涡叫耗散涡。含能涡为大尺度涡，在各向同性湍流中，可以认为大尺度涡体由它所包含的湍动总能量k ，以及向小尺度传递的能量ε决定。 长度尺度ε2/3k l =时间尺度εk t =速度尺度k u =积分尺度雷诺数1Re →>>=ν ul d 可知在含能尺度范围 内，惯性主宰湍流运动，因此含能尺度范围又 称惯性区。均匀湍流：统计上任何湍流的性质与空间位置无关，或者说，任何湍动量的平均 值及它们的空间导数，在坐标做任何位移下不 变。特征：不论哪个区域，湍流的随机特性是相同的，理论上说，这种湍流在无界的流场中 才可能存在。各向同性湍流：任何统计平均量与方向无关，或者说，任何湍动量在各个方向 都一样，不存在任何特殊地位的方向。任何统计平均湍动量与参考坐标轴的位移、旋转和反 射无关。特征：各向同性湍流，必然是均匀湍 流，因为湍流的任何不均匀性都会带来特殊的方向性。在实际中，只存在局部各向同性湍流 和近似各向同性湍流。各向同性下，雷诺应力 由9个量减为3个量。 了解时均动能方程、湍动能方程中各项的物理意义和特点，及能量平衡时均动能方程： 流体微团内平均动能变化率；外力的作功；平均压 力梯度所作的功； 雷诺应力所作功的扩散；雷诺应力所作的变形功；时均流粘性应力所作功 的扩散；时均流动粘性的耗散，即粘性应力的 变形功。 湍动能方程：

流体力学 大作业

一.选择题 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型曲线都是降水曲线，即，水深沿程减

计算流体力学大作业

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)

流体力学-大作业

一.选择题 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型曲线都是降水曲线，即，水深沿程减

湖北工业大学流体力学作业答案版

1、已知油的重度为7800N/m 3，求它的密度和比重。试求0.2m 3此种油的质量和重量各为多少 已已知知：：γ＝7800N/m 3；V ＝0.2m 3。 解析：(1) 油的密度为 3kg/m 79581 .97800===g γ ρ； 油的比重为 795.01000 795O H 2===ρρ S (2) 0.2m 3的油的质量和重量分别为 kg 1592.0795=?==V M ρ N 15602.07800=?==V G γ 2、体积为5m 3的水在温度不变的条件下，压力从1大气压增加到5大气压，体积减小了1L ，求水的体积压缩系数和弹性系数值。 已已知知：：V ＝5.0m 3，p 1＝×105Pa ，p 2＝×105Pa ，ΔV ＝1L 。 解析：由(1－9)和(1－10)式，得水的体积压缩系数及弹性系数值分别为 N /m 100.510 )0.10.5(0.5100.1d d 21053 p --?=?-??=-=p V V β 2910p m /N 100.2100.511?=?== -βE 3、一块可动平板与另一块不动平板同时浸在某种液体中，它们之间的距离为0.5mm ，可动板若以 0.25m/s 的速度移动，假设板间流体中的速度分布是线性的，为了维持这个速度需要单位面积上的作用力为2N/m 2，求这两块平板间流体的动力黏度μ。 解：板间流体的速度梯度可计算为 13du u 0.25500s dy y 0.510 --===? 由牛顿切应力定律d d u y τμ=，可得两块平板间流体的动力黏度为 3d 410Pa s d y u τμ-= =?? 4、流场的速度分布为2265375x y z u xy xt u y u xy zt =+=-=-，，，求流体在点（2，1，4）和时间t =3s 时的速度、加速度。

计算流体力学课程大作业

《计算流体力学》课程大作业 ——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟 张伊哲 航博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】的经典结论进行对比，评述本

2018《粘性流体力学》复习提纲

粘流复习大纲 1 卡门涡街、阻力危机和马格努斯效应等基本概念 2 流线、迹线、时间线和烟线的概念和物理含义（坐标系的影响） 3 涡量输运方程各项的物理意义，涡动力学亥姆霍兹三定理的内容、涵义及成立的条件，涡量以及流动‘有旋’或‘无旋’的定义，能判断简单流动是否有旋 4 推导N-S方程时所用到的Stokes三假设的内容 5 一些无量纲参数的定义和物理意义（Re, Ec, Pr），及其与速度边界层和温度边界层特性之间的内在关联，壁面恢复温度的概念 6 库特剪切流、突然起动平板流解的主要结论，库特剪切流的速度分布、温度分布，能够运用能量方程来分析库特剪切流的能量平衡 7 边界层的各种特征厚度及形状因子，边界层动量积分方程和计算，基于控制体积分方法分析边界层的流动 8 普朗特边界层理论，边界层微分方程的导出及主要结论，相似解的概念，布拉休斯解的主要结论 9 湍流的基本概念及主要特征（四个），湍流脉动与分子随机运动之间的差别 10 层流稳定性的基本思想，瑞利定理和费约托夫定理，中性稳定线，平板边界层稳定性研究得到的主要结果 11 猝发现象，能叙述边界层转捩的主要过程（典型流动现象） 12 影响转捩过程的主要因素以及控制边界层转捩的主要方法、判别转捩的试验方法 13 湍流的两种统计理论，能谱分析方法的主要结论，半经验理论中流场参数平均的三种方法 14 耗散涡、含能涡的尺度、特征与主要作用，及其特征尺度的描述参数 15 均匀剪切湍流、均匀湍流、各向同性湍流和局部平衡湍流的概念、特征和典型示例 16 不可压下的时均连续方程、动量方程，以及由此而来的方程组封闭性问题，雷诺应力的概念和物理意义 17时均动能方程、湍动能方程中各项的物理意义和特点，及能量平衡 18 目前，湍流的数值模拟的3个层次及各自的特点 19 湍流模型建立的基本法则和各项模化的一般方法 20 湍流模型的分类，涡粘模型的基本假设（布希内斯克的涡粘假定），普朗特混合长度理论，科尔莫果洛夫－普朗特理论，能量方程模型、k-e模型、k-w模型的湍流粘性系数的求法 21 湍流模型近壁区处理的几种方法及对计算网格的要求 22 ASM模型的优点和得出的基本假设 23湍流边界层的宏观结构和速度分布特性 湍流边界层内的湍动特性及能量平衡（包括时均动能和湍动能）

流体力学大作业

流体力学-大作业

————————————————————————————————作者: ————————————————————————————————日期：

一．选择题 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．所有

流体力学标准化作业答案

流体力学标准化作业（三） ——流体动力学 本次作业知识点总结 1.描述流体运动的两种方法 （1）拉格朗日法；（2）欧拉法。 2.流体流动的加速度、质点导数 流场的速度分布与空间坐标(,,)x y z 和时间t 有关，即 流体质点的加速度等于速度对时间的变化率，即 投影式为 或 ()du u a u u dt t ?==+??? 在欧拉法中质点的加速度du dt 由两部分组成， u t ??为固定空间点，由时间变化引起的加速度，称为当地加速度或时变加速度，由流场的不恒定性引起。()u u ??为同一时刻，由流场的空间位置变化引起的加速度，称为迁移加速度或位变加速度，由流场的不均匀性引起。 欧拉法描述流体运动，质点的物理量不论矢量还是标 量，对时间的变化率称为该物理量的质点导数或随体导数。例如不可压缩流体，密度的随体导数 3.流体流动的分类 （1）恒定流和非恒定流 （2）一维、二维和三维流动 （3）均匀流和非均匀流 4.流体流动的基本概念 （1）流线和迹线 流线微分方程 迹线微分方程 （2）流管、流束与总流 （3）过流断面、流量及断面平均流速

体积流量 3(/)A Q udA m s =? 质量流量 (/)m A Q udA kg s ρ=? 断面平均流速 A udA Q v A A ==? （4）渐变流与急变流 5. 连续性方程 （1）不可压缩流体连续性微分方程 （2）元流的连续性方程 （3）总流的连续性方程 6. 运动微分方程 （1）理想流体的运动微分方程（欧拉运动微分方程） 矢量表示式 （2）粘性流体运动微分方程（N-S 方程） 矢量表示式 21()u f p u u u t νρ?+?+?=+??? 7.理想流体的伯努利方 （1）理想流体元流的伯努利方程 （2）理想流体总流的伯努利方程 8.实际流体的伯努利方程 （1）实际流体元流的伯努利方程 （2）实际流体总流的伯努利方程 10.恒定总流的动量方程 投影分量形式

流体力学标准化作业答案第三章

流体力学标准化作业（三） ——流体动力学 本次作业知识点总结 1.描述流体运动的两种方法 （1）拉格朗日法；（2）欧拉法。 2.流体流动的加速度、质点导数 流场的速度分布与空间坐标(,,)x y z 和时间t 有关，即 (,,,)u u x y z t = 流体质点的加速度等于速度对时间的变化率，即 Du u u dx u dy u dz a Dt t x dt y dt z dt ????= =+++ ???? 投影式为 x x x x x x y z y y y y y x y z z z z z z x y z u u u u a u u u t x y z u u u u a u u u t x y z u u u u a u u u t x y z ?????=+++?????? ????? =+++???????????=+++?????? 或 ()du u a u u dt t ?==+??? 在欧拉法中质点的加速度du dt 由两部分组成， u t ??为固定空间点，由时间变化 引起的加速度，称为当地加速度或时变加速度，由流场的不恒定性引起。 ()u u ??v v 为同一时刻，由流场的空间位置变化引起的加速度，称为迁移加速度或位变加速度， 由流场的不均匀性引起。 欧拉法描述流体运动，质点的物理量不论矢量还是标量，对时间的变化率称为该物理量的质点导数或随体导数。例如不可压缩流体，密度的随体导数 D D u t t ρρ ρ?=+???（）

3.流体流动的分类 （1）恒定流和非恒定流 （2）一维、二维和三维流动 （3）均匀流和非均匀流 4.流体流动的基本概念 （1）流线和迹线 流线微分方程 x y z dx dy dz u u u == 迹线微分方程 x y z dx dy dz dt u u u === （2）流管、流束与总流 （3）过流断面、流量及断面平均流速 体积流量 3(/)A Q udA m s =? 质量流量 (/)m A Q udA kg s ρ=? 断面平均流速 A udA Q v A A == ? （4）渐变流与急变流 5. 连续性方程 （1）不可压缩流体连续性微分方程 0y x z u u u x y z ???++=??? （2）元流的连续性方程 12 1122 dQ dQ u dA u dA =?? =? （3）总流的连续性方程 1122u dA u dA = 6. 运动微分方程 （1）理想流体的运动微分方程（欧拉运动微分方程）

CFD 基 础(流体力学)解析

第1章 CFD 基 础 计算流体动力学(computational fluid dynamics ，CFD)是流体力学的一个分支，它通过计算机模拟获得某种流体在特定条件下的有关信息，实现了用计算机代替试验装置完成“计算试验”，为工程技术人员提供了实际工况模拟仿真的操作平台，已广泛应用于航空航天、 热能动力、土木水利、汽车工程、铁道、船舶工业、化学工程、流体机械、环境工程等 领域。 本章介绍CFD 一些重要的基础知识，帮助读者熟悉CFD 的基本理论和基本概念，为计算时设置边界条件、对计算结果进行分析与整理提供参考。 1.1 流体力学的基本概念 1.1.1 流体的连续介质模型 流体质点(fluid particle)：几何尺寸同流动空间相比是极小量，又含有大量分子的微 元体。 连续介质(continuum/continuous medium)：质点连续地充满所占空间的流体或固体。 连续介质模型(continuum/continuous medium model)：把流体视为没有间隙地充满它所占据的整个空间的一种连续介质，且其所有的物理量都是空间坐标和时间的连续函数的一种假设模型：u =u (t ,x ,y ,z )。 1.1.2 流体的性质 1. 惯性 惯性(fluid inertia)指流体不受外力作用时，保持其原有运动状态的属性。惯性与质量有关，质量越大，惯性就越大。单位体积流体的质量称为密度(density)，以r 表示，单位为kg/m 3。对于均质流体，设其体积为V ，质量为m ，则其密度为 m V ρ= (1-1) 对于非均质流体，密度随点而异。若取包含某点在内的体积V ?，其中质量m ?，则该点密度需要用极限方式表示，即 0lim V m V ρ?→?=? (1-2) 2. 压缩性 作用在流体上的压力变化可引起流体的体积变化或密度变化，这一现象称为流体的可压缩性。压缩性(compressibility)可用体积压缩率k 来量度

计算流体力学大作业报告(翼型空气动力分析)

课程综合作业课程名称：计算流体力学 专业班级：研究方向： 学生姓名：学号： 完成日期：

计算流体力学课程综合报告 1.简介 计算流体动力学（Computational Fluid Dynamics，简称CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。其基本思想为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看作是在流动基本方程（质量守恒方程、动量守恒方程、能量守恒方程）控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量（速度、压力、温度、浓度等）的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理星，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 2.计算流体动学的特点: ①流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和边界条件复杂，很难求得解析解，而用CFD方法则有可能找出满足工程需要的数值解。 ②可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较。 ③它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件。 ④数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差。 ⑤它不像物理模型实验一开始就能给出流动现象并定性地描述，往往需要由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模型进行验证。

流体力学基础

第 二章 流 体 力 学 基础 ξ 1 流体的主要性质 1.1 流体的主要物理性质 1、 流体的流动性 ——流体的易变形性 流体的基本属性 流体的力学定义：不能抵抗任何剪切力作用下的剪切变形趋势的物质 2、流体的连续性 ? 连续介质模型 ? 流体质点：含有大量分子的流体微团 3、质量和重力特性 – 密度（ kg/ m 3）、 比容： （ m 3/kg ）、比重：无量纲量 – 浓度：质量浓度（kg/m3）、摩尔浓度(mol/m3) 4、流体的可压缩性与热膨胀性 ? 等温压缩系数（ m 2/N ）： ? 热膨胀系数（1/K ） ? 液体的压缩系数和膨胀系数都很小 ? 压强和温度的变化对气体密度和体积的变化影响较大 ? 可压缩流体与不可压缩流体 5. 流体的粘滞性 (1) 粘滞性： ? ——由于相对运动而产生内摩擦力以反抗自身的相对运动的性质。 ? 一切流体都具有粘性，这是流体固有的特性。 ? 粘性的物理本质：分子间引力、分子的热运动,动量交换 (2) 牛顿粘性定律 ? 粘性力（内摩擦力）： ? ? 粘性切应力： （3） 粘性系数 动力粘滞系数或动力粘度(μ)（Pa·S ） 运动粘度 （m2/s ） : 理想流体（无粘性流体,μ＝0）与实际流体（粘性流体μ≠0） )/(2m N dy u d μτ-=) (N A dy u d F μ-=ρμν=

1.3 作用于流体上的力 1、 质量力 表征：单位质量力：----单位质量流体所受到的质量力 当质量力仅为重力：F bx =0，F by =0，F bz = -g 2 表面力 表征： 切向应力（剪切应力）：τ =T/（N/m 2） 法向应力（压应力）：p=P/A （N/m 2） ξ 2 流体运动的微分方程 1、 流体运动的描述 （1）描述流体运动的数学方法——拉格朗日法和欧拉法 拉格朗日法—— – 着眼于流体质点，设法描述每个流体质点自始至终的运动过程 – 描述流体质点的物理量表示为：f ＝f （a ，b ，c ，τ） – 欧拉法——空间描述 – 着眼于流体质点，设法描述每个流体质点自始至终的运动过程 – 描述流体物理量表示为：f ＝F （x ，y ，z ，τ） （2） 迹线与流线 迹线： – 同一流体质点在连续时间内的运动轨迹线 – 是拉格朗日法对流体运动的描述 流线： – 某一时刻流场中不同位置的连续流体质点的流动方向线 – 是欧拉法对流体运动的描述 – 流线的性质：流线不能相交，也不能是折线，流线只能是一条光滑的曲线；对于稳定流动，流线与迹线相重合； – (3)系统与控制体 系统 ——某一确定的流体质点集合的总体，与外界无质量交换 流体系统的描述是与拉格朗日描述相对应 控制体 ——流场中确定的空间区域 可与外界进行质量交换和能量交换 控制体描述则是与欧拉描述相对应 2、质量守恒定律——连续性方程 m F F m F F m F F z bz y by x bx ===,,

流体力学作业答案

例题1：如下图所示，一圆锥体绕自身轴线等速旋转，锥体与固定壁面间的距离为K ，空隙全部被动力粘滞系数为 μ的牛顿流体所充满。当旋转角速度为ω，锥体底部半径为R ， 高为H ，求作用于圆锥的阻力矩。 2α ω H K R 解： M= ????====K dh r K dA r K dAr K u dAr 3 22cos 2πμω απμωωμμτ= H K R α πμωcos 23 而 22cos R H H += α;故：M=2 23 2R H K R +?πμω 例题2：涵洞进口处，装有与水平线成600 倾角而边长为1m 的正方形平板闸门（AB=1m ），求闸门所受静水总压力的大小及作用点。

B A 60 1m 解：坐标只能建在水面上。 A B 600 x y 1m a A kp p 807.91807.9=?= a B kp p 300.18)2 3 1(807.9=+?=

KN p p P B A 050.14112 =??+= h h A y I y y C C C C C D 6.160sin 433.1121 60sin 433.1160 sin 121160sin 00 30=+=??+=+= 0=D x y B 1m A 60 x D P y D 矩形和圆形的 C y 和C I 值 矩形： 2h y C = 123 bh I C =

C y C y X h b 圆形： r y C = 4 4 r I C π= x y c r y C 例题3：一直立矩形闸门，用三根工字梁支撑，门高及上游水深H 均为3m,把此闸门所受静水压强分布图分为三等份，每根工字梁分别设在这三等份的重心，求三个工字梁的位置？ 解：设静水压力分布图的面积为A ，则每一等份为A/3

流体力学路线图

流体力学基础理论的学习历来被初学者视为畏途，每到学习结束要进入期末考试的时候，老师和学生一样心中难免忐忑，在流体力学这门课上挂科已经成为某种常态。即使是学习多年的老手也会在具体问题面前感到基础尚不完备，还不够扎实。这个问题的起源当然与流体运动规律本身的复杂性有关，这个复杂性导致流体力学与大家印象中的“学科”概念有一定的出入。比如我们在学习高等数学时，很容易发现，数学是一门“咬文嚼字”的学科，里面充满严格定义的概念，不论学习线性代数还是微积分，都是从一些基本公理出发，循着一条严格的逻辑路线，架构起整门课程。因为数学有这样逻辑严密的特点，所以虽然学起来也不容易，但大家一致认为数学是美的，而且不论谁写的数学书，比如微积分的书，内容都只有程度深浅的差异，而绝没有内容上的巨大差异。 流体力学则有所不同，流体的流动本身是一种连续不断的变形过程，经典的流体力学理论以连续介质假设为基础，将整个流体看作连续介质，同时将其运动看作连续运动。但是由于流体是复杂的，实际上至今还没有完全掌握其全貌，因此流体力学在建立了基本控制方程后，就开始转而从一些特殊的流动出发，采用根据流动特点进行简化的方式，先建立物理模型，再得到数学模型，进而得到我们在书中经常看到的很多“理论”，比如不可压无旋流、旋涡动力学、水波动力学、气体动力学等等，甚至理论中还包括理论，比如不可压无旋流中还有自由流线理论，等等。形成一个类似于俄罗斯套娃的学科结构，这种结构容易给人一种支离破碎的印象。特别是在各个理论之间联系比较薄弱的时候，更容易给人这种印象。似乎一门课中又包含了很多门“小课”，每门“小课”使用的数学工具也完全不同，甚至很多同行还进一步把自己分成是学气的（比如空气动力学），或者是学水的（比如学船舶的）等等。 就象旅行者要有一张地图才能更高效率地到达目的地一样，如果能有一张流体力学的地图，或者叫路线图（roadmap），应该对初学者有很大帮助。这张图就是这门学科的脉络，其中应包含流体力学的主要理论内容，扩展一步的话，还应该包括数学基础（先修课）和主要分支学科。先在这里做个记号，有时间的时候慢慢地先从流体力学基础理论入手，给出一个粗略的路线图，然后再逐渐给出分支学科的路线图，比如空气动力学、计算流体力学的路线图。希望能抛砖引玉，激发出同行们的兴趣，加入绘制路线图的工作。在想象中，这个路线图应该有学科的主要内容，同时应该有相关的参考书。这样初学者就可以按图索骥，沿着一