




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
热流问题的数值计算
NumericalSimulationsofThermal&FluidProblems第2章计算区域与控制方程离散主讲李炎锋2008年7月北京2.1空间区域的离散化(Domaindiscretization)2.1.1区域离散化的实质和内容
实质:用一组有限个离散的点代表原来的连续空间.实施过程:把所计算的区域划分成许多个互不重叠的子区域(sub-domain),确定每个子区域的节点位置及该节点所代表的控制容积(controlvolume)。4种几何要素:
节点(grid):需要求解的未知物理量的几何位置.控制容积(controlvolume):应用控制方程或守恒定律的最小几何单位.界面(surface):与各节点相对应的控制容积的分界面位置.网格线(gridline):沿坐标轴方向联结相邻两节点而形成的曲线簇。2.1.2两类设置节点的方法
1.外节点法(practiceA):先节点后界面法。2.内节点法(practiceB):先界面后节点法(见表2-1)用P,N,E,W,S表示所研究的节点及相邻的节点,用n,e,w,s表示相应的界面。两节点间的距离用x表示,x表示相邻两界面间的距离。2、非结构化网格(unstructuredgrid):生成复杂,对不规则区域适应性强;1、结构化网格(structuredgrid):生成简单,对不规则区域适应性差;两类网格比较:2.1.3两类节点设置方法的比较1、边界节点代表的控制容积不同。外节点法中,位于非角顶上的边界节点代表了半个控制容积;内节点法中,则应看成是厚度为零的控制容积的代表。2、当网格不均匀时,内节点法中节点永远位于控制容积的中心。外节点法则不然。
3、对一维问题,当节点不均匀布置时外节点法的离散误差比较小,但对二维流动问题用四边形网格来计算时,发现两种布置方法得出的结果几乎相同。2.1.4关于网格生成的进一步说明2、同一坐标方向上相邻两个子区域(或控制容积)的宽度应保持在一个合适的范围内。
1、对每个控制容积在不同方向的宽度应保持一个合适的比例。对椭圆型问题,应为接近1;对抛物型问题或某个方向的变化率明显大于另一个方向的椭圆型问题,变化剧烈的方向应取较小的宽度。对不均匀网格应注意:守恒型:非守恒型:2.2.1一维模型方程2.2建立离散方程的Taylor展开法及多项式拟合法2.2.2用Taylor展开法导出导数的差分方程式得:由:的截断误差小于等于来代替其中,(x)称为截断误差(truncationerror),表示随着x趋近于零,用(x2):中心差分(centraldifference)(x):向后差分(backwarddifference)(x):向前差分(forwarddifference)用一阶精度的差分表达式来代替,为:2.2.3一维模型方程的离散C-N格式:每一时层上空间导数的差分按每一层的中间时刻之值计算.隐式:每一时层上空间导数的差分按每一层的终了时刻之值计算;显式:每一时层上空间导数的差分按每一层的初始时刻之值计算;非稳态问题的三种格式:将一维模型方程的精确解在接点(i,n+1),(i+1,n)(i-1,n)上对节点(i,n)作Taylor展开,取时间的向前差分,空间的中心差分,可得:用近似值代替精确解(i,n),可得Taylor展开法导出的一维模型方程的显式离散格式为:2.2.4用多项式拟合法建立导数的差分表达式
设x0=0,x0为点(i,n)的x坐标,则有:设函数(x,t)在节点(i,n)附近对x的变化近似地成线性,则有:1.线性逼近原理:对未知函数的局部变化型线采用多项式逼近在点(i,n)的向前差分为:由上两式得,由此得,有:2.二次曲线拟合:一、二阶导数的差分逼近为:例题:设在如图所示的情形中,已知区域内部与边界节点的温度,物体导热系数为常数,试用多项式拟合法确定穿过壁面的热流密度。如果取温度分布为二次曲线,yTi,1Ti,5Ti,4Ti,3Ti,2xy解:设壁面附近温度T按线性变化,则
则有:该具有二阶精度.于是得:由此得:2.3建立离散方程的控制容积积分法及平衡法1.将守恒型的控制方程在任一控制容积及时间间隔内对空间和时间作积分;2.选定未知函数及其导数对时间及空间的局部分布曲线,即型线;3.对各个项按选定的型线作积分,并整理成关于节点上未知值的代数方程.2.3.1控制容积积分法的实施步骤及常用的型线主要步骤:分段线性分布(piecewiselinear)阶梯式分布(stepwise)常用的型线:(2)(1)将上式对控制容积P在t时间间隔内作积分,得:一维守恒型模型方程:2.3.2用控制容积积分法离散一维模型方程对对流项:选定随t变化的型线为阶梯式,得:对非稳态项:选定随x变化的型线为阶梯式,得:(3)(4)(8)(7)(6)若取随x呈分段线性变化,则式(4)(5)中界面的对流项(u)及扩散项可表示为:(5)对扩散项:选取一阶导数随时间作显式阶跃式变化,得:将式(3)-(9)代入(2)得:(9)对源项:假设S对t及x均呈阶梯式变化,则:把物理上的守恒定律直接用于所研究的控制容积,并把节点看成是控制容积的代表,导出节点上未知值间的代数关系式。2.3.3由控制容积平衡法导出离散方程注意:用控制容积积分法与用Taylor展开法导出的离散方程不一定相同,计算结果的精确度也不一样。
对一维问题,则有在t时间间隔内控制容积P中变量的增量,等于在同一时间间隔内由流及扩散作用进入该控制容积的的净值及源项所生之值的总和,有:2.3.4不同离散方法的比较Taylor展开法与多项式拟合法侧重于从数学角度进行推导,把控制方程中的各阶导数用相应的差分表示式代替;优点:易于对离散方程进行数学特性分析,缺点
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- Unit4 Section A (2a-2d) 教学设计2024-2025学年人教版英语八年级上册
- 圆的面积(一)(教案)2024-2025学年数学六年级上册-北师大版
- 2025年九江下载货运从业资格证模拟考试
- 2025年哈尔滨货运丛业资格证试题和答案
- 2025年泉州货运资格证模拟考试
- 七年级英语下册 Unit 8 The seasons and the Weather Topic 3 Lets celebrate Section C教学实录 (新版)仁爱版
- 2025年改性塑料粒子合作协议书
- 2025年异丙醚项目建议书
- 2025年节能技术服务合作协议书
- 教育热点话题分析
- 2024年陕西省中考数学试题含答案
- 2024年国网河南省电力公司招聘历年高频500题难、易错点模拟试题附带答案详解
- 2024年清理化粪池协议书
- 民营医院并购合同范本
- 中国高血压防治指南(2024年修订版)图文解读
- 2024年全国国家电网招聘之电网计算机考试历年考试题(附答案)
- (正式版)QB∕T 8049-2024 家用和类似用途微压富氧舱
- 传染病医院隔离手术操作流程
- 2024年广东省高三一模高考英语试卷试题答案祥解(含作文范文)
- 全新合伙协议三人模板下载
- 人体成分分析在健康管理中的应用
评论
0/150
提交评论