热物理过程的数值模拟-计算传热学1.doc
《热物理过程的数值模拟-计算传热学1.doc》由会员分享,可在线阅读,更多相关《热物理过程的数值模拟-计算传热学1.doc(44页珍藏版)》请在三一文库上搜索。
1、热物理过程的数值模拟Numerical Simulation of Thermophysics Process讲稿主讲:李隆键第一章 概 论11 流动与传热过程的予测方法及特点流动、传热、燃烧问题是热工类各专业和机械类动力机械专业所研究和解决的主要问题之一,燃烧问题实际上是有化学反应的流动与传热问题,推而广之,在所有热物理过程中,几乎都涉及到流动、传热问题。预测的重要性: 在规定设计参数的相应的结构下,热物理过程是否满足要求,达到预定的指标?要预测; 优化设计,不同方案的比较,要预测; 减少设计、生产、再设计和再生产的费用; 减少设计更改; 减少试验和测量次数。问题的核心:速度场、温度场(传热
2、量)、浓度场等。一、热物理问题的予测方法:理论分析法、实验测定、数值模拟1、理论分析以数学分析为基础,求解描述热物理过程的定解问题,获得函数形式的解,表示求解区域内物理量连续分布的场(速度场、温度场、浓度场)。控制方程+单值条件(数学模型)理论解(分析解,解析解)根据解的准确程度,又可再分为:(1)精确分析解(严格解)特点:函数形式的解;它在求解区域精确地满足定解问题。具体解法:直接积分法、分离变量法、积分变换法、热源法、映射法。(2)近似分析解法特点:函数形式的解,在求解区域上近似地满足定解问题(但在总量上满足相应的守恒原理,动量守恒、动量守恒、能量守恒、质量守恒)。具体解法:积分法(从积分
3、方程出发)变分近似解法摄动法(从微分方程出发)2、实验测定(1)纯实验法(2)相似理论实验法:同类相似,减少变量数目减少工作量,得到规律性结果,可直接应用。(3)实验类比法:异类相似物理现象不同,规律相同:微分方程形式相同,单值性条件类似电热类比,水热类比3、数值模拟以数值计算方法为基础,借助(利用)电子计算机求解物理过程的方法热物理过程的数值模拟,对传热过程称为传热的数值模拟、数值传热、计算传热。如前述,传热过程函盖了流动、燃烧,所以计算传热学实质上就代表了热物理过理过程的数值模拟。用电子计算机对热物理问题进行数值计算就象在实验室中对该现象进行实验测定一样,可称之为“数值实验”。随着高速、大
4、容量电子计算机的发展,特别是微型计算机的普及和推广,这种数值实验的方法越来越被更多的科技人员掌握和应用,成为解决热物理过程的一种重要方法。二、予测方法的比较与选择1、分析解法(1)精确预测了数学模型所控制的的热物理过程;(2)函数形式的解使得可以确定区域中任意位置物理量的大小;(3)以显函数的形式,展示各有关参量对该热物理过程的影响;(4)由于是函数形式的解,便于进一步的运算、处理,例如求导、积分。缺点:(1)获得分析解的可能较小;(2)即使能求得分析解,也常常是无究级数,特殊函数以及涉及特征值问题的超越函数,要得到具体的数值结果,也需要繁复的计算;(3)数学模型的结果也需要有实验检验。2、实
5、验方法(1)可以获得热物理过程可靠的数据资料;(2)全比例设备实验可予测由它完全复制的同类设备在相同条件下将如何运行和变化;(3)是研究一种新的基本现象的唯一方法;(4)是检验其它预测方法准确程度的标准。缺点:(1)全比例实验代价大(投资,物力,人力,周期);(2)缩小比例模型实验结果的外推受准则数实验范围的限制,有些在全比例设备上才能出现的特征在缩小比例模型上并非总是能模拟(例如流动的涡),降低了模型试验的效果;(3)测试困难及测量误差;(4)有些过程无法预先进行试验(航天,气象预报)。3、数值模拟(1)成本低:在大多数实际应用中,计算机运算的成本要比相应的实验研究成本低好几个数量级,对象愈
6、庞大,过程愈复杂,此优点愈突出;同时,与大多数物品价格不断上涨的趋势相反,计算成本还会降低;(2)速度快,周期短;不同方案的对比计算和优选,这对某些大型实验几乎是不可能的。(3)信息完整:能提供计算区域内所有各个位置上有关变量的值(速度、压力、温度、浓度等),而实验则不可能测出整个区域各点处所有变量的值。(4)具有模拟真实条件的能力:几何条件、边界条件、物性条件、初始条件很容易模拟真实条件,不需要采用缩小模型或冷态实验,无论大小、高位,低温、过程快、慢。(5)具有模拟理想条件的能力:对于研究物理现象而不是工程问题时,注意力集中几个基本参数而要设法消除所有无关的因素。几何条件(维数变化,尺寸)、
7、物性(常密度),BC(绝热表面),ic(特定的初始温度分布)。缺点:(1)数值模拟的对象是数学模型简化处理,结果的准确性有特价检验;(2)对一些十分复杂的问题(几何形状复杂,强烈非线性、物性变化大),数值解可能很难获得,或者即便可以获得,代价也是相当昂贵的,例如,对湍流问题,要想通过求解非稳态N-S方程来算出它们的全部与时间相关的结构,则仍然是计算所不能及的;(3)对解的唯一判断力较弱为了进一步讨论数值模拟的缺点,可以把所有的实际问题分成两大类:A类:有完整数学模型的一类问题,如热传导、层流问题、简单的湍流边界层问题;B类:迄今无完整数学模型的一类问题,如复杂湍流、某些非牛顿流体、某些两相流动
8、等,问题的分类还有一标准问题,即描述到什么样的程度可以认为是“够了”,“合适”的。A类缺点:对这类问题,用计算机求解的优越性远远大于实验研究。 有数值模拟缺点的(2)、(3)两条;在某些情况下也需要进行实验检验。对于此类问题,研究计算方法的目的在于使这些计算方法理加可靠、准确和有效,随着研究的进展,其缺点将被不断克服。B类缺点:A类的缺点B类全有,此外,必须进行实验检验。数学模型的研究不断地把B类问题转化为A类:试算与修正。先提出一个模型计算求解与实验结果进行比较修正模型,并不断完善、湍流模型的最新发展就是这种转换的一个典型例子,k-双方程模型最初建立在科尔莫戈洛夫(Kolmgorov)(19
9、42)及普朗特(Prandtl)(1945)的工作基础上的,但并未,也不可能付诸实现,只有到了20世纪70年代,当计算机和计算方法变得更加强有力的是候,该模型才逐步趋于完整并付诸实际应用。4、方法的选择三种方法或三种手段相辅相成,互为补充。(1)分析解可以为检验数值模拟结果的准确度提供比较依据;常常用有分析解的简单问题检验方法的准确度;(2)简单的解极解可以为发展数值方法中的某些算法提供理论依据,调和平均;(3)物理规律、数学模型的正确建立必须通过对现象的充分观察和测定,等;(4)出现在数学模型中的物性参数只有通过实验测定才能获得。数值模拟的对象是热物理过程的数学模型,所以其结果的准确度首先取
10、决于数学模型反映实际热物理过程的准确度(包括所用的特性参数),然后才是所采用的数值方法,计算机并不能创造信息,发现规律,它只能把人们所送入的信息,按照计算者安排的程式对信息进行加工、处理,从而得到相应的结果。但是,一旦确立了与实际物理过程相符合的物理模型、数学模型、数值模拟又可以发挥很大的作用,它可以减少实验工作量,拓宽实验研究的范围,实现对理想单值性条件的模拟;对那些耗资巨大,条件恶劣的实验,或者难于进行的实验来说,“数值实验”更是一种有吸引力的辅助或替代手段。理论分析,实验测定和数值模拟有机而协调地结合,是研究热物理过程理想而有效的方法。12 本课程的内容及安排一、内容两个组成部分1. 理
11、论部分: 基础理论(数学物理、数值方法)热物理过程的数值模拟:通用性,并以热传导、对流换热、流场的计算为例,推广到通用控制方程所描述的其它现象。论述方式的特点:(1)强调物理的概念和方法,而不过分倚重纯数学的推导。(2)以一维为基础,推广(扩展)到二维、三维。2、实践性环节:一个课程设计,程序设计,视情况而定。二、授课方式改变注入式,实行启发式,培养自学能力。三、教材、参考书1、SV帕坦卡著、张政译(郭宽良译),“传热与流体流动的数值计算”。(Numerical Heat Transfer and Fouid Flow),科学出版社,1984。2,陶文铨编著,数值传热学,西安交通学大出版社,1
12、988。第二章 物理现象的数学描述数值计算的对象过程的数学模型,核心是控制方程。(数学模型:控制方程+单值性条件,?单值性条件)21控制微分方程1、控制微分方程的意义控制微分方程是一定守恒原理的数学表达式,影响因变量的各因素之间必定存在某种联系。回忆导热方程:热力学第一定律,对任意控制容积V,导入控制容积的热流量+控制容积中的内发热量=控制容积中物质能量的增量。导入控制容积的净热流量,控制容积内的发热量=,控制容积中物质内能的增量=。 散度定理: 非稳态项 扩散项 源项比容,对常密度(不可压缩)物质运动(流动)时,增加:物体宏观运动带入的能量= 压缩机械动 = 粘性摩擦功 =流体的变形率张量第
13、二动力粘度,体积粘度对简单可压缩系统= 非稳 对流 扩散对理想气体及恒密度(固液), ,对固、液,忽略p变化,则各项则代表着各因素在单位容积时的作用效果. 单位容积焓的变化率单位时间、单位面积、传递的焓。对流流量密度JC 单位容积流出的净焓量扩散流量密度 源项2、化学组分方程的质量分量扩散流密:3、能量方程如果为常数,OK时h=0,则由4、动量方程:某方向上的动量变化率,单位体积质量x方向的动量,其它粘性力项,5、湍流的时间平均方程湍流:给定点处的物理量随时间而变化,随机性。工程上:关心的是运动状态的时间平均特性非稳定层流方程的时间平均方程,滞流流动的时间平均方程,假设湍流中存在相对于时均值的
14、脉动,平均化(时均化)运算附加项,雷诺应力(湍流热流),湍流扩散流量密度湍流模型,用平均性质来表示的附加项。引入湍流粘度(湍流扩散系数)滞流应力,流量必度时均方程层流方程,相应的层流交换系数用有效系数代替。模型6、通用方程(1)向量形式对于不同的,有特定的、S与之间相对应。大多数的传递过程,扩散流量密度Jd因变量的梯度确定,也有扩散流量密度不由梯度支配的情形,这时可将其(Jd)并入S中。连续性方程:(2)直角张量形式和法则:一项某下标重复,则该下标依次取1、2、3然后做和运算。如何把特定的控制微分方程改写成通用形式。把相关因变量的非稳态项,对流项及扩散项转化成标准形式;把扩散项内梯度前的系数取
15、为的表达式;把其余所有项之和置于方程右端,定义为源项。通用方程可以是有量纲形式,也可以是无量纲形式,、S也相应无量纲化。作用:通用方程通用数值方法公式通用计算机程度2-2 坐标的性质及控制方程的类型一、坐标的性质坐标自变量。坐标(系)或自变量的数目对问题的难易程度有很大影响,而对一定的问题而言,坐标(自变量)数目是可变的,既可以用这个坐标系来描述,又可以用别的坐标系来描述。1、自变量的作用:在数值计算中,将选择用来计算值的自变量值,所需计算值的位置的多少与自变量的数目有关,自变量数目所需计算的位置。因变量与自变量的相对性,适用于温度场是坐标的单调函数的情形。2、坐标的选择,原则:恰当,合适。自
16、变量数目最少的坐标系网格节点数少。(1)选择最简单的坐标系圆柱(简)中的轴对称导热:圆管内的轴对称流动:(2)运动坐标系准稳态的概念,(3)利用充分发展的概念:存在这样一个坐标,当过程发展到一定深度后,因变量的无量纲分布与该坐标无关。,充分发展后平面自由射流;但(4)相似变换:减少自变量数目的变换统称相似变换。半无限大物体(x0)在1B.C.下的非稳态导热:3、坐标的单、双向性采用单向坐标和双向坐标的概念,可以形象地描绘不同类型控制方程的物理作用上的区别。单向坐标:在一个坐标轴上,如果扰动(或影响)只能向一个方向传递,则称此坐标为单向坐标。何谓向一个方向传递?坐标上任意给定位置处因变量之值只受
17、该位置一侧条件变化的影响,且该点因变量之值也只对其一侧位置上的固变量值发生影响。(“抛物型”表示一种“单向作用”的概念)时间坐标是一个典型的单向坐标。高温固体的冷却,其在某瞬时的温度只受该瞬时以前的条件的影响。双向坐标:在一个坐标轴上扰动(或影响)可以向两侧传递,称为双向坐标。何谓向两个方向传递?坐标上任意给定位置处因变量之值要受该位置两侧条件变化的影响,且该处因变量之值也会对其两侧位置上的因变量值发生影响。空间坐标是典型的双向坐标。但在一定条件下,空间坐标也可以成为单向坐标;在流体流动时,如果在某一个坐标上有很强的单向流动,则重要的影响只能从上游传播到下游,某处状态也主要受其上游条件的影响,
18、受下游条件影响很小。总结:对流是一种单向过程,而扩散是一种双向过程:在流体流动时,二者同时存在,仅当流作用很强(流量很大)时,扩散作用可忽略不计,空间坐标才近似成为单向坐标。二、控制方程的类型讨论控制方程的类型对数值解的影响1、能量方程导热:对流换热:(已假定:)若常数,OK时h=0,则组成:非稳态项,对流项,扩散项,源项控制方程:微分容积中守恒原理的数学描述。2、能量方程的守恒性质数值计算中的控制容积无论多么小,总是有限容积,而不是微分容积。问题:对任意有限容积,控制方程是否也一定满足守恒原理。守恒型控制方程,对任意大小的有限容积能使守恒原理得到满足的控制方程。非 不能可以证明,在直角坐标中
19、,当对流项为散度形式时,控制方程是守恒的。将上面的控制方程对任一有限容积V作积分有限容积V上的能量平衡原理的表达式:单位时间内有限容积中物质内能的增加量对流带入的净能量率导入的净能量率生成的能量率若用非守恒型控制方程,有: 简证:取平行六面体 3、控制方程的分类(1)从数学上分类,限于二阶偏微分方程,对二元二阶线性偏微分方程,下标xy表示对该变量的导数,a、b、c、d、e、f、是x、y的函数,对求解区域R,方程的特性由系数之值决定,在区域内任意一点有两条实的特征线 双曲型,在区域内任意一点有一条实的特征线 抛物型,在区域内任意一点无实的特征线 椭圆型三类方程在数学上的主要区别:影响区域和依赖区
20、区域不相同。影响区域、依赖区域;R中任一点P的依赖区域,为了唯一定之道,必须给函数值的点B依赖区域的集合;P的影响区域:变化时,函数值发生变化的点的集合。椭圆型方程(即求解区域上每一点都是椭圆型的):对于能量方程无非稳态项,自变量与时间无关。特点:无特征线,任一点P的依赖区域是包围该点的区域封闭边界曲线,而P点的影响区域则是整个求解区域。边值问题,对应于物理学上的一类平衡问题,或稳态问题。求解区域内各点处的因变量值是相互影响的,与双向坐标相联系。结果:离散代数方程必须联立求解(直接解法或迭代解法),而不可能把区域中某一部分的值求得后再去确定其区域上的值。抛物型方程:因变量与时间有关,或问题中存
21、在类拟于时间的自变量,与单向坐标相联系,能量方程有非稳态项。对应于:物理学上一类非平衡问题,或非稳态问题,步进问题。特点:过区域中任一点P有一条实特征线,其方向与单向坐标相垂直,如图,P点的依赖区域和影响区域以特征线为分界线。对非稳态问题:某一瞬间物体中的温度(分布)取决于该瞬时以前的情况及边界条件,而与该瞬时以后将要发生的情况无关,反之,某一时刻的温度只影响此后的温度分布边界层类型的流动与换热:忽略主流方向的扩散作用,下游的物理量(u,v,t)取决于上游,上游 只会影响下游的物理量。结果:不必将单向坐标上所有位置处的离散方程联立求解!只需从某一初始值出发,结合边界条件,沿单向坐标一步步向前推
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 物理 过程 数值 模拟 计算 传热学
链接地址:https://www.31doc.com/p-2744372.html