流体力学-流体仿真介绍(共15页).doc
《流体力学-流体仿真介绍(共15页).doc》由会员分享,可在线阅读,更多相关《流体力学-流体仿真介绍(共15页).doc(15页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上流体力学研究流体(气体与液体)的宏观运动与平衡,它以流体宏观模型作为基本假说。 显然,流体的运动取决于每个粒子的运动,但若求解每个粒子的运动即不可能也无必要。对于宏观问题,必须在微观与宏观之间建立一座桥梁。 流体宏观模型认为流体是由无数流体元(或称流体微团)连续地组成的(即连续介质)。所谓流体元指的是这样的小块流体:它的大小与放置在流体中的实物比较是微不足道的,但比分子的平均自由程却要大得多,它包含足够多的分子,能施行统计平均求出宏观参量,少数分子出入于流体元不会影响稳定的平均值。 另一方面,对于进行统计平均的时间也应选得足够大,使得在这段时间内,微观的性质,例如分子
2、间的碰撞等已进行了许多次,在这段时间内进行统计平均能够得到稳定的数值。于是,从统计物理中得知,分子的物理量(质量、速度、动量和能量)经过统计平均后变成了流体元的质量,速度,压力和温度等宏观物理量,分子质量、动量和能量等输运过程,经过统计平均后表现为扩散,粘性,热传导等宏观性质。 上述微观上充分大、宏观上充分小的流体元称为流体质点,将流体运动的空间看作是由流体质点连续地无空隙地充满着的假设称为连续介质假设。应该指出,有了此假设才能把一个微观问题化成宏观问题,且数学上容易处理。实验和经验也表明在一般情况下这个假设总是成立的。 但是。在某些特殊问题中,连续介质的假设也可以不成立。例如在稀薄气体力学中
3、,分子间的距离很大,它能和物体的特征尺度比拟,这样虽然获得稳定平均值的流体元还是存在的,但是不能将它看成一个质点。又如考虑激波内的气体运动,激波的尺寸与分子平均自由程同阶,激波内的流体只能看成分子而不能当作连续介质来处理了。1 流体的基本性质1.1 易流动性 流体在静止时不能承受切向应力,不管多小的切向应力,都会引起其中各流体元彼此间的相对位移,而且取消力的作用后,流体元之间并不恢复其原有位置。正是流体的这一基本特性使它能同刚体和弹性体区别开来。刚体和弹性体也是连续介质,但是刚体中质点之间的相互距离不论其上作用的外力如何将保持不变;而在弹性体中,当作用力在数值上达到某一界限时,系统中各点间的相
4、互距离可以改变,但消除了力的作用之后,各点相互关系又恢复原有状态。相反地,流体能够有任意大的变形。因此流体在静止时只有法应力而没有切应力。流体的这个宏观性质称为易流动性。1.2 粘性 流体在静止时虽不能承受切应力,但在运动时,对相邻两层流体间的相对运动即相对滑动速度是有抵抗的,这种抵抗力称为粘性应力,流体所具有的这种抵抗两层流体相对滑动的性质称为粘性,粘性大小依赖于流体的性质,并显著地随温度而变化。实验表明,粘性应力的大小与粘性及相对速度成正比。 当流体的粘性较小,运动的相对速度也不大时,所产生的粘性应力比起其它类型的力(如惯性力)可忽略不计。此时,我们可以近似地把流体看成是无粘性的,这样的流
5、体称为理想流体。十分明显,理想流体对于切向变形没有任何抗拒能力。这样对于粘性而言,我们可以将流体分成理想流体和粘性流体两大类。应该强调指出,真正的理想流体在客观实际中是不存在的。它只是客观流体在某种条件下的一种近似模型。 除了粘性外,流体还有热传导及扩散等性质。 流体的宏观性质,扩散,粘性,热传导等是分子输运性质的统计平均。由于分子的不规则运动,在各层流体间将交换着质量,动量和能量,使不同流体层内的平均物理量均匀化,这种性质称为分子运动的输运性质。质量输运在宏观上表现为扩散现象,动量输运表现为粘性现象,能量输运则表现为热传导现象。1.3 压缩性 流体质点的体积或密度在受到一定压力或温度差的条件
6、下可以改变,这个性质称为压缩性。真实流体都是可以压缩的。它的压缩程度依赖子流体的性质及外界的条件。液体在通常的压力或温度下,压缩性很小。因此在一般情形下液体可以近似地看成是不可压缩的。2 描写流体运动的两种方法2.1 拉格朗日方法(随体法) 在拉格朗日方法中,注意的中心即着眼点是流体质点,确定所有流体质点的运动规律,即它们的位置随时间变化的规律。十分明显,如果知道了所有流体质点的运动规律,那么整个流体运动的状况也就清楚了。 现在我们将描写运动的观点和方法用数学式子表达出来,为此首先必须用某种数学方法区别不同的流体质点。通常利用初始时刻流体质点的坐标作为区分不同流体质点的标志。设初始时刻 t =
7、 t0 时,流体质点的坐标是 a,b,c,它可以是曲线坐标,也可以是直角坐标,重要的是给流体质点以标号而不在于采取什么具体的方式。 我们约定采用 a,b,c 三个数的组合来区别流体质点,不同的 a,b,c 代表不同的质点,于是流体质点的运动规律可表为下列矢量形式: r=r(t,a,b,c) 其中 r 是流体质点的矢径。在直角坐标系中,有分量式: x=x(t,a,b,c) y=y(t,a,b,c) z=z(t,a,b,c) 变数 t; a,b,c 称为拉格朗日变数。2.2 欧拉方法(当地法) 欧拉方法不直接考虑个别流体质点如何运动,而是用场的观点研究流体运动。它只集中注意力于那些发生在空间给定点
8、的流动情况;对于流体质点从什么地方和如何在给定时刻达到这一点,经过这点以后又会运行到别的什么地方和怎样运行到那些地方的,这一切问题从欧拉方法观点看来并不是基本的。这样,欧拉方法是把空间某一固定点 (x, y, z) 的流体质点的速度当作时间的函数来研究的;显然,这个速度也是坐标 (x, y, z) 的函数。因此, 其分量为: ux=ux(t,x,y,z) uy=uy(t,x,y,z) uz=uz(t,x,y,z) 变数 t; x, y, z 称为欧拉变数。如果在上式中把 t 当作可变的,而把 x, y, z 当作常数,则对不同的 t 我们得到不同时刻经过空间中确定点的不同流体质点的速度;而如把
9、 t 当作常数, x, y, z 当作变数,则可得到对于确定时刻空间中流体质点的速度分布。 由于上式确定的速度函数是定义在空间点上的,它们是空间点坐标 x, y, z 的函数,所以我们研究的是场,如速度场等。因此当我们采用欧拉观点描述运动时,就可以利用场论的知识。若场内函数不依赖矢径 r 则称之为均匀场,反之称之为非均匀场;若场内函数不依赖时间则称为定常场,反之称为非定常场。 描述场的几何方法是引入所谓的场线,就像静电场中引入电力线,磁场中引入磁力线一样,在流速场中可以引入流线。流线是这样规定的: 流线为流体内的一条连续的有向曲线,流线上每一点的切线方向代表流体内微粒经过该点时的速度方向。 一
10、般情况下空间各点的流速随时间 t 变化,因此流线也是随时间变化的。由于流线分布与一定的瞬时相对应,所以在一般情况下,流线并不代表流体中微粒运动的轨迹。 只有在稳定流动中,流线不随时间变化,此时流线才表示流体中微粒实际经过的轨迹。只有此时流线才与迹线重合。 另外,由于流线的切线表示流体内微粒运动的方向,所以流线永远不会相交,因为如果流线在空间某处相交就表示流体中的微粒经过该点时同时具有两个不同的速度,这当然是不可能的。 在流体内部取一微小的封闭曲线,通过曲线上各点的流线所围成的细管就称为流管。 由于流线不会相交,因此流管内、外的流体都不具有穿过流管的速度,也就是说流管内部的流体不能流到流管外面,
11、流管外的流体也不能流入流管内。计算流体动力学(Computational Fluid Dynamics,简称CFD)是通过计算机数值计算和图像显示,对包含有流体流动和热传导等相关物理现象的系统所做的分析。CFD的基本思想可以归结为:把原来在时间域及空间域上连续的物理量的场,如速度场和压力场,用一系列有限个离散点上的变量值的集合来代替,通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组,然后求解代数方程组获得场变量的近似值CFD可以看做是在流动基本方程(质量守恒方程飞动量守恒方程、能量守恒方程)控制下对流动的数值模拟。通过这种数值模拟,我们可以得到极其复杂问题的流场内各个位置上
12、的基本物理量(如速度、压力、温度、浓度等)的分布,以及这些物理量随时间的变化情况,确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理量,如旋转式流体机械的转矩、水力损失和效率等。此外,与CAD联合,还可进行结构优化设计等。CFD方法与传统的理论分析方法、实验测量方法组成了研究流体流动问题的完整体系,图1给出了表征三者之间关系的“三维”流体力学示意图理论分析方法的优点在于所得结果具有普遍性,各种影响因素清晰可见,是指导实验研究和验证新的数值计算方法的理论基础。但是,它往往要求对计算对象进行抽象和简化,才有可能得出理论解。对于非线性情况,只有少数流动才能给出解析结果。 实验测量方法所
13、得到的实验结果真实可信,它是理论分析和数值方法的基础,其重要性不容低估。然而,实验往往受到模型尺寸、流场扰动、人身安全和测量精度的限制,有时可能很难通过试验力一法得到结果。此外,实验还会遇到经费投入、人力和物力的巨大耗费及周期长等许多困难。 而CFD方法恰好克服了前面两种方法的弱点,在计算机上实现一个特定的计算。就好像在计算机上做一次物理实验。例如,机翼的绕流,通过计算并将其结果在屏幕上显示,就可以看到流场的各种细节:如激波的运动、强度,涡的生成与传播,流动的分离、表面的压力分布、受力大小及其随时间的变化等。数值模拟可以形象地再现流动情景,与做实验没有什么区别。2 计算流体动力学的特点 CFD
14、的长处是适应性强、应用面广。首先,流动问题的控制方程,般是非线性的,自变量多,计算域的几何形状和边界条件复杂,很难求得解析解,而用CFD方法则有可能找出满足工程需要的数值解;其次,可利用计算机进行各种数值试验,例如,选择不同流动参数进行物理方程中各项有效性和敏感性试验,从而进行方案比较。再者,它不受物理模型和实验模型的限制,省钱省时,有较多的灵活性,能给出详细和完整的资料,很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件。CFD也存在一定的局限性。首先,数值解法是一种离散近似的计算方法,依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限数学模型
15、,且最终结果不能提供任何形式的解析表达式,只是有限个离散点上的数值解,并有一定的计算误差;第二,它不像物理模型实验一开始就能给出流动现象并定性地描述,往往需要由原体观测或物理模型试验提供某些流动参数,并需要对建立的数学模型进行验证;第三,程序的编制及资料的收集、繁理与正确利用,在很大程度上依赖于经验与技巧。此外,因数值处理方法等原因有可能导致计算结果的不真实,例如产生数值粘性和频散等伪物理效应。当然,某些缺点或局限性可通过某种方式克服或弥补,这在本书中会有相应介绍。此外,CFD囚涉及大量数值计算,因此,常需要较高的计算机软硬件配置。 CFD有自已的原理、方法和特点,数值计算与理论分析、实验观测
16、相互联系、相互促进,但不能完全替代,三者各有各的适用场合。在实际工作中,需要注意三者有机的结合,争取做到取长补短。3 计算流体动力学的应用领域 近十多年来,CFD有了很大的发展,替代了经典流体力学中的一些近似计算法和图解法:过去的一些典型教学实验,如Reynolds实验,现在完全可以借助CFD手段在计算机上实现。所有涉及流体流动、热交换、分子输运等现象的问题,凡乎都可以通过计算流体力学的方法进行分析和模拟。CFD不仅作为一个研究工具,而且还作为设计工具在水利工程、土木工程、环境工程、食品工程、海洋结构工程、工业制造等领域发挥作用。典型的应用场合及相关的工程问题包括: 水轮机、风机和泵等流体机械
17、内部的流体流动 飞机和航天飞机等飞行器的设计 汽车流线外型对性能的影响 洪水波及河口潮流计算 风载荷对高层建筑物稳定性及结构性能的影响 温室及室内的空气流动及环境分析 电子元器件的冷却 换热器性能分析及换热器片形状的选取 河流中污染物的扩散 汽车尾气对街道环境的污染 食品中细菌的运移 对这些问题的处理,过去主要借助于基本的理论分析和大量的物理模型实验,而现在大多采用CFD的方式加以分析和解决,CFD技术现己发展到完全可以分析三维粘性湍流及旋涡运动等复杂问题的程度。4 计算流体动力学的分支 经过四一十多年的发展,CFD出现了多种数值解法。这些方法之间的上要区别在于对控制方程的离散方式。根据离散的
18、原理不同,CFD大体上可分为三个分支: 有限差分法(Finite Difference Method,FDM) 有限元法(Finite Element Method,FEM) 有限体积法(Finite Volume Method,FVM) 有限差分法是应用最早、最经典的CFD方法,它将求解域划分为差分网格,用有限个网格节点代替连续的求解域,然后将偏微分方程的导数用差商代替,推导出含有离散点上有限个未知数的差分方程组。求出差分方程组的解,就是微分方程定解问题的数值近似解。它是一种直接将微分问题变为代数问题的近似数值解法。这种方法发展较早,比较成熟,较多地用于求解双曲型和抛物型问题。在此基础上发展
19、起来的方法有PIC(Particle-in-Cell)法、MAC(Marker-and-cell)法,以及由美籍华人学者陈景仁提出的有限分析法(Finite Analytic Method)等有限元法是20世纪80年BANNED始应用的一种数值解法,它吸收了有限差分法中离散处理的内核,又采用了变分计算中选择逼近函数对区域进行积分的合理方法。有限元法因求解速度较有限差分法和有限体积法慢,因此应用不是特别广泛。在有限元法的基础上,英国C.A.Brebbia等提出了边界元法和混合元法等方法。 有限体积法是将计算区域划分为一系列控制体积,将待解微分方程对每一个控制体积积分得出离散方程。有限体积法的关键
20、是在导出离散方程过程中,需要对界面上的被求函数本身及其导数的分布作出某种形式的假定,用有限体积法导出的离散方程可以保证具有守恒特性,而且离散方程系数物理意义明确,计算量相对较小。1980年,S.V.Patanker在其专著Numerical Heat Transfer and FluidFlow中对有限体积法作了全面的阐述。此后,该方法得到了广泛应用,是目前CFD应用最广的一种方法。当然,对这种方法的研究和扩展也在不断进行,如P.Chow提出了适用于任意多边形非结构网格的扩展有限体积法等。9 显示和输出计算结果 通过上述求解过程得出了各计算节点上的解后,需要通过适当的手段将整个计算域上的结果表
21、示出来。这时,我们可采用线值图、矢量图、等值线图、流线图、云图等方式对计算结果进行表示。 所谓线值图,是指在二维或二维空间上,将横坐标取为空间长度或时间历程,将纵坐标取为某一物理量,然后用光滑曲线或曲面在坐标系内绘制出某一物理量沿空间或时间的变化情况。矢量图是直接给出二维或三维空间里矢量(如速度)的方向及大小,一般用不同颜色和长度的箭头表示速度矢量。矢量图可以比较容易地让用户发现其中存在的漩涡区。等值线图是用不同颜色的线条表示相等物理量(如温度)的一条线。流线图是用不同颜色线条表示质点运动轨迹。云图是使用渲染的方式,将流场某个截面上的物理量(如压力或温度)用连续变化的颜色块表示其分布。 现在的
22、商用CFD软件均提供了上述各表示方式。用户也可以自己编写后处理程序进行结果显示。CFD的软件结构为方便用户使用CFD软件处理不同类型的工程问题,一般的CFD商用软件往往将复杂的CFD过程集成,通过一定的接口,让用户快速地输入问题的有关参数。所有的商用CFD软件均包括三个基本环节:前处理、求解和后处理。与之对应的程序模块常简称前处理器、求解器、后处理器。以下简要介绍这三个程序模块。1 前处理器 前处理器(preprocessor)用于完成前处理工作。前处理环节是向CFD软件输入所求问题的相关数据,该过程一般是借助与求解器相对应的对话框等图形界面来完成的。在前处理阶段需要用户进行以下工作: 定义所
23、求问题的几何计算域 将计算域划分成多个互不重叠的子区域,形成由单元组成的网格 对所要研究的物理和化学现象进行抽象,选择相应的控制方程 定义流体的属性参数 为计算域边界处的单元指定边界条件 对于瞬态问题,指定初始条件 流动问题的解是在单元内部的节点上定义的,解的精度由网格中单元的数量所决定。一般来讲,单元越多、尺寸越小,所得到的解的精度越高,但所需要的计算机内存资源及CPU时问也相应增加。为了提高计算精度,在物理量梯度较大的区域,以及我们感兴趣的区域,往往要加密计算网格。在前处理阶段生成计算网格时,关键是要把握好计算精度与计算成本之间的平衡。 目前在使用商用CFD软件进行CFD计算时,有超过50
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 流体力学 流体 仿真 介绍 15
限制150内