第七章 不可压缩流体动力学基础.

Slides:



Advertisements
Similar presentations
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
Advertisements

第五节 全微分方程 一、全微分方程及其求法 二、积分因子法 三、一阶微分方程小结. 例如 所以是全微分方程. 定义 : 则 若有全微分形式 一、全微分方程及其求法.
第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
2.6 隐函数微分法 第二章 第二章 二、高阶导数 一、隐式定义的函数 三、可微函数的有理幂. 一、隐函数的导数 若由方程 可确定 y 是 x 的函数, 由 表示的函数, 称为显函数. 例如, 可确定显函数 可确定 y 是 x 的函数, 但此隐函数不能显化. 函数为隐函数. 则称此 隐函数求导方法.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
信号与系统 第三章 傅里叶变换 东北大学 2017/2/27.
§3.4 空间直线的方程.
《解析几何》 -Chapter 3 §7 空间两直线的相关位置.
3.4 空间直线的方程.
1.非线性振动和线性振动的根本区别 §4-2 一维非线性振动及其微分方程的近似解法 方程
碰撞 两物体互相接触时间极短而互作用力较大
碰撞分类 一般情况碰撞 1 完全弹性碰撞 动量和机械能均守恒 2 非弹性碰撞 动量守恒,机械能不守恒.
《解析几何》 乐山师范学院 0 引言 §1 二次曲线与直线的相关位置.
一、二阶行列式的引入 用消元法解二元线性方程组. 一、二阶行列式的引入 用消元法解二元线性方程组.
粘性流体力学 阎超 北京航空航天大学.
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
第五节 微积分基本公式 、变速直线运动中位置函数与速度 函数的联系 二、积分上限函数及其导数 三、牛顿—莱布尼茨公式.
一、原函数与不定积分 二、不定积分的几何意义 三、基本积分公式及积分法则 四、牛顿—莱布尼兹公式 五、小结
第二节 微积分基本公式 1、问题的提出 2、积分上限函数及其导数 3、牛顿—莱布尼茨公式 4、小结.
数 学 分 析 第九章 定积分 第二节 微积分学基本公式 主讲:师建国.
第四章 函数的积分学 第六节 微积分的基本公式 一、变上限定积分 二、微积分的基本公式.
系统 控制体 输运公式 1. 系统(system)——由确定的流体质点组成的流体团或流体体积V(t)。
2.3 液体动力学基础 本节主要讨论液体的流动状态、运动规律、能量转换以及流动液体与固体壁面的相互作用力等问题。
第四节 一阶线性微分方程 线性微分方程 伯努利方程 小结、作业 1/17.
第三节 格林公式及其应用(2) 一、曲线积分与路径无关的定义 二、曲线积分与路径无关的条件 三、二元函数的全微分的求积 四、小结.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
1.5 场函数的高阶微分运算 1、场函数的三种基本微分运算 标量场的梯度f ,矢量场的散度F 和F 旋度简称 “三度” 运算。
全 微 分 欧阳顺湘 北京师范大学珠海分校
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
习题六 1. 判断下列流场是否有旋?并分别求出其流线、计算oxy平面的单位圆周上的速度环量。 柱坐标 [解] 计算旋度 计算流线 速度环量
§3.7 热力学基本方程及麦克斯韦关系式 热力学状态函数 H, A, G 组合辅助函数 U, H → 能量计算
全国高校数学微课程教学设计竞赛 知识点名称: 导数的定义.
太沙基一维固结理论 饱和土体在外力作用下其孔隙水压力逐渐消散、有效应力同步增长、体积压缩。
第8章 静电场 图为1930年E.O.劳伦斯制成的世界上第一台回旋加速器.
2.1.2 空间中直线与直线 之间的位置关系.
流体力学基础 流体静力学 连续性原理 伯努利方程.
第一章 函数与极限.
第二十二章 曲面积分 §1 第一型曲面积分 §2 第二型曲面积分 §3 高斯公式与斯托克斯公式.
实数与向量的积.
正方形 ——计成保.
第七章 理想不可压缩流体的有旋流动和无旋流动
⑴当∠MBN绕点B旋转到AE=CF时(如图1),比较AE+CF与EF的大小关系,并证明你的结论。
§1体积求法 一、旋转体的体积 二、平行截面面积为已知的立体的体积 三、小结.
第五节 对坐标的曲面积分 一、 对坐标的曲面积分的概念与性质 二、对坐标的曲面积分的计算法 三、两类曲面积分的联系.
作业 P158 习题 2 1(2)(4) (5). 2(1). 预习 P156— /5/2.
作业 P152 习题 复习:P 预习:P /5/2.
人工压缩性方法 郭 红.
(七)不可压缩流的数值方法 7.1 MAC方法 7.2 投影法 7.3 人工压缩性方法 7.4 SIMPLE方法 7.5 其他方法:
流体佯谬 由于牛顿力学的巨大成功,人们对牛顿确定的三大定律深信不疑,奉之为金科玉律,然而在生活中,我们常常会惊异的发现流体表现出一些意想不到的效应。例如:
复习: 若A(x1,y1,z1) , B(x2,y2,z2), 则 AB = OB - OA=(x2-x1 , y2-y1 , z2-z1)
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
3.1 变化率与导数   3.1.1 变化率问题 3.1.2 导数的概念.
3.1.2 空间向量的数量积运算 1.了解空间向量夹角的概念及表示方法. 2.掌握空间向量数量积的计算方法及应用.
相关与回归 非确定关系 在宏观上存在关系,但并未精确到可以用函数关系来表达。青少年身高与年龄,体重与体表面积 非确定关系:
第三章 函数的微分学 第二节 导数的四则运算法则 一、导数的四则运算 二、偏导数的求法.
高中数学必修 平面向量的基本定理.
直线的倾斜角与斜率.
第四节 向量的乘积 一、两向量的数量积 二、两向量的向量积.
2.2.1质点的动量及动量定理 2.2 动量 动量守恒定律 1. 冲量 力在时间上的积累,即冲量。 恒力的冲量 (t1 → t2): z
3.2 平面向量基本定理.
制作者:王翠艳 李晓荣 o.
第三节 数量积 向量积 混合积 一、向量的数量积 二、向量的向量积 三、向量的混合积 四、小结 思考题.
§2.高斯定理(Gauss theorem) 一.电通量(electric flux) 1.定义:通过电场中某一个面的电力线条数。
《偏微分方程》第一章 绪论 第一章 绪论 1.1.
第三章 图形的平移与旋转.
3.3.2 两点间的距离 山东省临沂第一中学.
Presentation transcript:

第七章 不可压缩流体动力学基础

本章主要讨论三元流动问题,即讨论有关流动问题的基本概念和基本原理以及描述不可压缩流体运动的基本方程和定解条件。 本章主要内容 本章主要讨论三元流动问题,即讨论有关流动问题的基本概念和基本原理以及描述不可压缩流体运动的基本方程和定解条件。 重点、难点内容 流体微团运动的分析 有旋流动、无旋流动 理想流体运动微分方程 涡线、涡管以及斯托克斯定理

第一节 流体微团运动的分析 分析流场中任意流体微团运动是研究整个流场运动的基础。 第一节 流体微团运动的分析 分析流场中任意流体微团运动是研究整个流场运动的基础。 流体运动要比刚体运动复杂得多,流体微团基本运动形式有平移运动,旋转运动和变形运动等,而变形运动又包括线变形和角变形两种。

平面流动 平移 转动 线变形 角变形

平移运动、旋转运动、线变形运动和角变形运动 右图为任意t时刻在平面流场中所取的一个正方形流体微团。由于流体微团上各点的运动速度不一致,经过微小的时间间隔后,该流体微团的形状和大小会发生变化,变成了斜四边形。

流体微团的运动形式与微团内各点速度的变化有关。设方形流体微团中心 M 的流速分量为 ux 和 uy (图 7-1 ) ,则微团各侧边的中点 A 、 B 、 C 、 D 的流速分量分别为:

可见,微团上每一点的速度都包含中心点的速度以及由于坐标位置不同所引起的速度增量两个组成部分。

平移运动速度 微团上各点公有的分速度 ux 和uy ,使它们在 dt 时间内均沿 x 方向移动一距离 uxdt , 沿 y 方向移动一距离 uydt 。因而,我们把中心点 M 的速度 ux和 uy ,定义为流体微团的平移运动速度。

线变形运动 微团左、右两侧的 A 点和 C 点沿 x 方向的速度差为 当这速度差值为正时,微团沿 x 方向发生伸长变形;当它为负时,微团沿 x 方向发生缩短变形。

线变形速度 单位时间,单位长度的线变形称为线变形速度。 以θx表示流体微团沿 x 方向的线变形速度,则:

三元流动线变形速度

微团的旋转和角变形 流体微团的运动过程可以看作下述两种基本运动的组合: 流体微团绕M点作无角变形的旋转运动,ABCD-A’B’C’D’;

旋转角速度 设沿逆时针方向旋转为正,则AMC线的旋转角速度为: BMD线的旋转角速度为:

把对角线的旋转角速度定义为整个流体微团在平面上的旋转角速度。 因而,角速度矢量为: 角速度的大小为:

角变形速度 直角边 AMC (或BMD)与对角线 EMF 的夹角的变形速度定义为流体微团的角变形速度。

亥姆霍兹速度分解定理

亥姆霍兹速度分解定理

第二节 有旋流动 流体微团的旋转角速度在流场内不完全为零的流动称为有旋流动。 自然界和工程中出现的流动大多数是有旋流动,例如: 龙卷风 第二节 有旋流动 流体微团的旋转角速度在流场内不完全为零的流动称为有旋流动。 自然界和工程中出现的流动大多数是有旋流动,例如: 龙卷风 管道流体运动 绕流物体表面的边界层及其尾部后面的流动。

有旋流动与无旋流动 无旋流动 有旋流动

有旋流动与无旋流动 无旋流动 有旋流动

涡量

涡量连续性微分方程

涡线 在涡量场中可以画出表征某一瞬时流体质点的旋转角速度向量方向的曲线,称为涡线。 在给定的瞬时,涡线上各点的角速度向量在该点处与涡线相切。

涡线微分方程 沿涡线取一微小线段ds,由于涡线与角速度向量的方向一致,所以,沿三个坐标轴方向的分量dx, dy, dz必然与角速度向量的三个分量成正比,即:

涡管 在涡量场中任意画一封闭曲线,通过这条曲线上的每一点所作出的涡线构成一管状的曲面,称为涡管。

涡通量

涡管强度守恒定理 可以证明: 对于微元涡管,可以近似认为截面上各点的涡量为常数,则由上式: 涡管截面愈小的地方,流体的旋转角速度愈大。

涡量场中的涡线、涡管、涡通量等概念分别与流速场中的流线、流管、流量等概念相对应,而涡线方程与涡管的涡通量方程则分别与流线方程和元流连续性方程相对应。

速度环量 通常,涡通量是利用速度环量这个概念来计算的。 在流场中任取一封闭曲线s,则流速沿曲线s的积分称为曲线s上的速度环量。并规定积分沿s逆时针方向绕行为s的正方向。

斯托克斯定理 沿任意封闭曲线s的速度环量等于通过以该曲线为边界的曲面 A 的涡通量。

第三节 不可压缩流体连续性微分方程 - = 直角坐标系中的连续性方程 质量守恒 输入微元体 的质量流量 输出微元体 的质量流量 微元体内的 第三节 不可压缩流体连续性微分方程 直角坐标系中的连续性方程 质量守恒 y x z dz dx dy 输入微元体 的质量流量 输出微元体 的质量流量 - 微元体内的 质量变化率 = 微元体及其表面的质量通量

连续性方程 1、x方向:dt时间内沿从六面体 x 处与 x+dx 处输入与输出的质量差: y方向: z方向:

连续性方程 3、微元体内的质量变化: 从而有: 或: 矢量形式: (适用于层流、湍流、牛顿、非牛顿流体) 矢量形式: 连续方程物理意义:流体在单位时间内流经单位体积空间输出与输入的质量差与其内部质量变化的代数和为零。

若流体不可压缩: 上式表明,对于不可压缩液体,单位时间单位体积空间内流入与流出的液体体积之差等于零,即液体体积守恒。 恒定流或非恒定流; 理想液体或实际液体。 适用范围: 一维流动的连续方程 连续性方程是流体流动微分方程最基本的方程之一。任何流体的连续运动均必须满足。

例:已知不可压流体速度, 估算w。 解: 不可压流体

第四节 以应力表示的粘性流体运动微分方程 应力状态及切应力互等定律 应力状态: 切应力互等定律 y x z 应力状态: 粘性流场中任意一点的应力有9个分量,包括3个正应力分量和6个切应力分量: 切应力互等定律 微元体上X和Z方向的表面力 在6个切应力分量中,互换下标的每一对切应力是相等的。

微元体表面力的总力分量 x方向的表面力: y方向的表面力: z方向的表面力:

x = 动量流量及动量变化率 动量流量 动量通量 流通面积 动量流量 dy z dx dz y x 动量在微元体表面的输入与输出 图中标注的是动量的输入或输出方向,而动量或其通量本身的方向均指向x方向,即分速度vx的方向。 动量在微元体表面的输入与输出

输入输出微元体的动量流量 微元体内的动量变化率 x方向: y方向: z方向: x方向: y方向: z方向: 流体的瞬时质量为

以应力表示的运动方程 x方向的运动方程: y方向的运动方程: z方向的运动方程: 注:上式就是以应力表示的粘性流体的运动方程,   适用于层流、湍流、牛顿、非牛顿流体。

方程的物理意义: 方程可简略表示成: 这就是以单位体积的流体质量为基准的牛顿第二运动定律 方程左边是:任意时刻t通过考察点A的流体质点加速度的三个分量; 方程右边是:作用在单位体积流体上的表面力和体积力在各坐标上的分量。 方程可简略表示成: 这就是以单位体积的流体质量为基准的牛顿第二运动定律

第五节 应力和变形速度的关系

粘性流体运动微分方程 Navier-Stokes方程 对一维流动问题: 对粘性流体流动问题: 关键:寻求流体应力与变形速率之间的关系 以应力表示的运动方程,需补充方程才能求解。 对一维流动问题: 补充方程:牛顿剪切定律 对粘性流体流动问题: 补充方程:广义的牛顿剪切定律 即:牛顿流体本构方程 目的 关键:寻求流体应力与变形速率之间的关系 将应力从运动方程中消去,得到由速度分量和压力表示的粘性流体运动微分方程,即N-S方程。

N-S方程 牛顿流体的本构方程 引入的基本假设: 应力与变形速率成线性关系; 应力与变形速率之间的关系各向同性; 为了寻求流体应力与变形速率之间的关系,Stokes提出三个基本假设: 应力与变形速率成线性关系; 应力与变形速率之间的关系各向同性; 静止流场中,切应力为零,各正应力均等于静压力

牛顿流体的本构方程:

本构方程的讨论: 正应力与线变形速率: 线变形率与流体流动: 正应力中的粘性应力: 流体正应力与三个速度偏导数有关 (即:线变形率),同固体力学中的虎克定律。 线变形率与流体流动: 从流体流动角度看,线变形率的正负反映了流体的流动是加速还是减速;体变形率的正负反映了流动过程中流体体积是增加还是减少。 正应力中的粘性应力: 附加粘性正应力 附加粘性正应力的产生是速度沿流动方向的变化所导致的。

牛顿流体本构方程反映了流体应力与变形速率之间的关系,是流体力学的虎克定律(反映应力和应变的关系)。 正应力与压力: 由于粘性正应力的存在,流动流体的压力在数值上一般不等于正应力值。但有: 这说明:三个正压力在数值上一般不等于压力,但它们的平均值却总是与压力大小相等。 切应力与角边形率: 流体切应力与角变形率相关。 牛顿流体本构方程反映了流体应力与变形速率之间的关系,是流体力学的虎克定律(反映应力和应变的关系)。

适用于牛顿流体 第六节 N-S方程 流体运动微分方程——Navier-Stokes方程

适用于牛顿流体 常见条件下N-S方程的表达形式: 常粘度条件下N-S方程: 矢量形式:

适用于牛顿流体 不可压缩流体的N-S方程: 矢量形式:

适用于牛顿流体 常粘度条件下不可压缩流体的N-S方程: 矢量形式: 扩散项(粘性力项) 对静止或理想流体为0 高速非边界层问题≈0 非定常项 定常流动为0 静止流场为0 对流项 静止流场为0 蠕变流时≈ 0 单位质量流体 的体积力 单位质量流体 的压力差

第八节 流体流动的初始条件和边界条件 N-S方程应用概述 连续方程和N-S方程是粘性流体流动应遵循的质量守恒和 动量守恒的数学表达式。 封闭条件:理论上方程是封闭的,但若要考虑到物性参数的变化,应将物性变化的关系作为补充方程。 应用条件:只适用于牛顿流体 方程求解:N-S方程无普遍解;特殊条件下,有可能获得准确或近似的分析解;通常通过数值计算获得离散解。

流动微分方程的应用求解步骤 根据问题特点对一般形式的运动方程进行简化,获得针对具体问题的微分方程或方程组。 提出相关的初始条件和边界条件。   初始条件:非稳态问题 流体具有粘性,在与壁面接触处流体速度为零。 固壁-流体边界: 边界条件 对非高速流,气液界面上,液相速度梯度为零。 液体-气体边界: 液体-液体边界: 液液界面两侧的速度或切应力相等。

N-S方程应用举例: 例: 圆管内的一维稳态流动分析。 例: 圆管内的一维稳态流动分析。 不可压缩流体在水平 圆管内作一维稳态层流流动。试写出该条件下的连续性方程和运动微分方程。并证明管道截面上任一点的总势能和轴向压力梯度为常数。

例题

表明vz只是r的函数; 只是z的函数

作业 7-2、7-3、7-4、7-5