首页/文章/ 详情

EnergyPlus辐射制冷屋顶节能仿真

49分钟前浏览0

关键词:EnergyPlus;辐射制冷;PDRC;屋顶模型;建筑节能

一、文章简要介绍

被动日间辐射制冷(PDRC)材料能把屋顶表面温度降到环境温度以下,是建筑节能的热门方案。但主流建筑能耗模拟软件EnergyPlus默认用恒定发射率,模拟不了PDRC在8到13微米大气窗口的波长选择性发射特性。这篇发表在Heliyon 2024(DOI: 10.1016/j.heliyon.2024.e26428)的研究,来自丹麦奥胡斯大学,提出用RC屋顶模型把波长相关的PDRC模型桥接进EnergyPlus,再系统评估不同气候、不同建造年代、不同高度的建筑应用PDRC屋顶的节能潜力。

图1 EnergyPlus与波长相关PDRC模型耦合流程图:RC屋顶模型作为桥梁,迭代交换内外表面温度边界条件直至收敛(温差小于0.1K)

图1是整体耦合流程:先在EnergyPlus里建无PDRC的建筑模型,拿到室内热边界条件(内墙和地板表面温度、室内空气温度、室内短波辐射);再把这些条件连同TMY气象文件的室外气候数据,输入耦合了已验证PDRC模型的RC屋顶模型,算出带PDRC屋顶的外表面温度;最后把外表面温度作为EnergyPlus的热边界条件回代,直到RC模型算的吊顶温度和EnergyPlus算的吊顶温度差小于0.1K才收敛。

二、仿真步骤

步骤一:RC屋顶模型搭建。典型多层屋顶结构类比电路,等效成两个热阻(内热阻Ri、外热阻Ro)加一个总热容CTotal,瞬态传热方程是CTotal乘平均温度对时间导数,等于外表面到平均温度的热流加吊顶到平均温度的热流。热阻分配用电路里常用的Elmore延迟模型,按各材料层的厚度、密度、比热容算热容,按响应时间分热阻。

图2 屋顶结构等效RC电路模型:(a)典型多层屋顶结构,(b)简化为两个热阻(Ri、Ro)加一个热容(CTotal)

图2是RC等效电路示意:左边是真实多层结构,右边是简化后的2R1C网络。这套模型的好处是既能保留屋顶的瞬态热响应,又不需要EnergyPlus原生支持波长相关发射率。

步骤二:外表面能量守恒建模。波长相关模型下,屋顶外表面释放的辐射用普朗克黑体辐射公式对发射率谱积分(式8),吸收的太阳辐射和大气辐射也按发射率谱和角度积分(式9、10),大气发射率用MODTRAN透过率和云量修正(式12到15),对流换热系数按风速线性关联(hc=3.8v+5.7)。EnergyPlus默认的恒定发射率模型则简单得多:净长波辐射等于发射率乘斯蒂芬-玻尔兹曼常数乘(表面温度四次方减有效天空温度四次方)。

图3 PDRC材料发射率谱(Raman等材料):(a)太阳短波波段低吸收,(b)大气窗口8-13微米高选择性发射率

图3是论文采用的Raman等PDRC材料发射率谱:太阳波段吸收极低(减少太阳得热),大气窗口8到13微米发射率接近1(把热量直接辐射到外太空)——这正是恒定发射率模型模拟不了的关键特性。

步骤三:EnergyPlus建筑模型与验证。200平方米足迹的典型住宅,按ASHRAE标准分5个热区(4个周边区加1个核心区);参数设置参考中国标准:制冷26摄氏度、制热18摄氏度,内热源按JGJ规范设定,人员作息用JGJ/T 449-2018的日程。三组验证全部通过:PDRC模型此前已验证;RC屋顶模型按EN ISO 52016参考建筑在广州8月与EnergyPlus对比,R平方0.98;建筑能耗模型与Deng等实测数据及GB/T 51161-2016限值对比吻合。

图4 RC屋顶模型验证:与传统屋顶(反射率0.3、恒定发射率0.9)在广州8月的表面温度对比,与EnergyPlus计算结果吻合

图4是RC屋顶模型验证结果:用反射率0.3、恒定发射率0.9的传统屋顶,把RC模型算的外表面温度和EnergyPlus直接算的结果叠在一起,两条曲线贴合得很好。

图5 RC模型与EnergyPlus模拟表面温度相关性:R平方达0.98,验证RC模型可准确捕捉屋顶温度响应

图5是相关性验证:RC模型与EnergyPlus的逐时表面温度散点高度线性相关,R平方0.98,说明简化模型能准确捕捉屋顶对环境的响应,可以放心作为PDRC与EnergyPlus的桥梁。

三、关键结果与发现

第一个发现:气候决定成败。在pre-2001低层建筑上,广州制冷需求省2.92千瓦时每平方米,上海全年总需求省1.31,北京省0.70;但哈尔滨因为制热惩罚(制热需求增加1.97)反而全年多耗1.15千瓦时每平方米,昆明也多耗0.51。图6给出五城市制热制冷需求全景——辐射制冷省的是冷,严寒地区冬季"冷源"反而成了负担。

图6 五个气候代表城市(哈尔滨/北京/上海/广州/昆明)pre-2001低层建筑(a)制热与(b)制冷需求:PDRC屋顶与传统屋顶对比

图6是五城市需求对比:左图制热、右图制冷。PDRC屋顶把外表面温度压得比环境还低,夏季是天然冷源;但冬季它也在持续向太空散热,抵消了被动得热,所以制热需求几乎必然上升。

第二个发现:建筑高度越矮越划算。上海低层(3层)建筑制冷需求省13.81%,高层(18层)只有1.77%——因为屋顶面积占外围护比例随高度下降。制热惩罚也从低层的1.36%降到高层的0.26%。结论很明确:PDRC屋顶适合低层建筑,高层建筑收益可以忽略。

第三个发现:老建筑收益大,新建筑收益小。上海3层建筑:pre-2001年建成的省1.31千瓦时每平方米、2002到2009年建成的省0.84、2010年后建成的只省0.45。原因是新建筑外墙U值从1.96降到0.8、屋顶从1.66降到0.5、制热COP从1.0提到2.2,围护结构和设备都更高效,PDRC的边际贡献自然变小。

第四个发现:保温厚度是把双刃剑。PDRC屋顶加保温,制冷需求反而上升(保温阻碍了屋顶冷量向室内传递),但制热需求大幅下降——70毫米保温下哈尔滨制热省19.13%、北京20.57%、昆明33.71%。墙体保温比屋顶保温影响更大(墙占外围护面积73.27%),但PDRC上立面有眩光问题,需谨慎。

图7 广州PDRC屋顶与传统屋顶外表面温度对比:(a)全年变化,(b)7月30日逐时——白天最大温差达21.38摄氏度

图7是广州的逐时对比:白天太阳辐射最强时表面温差最大,峰值出现在13点,最高差21.38摄氏度;夜里没有太阳,两种屋顶温度基本重合。这套"RC屋顶桥接"方法的价值在于:把波长相关发射率、大气发射率、屋顶逐层导热全部纳入EnergyPlus框架,比等效恒定发射率法更接近真实物理。做建筑节能、辐射制冷、围护结构仿真的同学可以直接复用这个耦合思路。

四、我们提供的服务

我们可提供以下仿真技术服务:①EnergyPlus建筑能耗模拟(全年逐时负荷、多气候区对比);②辐射制冷(PDRC)屋顶建模与节能评估(RC模型耦合、波长相关发射率);③建筑围护结构热工分析(保温层优化、U值敏感性);④多城市气候工况批量仿真(TMY气象数据处理);⑤建筑节能改造方案评估(不同建造年代建筑对比)。如果你也在做建筑节能或辐射制冷仿真,这篇论文的耦合建模思路很值得拆一遍。

来源:320科技工作室
电路建筑材料气象EnergyPlus
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-09-30
最近编辑:49分钟前
320科技工作室
硕士 | 结构工程师 微❤️ CAE320
获赞 243粉丝 447文章 493课程 0
点赞
收藏
作者推荐

STAR-CCM+双体船浅水阻力仿真

关键词:STAR-CCM+;浅水阻力;双体船;RANS;网格验证一、文章简要介绍太阳能双体船SolarCat跑在浅水区,阻力会涨多少?这个问题直接关系到航行安全和能耗。克罗地亚萨格勒布大学的Martic等人用STAR-CCM+对SolarCat做了系统的阻力数值评估,发表于Journal of Marine Science and Engineering 2023年11卷。计算基于RANS方程加k-omega SST湍流模型,自由液面用VOF捕捉,浅水工况还启用了网格变形算法。四个水深(无限深、h/T=7.6、4、2)两个航速(4节、5.5节)共七组算例,外加一套完整的网格无关性验证。这篇的工作量很扎实,对做船舶水动力仿真的同行有参考价值。SolarCat是太阳能电动客渡双体船,模型左右对称,仿真只取一半计算域。二、STAR-CCM+仿真设置计算域和边界条件是第一步。双体船左右对称,只建一半域。深水工况:入口和顶部离船1.5倍LPP,出口在船后3倍LPP,底部在船下2.5倍LPP,侧面离对称面2倍LPP。入口、顶部、底部用速度入口,出口用压力出口,对称面设对称边界,船体是无滑移壁面。浅水工况的底部位置按h/T比例放:h/T=7.6、4、2分别对应龙骨线下0.3、0.14、0.05倍LPP,底部也加无滑移壁面,相对速度与船速等值反向。深水工况边界:速度入口、压力出口、对称面、无滑移船体壁面,域尺寸按LPP倍数布置。有限水深工况底部按h/T比例放置并施加无滑移壁面,模拟水底对流动的约束。网格策略上,自由液面区域加密以捕捉Kelvin尾流,两个片体之间以及片体与底部之间也做了局部加密。浅水工况底部有边界层,用棱柱层网格细化。y+按ITTC建议控制:船体用壁面函数,y+落在30到300之间;底部y+控制在1以下,直接解析黏性底层。时间步长按LPP/航速除以系数设定,细、中、粗三档分别是0.01325秒、0.0265秒、0.053秒,速度亚松弛0.7、压力0.4,内迭代5次。船体表面棱柱层网格加自由液面加密,浅水工况底部再叠一层棱柱层。船体y+分布在30到300区间,满足壁面函数的适用条件。验证研究用的是GCI方法,网格加密比取根号2、时间步加密比取2。深水工况总阻力在细网格和细时间步下GCI都低于1%,中网格中时间步低于1.5%;h/T=2浅水工况网格不确定度变大,细网格GCI到3.79%,中网格4.87%,时间步GCI仍在1%上下。下沉量各工况GCI基本都在2%以内。结论是中等网格加中等时间步的配置精度可接受,后续算例统一用这套配置。h/T=2的网格截面,可见底部间隙内网格明显加密。三、关键结果与发现总阻力随水深变浅明显上升。5.5节时,h/T=7.6工况总阻力710.57牛,h/T=4是749.54牛,h/T=2涨到989.51牛,比h/T=7.6高了约40%。4节航速下h/T=2比h/T=7.6高约14.6%。水深变浅对压阻力影响远大于摩擦阻力:5.5节时h/T=4压阻力增13%,h/T=2几乎翻倍;摩擦阻力占比从h/T=4的约70%降到h/T=2的约60%,说明浅水里压阻成分在变大。下沉量方面,h/T=2在5.5节时出现明显的蹲效应,下沉约10厘米,湿表面积随之增加,摩擦阻力也往上走。4节时各水深下沉量都小很多。纵倾角基本不变,浅水对纵倾的影响可以忽略。流场细节上,水深变浅后对称面速度幅值增大,自由液面波型在片体内侧变化明显,纵向波切面的对比能直接看出浅水对兴波的影响。这套仿真流程从建模、网格到验证闭环完整,算例设置和不确定度评估方法可以直接借鉴。水深从7.6缩到2,对称面速度幅值逐渐增大,流动加速更明显。不同水深的纵向波切面直接对比,浅水兴波特性差异一目了然。来源:320科技工作室

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈