TM01圆柱体稳态热传导对比分析
写在之前
ANSYS中有Steady-State Thermal(SST)分析系统与Fluent分析系统,都可以做传热分析,那到底有什么不同呢?下面详细展开……
- 在SST中,传热分析基于Fourier’s law的传导方程,将热对流Convection处理成简单的边界条件
- SST计算其实调用Mechanical的内核,只能用于固体传热
- 稳态传热分析往往用作瞬态热分析的第一步或最后一步,仿真只需要材料的导热系数lamda
- Fluent同时具备流体与固体传热计算的能力,可以共轭传热/流动问题、流固耦合传热仿真,能力更强
对流换热系数
/表面传热系数
/换热系数
/膜系数
有可能与温度相关,SST和Fluent中均可通过函数或用户数据定义- Fluent通过PROFILEs或UDF来定义复杂边界条件
工况①
描述
直径为1m,长为5m的圆柱体在一定热边界条件下分别在Steady-State和Fluent里仿真,对比圆柱中心轴线上的温度曲线,具体分为以下内容:
SST添加探针提取各热边界的热流计算结果,看整体计算是否平衡;Fluent通过flux reports获取各热边界的通量结果,看整体计算是否平衡;
整体温度分布云图和热流密度分布图,单位W/m2
SST添加path的方式获取轴线温度曲线,Fluent通过添加Line/Rake的方式获取轴线温度曲线,并导出到Matlab里面画图对比
工程原理图如下所示:
注:Fluent仿真时勾选双精度,不然误差会很大
材料参数
我们研究的是稳态传热问题,只与导热系数相关,如果是瞬态热分析,则还需要比热容Cp
默认选钢,导热系数:60.5w/(m·℃)
热边界条件
- Left:温度边界50℃
- Right:温度边界0℃
- Side:温度边界0℃
网格划分
使用自动网格划分,网格size设置为0.05m,并创建各面名称Left、Right、Side
网格信息 | SST | Fluent |
---|---|---|
Nodes | 307074 | 88944 |
Elements | 74033 | 85042 |
注:节点和元素数量不同证明两种仿真进行网格操作的内核不同
结果分析
三个边界条件热流量
边界名称 | SST仿真热流结果 | Fluent仿真热流结果 |
---|---|---|
Left | 17451.250841187943W | 20203.408W |
right | -5.6861427486147756E-07W | -5.7704047e-07W |
Side | -22395.632164042483W | -20203.408W |
SST热流相加结果等于-4.944381323423157e+03,不守恒?我们暂且不管,先看看其他的对比效果
Fluent热流相加等于2.582965e-10,基本等于0,小于通过边界最小通量的1%,满足热平衡要求,计算收敛
温度云图
SST温度云图 | Fluent温度云图 |
热流云图
SST的total Heat Flux温度云图 | Fluent的Total Surface Heat Flux |
注:图中箭头方向表示热流方向,箭头长短表示数值大小,可以看到两者完全不一样,需要后续进一步查明两者具体含义
中轴线温度分布
总结
从以上结果看,虽然SST的热流量报告显示不平衡,但其余计算结果均相同,后查阅官方文件发现,在SST中,当多个热边界条或者热约束共同作用一个边或面时,reaction heat flow rate
的报告就会变得不准,但计算结果依然正确,为了验证这一说法,又进行了工况②的仿真
工况②
处边界条件不同外其余均相同
边界条件
- Left:温度边界50℃
- Right:温度边界0℃
- Side:对流边界5W/m2,环境温度0℃
结果分析
三个边界条件热流量
边界名称 | SST仿真热流结果 | Fluent仿真热流结果 | 计算差值 |
---|---|---|---|
Left | 1368.6 W | 1367.0047W | 0.12% |
right | -156.16 W | -155.73317W | 0.27% |
Side | -1212.5W | -1211.2715W | 0.10% |
SST热流相加结果等于-0.06,基本等于0,满足热平衡要求
Fluent热流相加等于-3.03703e-09,基本等于0,小于通过边界最小通量的1%,满足热平衡要求,计算收敛
可以看到,两者计算几乎没有差别,均可以实现对工况传热的准确模拟
温度云图
SST温度云图 | Fluent温度云图 |
热流云图
SST的total Heat Flux温度云图 | Fluent的Total Surface Heat Flux |
注:图中箭头方向表示热流方向,箭头长短表示数值大小,可以看到两者完全不一样,需要后续进一步查明两者具体含义
中轴线温度分布
总结
从以上结果看,SST与Fluent在做固体传热时,只要计算收敛,计算结果基本没有差距。
写在之后
可以见得,SST在复杂热边界条件可能出现热通量报告不准,但计算仍然有效的情况;反观Fluent在两种工况里都体现里良好的热平衡求解性,但还要再去深究的是热流云图,两者的具体含义,以及Fluent的热流是否可以类似的由向量箭图的表示方式
TM01圆柱体稳态热传导对比分析相关推荐
- python有限元传热求解_二维稳态热传导基本方程的有限元求解(2)
四节点矩形单元 在二维稳态热传导基本方程的有限元求解(1)这篇文章中,我们仅仅给出了有限元单元方程的一种比较标准的推导步骤,并未涉及某种具体的单元.且在式(20)中,单元 上温度 的近似函数表示成节点 ...
- 二维稳态热传导基本方程的有限元方程
二维稳态热传导基本方程的有限元方程 - 知乎 二维稳态热传导基本方程的有限元求解(2) - 知乎 关于学习拉格朗日矩形单元和serendipity四边形单元形函数的构造方法 - 知乎 CFD理论学习- ...
- 【Comsol学习】二维稳态热传导问题
问题分析 根据comsol例程,问题给出条件如下 稳态建模仿真 选择二维稳态研究,选择固体传热物理场,开始设定各类条件 几何条件 根据问题描述构建几何模型 物理条件 按照问题描述,为矩形设置其热导率等 ...
- 二维稳态热传导 代码实现_常用振动激励方式:稳态正弦 、瞬态和随机激振
在测量机械设备或结构的振动力学参量或动态性能,如固有频率.阻尼.刚度.响应和模态等时,需要对被测对象施加一定的外力,让其作受迫振动或自由振动,以便获得相应的激励及其响应.激励方式通常可以分为稳态正弦激 ...
- matlab电力系统潮流计算软件matpower_衡真课堂 | 电力系统稳态与暂态分析
电力系统稳态分析 电力系统稳态运行指电力系统正常的.相对静止的运行状态. 电力系统稳态分析包括三个方面 ①电力系统的基本知识和等值网络: ②电力系统正常运行状况的分析和计算: ③电力系统有功功率-频率 ...
- MOOSE多物理场耦合平台入门学习记录一(稳态热传导程序实例)
MOOSE多物理场耦合平台入门学习记录 MOOSE的简介 MOOSE的安装 Linux和Mac Windows MOOSE程序的一般开发流程-以导热微分方程为例 简单问题的有限元处理 MOOSE程序的 ...
- 二维稳态热传导 代码实现_博世壁挂炉发生故障代码及处理办法
博世bosch是来自德国的知名壁挂炉品牌,由于其出色的工艺技术,故障在同类产品中属于较低一类,国内业主普遍反映较好.那么,假如我们家里的博世燃气壁挂炉出现故障代码时,该如何解决呢?这里为大家就博世壁挂 ...
- 电路复习——正弦稳态交流电路分析
正弦稳态交流电路分析 作者:Elwin [!] 以下笔记内容可能会出现部分错误的地方,恳请各位师生批评指正,谢谢! 主要内容 正弦量 正弦量的相量表示 基本元件伏安关系的相量模型 基尔霍夫定律的相量模 ...
- 用MATLAB编程正弦稳态相量图,matlab课程设计--利用MATLAB对线性电路正弦稳态特性分析...
matlab课程设计--利用MATLAB对线性电路正弦稳态特性分析 课程设计任务书 学生姓名: 专业班级: 指导教师: 刘 新 华 工作单位:信息工程学院 题 目: 利用MATLAB对线性电路正弦稳态 ...
最新文章
- Service的基本组成
- ZOJ-1654 Place the Robots---二分图最小点覆盖+构图
- sts可以安装stylus插件吗_stylus 介绍 , 安装 , 使用
- PPRevealSideViewController
- Andros系列排爆机器人原理_中国製造2025系列M之二:高档数控机床和机器人
- android开发之后端云bmob的使用
- 直播电商只能卖便宜货吗?
- java scanner 回车_Java Scanner类用法及nextLine()产生的换行符问题实例分析
- bzoj 1726: [Usaco2006 Nov]Roadblocks第二短路(A*第k短路)
- 超级素数幂--全国模拟(一)
- OOA、OOD、OOP 区别与思想
- 计算机概论综述,计算机组成原理唐朔飞PPT 第1章 计算机系统概论综述.ppt
- 日日日日日日日日日日
- 【高级篇 / SDWAN】(7.0) ❀ 03. SD-WAN 链路负载均衡的模式 ❀ FortiGate 防火墙
- 微信公众号数据2019_微信公众号榜单排名,2020微信公众号排名
- 利用python搭建“5433小游戏集成平台”
- python pdf报告_利用python设计PDF报告,jinja2,whtmltopdf,matplotlib,pandas
- 利用飞信自动发短信脚本
- 小红书血洗1w+大v,他们该何去何从?
- Java数据结构—图(Graph)
热门文章
- mysql server nchar_SQL Server NChar数据类型
- 政务服务中心工作人员是公务员吗?
- Excel制作九九乘法表
- xtrabackup出现‘FLUSH NO_WRITE_TO_BINLOG TABLES‘: 1205 (HY000) Lock wait timeout exceeded;
- P value校正思想与实现
- 傅里叶周期分析软件专业版
- Excel 用D列的值去A列中查找,将A列所对应的B列值填入E列
- AI绘图生成的PHP吉祥物素材图
- autocad常用的快捷键
- 技嘉主板开机代码15_技嘉AORUS LIQUID COOLER 360评测:LCD屏幕,信仰加倍