- 1、本文档共20页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
传热学热传导问题的数值解珐培训教案
传热学 第 4 章 热传导问题的数值解法 * * 数值解法: 有限差分法,有限元法及边界元法。 分析解法: (1)对导热微分方程式在规定的边界条件和 初始条件下积分求解。 (2)求解结果能清楚显示各种因素对温度分布 的影响,但仅适用于简单的导热问题,同时解 的形式复杂。 导热问题的求解方法: 实验方法: 相似原理指导下的实验方法。 4.1 导热问题数值求解的基本思想 理论基础:离散数学。 基本思想:把原来在时间,空间坐标系中连续的 物理量的场(如导热物体的温度场)用有限个 离散点上的数值的集合来代替,通过求解按一定 方法建立起来的关于这些值的代数方程,来获得 离散点上被求物理量的值,这些离散点上被求物 理量值的集合称为该物理量的数值解。 基本方法: 有限差分法 (finite-difference) 有限元法 (finite-element) 边界元法 (boundary-element) 4.1.2 导热问题数值求解基本思路 分析;二维矩形区域内的稳态,无内热源, 常物性的导热问题。 数值求解的步骤; 1. 建立控制方程及定解条件 控制方程:描述该导热问题的导热微分方程。 定解条件:边界条件。 二维,稳态,无内热源,常物性的导热: 举例 导热微分方程: t0 W H h1 tf h2 tf h3 tf x y 2. 区域离散化 (discretization) 沿x方向和y方向分别以Δx,Δy为间隔把 求解区域划分成很多个小的子区域。 步长:相邻两节点间的距离Δx, Δy。 节点:网格线(边界线)的交点。 节点表示: n N M Δx h3 tf m Δy m,n m+1,n m-1,n m,n+1 m,n-1 控制容积 (control volum): 以节点为中心, 边长等于Δx, Δy的小区域。 节点温度代表以它为中心的控制容积的平均温度。 3. 建立节点物理量的代数方程 离散方程:节点上物理量的代数方程。 二维稳态导热均匀步长时内部节点温度 的代数方程: 举例 4. 设立迭代初场 迭代初场:采用迭代法求解时,需要对被求解 的温度场预先假定一个解,称为初场。 5. 求解代数方程组 6. 解的分析 可进一步计算热流量,热应力及热变形等。 4.2 内节点离散方程的建立方法 数值计算过程的核心内容 . 两种方法: 泰勒级数展开法; 控制容积热平衡法 . 4.2.1 泰勒级数展开法 根据泰勒级数,导出节点(m,n)处二阶偏导数 的差分表达式: 上两式相加得: 整理得: 其中 称为截断误差。 整理成二阶导数的近似代数关系式: ——二阶导数的中心差分表达式 (a) 同理可得: (b) 将(a), (b)代入导热微分方程式中得离散方程: 设步长Δx=Δy,有: —— 二维稳态导热均匀步长时内部节点 温度差分方程 说明 物体内每一个节点的温度都等于它周围相邻 4 个节点温度的算术平均值。 4.2.2 热平衡法 对控制容积写出能量守恒的表达式: n Δx w Δy m,n m+1,n m-1,n m,n+1 m,n-1 s e Φw Φe Φn Φs 根据傅立叶定律: 设步长Δx=Δy,有: 垂直纸面单位宽度 说明 稳态时流向任何节点的热量总和必等于零。 4.3 边界节点离散方程的建立及代数方程的求解 4.3.1 边界节点离散方程的建立 第一类边界条件 — 边界温度已知 第二类边界条件 第三类边界条件 需建立边界节点温度 的差分方程 1. 位于平直边界上的节点 若步长Δx=Δy,有: n Δx m,n m m,n+1 m,n-1 q w Δy m-1,n 2. 外部角点 控制容积的热平衡为: 若步长Δx=Δy,有: Δx Δy m+1,n m-1,n m,n+1 m,n-1 A C B E D F m,n m,n-1 q w m-1,n q w 3. 内部角点 控制容积的热平衡为: 若步长Δx=Δy,有: Δx Δy m+1,n m-1,n m,n+1 m,n-1 A C B E D F m,n m,n-1 q w m-1,n q w 4. 边界热流密度的三种情况 (1)绝热边界: (2) 值不为零:代入给定的 值。 (3)对流边界: 平直边界节点: 外部角点 内部角点 式中:
您可能关注的文档
最近下载
- 管理学-PPT课件(全).pptx
- 高峰数字(四)_原创精品文档.pdf VIP
- “深圳工匠杯”深圳市第十四届职工技术创新运动会暨2024年深圳技能大赛碳排放管理员技能竞赛理论题库及答案(300题).docx VIP
- 粉针剂包装技术与设备.pptx
- 全文解读新修订的2021《中国共产党统一战线工作条例》ppt精品.pptx VIP
- 白蚁防治服务合同协议书范本.docx VIP
- (2022年新版)义务教育初中英语课程标准详细解读 (新课标解读讲座培训课件).pptx VIP
- 湖南《住宅工程质量常见问题防治技术规程》.docx
- 葡萄酿酒有限公司生产安全事故应急预案.docx VIP
- 路基爆破设计及施工组织方案.doc
文档评论(0)