第二章计算流体力学的基本知识

合集下载

第二章计算流体力学的基本知识

第二章计算流体力学的基本知识

第二章计算流体力‎学的基本知‎识流体流动现‎象大量存在‎于自然界及‎多种工程领‎域中,所有这些工‎程都受质量‎守恒、动量守恒和‎能量守恒等‎基本物理定‎律的支配。

这章将首先‎介绍流体动‎力学的发展‎和流体力学‎中几个重要‎守恒定律及‎其数学表达‎式,最后介绍几‎种常用的商‎业软件。

2.1 计算流体力‎学简介2.1.1计算流体‎力学的发展‎流体力学的‎基本方程组‎非常复杂,在考虑粘性‎作用时更是‎如此,如果不靠计‎算机,就只能对比‎较简单的情‎形或简化后‎的欧拉方程‎或N-S方程进行‎计算。

20世纪3‎0~40年代,对于复杂而‎又特别重要‎的流体力学‎问题,曾组织过人‎力用几个月‎甚至几年的‎时间做数值‎计算,比如圆锥做‎超声速飞行‎时周围的无‎粘流场就从‎1943年‎一直算到1‎947年。

数学的发展‎,计算机的不‎断进步,以及流体力‎学各种计算‎方法的发明‎,使许多原来‎无法用理论‎分析求解的‎复杂流体力‎学问题有了‎求得数值解‎的可能性,这又促进了‎流体力学计‎算方法的发‎展,并形成了"计算流体力‎学"。

从20世纪‎60年代起‎,在飞行器和‎其他涉及流‎体运动的课‎题中,经常采用电‎子计算机做‎数值模拟,这可以和物‎理实验相辅‎相成。

数值模拟和‎实验模拟相‎互配合,使科学技术‎的研究和工‎程设计的速‎度加快,并节省开支‎。

数值计算方‎法最近发展‎很快,其重要性与‎日俱增。

自然界存在‎着大量复杂‎的流动现象‎,随着人类认‎识的深入,人们开始利‎用流动规律‎来改造自然‎界。

最典型的例‎子是人类利‎用空气对运‎动中的机翼‎产生升力的‎机理发明了‎飞机。

航空技术的‎发展强烈推‎动了流体力‎学的迅速发‎展。

流体运动的‎规律由一组‎控制方程描‎述。

计算机没有‎发明前,流体力学家‎们在对方程‎经过大量简‎化后能够得‎到一些线形‎问题解析解‎。

但实际的流‎动问题大都‎是复杂的强‎非线形问题‎,无法求得精‎确的解析解<。

第二讲 流体力学的基本知识

第二讲 流体力学的基本知识
◆ 液压传动中的压力就是我们物理学中压强的概 念。是指当液体相对静止时,液体单位面积上所受 的法向力。 ◆ P=F/A ◆ 单位是帕。1MPa=_____Pa ◆ 2.压力是怎样产生的? ◆ 我们可以这么认为,液体或气体受到挤压会膨胀产 生压力。
密闭容器内的流体的特点 密闭容器中液体各点的压力是相等的。 密闭液体可以用于管路中向各个方向传递动力。
2.如题图所示连通器,中间有一活动隔板T,已知活塞面 积A1=1×10-3 m2, A2=5×10-3 m2,F1=200N,G=2500N, 活塞自重不计,问: (1)当中间用隔板T隔断时,连通器两腔压力P1、P2各是 多少? (2)当把中间隔板抽去,使连通器连通时,两腔压力P1、 P2各是多少?力F1能否举起重物G? (3)当抽去中间隔板T后若要使两活塞保持平衡,F1应是 多少?
液压油
学习目标 1.理解掌握液压油的性质 2.掌握液压油的类型 3.能够选用正确的液压油 4.学会分析液压油的故障
液体是液压传动的工作介质。最常用的工作介 质是液压油。 ◆ 1.液压油的性质 ◆ 1)密度 M M-液体的质量; V V-液体的体积。 一般液压油的密度为:900kg/m3 ◆ 2)可压缩性 ◆ 指液体在外力作用下体积减小的特性; ◆ 一般认为油液是不可压缩的。
料脱落的颗粒和纤维剥落的油漆、碎渣等。
• • • •
2.油液污染的危害 污染物包括:金属材料75%、尘埃15%、其它10% 1)对油泵的危害:使油泵润滑部分磨损加剧。 2)对液压阀的危害:使阀心移动困难或卡住阀口 密封不严,使阀失去控制性能。 • 3)对油缸危害:加速密封的损坏,油缸内表面拉 伤,内外泄露增加。 • 4)对过滤器的危害:会使滤网阻塞,油泵吸油困 难,回油不畅。严重时击穿滤心。 • 5)油液变质降低油液原有的特性和使用期。

(完整版)流体力学重点概念总结

(完整版)流体力学重点概念总结

第一章绪论表面力:又称面积力,是毗邻流体或其它物体,作用在隔离体表面上的直接施加的接触力。

它的大小与作用面积成比例。

剪力、拉力、压力质量力:是指作用于隔离体内每一流体质点上的力,它的大小与质量成正比。

重力、惯性力流体的平衡或机械运动取决于:1.流体本身的物理性质(内因)2.作用在流体上的力(外因)流体的主要物理性质:密度:是指单位体积流体的质量。

单位:kg/m3 。

重度:指单位体积流体的重量。

单位: N/m3 。

流体的密度、重度均随压力和温度而变化。

流体的流动性:流体具有易流动性,不能维持自身的形状,即流体的形状就是容器的形状。

静止流体几乎不能抵抗任何微小的拉力和剪切力,仅能抵抗压力。

流体的粘滞性:即在运动的状态下,流体所产生的阻抗剪切变形的能力。

流体的流动性是受粘滞性制约的,流体的粘滞性越强,易流动性就越差。

任何一种流体都具有粘滞性。

牛顿通过著名的平板实验,说明了流体的粘滞性,提出了牛顿内摩擦定律。

τ=μ(du/dy)τ只与流体的性质有关,与接触面上的压力无关。

动力粘度μ:反映流体粘滞性大小的系数,单位:N•s/m2运动粘度ν:ν=μ/ρ第二章流体静力学流体静压强具有特性1.流体静压强既然是一个压应力,它的方向必然总是沿着作用面的内法线方向,即垂直于作用面,并指向作用面。

2.静止流体中任一点上流体静压强的大小与其作用面的方位无关,即同一点上各方向的静压强大小均相等。

静力学基本方程: P=Po+pgh等压面:压强相等的空间点构成的面绝对压强:以无气体分子存在的完全真空为基准起算的压强 Pabs相对压强:以当地大气压为基准起算的压强 PP=Pabs—Pa(当地大气压)真空度:绝对压强不足当地大气压的差值,即相对压强的负值 PvPv=Pa-Pabs= -P测压管水头:是单位重量液体具有的总势能基本问题:1、求流体内某点的压强值:p = p0 +γh;2、求压强差:p – p0 = γh ;3、求液位高:h = (p - p0)/γ平面上的净水总压力:潜没于液体中的任意形状平面的总静水压力P,大小等于受压面面积A与其形心点的静压强pc之积。

流体力学基础知识

流体力学基础知识

第一章,绪论1、质量力:质量力是作用在流体的每一个质点上的力。

其单位是牛顿,N。

单位质量力:没在流体中M点附近取质量为d m的微团,其体积为d v,作用于该微团的质量力为dF,则称极限lim(dv→M)dF/dm=f,为作用于M点的单位质量的质量力,简称单位质量力。

其单位是N/kg。

2、表面力:表面力是作用在所考虑的或大或小得流体系统(或称分离体)表面上的力。

3、容重:密度ρ和重力加速度g的乘积ρg称容重,用符号γ表示。

4、动力黏度μ:它表示单位速度梯度作用下的切应力,反映了黏滞性的动力性质。

其单位为N/(㎡·s),以符号Pa·s表示。

运动黏度ν:是单位速度梯度作用下的切应力对单位体积质量作用产生的阻力加速度。

国际单位制单位㎡/s。

动力黏度μ与运动黏度ν的关系:μ=ν·ρ。

5、表面张力:由于分子间的吸引力,在液体的自由表面上能够承受的极其微小的张力称为表面张力。

毛细管现象:由于表面张力的作用,如果把两端开口的玻璃细管竖立在液体中,液体就会在细管中上升或下降h高度的现象称为毛细管现象。

6、流体的三个力学模型:①“连续介质”模型;②无黏性流体模型;③不可压缩流体模型。

(P12,还需看看书,了解什么是以上三种模型!)。

第二章、流体静力学1、流体静压强的两个特性:①其方向必然是沿着作用面的内法线方向;②其大小只与位置有关,与方向无关。

2、a流体静压强的基本方程式:①P=Po+rh,式中P指液体内某点的压强,Pa(N/㎡);Po指液面气体压强,Pa(N/㎡);r指液体的容重,N/m³;h指某点在液面下的深度,m;②Z+P/r=C(常数),式中Z指某点位置相对于基准面的高度,称位置水头;P/r指某点在压强作用下沿测压管所能上升的高度,称压强水头。

两水头中的压强P必须采用相对压强表示。

b流体静压强的分布规律的适用条件:只适用于静止、同种、连续液体。

3、静止均质流体的水平面是等压面;静止非均质流体(各种密度不完全相同的流体——非均质流体)的水平面是等压面,等密度和等温面。

(完整版)流体力学知识点总结汇总

(完整版)流体力学知识点总结汇总

流体力学知识点总结 第一章 绪论1 液体和气体统称为流体,流体的基本特性是具有流动性,只要剪应力存在流动就持续进行,流体在静止时不能承受剪应力。

2 流体连续介质假设:把流体当做是由密集质点构成的,内部无空隙的连续体来研究。

3 流体力学的研究方法:理论、数值、实验。

4 作用于流体上面的力(1)表面力:通过直接接触,作用于所取流体表面的力。

作用于A 上的平均压应力作用于A 上的平均剪应力应力法向应力切向应力(2)质量力:作用在所取流体体积内每个质点上的力,力的大小与流体的质量成比例。

(常见的质量力:重力、惯性力、非惯性力、离心力)单位为5 流体的主要物理性质 (1) 惯性:物体保持原有运动状态的性质。

质量越大,惯性越大,运动状态越难改变。

常见的密度(在一个标准大气压下): 4℃时的水20℃时的空气(2) 粘性ΔFΔPΔTAΔAVτ法向应力周围流体作用的表面力切向应力A P p ∆∆=A T ∆∆=τAF A ∆∆=→∆lim 0δAPp A A ∆∆=→∆lim 0为A 点压应力,即A 点的压强 ATA ∆∆=→∆lim 0τ 为A 点的剪应力应力的单位是帕斯卡(pa ),1pa=1N/㎡,表面力具有传递性。

B Ff m =2m s 3/1000mkg =ρ3/2.1mkg =ρ牛顿内摩擦定律: 流体运动时,相邻流层间所产生的切应力与剪切变形的速率成正比。

即以应力表示τ—粘性切应力,是单位面积上的内摩擦力。

由图可知—— 速度梯度,剪切应变率(剪切变形速度) 粘度μ是比例系数,称为动力黏度,单位“pa ·s ”。

动力黏度是流体黏性大小的度量,μ值越大,流体越粘,流动性越差。

运动粘度 单位:m2/s 同加速度的单位说明:1)气体的粘度不受压强影响,液体的粘度受压强影响也很小。

2)液体 T ↑ μ↓ 气体 T ↑ μ↑ 无黏性流体无粘性流体,是指无粘性即μ=0的液体。

无粘性液体实际上是不存在的,它只是一种对物性简化的力学模型。

气体动力学讲义吴子牛第二章流体力学问题模拟的基础知识

气体动力学讲义吴子牛第二章流体力学问题模拟的基础知识

2 V
同理,紊动阻力也必须满足牛顿相似定律,于是得到:
于是得到:
L2
2 V
L2
2 V
1
2.2模型相似准则
上式表明,如果两个几何相似的水流是在紊动阻力 作用下达成动力相似,则它们的沿程阻力系数相等;反 之,如果两个水流的沿程阻力系数相等,则这两个水流 一定是早紊动阻力作用下动力相似的。对于气流流动, 结论也是相同的。这就是紊动阻力相似准则。
2.2模型相似准则
一、重力相似准则(佛汝德相似准则)
如果两个相似水流中起主导作用的是重力,作用于原型和模
型的相应部分的重力分别以Gy和Gm表示,由于 G gW,因 而重力的比尺为:
G
Gy Gm
y g yWy m g mWm
g
3 L
要使在重力作用下原型与模型相似,同样必须满足动力相似的
一般规律,而由于作用力F中仅考虑重力G,因而F=G,即。
2.2模型相似准则
在水力学中,根据重力相似准则设计的模型实验有堰流、 孔口管嘴泄流以及流过水工建筑物等以重力为主导的水流 运动。对于波浪运动、船舶和水上飞机浮筒等水上运动物 体的波浪阻力问题,也必须根据重力相似准则设计模型。 在空气动力学中,对于具有加速度的运动物体,例如飞机 的机动飞行,设计模型实验也要满足重力相似准则,才能 保证机动飞行姿态相似。
const.
W
Wy Wm
3 L
const.
2.1 流动的力学相似
从理论上讲,最好能做到所有模型尺寸全按一个比尺 缩小或放大,这种长、宽、高比尺均一致的模型称为正 态模型。在流体力学模型实验中,通常遇到的是这类模 型。在有的情况下,不能做到这一点,例如进行天然河 道流动的模型实验,由于天然河道的长度比宽度和水深 要大得多,如果按照同一比例尺缩制模型,势必造成水 深太小甚至改变了模型中水流的性质。对于这种情况, 就要分别采用不同的长度比尺、宽度比尺和高度比尺, 因而这种模型就改变了原有的形状。这种比尺不一样的 模型称为变态模型。

流体力学基本知识

流体力学基本知识
粘滞性的大小随流体的种类及所处的外界 条件而不同。例如流体中的水与重油, 温度的高低都影响其粘滞性的变化。
牛顿试验研究提出与粘滞性有关的内摩擦 定律为
李峥嵘 博士
9
2019/9/22
李峥嵘 博士
10
2019/9/22
李峥嵘 博士
11
பைடு நூலகம்
2019/9/22
四、计量单位
1、国际单位
1)基本单位:长度、质量、时间、热力
李峥嵘 博士
43
2019/9/22
3、沿程损失和局部损失 1)沿程损失
流体流动中为克服摩擦阻力而损耗的能量
称为沿程损失。沿程阻力损失与长度、粗糙 度及流速的平方成正比,而与管径成反比, 通常采用达西一维斯巴赫公式计算:
李峥嵘 博士
44
2019/9/22
2)局部损失 流体运动过程中,通过断面变化处、转向 处、分支或其他使流体流动情况改变时,
阻力损失值视流体的流行形态而 不同,因此计算流体的阻力损
失.应了解水流的形态。
李峥嵘 博士
40
2019/9/22
在一一端、装有流阀体门的的流长玻动璃形态
管放色中水水流充,流种体满 并 ,不在水 由 则同, 小 可流的稍 管 见动流开 注 管中动启 入 内,阀有颜形由门颜色态于—流—速层不流同和而紊呈流现。出两
传递的能量,是物体间(内)通过分子 运动相传递的能量。给物体加热,实
际就是增加使物体分子运动的能量,物 体的温度就将升高,反之使物体散热减 小分子运动的能量,物体温度降低。
李峥嵘 博士
49
3、温度
2019/9/22
李峥嵘 博士
50
4、热膨胀
2019/9/22

第二章--计算流体力学的基本知识

第二章--计算流体力学的基本知识

第二章--计算流体力学的基本知识第二章计算流体力学的基本知识流体流动现象大量存在于自然界及多种工程领域中,所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。

这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式,最后介绍几种常用的商业软件。

2.1计算流体力学简介2.1.1计算流体力学的发展流体力学的基本方程组非常复杂,在考虑粘性作用时更是如此,如果不靠计算机,就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。

20世纪30~40年代,对于复杂而又特别重要的流体力学问题,曾组织过人力用几个月甚至几年的时间做数值计算,比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。

数学的发展,计算机的不断进步,以及流体力学各种计算方法的发明,使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性,这又促进了流体力学计算方法的发展,并形成了"计算流体力学"。

从20世纪60年代起,在飞行器和其他涉及流体运动的课题中,经常采用电子计算机做数值模拟,这可以和物理实验相辅相成。

数值模拟和实验模拟相互配合,使科学技术的研究和工程设计的速度加快,并节省开支。

数值计算方法最近发展很快,其重要性与日俱增。

自然界存在着大量复杂的流动现象,随着人类认识的深入,人们开始利用流动规律来改造自然界。

最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。

航空技术的发展强烈推动了流体力学的迅速发展。

*流体运动的规律由一组控制方程描述。

计算机没有发明前,流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。

但实际的流动问题大都是复杂的强非线形问题,无法求得精确的解读解。

计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现,从而催生了计算流体力学这门交叉学科。

计算流体力学是一门用数值计算方法直接求解流动主控方程(Euler 或Navier-Stokes方程)以发现各种流动现象规律的学科。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第二章计算流体力学的基本知识流体流动现象大量存在于自然界及多种工程领域中,所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。

这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式,最后介绍几种常用的商业软件。

2.1 计算流体力学简介2.1.1计算流体力学的发展流体力学的基本方程组非常复杂,在考虑粘性作用时更是如此,如果不靠计算机,就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。

20世纪30~40年代,对于复杂而又特别重要的流体力学问题,曾组织过人力用几个月甚至几年的时间做数值计算,比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。

数学的发展,计算机的不断进步,以及流体力学各种计算方法的发明,使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性,这又促进了流体力学计算方法的发展,并形成了"计算流体力学"。

从20世纪60年代起,在飞行器和其他涉及流体运动的课题中,经常采用电子计算机做数值模拟,这可以和物理实验相辅相成。

数值模拟和实验模拟相互配合,使科学技术的研究和工程设计的速度加快,并节省开支。

数值计算方法最近发展很快,其重要性与日俱增。

自然界存在着大量复杂的流动现象,随着人类认识的深入,人们开始利用流动规律来改造自然界。

最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。

航空技术的发展强烈推动了流体力学的迅速发展。

流体运动的规律由一组控制方程描述。

计算机没有发明前,流体力学家们在对方程经过大量简化后能够得到一些线形问题解析解。

但实际的流动问题大都是复杂的强非线形问题,无法求得精确的解析解。

计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现,从而催生了计算流体力学这门交叉学科。

计算流体力学是一门用数值计算方法直接求解流动主控方程(Euler或Navier-Stokes方程)以发现各种流动现象规律的学科。

它综合了计算数学、计算机科学、流体力学、科学可视化等多种学科。

广义的CFD包括计算水动力学、计算空气动力学、计算燃烧学、计算传热学、计算化学反应流动,甚至数值天气预报也可列入其中。

自20世纪60年代以来,CFD技术得到飞速发展,其原动力是不断增长的工业需求,而航空航天工业自始至终是最强大的推动力。

传统飞行器设计方法实验昂贵、费时,所获信息有限,迫使人们需要用先进的计算机仿真手段指导设计,大量减少原型机实验,缩短研发周期,节约研究经费。

四十年来,CFD在湍流模型、网格技术、数值算法、可视化、并行计算等方面取得飞速发展,并给工业界带来了革命性的变化。

如在汽车工业中,CFD和其它计算机辅助工程(CAE)工具一起,使原来新车研发需要上百辆样车减少为目前的十几辆车;国外飞机厂商用CFD取代大量实物实验,如美国战斗机YF-23采用CFD进行气动设计后比前一代YF-17减少了60%的风洞实验量。

目前在航空、航天、汽车等工业领域,利用CFD进行的反复设计、分析、优化己成为标准的必经步骤和手段。

当前CFD问题的规模为:机理研究方面如湍流直接模拟,网格数达到了109(十亿)量级,在工业应用方面,网格数最多达到了107(千万)量级。

与实验研究相比,理论计算具有花费少、速度快、信息完整、模拟能力强等优点,特别是大量的计算流体力学软件的出现,大大减少了计算流体力学研究的工作量,从而扩大了计算流体力学的应用范围,推动了流体力学更深入的发展。

计算流体力学还不是一项很成熟的技术,在用计算流体力学对流动现象进行预测的时候,需要对复杂的流动现象进行处理,然后用数学模型来描述它,计算的结果既取决于计算方法,也取决于数学模型本身,如果数学模型的描述不够精确,甚至不恰当,其计算结果也就没有任何价值可言。

尽管作为一门新兴的学科,计算流体力学还有缺陷,但它会随着技术的进步和发展而日趋成熟,并将在化工领域得到广泛的应用。

一个完整的计算流体力学模型应包含如下几个方面的内容: 本构方程,即流体力学基本方程:连续性方程(质量方程)、动量方程、能量方程、状态方程等。

湍流模型,不同于层流,必须考虑流体单元的脉动速度,脉动是湍流流动的基本特征。

从模型的建立及求解过程可以看出,其实质是寻找出由于脉动而起的运动粘度的表达式。

多相流模型,对于多相流模拟计算来说,基本的湍流模型还不够用,需要进一步寻找各相运动规律及相间作用力规律。

模型的求解数值方法,对模型进行计算时,需要选择好的差分格式、松弛因子、时间步长等,以使结果收敛尽量减少CPU运算时间。

2.1.2 计算流体力学的定义计算流体动力学(Computational Fluid Dynamics ,简称CFD)是建立在经典流体力学与数值计算方法基础上的新型独立的学科,通过计算机数值计算和图像显示的方法,在时间和空间上定量描述流场的数值解,从而达到对物理问题研究的目的。

它兼有理论性和实践性的双重特点,建立了理论和方法,为现代科学中许多复杂流动和传热问题提供了有效的计算技术。

计算流体动力学(CFD)是通过计算机数值计算和图像显示,对包含有流体流动和热传导等相关物理现象的系统所做的分析。

它的基本思想是:把原来在时间域及空间域上连续的物理量的场,如速度场和压力场,用一系列有限个离散点上的变量值的几何来代替,通过一定的原则和方式建立起来的关于这些离散点上场变量之间关系的代数方程组,然后代数方程组获得场变量的近似值[5]。

CFD方法和传统的理论分析方法、实验测量方法组成了研究流体流动问题的完整体系。

理论分析方法的优点在于所得结果具有普遍性,各种影响因素清晰可见,是指导实验研究和验证数值计算方法的理论基础,但是它往往要求对计算进行抽象和简化,才可能得出理论解。

对于非线性情况,只有少数流动才能给出解析结果。

实验测量方法所得到的实验结果真实可信,它是理论分析和数值方法的基础,其重要性不容低估。

然而,实验往往受到模型尺寸、流场流动、人身安全和测量精度的限制,有时可能很难通过试验的方法得到满意的结果。

而CFD方法恰好克服了前面两种方法的弱点,在计算机上实现一个特定的计算,就好像在计算机上做一个物理实验。

例如,机翼的绕流,通过计算机并将其结果在屏幕上显示,就可以看到流场的各种细节:如激波的运动、强度,涡的生成与传播,流动的分离、表面的压力分布、受力大小及其随时间的变化等。

数值模拟可以形象地再现流动情景,与做实验没有什么区别。

2.1.3 计算流体力学的计算步骤采用CFD的方法对流体流动进行数值模拟,通常包括如下步骤:(1)建立反映工程问题或物理问题本质的数学模型。

具体的说就是要建立反映问题各个量之间关系的微分方程及相应的定解条件,这是数学模型的出发点。

没有正确完善的数学模型,数值模拟就没有任何意义。

流体的基本控制方程通常包括质量守恒方程、动量守恒方程、能量守恒方程,以及这些方程相应的定解条件。

(2)寻求高效率、高准确度的计算方法,即建立针对控制方程的数值离散化方法,如有限差分法、有限元法、有限体积法等。

这里的计算方法不仅包括微分方程的离散化方法及求解条件,还包括体坐标的建立,边界条件的处理等。

这些内容可以说是CFD的核心。

(3)编制程序和进行计算。

这部分工作包括计算网格划分、初始条件和边界条件的输入,控制参数的设定等。

这是整个工作中花时间最多的部分。

由于求解的问题比较的复杂,比如Navier-Stokes方程就是一个十分复杂的非线性方程,数值求解方法在理论上不是绝对完善的,所以需要通过实验加以验证。

正是从这个意义上讲,数值模拟又叫数值实验。

(4)显示实验的结果,计算结果一般通过图表等方式显示,这对检查和判断分析质量和结果有重要的意义。

2.1.4 计算流体力学的局限性虽然CFD具有许多的优点,但是也存在一定的局限性。

首先,数值解法是一种离散近似的计算方法,依赖于物理上合理、数学上适用、适合在计算机上进行计算的离散的数学模型,且最终结果不能提供任何形式的解析表达式,只是有限个离散点上的数值解,并有一定的计算误差;第二,它不像物理模型试验一开始就能给出流动现象并定性的描述,往往需要由原体观测或物理模型试验提供某些流动参数,并需要对建立的数学模型进行验证;第三,程序的编制及资料的收集、整理与正确利用,在很大程度上取决于经验和技巧。

此外,因数值处理方法等原因有可能导致计算结果的不真实,例如产生数值粘性和频散等伪物理效应。

当然,某些缺点或局限性可以通过某种方式克服或弥补。

最后,CFD因涉及大量的数值计算,因此,需要较高的计算机软硬件配置。

2.1.5 几种数值解法经过四十多年的发展,CFD出现了多种数值解法。

这些方法之间的主要区别在于对控制方程的离散方式。

根据离散的原理不同,大体上可以分为三个分支:有限差分法、有限元法、有限体积法。

有限差分法是运用最早、最经典的CFD方法,它将求解域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。

求出差分方程组的解,就是微分方程定解问题的数值近似解。

它是一种直接将微分问题变成代数问题的近似数值解法。

这种方法发展较早,比较成熟,较多的用于求解双曲型和抛物型问题。

在此基础上发展起来的方法有PIC(Particle-in-Cell)法、MAC(Marker-and-Cell)法,以及由美籍华人学者陈景仁提出的有限分析法(Finite-Analytic-Method)等[6]。

有限元法是20世纪80年代开始应用的一种数值解法,它吸收了有限差分法中离散处理的内核,又采用了变分计算中选择逼近函数对区域进行积分的合理方法。

有限元法因求解速度较有限差分法和有限体积法慢,因此应用不是很广泛。

在有限元法的基础上,英国C.A.Brebbia等提出了边界元法和混合元法等方法。

有限体积法是将计算区域划分为一系列控制体积,将待解微分方程对每一个控制体积进行积分,得出离散方程。

有限体积法的关键是在导出离散方程过程中,需要对界面上的被求函数本身及导数的分布做出某种形式的假定。

用有限体积法导出的离散方程可以保证具有守恒特性,而且离散方程系数物理意义明确,计算量相对较小。

它是目前CFD应用最广的一种方法。

当然这种方法的研究和扩展也在不断的进行,有的学者提出了适用于任意多边形非结构网格的扩展有限体积法[7]。

2.2 流体动力学控制方程2.2.1 流体的质量守恒方程任何流体问题都必须满足质量守恒定律。

该定律可表达为:单位时间内流体微元体中质量的增加,等于同一时间间隔内流入该微元的净质量。

按照这一定律,可以得出质量守恒方程(mass conservation equation )[9]:0)()()(=∂∂+∂∂+∂∂+∂∂tw t v t u t ρρρρ (2.2) 引入矢量符号div(a)=za y a x a z y x ∂∂+∂∂+∂∂,则上式写成: 0)(=+∂∂u div tρρ (2.3) 有的文献使用符号∇表示散度,即∇•a= div(a)=za y a x a z y x ∂∂+∂∂+∂∂,这样,上式又可以写成:0)(=∇+∂∂u tρρ (2.4)上式中:ρ是密度,t 是时间,u 是速度矢量,u 、v 、w 是速度矢量在x 、y 、z 方向的分量。

相关文档
最新文档