基于PHOENICS3.3的辐射、导热、对流联合模拟方法的研究.pdfVIP

基于PHOENICS3．3的辐射、导热、对流联合模拟方法研究 杜敬三1，李晓锋1，朱颖心1，金招芬i， (1清华大学建筑技术科学系，北京100084) 摘要： 在建筑热环境领域，为了预测和评价受太阳辐射影响比较大的建筑(例如中庭内)的热环 境，需要将辐射、对流和导热进行联合模拟计算。不少研究人员在这方面做了大量的工作， YoshllchlOzekl、Shuzo Murakami、Shinsuke Kato等人利用该算法对建筑中庭热环境进行了研 究。但是由于该方法复杂，编制通用程序困难，使得该方法至今只被极少数学者所掌握，而且 由于几种主要的商业CFD软件都没有提供适合建筑热环境研究领域的联合模拟功能，使得该 算法不能得到推广应用。而对于PHOENlCS来说，由于它本身提供易于使用的用户程序接13， 使得在CFD基础上进行联合模拟计算程序的二次开发变得相对简单。本文作者在 PHOENICS3．3的基础上，利用其程序接口，成功开发出联合模拟计算程序，既具有通用性和 整体性，又相对简单，可广泛应用于建筑热环境领域各种问题尤其是长波、短波辐射问题的研 究。并将模拟结果和典型实验结果进行比较，验证了该程序的准确性和可靠性。 关键字：中庭，太阳短波辐射，长波辐射，CFD，联合模拟 ‘。，( ， _|( 0 引言 CFD方法在建筑热环境研究领域的应用非常广泛，很多人用它来预测或者评价室内的温 在建筑热环境领域的学术研究或工程实例中，应用比较普遍。 然而对于受太阳短波辐射影响比较大的建筑，比如中庭，室内各壁面受到的太阳短波辐 射很不均匀。在太阳短波辐射、壁面间的长波辐射、墙体导热以及壁面与空气的对流换热共同 作用下，导致不同墙体壁面温度相差比较大，从而使得室内流场和温度场比较复杂。 为了能对该类建筑的热环境进行预测和评价，需要将辐射、对流和导热进行联合模拟计 算，世界上很多学者做了大量的工作。日本学者Toshiaki Omori等人利用射线追踪法计算具有 复杂形状的室内太阳得热量，为室内太阳辐射得热的计算开辟了新途径，并对位于东京的一间 在南面和西面带有大面积玻璃窗的房间进行了联合模拟计算，得到较好的结果11】：Shinsuke Kato、Shuzo 对不同的工况进行了比较【2l。联合模拟计算方法虽然在建筑热环境领域获得了成功的应用， 但是其算法复杂，编制通用程序困难，使得该方法至今只被极少数学者所掌握，为了使得该方 法获得广泛的应用，需要将该算法编制成比较通用的程序或软件。 现有的各大商业CFD软件，尽管考虑了长波辐射与对流换热的耦合作用，但是由于无法 计算复杂建筑室内的日射短波辐射得热，也不能方便地定义非均匀表面温度的边界条件，至今 都没有适合于建筑热环境领域联合模拟计算的模块或功能．而对于PHOENICS来说，由于它 本身提供易于使用的用户程序接口，使得在CFD基础上进行联合模拟计算程序的二次开发变 得相对简单。 基于以上原因，本文作者在PHoENIcS3-3的基础上，利用其程序接口，成功开发出了适 1053 合于建筑热环境领域的辐射、导热、对流联合模拟计算程序，既具有通用性和整体性，又相对 简单，可广泛应用于建筑热环境领域各种问题尤其是长波、短波辐射问题的研究。并将模拟结 果和典型实验结果进行比较，验证了该程序的准确性和可靠性。 1 联合模拟计算原理 对于中庭，其围护结构的壁面可以分为两种：(1)透明材料壁面：包括玻璃窗、玻璃屋顶、 玻璃幕墙等；(2)非透明材料壁面：屋顶、地板、墙面、家具等。 对于由透明材料构成的围护结构，其传热过程可分为三个部分： (1)辐射换热：包括太阳短波辐射换热和长波辐射换热。本文中太阳短波辐射包括：太阳直射 辐射、天空散射辐射以及其他建筑和地面反射辐射。太阳辐射到达玻璃表面后，其能量被 分为三部分：一部分被玻璃反射，一部分直接透过玻璃到达室内，另一部分被玻璃吸收。 而在玻璃内表面，接收其它表面反射和发射过来的辐射(既有短波辐射，也有长波辐射)， 并向其它表面发射长波辐射。 (2)导热：由于玻璃内外表面存在温度差，玻璃本身存在导热传热过