流体分析 CFD 基础简介.ppt

CFD基础简介,概要,什么是CFD?CFD的历史CFD基本概念CFX-5.7功能简介CFX-5主要应用CFX-5求解问题的主要过程CFX-5求解器的能力CFX-5用户界面,什么是CFD?,CFD是计算流体动力学的缩写，即ComputationalFluidDynamics。CFD是模拟包含流体流动、传热以及其它附加的物理化学现象的工具。,Coalcombustioninamulti-burnerfurnace.CourtesyofInternationalCombustionLimited,CFD的历史,1970s中期开始理解CFD的本构方程，发展通用计算程序.1980s早期-第一代通用的CFD程序诞生了,需要在大型计算机上运行，使用者需要具有专业的计算流体动力学知识。仅在研究领域使用。目前-计算机性能的极大提高，计算图形学的广泛使用以及稳健的CFD求解器的开发使得CFD成为大多数工程师喜爱的分析工具。,CFD基本概念,质量守衡定律控制体牛顿运动定律能量守衡线性和非线形的偏微分方程离散原理非线性方程边界条件,守衡,所有设备中的流体流动都可以用质量,动量和能量守衡方程描述质量守衡指物质既不会产生也不会灭亡,只会保持一致动量守衡方程是牛顿运动定律的表述能量守衡方程是热力学第一定律的表述这三个守衡方程是CFD的理论基石！,控制体指空间的一个区域该区域被控制面所包围其大小和形状是任意的,可以是管,设备,旋转叶片等适用于有限的和极小尺寸的区域基本的守衡定律适用于这种定义的区域,控制体,质量流入,质量流出,质量总量,控制体质量增加率,=,流入的质量-流出的质量,表述,积分方程,控制体,控制面,例:,均布,1D,稳态流动,质量守衡,动量流入,动量流出,动量总量,动量增加率,=,+流出的动量-流入的动量,表述,积分方程,净力,净力,例:均布,1D,稳态流动,动量守恒,流入的总焓,流出的总焓,能量总量,内部能量变化率,流入热量-输出功=,表述,积分方程,流入热量,输出功,+流出的焓-流入的焓,例子:,1D稳态流动,能量守衡,偏微分方程,守衡的表达式,以不同阶次偏微分方程形式描述以微分形式的方程而不是以积分方程开始,是因为微分方程用起来更方便理论解通常只有对一些简单情况,具有简单的边界条件时才能获得例如:具有温度边界条件的矩形空间的线性导热方程,方程特性,流动方程通常是复杂和非线性的求解的是一组耦合(或相关)的方程而不仅仅是一个方程方程具有强非线性无法获得简单的理论结果必须借助于近似的离散方法,如有限元,有限体或有限插分原理CFX5采用基于有限元的有限体积法,Example:2DHeatConduction,T=50,T=100,T=0,q=0,线性和非线性方程,线性方程容易求解解有保证例如:简单的导热(导热系数为常数),势流方程,和弹性应力分析,离散原理,步骤:确定分析区域将控制区划分为离散区域(网格)对每一个网格应用守衡准则:对每个网格上的微分方程积分应用Gauss原理将体积分转换为面积分插值所有网格和网格表面的热流相如需要,线性化应用边界条件求解一组代数方程的结果如果是非线性的,迭代直至收敛,非线性方程,非线性意味着:系数取决于结果,因此:需要迭代,来更新非线性无法用代数法直接求解离散方程求解非线方程的步骤线性化方程已假定的初场开始求解代数方程更新非线性系数迭代至收敛,边界条件,完全确定问题必须有边界条件边界条件必须可靠且连续例如包括:在壁面,流体速度相对于壁面为零进口总压或速度出口静压或质量流量,CFX-5.7功能简介,稳态/瞬态流动层流/湍流流动单相/多相流动(欧拉多相流，拉格朗日多相流)不可压流动/可压缩流动传热/流固共扼传热浮力驱动流非牛顿流多组分流动燃烧和辐射标量传输方程旋转坐标系多重并行计算通用网格界面（GGI)二次开发工具,CFX5应用举例,CFX-5在以下领域被广泛应用化学过程工业旋转机械工业汽车工业航空航天工业国防工业能源工业环境工业.,混合器,混合器中的等值面和流线.,静态混合器,静态混合器中的流线.采用CFXExpressionLanguage(CEL)定义非牛顿流.,旋风分离器,采用雷诺应力湍流模型的结果.从流线的曲率上可以看出雷诺应力模型的计算结果更精确,采用k-e湍流模型的结果。,催化净化器,压力和速度矢量分布.采用了多孔介质模型.多孔介质引起了很大的压差，并使进口突扩段的回流区光滑。,空冷换热器,空冷换热器中的速度矢量和空气温度分布,径流泵,径流泵叶片表面的压力分布。CFX-5的多重参考坐标系功能在许多种的旋转机械上都有应用。,F22飞机,飞机表面的压力分布.采用了CFX-5强大的并行功能。来流马赫数为0.9，攻角为5.。,GeometryandgridprovidedbyDaimlerChryslerAerospaceAG,楼道内的烟雾流动,详细的几何特征（例如楼梯扶手)。求解结果可以输给其它后处理软件如Ensight。,CFX5.7分析问题的主要过程,StepI:几何模型和网格划分定义几何定义关心的区域（2D或3D)创建网格StepII:设置物理问题和条件-CFX-Pre定义域定义边界条件和初场定义求解参数StepIII:求解-CFX-Solve求解控制方程组StepIV:后处理-CFX-Post分析结果,定义流体域.创建流体域(流体,固体,导热固体,多孔介质)设置流体属性(粘性,导热系数,比热,等)选择物理模型(湍流模型,可压缩流,浮力流,多相流,等)定义边界条件几何体的所有外部表面都要定义。定义求解参数和输出选项,II.设置物理问题和条件,CFX-Pre,III.求解过程,CFX-5求解器求解控制方程组监视收敛过程,IV:后处理,在后处理模块CFX-Post中分析计算的结果数值的结果图形的结果批处理交互式,CFX-PreGUI,ToolsMenuBar,DefineMenuBar,ObjectSelector,PhysicsErrorSummary,Viewer,CFX-5.7文件结构,CFX-Post-.res-.def-.gtm,CFX-Prename.cfxCfx.#.pre,ThirdPartyMeshes-CFX-4-TASCflow-ICEM-.,CFX5.7文件类型,.cfxCFX-Pre产生的case文件.preCFX-Presession文件.defCFX-Pre产生的包括物理问题及网格的定义文件，对求解器来讲是输入文件。可以用CFX-Post打开来检查网格。.res求解器产生的结果文件，包括网格、计算结果和.def文件中的其它信息.out求解器产生的文本文件，记录了求解过程。,压力p、温度T、比容v一、压力(pressure)p物理中压强，单位:Pa,N/m2,流体分析基本状态参数,常用单位：1bar=105Pa1MPa=106Pa1atm=760mmHg=1.013105Pa1mmHg=133.3Pa,(1)ppbp=pb+pe;(2)ppbp=pb-pvpe:表压力pv:真空度,表压、真空度、绝对压力的关系,p,p,pe,pb,p,pb,pv,p,(1),(2),关于压力的几点说明,液柱高度表示的压力单位与压力的关系：p=gh；真空状态下计算绝对压力时，大气压的取值对计算结果影响较大，不能简单取pb=0.1MPa;表压力不是状态参数。,定义：vV/mm3/kg物理意义：工质聚集的疏密程度物理上常用密度m/Vkg/m3v=1,二、比体积(specificvolume)v,定义传统：冷热程度的度量。感觉，导热，热容量微观：分子平均动能的量度T0.5mw21)气体分子运动论结论2)0.5mw20,T0,1951年核磁共振法对氟化锂晶体的实验发现负的开尔文温度。3)T=00.5mw2=0分子一切运动停止，零点能。,三、温度(temperature)T,T/K=273.15+t/t/=5(t/-32)/9T/0R=t/+459.67T/0R=9/5T/K,4.常用温标,摄氏度K开尔文度华氏度0R兰金温标度,出口质量m2、流速cf2、标高z2,开口系统的稳定流动能量方程,在时间内,进口质量m1、流速cf1、标高z1,稳定流动：,在时间内进入系统的能量,在时间内离开系统的能量,根据热力学第一定律可得,令,，h称为比焓。,比焓的物理意义：,比焓是状态参数；对于流动工质，比焓表示每千克工质沿流动方向向前传递的总能量中取决于热力状态的部分。,上式可整理成,令,，上式改成,以上两式称为开口系统的稳定流动能量方程。,对于单位质量工质，,对于微元过程，稳定流动能量方程写成,注意：,（1）无论对于流动工质还是不流动工质，比焓都是状态参数；,（2）对于流动工质，流动功等于pv，比焓表示单位质量工质沿流动方向向前传递的总能量中取决于热力状态的部分；,（3）对于不流动工质，不存在流动功，比焓也不表示能量，仅是状态参数。,（4）工程上一般只需要计算工质经历某一过程后焓的变化量，而不是其绝对值，所以焓值的零点可人为地规定。,绝热节流,ws=0c2=0z=0q=0h=0h1=h2,节流前后h不变，不是等h过程。,喷管与扩压管,喷管：获得动能扩压管：获得压力能,动能参与转换，不能忽略,z=0ws=0q=0,h=-c2/2,动能与焓变相互转换。,伯努利方程,对于不可压缩均质流体，即=常数；如z1=z2，则：,即：,湿饱和蒸气vvv”hhh”sss”,未饱和水(过冷水)hh”ss”,水蒸气的定压发生过程,水蒸气定压加热过程在常用的p-v图，T-