数值传热学陶文铨第三章
数值传热学(陶文铨)第二章习题答案
p1 p2
（1）
p 2 p x2 x 2 x i 2 x i 2 2!
p 2 p x2 p3 p2 x 2 x i 2 x i 2 2!
（2） 式（2）-（1）得
p p1 p 3 x i 2 2 x
2-11.解：对于均分网格用泰勒级数展开法，用 2 点分别表示 1,3,4 处热流量得
（i+1,n）= （i,n）+
(1)
x
( x)
i ,n
2 x 2
i ,n
(( x) )2 2!
（i-1,n）= （i,n）+
(2) 式（1）-（2）得：
x
( x)
i ,n
2 x 2
i ,n
(( x) )2 2!
（i+1,n）- （i-1,n）=
即
krp k krp k 2k rpTp ( )Te ( )T r s r r r 2 r 2 w p
化简后得
2krp
r
Tp
kre kr S Te w Tw rp r r r 2
计算结果与控制容积积分法一致。 2-9.解：对于均分网格用泰勒级数展开法分别表示 1 和 3 点处的压力值
（1） 对 T 随 r 由 Tw 变到 Te 的过程进行积分
（2） 可化为
rk
dT dr
e
w
S 2 r 2 w
e
（3） 取 T 随 r 呈分段线性的变化，则（3）式中
Tw
Tp Tw 2 Tp Te 2
（4）
Te
（5）
dT Te Tp r dr e
（6）
数值传热学(课件)-1
热流问题的数值计算Numerical Simulations of Thermal & Fluid Problems第一章 绪论主讲 陶文铨西安交通大学能源与动力工程学院 热流中心 CFD-NHT-EHT CENTER 2007年10月16日, 西安1/88物理问题数值解的基本思想 把原来在空间与时间坐标中连续的物理量的场 (如速度场,温度场,浓度场等),用一系列有限 个离散点(称为节点,node)上的值的集合来代替; 通过一定的原则建立起这些离散点上变量值之间关 系的代数方程 (称为离散方程,discretizationequation);求解所建立起来的代数方程以获得所求解变量的近似解.2/88大规模科学计算的重要性 传热与流动问题数值计算是应用计算机求解热量传 递过程中的速度场,温度场等的分支学科,是大规模 科学计算的重要组成部分,其重要性不言而喻. 2005年美国总统顾问委员会向美国总统提出要大 力发展计算科学以确保美国在世界上的竞争能力. 波音公司实现了对航空发动机的网格数达10亿量 级的直接数值模拟,以研究所设计发动机的性能.3/88现代科学研究的三大基本方法及其关系理论分析Analytical实验研究Experimental数值模拟Numerical4/88课程简介1. 学时- 30学时理论教学;6学时计算机作业 2. 考核- 平时作业/计算机大作业/考试: 20/30/50 3. 方法- 理解,参与,应用 努力将与数学处理相对应的物理背景联系起来理解. 4. 助手- 于乐 5. 参考教材-《计算流体力学与传热学》,中国建筑 工业出版社,19915/88学习方法建议1. 善于从物理过程基本特性来掌握理解数值方法; 2. 对数值方法-明其全而析其微:明其全-了解基本原理;析其微-掌握实施细节;3. 努力上机实践; 4. 学会分析计算结果: 合理性,规律性; 5. 应用商业软件与自编程序相结合.6/88《热流问题的数值计算》 主要教学内容第一章 绪论(物理与数学基础) 第二章 一维导热问题的数值解 第三章 多维导热问题的数值解 第四章 势流及管道内充分发展流动与换热的数值解 第五章 有回流的动与换热问题的数值解 第六章 二维涡量-流函数法通用程序介绍 第七章 原始变量法与湍流数值模拟简介7/88绪论1.1 流动与传热问题控制方程的基本类型 1.2 流动与传热问题数值计算的基本步骤 1.3 建立离散方程的方法 1.4 离散方程数学与物理特性分析简介8/881.1 流动与传热问题控制方程的基本类型1.1.1 流动与传热问题完整的数学描写 1.1.2 控制方程 1. 质量守恒方程 3. 能量守恒方程 1.1.3 单值性条件 1.1.4 建立数学描写举例 1.1.5 控制方程式的分类9/882. 动量守恒方程1.1 流动与传热问题控制方程的基本类型1.1.1 流动与传热问题完整的数学描写 1. 有关的守恒定律的偏微分方程(控制方程)一切宏观的流动与传热问题都由三个守恒定律所 支配:质量,动量与能量守恒(conservation law).2. 与表述守恒定律的偏微分方程相关的单值性条件.不同问题的区别主要在于单值性条件 (conditions for unique solution) 的不同:初始条件以,边界条件 以及物性数据.10/881.1.2 控制方程(Governing equations) Mass conservation1. 质量守恒方程r ( r u ) ( r v) ( r w) + + + =0 t x y z单位时间 内质量的 增加 单位时间内流 进微元体的净 质量物理意义:单位时间内空 间某一微元容积质量的增 加等于流入该微元容积的 净质量.11/88对不可压缩流体: r = const 对二维不可压缩流体:u v + =0 x yu v w + + =0 x y z对二维问题,速度矢量:ur u v 数学上称: + = div(U ) x yur r ur U =ui+v j为速度矢量的散度,因此对二维不可压流体有:ur div(U ) = 0下面只讨论不可压缩流体(incompressible flow).12/882. 动量守恒方程(Momentum conservation)对上图所示的微元体分别在三个坐标方向上应用 Newton第2定律(F=ma)在流体中的表现形式: [微元体内动量的增加率]=[作用在微元体上各种力之和] 可得出三个坐标方向的动量方程:u uu uv uw 1 p 2u 2u 2u + + + =+ n ( 2 + 2 + 2 ) + Fx t x y z r x x y z 1 p v vu vv vw 2v 2v 2v + + + =+ n ( 2 + 2 + 2 ) + Fy t x y z r y x y z 1 p w wu wv ww 2 w 2 w 2 w + + + =+ n ( 2 + 2 + 2 ) + Fz t x y z r z x y z微元体内动 量的增加率压力粘性力体积力13/883. 能量守恒方程(Energy conservation)[微元体内热力学能的增加率]=[通过流动与导热进入 微元体内的净热流量]+[体积力与表面力对微元体所做 的功率] 引入导热Fourier定律,假定热物性为常数,可得T (uT ) (vT ) ( wT ) 2T 2T 2T rcp[ + + + ] = l( 2 + 2 + 2 ) + S t x y z x y z微元体 内能增 加率 由于流动被带出 微元体的净功率 由于导热而进入 源项 微元体的净功率 生成 热14/88l =a rcp流体的热扩散率(thermal diffusivity)4. 对于二维稳态对流换热问题控制方程汇总u v + =0 x yuu uv 2u 2u 1 p + =+ n ( 2 + 2 ) + Fx y z r x x yvu vv 2v 2v 1 p + =+ n ( 2 + 2 ) + Fy y z r y x y(uT ) (vT ) 2T 2T + = a( 2 + 2 ) + ST x y x y对流项扩散项源项数值计算中常用的术语.15/88不同的二维,稳态求解问题之间的区别在于: (1)边界条件不同; (2)源项与扩散系数不同.5. 二点说明1. 所导出的三维非稳态Navier-Stokes方程,无论对 层流或是湍流都是适用的. 2. 辐射换热需要用积分方程来描述,课程中将不涉及 这类问题.16/881.1.3 单值性条件 1. 初始条件 2. 边界条件 (1) 第一类 (Dirichlet):t = 0, T = f ( x, y, z )TB = Tgiven(2) 第二类 (Neumann): qB = -l (T ) B = qgiven n(3) 第三类 (Rubin):规定了边界上被求函数的一阶导数与函数之间的关系: -l ( T ) B = h(TB - T f )n数值计算中计算区域的出口边界条件常常最难 确定,要做近似处理.17/881.1.4 建立数学描写举例 1. 问题与假设条件突扩区域中的对流传热:二维,稳态,不可压缩, 常物性,不计重力与黏性耗散.18/882. 控制方程u v + =0 x y1 p u u u u u +v =+n ( 2 + 2 ) r x x y x y 2 2 v v 1 p v v u +v =+n ( 2 + 2 ) x y r y x y2 2T T T T u +v = a( 2 + 2 ) x y x y2 219/883. 边界条件 (1)进口边界条件:给定u,v,T随y 的分布; (3)中心线: u = T = 0; v = 0 y y(4)出口边y x界:数学上要 求给定u,v,T 或其导数随y 的分布;实际 上做不到;数 值上近似处理20/88(2)固体边界条件:速度无滑移,温度无跳跃1.1.5 传热与流动问题的数学描写的分类 1. 从数学角度分类-椭圆型与抛物型椭圆型 (Elliptic)椭圆型方程数学上的特点是:所求解的因变量对每个 空间自变量均存在二阶导数项: 导热方程-所求解的因变量为温度T ,空间自变量x,y; 动量方程-所求解的因变量为速度u ,空间自变量x,y.21/88抛物型(Parabolic)抛物型方程数学上的特点是:所求解的因变量对某个 个自变量只存在一阶导数项: 非稳态导热方程-因变量T 对时间t仅有一阶导数; 边界层动量方程-u对空间自变量x仅有一阶导数. 仅存在一阶导数的自变量在物理过程上的重要特 点:过程只能沿该坐标的单个方向进行而不能逆向进 行.22/88抛物型与椭圆型流动的例子椭圆型方程的求解必须全场联立进行,而抛物性 方程的求解可以沿坐标正向逐步推进, 大大节省时间.23/88(1)椭圆型问题: 流动有回流,必须 全场同时求解; (2)抛物型问题:流动无回流,可以沿主流方向步 步逼进,不必全场同时求解,大大节省时间.Marching method24/882. 从物理角度分类-守恒型与非守恒型守恒型( Conservative)-对任意大小容积守恒特性 都能得到满足的方程; 凡对流项表示成散度形式的方程具有守恒性 . 非守恒型方程+u v v u u v u ++ u = 0= 0 u ( + ) = 0 x x y y x y (uu ) (uv) 1 p 2u 2 v =+n ( 2 + 2 ) + r x x x y x守恒型方程凡是从守恒型控制方程推导得到的用于数值求解 的代数方程也具有守恒特性.25/881.2 流动与传热问题数值求解的基本步骤1.2.1 流动与传热问题数值求解步骤 1. 建立数理模型 3. 方程的离散化 5.代数方程求解 1.2.2 区域离散化方法 2.区域的离散化 4. 边界条件离散 6. 求解结果分析1.区域离散化的任务 2. 区域离散方法1.2.3 网格系统标记方法26/881) 外节点法2. 内节点法1.2.1 流动与传热问题数值求解步骤把原来在空间与时间坐标中连续的物理量的场 (如速度场,温度场,浓度场等),用一系列有限个 离散点(称为节点,node)上的值的集合来代替;通过 一定的原则建立起这些离散点上变量值之间关系的代 数方程(称为离散方程,discretization equation);求 解所建立起来的代数方程以获得所求解变量的近似解.27/88(1) 区域离散 (2) (3) (4) (5) 代数求解 (6)28/88方程离散结果分析1.2.2 区域离散化1.区域离散化的任务将所计算的区域分割成许多不重叠的子区域,确 定每个子区域中节点的位置以及所代表的控制容积. 离散结果得出四种几何要素: (1) 节点(node):所求解未知量的位置; (2) 控制容积(control volume):实施守恒定律的最 小几何单位; (3) 界面(interface):控制容积的分界位置; (4) 网格线(grid lines):沿坐标方向相邻节点连接 成的曲线簇.29/882. 区域离散方法 (a) 外节点法:节点位于子区域的角顶;控制容积界 面位于两节点之间;生成过程:先节点后界面;又 称 Practice A.子区域控制容积30/88YPractice A-外节点法 x31/88(b) 内节点法:节点位于子区域的中心;子区域即为 控制容积;生成过程:先界面,后节点,又称 Practice B.子区域即为控制容积32/88YPractice B-内节点法 x33/88 1.2.3 内接点与外节点法的比较 (a)边界节点所代表的控制容积不同 方法A 边界节点代表半个CV方法B 边界节点代表零个CV(b)网格非均分时,节点作为控制容积的代表方法B 更合理 方法A 方法B34/881.2.3 网格系统表示方法 网格线-节点间连线,用实线表示;界面为虚线; 节点间距离-dx;界面间距离-Dx .35/881.2.4 网格独立解 当网格足够细密以至于再进一步加密网格已对 数值计算结果基本上没有影响时所得到的数值解称 为网格独立解(grid-independent solution).Int. Journal Numerical Methods in Fluids, 1998, 28: 1371-1387.36/881.3 建立离散方程的方法 1.3.1 一维模型方程( 1-D model equation ) 1.3.2 由Taylor 展开法导出导数的差分表示式 1.3.3 控制容积积分法导出导数的差分表示式 1.3.4 讨论37/881.3 建立离散方程的方法 1.3.1 一维模型方程( 1-D model equation ) 一维模型方程是一维非稳态有源项的对流-扩 散方程,具有四个特征项,便于离散方法的研讨. 非守恒型 守恒型 ( rf ) f f + ru = (G ) + Sf t t x xFDM采用 ( rf ) ( r uf ) f + = (G ) + Sf FVM采用 t t x x 瞬态 对流 扩散 源项38/88"麻雀虽小,五脏俱全!"1.3.2 由Taylor 展开法导出导数的差分表示式 1. 一阶导数的差分表达式的导出 将函数f ( x, t ) 在(i+1,n)的值对(i,n)点做Taylor展开:f 2f Dx 2 2 f (i + 1, n) = f (i, n) + )i ,n Dx + 2 )i ,n Dx + ..... x x 2!f f (i + 1, n) - f (i, n) Dx 2f ) i ,n = - ( 2 )i ,n + ... x Dx 2 x39/88O ( Dx ) 称为截断误差, truncation error,表示:随 Dx 的趋于零,用 f (i + 1, n) - f (i, n) 代替 f )i ,n 的误差 x Dxf f (i + 1, n) - f (i, n) )i ,n = + O(Dx) x Dx KD x, K 与 Dx 无关.D x 的方次称为截差的阶数(order of TE).用数值计算的近似解 fin 代替精确解 f (i, n)fin 1 - fin f )i ,n @ + , O(Dx) 得向前差分: x Dx40/88f -f f )i ,n @ 向后差分: x Dxn in i -1, O (Dx )fin 1 - fin 1 f )i , n @ + , O(Dx 2 ) 中心差分: x 2Dx2. 一,二阶导数的各种差分表达式. 表达差分结构的格式图案o构筑差分表达式的位置; 构筑差分表达式所用到的节点.41/88一阶导数的 常用差分表达式42/88二阶导数的常用差分表达式定性判别导数的差分表达式正确与否的方法: (1)量纲是否正确-与导数本身一致; (2)均匀场的各阶导数应为零.43/883. 一维模型方程的有限差分显式离散表示式 微分方程形式: 假设 ( rf ) f f + ru = (G ) t t x xr , u, G均为常数,显式差分表达式:fin +1 - fin fin 1 - fin 1 r + ru + = Dt 2Dx fin 1 - 2fin + fin 1 G + , O (Dt , Dx 2 ) Dx 2差分方程 截断误差44/88显式(Explicit)-空间导数均以初 始时刻之值计算.1.3.3 控制容积积分法导出导数的差分表示式 1. 控制容积积分法实施步骤 1. 将守恒型的方程对控制容积做积分; 2. 选定被求函数及其一阶导数对时间,空间的变化 曲线-型线; 3. 完成积分,整理成相邻节点间未知量的代数方程. 2. 两种常用型线 型线-被求函数随自变量的局部变化方式,本是 所求内容,近似求解需先假定.45/88随空间自变量的变化型线 型线 型线分段线性阶梯逼近46/88piece-wise linear step-wise approximation随时间自变量的变化型线分段线性 piece-wise linear阶梯逼近 step-wise approximation47/883. 一维模型方程的控制容积积分法离散 将守恒型控制方程对控制容积P 在[t, t+ Dt ]内 做积分, ( rf ) ( r uf ) ft立即可得e+xt +Dt t=xe(Gx)r ò (ft +Dt -ft )dx +rwò [(uf)òt- (uf)w ]dt =t +Dt=Gf f [( )e - ( ) w ]dt x xf 以及 x48/88继续积分,需要知道:f对空间与时间的变化型线.1. 非稳态项假设 f 对空间呈阶梯型变化:t t r ò (f t +Dt - f t )dx = r (f P+Dt - f P )Dx w e2. 对流项假设 f 对时间呈显示阶梯型变化:rt +Dtòt[(uf )e - (uf ) w ]dt = r[(uf )te - (uf )tw ]Dt49/88假设 f 对空间呈分段线性变化:fE + fP fP + fW fE - fW r[(uf ) - (uf ) ]Dt = r uDt ( ) = r uDt 2 2 2t e t w均分网格3. 扩散项f 假设 对时间呈显式阶梯型变化: xt +DtGòtf f f t f t [( )e - ( ) w ]dt = G[( )e - ( ) w ]Dt x x x x50/88假设 f 对空间呈分段线性变化:。
第一章数值传热学
(uT ) (vT ) T T a( 2 2 ) x y x y
2 2
19/80
MOE KLTFSE
3. 边界条件
定u,v,T随 y 的分布；
（1）进口边界条件：给
u T （3）中心线： 0; v 0 y y
y x
界：数学上要 求给定u,v,T或 其导数随 y 的 分布；实际上 做不到；数值 上近似处理。
cp
c p
( ) c p
Pr
12/80
MOE KLTFSE
4. 通用控制方程
( ) * * div( U ) div( grad ) S t
瞬态项 对流项 扩散项 广义源项 不同求解变量之间的区别： （1）边界条件与初始条件不同； （2）广义源项表达式不同； （3）广义扩散系数不同。 文献中常以表格形式给出所求解变量的源项与 广义扩散系数的表达式。
常物性不可压缩流体动量方程源项中显含速度部分 为零。
11/80
MOE KLTFSE
3. 能量守恒方程
[微元体内热力学能的增加率]=[进入微元体内的净热 流量]+[体积力与表面力对微元体所做的功] 引入导热Fourier定律，忽略力所作的功， 设hc
pT ;
c p 为常数
( T ) div( T U ) div( gradT ) ST cp t
4/80
MOE KLTFSE
绪论教学目录
1.1 传热与流动问题的数学描写 1.2 传热与流动问题数值计算的基本思想及应 用举例 1.3 传热与流动问题的数学描写的分类及其对 数值解的影响 1.4 传热与流动问题的数值计算的近代发展
5/80
数值传热学3
故得
u
n i 1
2x
n i 1
2 u ( ) i ,n O ( x ) x
假设源项不存在截断误差，则： 1-D模型方程的FTCS格式的截断误差即为：
O(t, x )
2
截断误差的数学含义是： 存在着两个正的常数，K1，K2，当 t 0, x 0
2/61
3.1 离散方程的相容性、收敛性与稳定性 3.1.1 截断误差及相容性 3.1.2 离散误差与收敛性 3.1.3 舍入误差与初值问题的稳定性 3.1.4 数值特性分析举例
3/61
3.1 离散方程的相容性、收敛性与稳定性 3.1.1 截断误差及相容性（consistency)
1. 离散方程的精确解
18/61
at r , 离散方程可化为 记 2 x
n n Ti n1 Ti n (1 2r ) r (Ti T 1 i 1 ), i 1, 2,....( I 1)
对确定的时层n,上式可展开写出为
i 1, i 2, i 3,
T T T
n 1 1 n 1 2 n 1 3
n1 n T AT g 0 0 T F
设包含误差的解为 T ， 则有
(b)－(a)
n1 n T AT g (b) 0 T F n1 n1 n n T T A(T T ) 0 0 0 T T
d 2 d 2 0 2 dx dx (0) 0; (4) 1
d 2 d 解：采用有限差分法，用差分式代替 dx 2 , dx 。
采用二阶截差时，可对节点2，3，4建立离散方程； 另外可对节点3采用四阶截差，对节点2，4二阶截 差。
沈阳建筑大学加试《传热学》2020年考研专业课复试大纲
2020年硕士研究生考试
同等学力加试传热学科目考试大纲
一、考查目标
按全国硕士研究生入学考试要求为沈阳建筑大学招收建筑设备与环境、供暖通风与空调专业硕士研究生而设置的专业课程考试科目。
其中，传热学是属招生学校自行命题的性质。
它的考查目标是高等学校优秀本科毕业生能达到的及格或及格以上水平，以保证被录取者具有基本的传热理论知识并有利于招生学校在专业上择优选拔。
传热学考试的目标在于考查考生对传热学基本概念、基本理论的掌握和分析求解基本问题的能力。
考生应能：
1. 准确地把握定义的物理量以及它们的量纲；
2. 正确理解基本概念和基本规律；
3. 正确应用基本理论知识分析和处理实际传热问题；
4. 掌握基本计算方法，准确完成传热问题的定量计算。
二、考试形式与试卷结构
（一）试卷满分及考试时间
传热学满分为100分，考试时间为2小时。
（二）答题方式
答题方式为闭卷、笔试。
（三）试卷内容结构。
数值传热学习题答案(汇总版)
2-4-9
= rP rS
式（2-4-9）也可以写成 a PTP = a E TE + aW TW + b 的形式。而且两种结果是一致的。
2—6：
n n TE −TW dT P , n = 解：将 ， dx 2x n n TE −2TPn + TW d 2T P , n = , dx2 x 2
dk = f (x ) 代入原方程，得： dx
令
2-4-4
rk rk a E = ， aW = ， a P = a E + aW ， b x w x e
= SrP r ，
式（2-4-4）可以写成 a PTP = a E TE + aW TW + b 的形式。 2. 再用 Taylor 展开法导出 k
2 2 uE + uP u = ， 2 2 e
2 2 uW + uP u = 2 2 w
t u ut N − uP y = (y ) ， n n
t
t ut u p − uS y = (y ) 。 s s
t
(y ) n = (y ) s = y
n n n n TE −TW TE −2TPn + TW k + f (x ) +S=0 整理得： 2x x 2
4kT P= 2k + xf ( x)T E+2k − xf ( x)T W +2x 2 S
− 2k 时， a E 会成为负值， x 2k 当 f(x）> 时， aW 会成为负值。 x
rk dr = rk r r dr dr dr
w
e
1 d
陶文铨老师_数值传热学_热流问题的数值计算02
Numerical Simulations of Thermal & Fluid Problems
第二章 一维导热问题的数值解
主讲 陶文铨
西安交通大学能源与动力工程学院 热流中心 CFD-NHT-EHT CENTER 2007年10月23日, 西安
1/62
第二章 一维导热问题的数值解 2.1 一维稳态导热 2.2 边界条件及源项的处理 2.3 一维非稳态导热 2.4 代数方程的求解方法 2.5 一维非稳态导热通用程序介绍 2.6 数值计算举例
A= B=
δ xiΔxi S hiδ x )T f TM 1−1 + +( q －第二类条件 T = λ λ M1 hiδ x +1 λ hT f －第三类条件
δ xiΔxi S δ x TM 1−1 + + A λ λ h －第三类条件 TM 1 = B iδ x +1 λ 此式适合于两种边界条件
0
－第二类条件
22/62
（3）区域离散方法B 区域离散方法B的边界节点的控制容积为零，只要令 以上公式中的
Δx = 0 即得方法B中的相应公式：
0 TM1−1 −TM1 qB +λ + ΔxiS = 0 δx qB iδ x 第二类边界－ TM 1 = TM 1−1 + λ hiδ x )T f TM 1−1 + ( λ TM 1 = 第三类边界－ 控制容积为零 hiδ x 1+ λ 离散方法B公式具有二阶截差的精度。
边界节点不同截差离散对数值计算结果的影响
26/62
离散方法B，内部取3点， 内部结点 T2 , T3 , T4 用二阶截差离散，边界节点 T5 按前述公式 TM 1 = TM 1−1 +
数值传热学陶文铨课后答案第三章
数值传热学陶文铨课后答案第三章A.陶文铨,B.周伟俊D.黄仁龙，黄国生，周伟俊，周国生，周大白。
在计算流体力学中，流场是一个基本问题。
计算流体力学中常用到的主要方法是()。
一、对于固体，流场和速度的分布问题的计算是通过计算流体的物理量来进行的。
在计算流体力学中，流场和速度的分布问题的计算是通过计算流体的物理量来进行的。
A.质量流场 F:计算流体力学中流场和速度分布问题方法，其目的是用物理量来表示（质量）。
D.速度分布方程 F:计算流体力学中非稳态速度分布方程。
二、对于固体流动，常用的方法有流线理论法、流体压力分析法、边界层理论法、粒子质量守恒定律法、自由流动技术法、流动力学理论以及流体力学仿真技术。
解析：本题考查的是流线理论中的流场理论。
固体流动方向上分为轴流和向后流场。
在轴流流动中，固体的扩散与温度成正比；向后流一般不流动。
三、流体粘性为固液相物之间的粘滞系数，且在固态中较为稳定，因此研究固液相物之间结合力和粘滞系数问题的方法通常称为粘性流场方法。
【参考答案】 C解析：流体粘性指流体相物之间形成的粘性，根据粘性对固液相物有粘滞力或阻力，流体粘性小于固体流体之间粘滞力或阻力时固体粘滞力>固体之间粘性；流体粘性大于固体之间粘性时固体粘滞力>固体之间粘性。
固液相混是指固液相混在一起。
流体特性参数是用于定量计算和评价流体性能的重要参数。
流体特性参数如流速、温度、压力、密度、粘度等。
A选项; D选项四、流体流动分为粘性流动和非粘性流动。
A.粘性流动是指流体在流动过程中具有一定的粘性而非粘性，其流动过程主要分为四个阶段：B.非粘性流动是指流体在流动过程中不具有任何粘性，其流场主要分为湍流、不稳定、流动速度等多种形式。
C.湍流是指流体在湍流状态下不能充分发挥其流动能力，其流场主要分为两种形式：湍流和不稳定。
D.不稳定是指流体在流动过程中没有任何不稳定状态，其流场主要分为两种形式：一种是无规则流动（在流动过程中无规则流动),另一种是不规则流动时无规则流动）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3-7证明对流项的背风差分总使扰动逆流而传递。
证明：Taylor 展开法中逆风差分的构造法：
1,i i i x x φφφ+-∂=∂∆ u>0 1,i i i x x
φφφ--∂=∂∆ u<0 下面以u>0的情形来分析.对于节点i+1,在n 时层产生在节点i 的扰动对i+1的影响由下式确定: 11112n n n n i i i i u t x φφφφ+++++--=-∆∆ (1n i φ+=0,2n i φ+=0)
由此得 11n i φ++=0
而i-1处则有 1111n n n n i i i i u t x φφφφ+-----=-∆∆ （1n i φ-=0）
得 11n i u t x φε+-∆⎛⎫= ⎪∆⎝⎭
因此可知对流项的背风差分总使扰动逆流而传递。
3-10一阶导数的而二阶差分格式称为二阶迎风格式（在来流方向区节点构成差分格式）。
试分析其迁移性。
解：经查表2-1可知在来流方向区节点的一阶导数二阶迎风格式为：
n n n i i-1i-2i n 34=x 2x φφφφ
-+∂∂∆， u>0 下面以u>0的情形来分析.对于节点i+1,在n 时层产生在节点i 的扰动对i+1的影响由下式确定: n+1n n n n i+1i+1
i+1i i-1n+1i+134=-u t 2x 2u t =x φφφφφφε--+∆∆∆⎛⎫ ⎪∆⎝⎭ （n i+1φ=0，n i-1φ=0） 得 n+1i+12u t =x φε∆⎛⎫ ⎪∆⎝⎭
而i-1处则有 n n n n+1n i-1i-2i-3i-1i-134=-u t
2x φφφφφ-+-∆∆ （n i-1φ=0，n i-2φ=0，n i-3φ=0） 因此得
n+1i-1φ=0
因此可知一阶导数的二阶迎风格式（在来流方向区节点构成差分格式）具有迁移性。
扰动只向后传动！！！。