Fluent模拟太阳辐射的问题
我想做一个模型来模拟地物在自然条件下的温度分布,需要考虑太阳辐射、大气对流传导和地面辐射等因素,用Do模型能计算吗?
我自己试了一下,有几个问题没有搞明白,望高手指点:
1、如果用稳态计算出来的结果是在初始状态下太阳辐射条件(方向矢量、热流通量等)下达到的平衡状态吗?如果需要知道不同时刻的温度场,是不是应该用瞬态模型计算呢,但是用瞬态模型计算总是不收敛,这样的结果可信吗?
2、取一个比目标大几倍的空间作为计算域,空间的边界条件怎么设置,除了入口与出口,其他的边不能设置成壁面吧,设置成对称的边界行吗?
3、地面作为大空间的地面,设置成对流传热的壁面,为什么温度总是过高呢(390K)?应该怎么设置壁面的边界条件呢?
1 是平衡状态的。是应该用瞬态模型,不收敛,如果发散,当然不行。但是如果你是指未能残差降低一定的幅度,那么应该了解这种判据未必合理,应该用温度、速度或通量等更直接的判据。
2 流动是自然对流吧?是壁面可以设置成壁面,不是当然不能。你提到“边”,难道是二维模型?那么怎么模拟太阳辐射呢?
3 你提到“地面设置成对流传热的壁面”,就是说你认为地面的传热边界条件是预设的,那么有什么必要用CFD呢?
你把问题说得再明白一些。
目的是计算在太阳辐射条件下物体与地面的温度分布,最好能将地面对物体的辐射也考虑进去。
模型介绍:
地面取7*7平方米,深1米,物体为一个边长为1米的立方体,放在地面中心,空间高度取4米,大小与地面相同。空气采用材料库中的air材料,地面采用材料库中的Caco3(碳酸钙)材料,物体为自定义的cube材料。
边界条件设置:
设置地面以下1米处为恒温壁面,温度为288K,四个侧面为绝热壁面。空气与地面的接触面为耦合壁面,物体与地面接触面conduct-face为导热壁面(假设传热系数为20),物体其余5个与空气接触的面设置成辐射传热。空气设置速度入口与outflow出口,其余三个面设置为symmetry。
用Do模型自带的Solar Calculator计算阳光辐射。
以上的模型建立以及设置是否合理,望不吝赐教,谢谢!
我用这个模型计算,在稳态条件下,迭代了140次收敛(亚松弛因子等都是用的默认设置)。为什么物体的温度始终都和空气温度一样呢,而且地面温度也太高了,什么地方出了问题呢?
是漏掉了与空气的对流换热吧,我把地面、物体与空气接触的面设置成混合的边界条件,计算结果还是不理想啊,而且对流换热系数不知道,只能假设,这样也会影响计算结果啊! 另外怎么考虑地面对物体的辐射传热呢?
而且我把太阳辐射的方向角调得比较大(如下午5点)时,按理说物体各个表面的温度应该不同,但是为什么出不来这样的模拟结果呢?
分析一下传热过程和求解模型:
空气:热对流,求解N-S方程和能量方程
地面土壤、物体:导热,求解能量方程
地面与物体之间:接触传热,用薄壁热阻模型
地面或物体与空气之间:对流换热,Coupled传热边界条件
太阳对地面或物体表面的辐射:投射辐射热流,太阳辐射模型
地面或物体与周遭物体之间:辐射换热,由于处于计算域内部,必须用辐射模型。这部分传热过程是最容易被忽视的,但很多时候是不可忽略的。一般情况比较复杂。在室外,周遭物体包括地面上的物体和天空,它们的辐射性质和温度又不一样,需要适当简化与近似。 有辐射换热时,流体区边界条件不能随便用symmetry边界条件。
不好意思,前两天有点忙,没时间上网,让帖子沉下去了,呵呵
由于是新手,关于版主的解答,我还是有些地方不太明白,希望能说得更详细些。
1、建模时空气区域是否应该为一个六面体抠去一个小立方体?
2、空气区域除了入口、出口以及与地面、物体耦合的面,其余三个面如果不能设置成Symmtry,应该如何设置边界条件?
3、耦合的边界条件是这样设置吗:我发现在Gambit中将两个体共有的面合并之后,生成的网格导入到Fluent中才会有耦合选项(即自动生成一个shadow面),但是如果空气区域减去物体(立方体)后,地面就不再是正方形(中间抠去了一个小正方形),这样怎么与仍是正方形的地表耦合呢?用split face命令将地表面用立方体的地面分割成两个面,可以吗?
4、薄壁热阻模型具体在Fluent软件中怎么设置边界条件呢?本例中物体立方体的底面与地面是大小不同的面,两个面应如何设置边界条件呢?
5、版主列出了各个接触面之间的方程,分析得很全面,但是这些控制方程在Fluent软件中怎么计算呢,如“空气:热对流,求解N-S方程和能量方程; 地面土壤、物体:导热,求解能量方程; 地面与物体之间:接触传热,用薄壁热阻模型”,这些方程是分别计算还是通过Fluent软件的设置,一次就计算完成了呢?
6、投射辐射热流怎么计算啊,这个模型在Fluent中什么地方设置呢?
感谢版主耐心的支持,由于我没有学过流体力学,只是需要用Fluent计算才开始接触的这个软件,所以有许多问题比较低级,还请版主不吝赐教,感激不尽!
1、建模时空气区域是否应该为一个六面体抠去一个小立方体?
最好不是抠去,而是割(split)出一个小立方体,代表物体,并保持分割面处connected。如果抠去,那么要保留小立方体,并将其与空气区域的对应面联结(connect)起来,然后再划分网格。否则网格就是non-conformal的,要将两区域的对应面设成interface类型的边界。土壤与空气区域也最好用分割的方法。
2、空气区域除了入口、出口以及与地面、物体耦合的面,其余三个面如果不能设置成Symmtry,应该如何设置边界条件?
四周面不妨都设成流动入口或出口。设成透明的壁也行,因为你要模拟的是一个大空间,只要空间足够大。但不能设成对称,如果设成对称面,就是一面镜子,你是学光学的,对吧?
3、耦合的边界条件是这样设置吗:我发现在Gambit中将两个体共有的面合并之后,生成的网格导入到Fluent中才会有耦合选项(即自动生成一个shadow面),但是如果空气区域减去物体(立方体)后,地面就不再是正方形(中间抠去了一个小正方形),这样怎么与仍是正方形的地表耦合呢?用split face命令将地表面用立方体的地面分割成两个面,可以吗? 参见1。地面与物体,以及地面与空气区域之间的面均设为wall边界,这应该是两个面。
4、薄壁热阻模型具体在Fluent软件中怎么设置边界条件呢?本例中物体立方体的底面与地面是大小不同的面,两个面应如何设置边界条件呢?
将对应wall边界的thickness设定为非零值,并选择相应的material,以模拟热阻。
5、版主列出了各个接触面之间的方程,分析得很全面,但是这些控制方程在Fluent软件中怎么计算呢,如“空气:热对流,求解N-S方程和能量方程; 地面土壤、物体:导热,求解能量方程; 地面与物体之间:接触传热,用薄壁热阻模型”,这些方程是分别计算还是通过Fluent软件的设置,一次就计算完成了呢?
我列出的是你的问题中涉及的主要物理过程和相应的物理模型。对于传热问题,fluent求解的基本控制方程是N-S方程和能量方程,这也就是最基本的物理模型,其它物理模型根据问题的具体情况是可选的。一旦在fluent中选择或启用了某些物理模型,也就定义了计算模型,就意味着让fluent求解相应的控制方程。
6、投射辐射热流怎么计算啊,这个模型在Fluent中什么地方设置呢?
太阳辐射可以用fluent的太阳辐射模型。如果还需要考虑其它物体之间的辐射则需要选择一种辐射模型。如果仅需考虑计算域外边界壁面与计算域外部具有均一温度且壁面对其辐射角系数为1的黑体之间的热辐射传热,且假定壁面面向外部的面为漫灰表面,只需要定义壁面辐射边界条件。
谢谢版主的解答!
对于第二个问题,如果四周面都设成流动入口或出口,计算时会出现reversed flow,这样很难收敛,对计算结果应该也有影响吧
对于第六个问题,“还需要考虑其它物体之间的辐射则需要选择一种辐射模型”,此处说的
辐射模型是指Fluent里自带的P1、Do等模型吗,我用的就是Do模型,是不是这样就默认为考虑了地面与物体之间的辐射呢?
前几天怎么论坛打不开啊,终于能进来了,呵呵
这两天计算了一下,空气设置为透明壁面挺好的。但是出现另外一个问题,最高温度总是出现在地面以下约10厘米处,我分析了一下,主要原因一方面是地面受到太阳辐射后热量向地下传导,另一方面是空气对流对地面有冷却作用。向下传导是正常的,但是实际情况下,空气对地面的冷却作用不应该有这么明显,而且地面附近的空气温度应该很高才对。而我在速度入口的边界条件处设置的温度是常数,这个问题应该怎么解决呢?是不是要通过udf文件来定义入口处空气温度随高度变化呢?
Fluent模拟太阳辐射的问题
最近想用Fluent做一个模拟太阳辐射的模型,但自己对Fluent所知甚少,在版主xrs333不厌其烦的指导下,终于建立起了模型,在此向xrs333致敬!这是我自己建立模型后的计算结果:
参数设置:Do模型,solar caculator地点设置为北京的经纬度,Gmt=+8,时间为6月21日上午10点;入口风速1m/s,温度288K,RNG湍流模型;物体与地面外表面都为耦合传热,空气外壁设为透明。
材料性质:
空气:吸收系数1,散射系数0;物体吸收系数0.1,散射系数0;地面:吸收系数0.2,散射系数0
(关于材料的吸收系数与散射系数的设置有点不明白,我看Fluent User's Guide中将吸收系数定义为每单位长度上辐射强度的变化,并给出了朗伯吸收定律I=I0*exp(-ax),(原文:The absorption coefficient is requested in units of 1/length. Along with the scattering coefficient, it describes the change in radiation intensity per unit length along the path through the fluid medium.)按照给出的公式来看,吸收系数a在负指数项上,应该是a越大,辐射强度减小得越快,但是在计算中我发现,将空气的吸收系数设置得越小(如0.001),到达地面的热量几乎没有,设置得越大(如100),地面温度越高,这和公式上的正好相反,是不是我对公式的理解有误呢?
分析稳态计算结果,我发现最高温度不是出现在地面或物体的表面上,而是在地面以下约10cm处(如图),这与实际情况也不符合,这个问题应该如何解决呢?
相关文章
- 石蜡相变蓄热过程数值模拟
- CFD仿真
- 教学楼中庭的自然通风与热环境分析
- 毕业设计说明书(论文)
- 对圆柱体网格划分的一些经验总结
- 随着电子科技水平的不断进步
- Fluent理论手册(2)-旋转坐标系
- 管内湍流的数值模拟
- 第01次课教案-绪论及导热理论基础
第24卷第2期 2008年6月 北京建筑工程学院学报 JournalofBeijingUniversityofCivilEngineeringandArchitecture V01.24NO.2 Jun.2008 文章编号:1004-601 ...
3.1气体泄漏扩散的模拟方法 目前在研究气体扩散领域应用较多的模拟方法主要有三种,即:物理模拟方法. 数学模拟方法和CFD 数值模拟方法.当然在实际的模拟仿真过程中,经常是两种或是三种方法同时使用,以此来验证模拟的准确性. 3.1.1物理模 ...
2013年第2期(总第41卷第264期) No.2in 2013(TotalNo.264,Vol.41) doi :10.3969/j.issn.1673-7237.2013.02.001 建筑节能 ■暖通与空调 HEATING, VENT ...
武汉理工大学 毕业设计(论文) GP车型排气系统 主消声器内流场分析 学院(系): 汽车工程学院 专业班级: 热动0602班 学生姓名: X X X 指导教师: X X X 武汉理工大学学士论文原创性声明 本人郑重声明:所呈交的学位论文,是 ...
对圆柱体网格划分的一些经验总结(转) jiaogl - BY - 2007-6-20 15:15:00 最近一段时间在做辐射废锅换热的研究.在对其流场进行数值模拟的过程中,在Gambit 中试验了一些关于圆柱体的网格划分方法,并将其导入Fl ...
因此,为适应现代电子设备的冷却需要而迅速发展起来的热设计及热分析技术,受到了 广泛的重视.电子设备热技术主要包括热设计.热分析和热测试. 热分析是指用理论或数值模拟的方法对将要设计的电子设备进行分析,它有两条重要的 要求:[6] (1) 预 ...
2 旋转参考系流动 本节主要描述旋转参考系流动模拟的数学背景.以下章节包含的内容包括: ● 2.1:介绍 ● 2.2:旋转参考系流动 ● 2.3 多旋转参考系流动 2.1 简介 ANSYS FLUENT求解流体流动及热传递方程,默认情况下求 ...
管内湍流的数值模拟 摘 要:当Reynolds数大于临界值时,平滑流动会出现一系列复杂的变化,最终会导致流动特征的本质变化,流动呈无序的混乱状态,这种状态称为湍流.计算流体力学是通过计算机数值计算和图像显示,对包含有流体流动和热传导等相关物 ...
课程内容及任务 热传导.对流换热.辐射换热基本概念.基本理论.分析解决问题的基本方法和思路.计算传热学基础. 通过经典传热现象的分析,加强其研究思路.方法的学习,为今后实际问题分析研究打下理论,尤其是研究方法的基础. 数值计算:有限元计算软 ...