Establishing Boundary Conditions Considering Influence Factors of the Room Equipped with a Ceiling Radiant Cooling Panel

The objective of this study is to establish boundary conditions to evaluate the cooling capacity of the Cooling Radiant Ceiling Panel (CRCP) considering the environment of a room equipped with the CRCP. The current study investigated the boundary conditions and derivation techniques from previous studies. Based on the results of the analysis, a heat transfer model was derived for a room fitted with CRCP. In addition, the heat transfer model was used to derive the factors affecting the cooling capacity of the CRCP and each factor was simulated and verified through this model. The effects of these factors on the capacity of the CRCP was established by using various boundary conditions. To verify the validity of the simulation model, the experimental results were compared with the cooling capacity for a specific case. As a result, it was established that even for the same panel, there was a variance in the cooling capacity of the CRCP based on the boundary conditions and that the influence of the surface exposed to the outdoors had more implications. Consequently, this study presents the influence factors to be considered when designing CRCP.


Introduction
Cooling Radiant Ceiling Panel (CRCP) is a system that uses cold water to reduce the thermal load and is installed on the ceiling of buildings [1]. CRCP can be applied not only to homes but also to spaces with large ceiling areas such as classrooms [2]. In order to eliminate the thermal load efficiently, it is crucial to estimate the cooling capacity of the CRCP [3]. For accurate evaluation, it is necessary to evaluate the cooling capacity on the boundary condition when the panels are installed in the actual building ( Figure 1). In the actual building space, the CRCP uses radiation and convection heat transfer to cool the air. Therefore, it is necessary to implement the environment in which the panel can radiate and transfer convective heat as much as in the actual building space.
To evaluate the amount of heat transfer in the actual building space, the air and surface temperature should be set at the same temperature as the actual building during the evaluation ( Figure 2). However, the techniques used to evaluate the cooling capacity of the panel estimates the amount of heat that a panel removes from a thermal load generated by a cooling load simulator fitted in the test room [4]. In this case, the surface temperature and the air temperature are kept constant, in order to remove only the thermal load generated by a cooling load simulator. The same surface temperature is not applied since the variations due to the influencing factors in the surrounding environment are unknown. The aim of this study is to clarify the relationship between the influencing factors and the surface temperature, and to establish the same boundary conditions in order to evaluate the cooling capacity of the CRCP. It is expected that the actual amount of cooling capacity of the CRCP could be measured in the boundary conditions when the CRCP is installed in the actual building space for the cooling system design.

Types of CRCP
The cooling radiant ceiling panel consists of pipes that circulate cold water to the panel in order to remove the thermal load generated in the room and maintain the room in a thermally comfortable state. In order to derive the boundary condition for evaluating the cooling capacity of CRCP, it is necessary to describe the concept of the CRCP and the cooling capacity. In this study, the type of ceiling panel is defined to be evaluated based on the distinction by the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) and Representatives of European Heating and Ventilation Associations (REHVA) [4][5][6]. In both ASHRAE and REHVA, CRCPs are defined as In the actual building space, the CRCP uses radiation and convection heat transfer to cool the air. Therefore, it is necessary to implement the environment in which the panel can radiate and transfer convective heat as much as in the actual building space.
To evaluate the amount of heat transfer in the actual building space, the air and surface temperature should be set at the same temperature as the actual building during the evaluation ( Figure 2). However, the techniques used to evaluate the cooling capacity of the panel estimates the amount of heat that a panel removes from a thermal load generated by a cooling load simulator fitted in the test room [4]. In this case, the surface temperature and the air temperature are kept constant, in order to remove only the thermal load generated by a cooling load simulator. The same surface temperature is not applied since the variations due to the influencing factors in the surrounding environment are unknown. The aim of this study is to clarify the relationship between the influencing factors and the surface temperature, and to establish the same boundary conditions in order to evaluate the cooling capacity of the CRCP. It is expected that the actual amount of cooling capacity of the CRCP could be measured in the boundary conditions when the CRCP is installed in the actual building space for the cooling system design. In the actual building space, the CRCP uses radiation and convection heat transfer to cool the air. Therefore, it is necessary to implement the environment in which the panel can radiate and transfer convective heat as much as in the actual building space.
To evaluate the amount of heat transfer in the actual building space, the air and surface temperature should be set at the same temperature as the actual building during the evaluation ( Figure 2). However, the techniques used to evaluate the cooling capacity of the panel estimates the amount of heat that a panel removes from a thermal load generated by a cooling load simulator fitted in the test room [4]. In this case, the surface temperature and the air temperature are kept constant, in order to remove only the thermal load generated by a cooling load simulator. The same surface temperature is not applied since the variations due to the influencing factors in the surrounding environment are unknown. The aim of this study is to clarify the relationship between the influencing factors and the surface temperature, and to establish the same boundary conditions in order to evaluate the cooling capacity of the CRCP. It is expected that the actual amount of cooling capacity of the CRCP could be measured in the boundary conditions when the CRCP is installed in the actual building space for the cooling system design.

Types of CRCP
The cooling radiant ceiling panel consists of pipes that circulate cold water to the panel in order to remove the thermal load generated in the room and maintain the room in a thermally comfortable state. In order to derive the boundary condition for evaluating the cooling capacity of CRCP, it is necessary to describe the concept of the CRCP and the cooling capacity. In this study, the type of ceiling panel is defined to be evaluated based on the distinction by the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) and Representatives of European Heating and Ventilation Associations (REHVA) [4][5][6]. In both ASHRAE and REHVA, CRCPs are defined as

Types of CRCP
The cooling radiant ceiling panel consists of pipes that circulate cold water to the panel in order to remove the thermal load generated in the room and maintain the room in a thermally comfortable state. In order to derive the boundary condition for evaluating the cooling capacity of CRCP, it is necessary to describe the concept of the CRCP and the cooling capacity. In this study, the type of ceiling panel is defined to be evaluated based on the distinction by the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) and Representatives of European Heating and Ventilation Associations (REHVA) [4][5][6]. In both ASHRAE and REHVA, CRCPs are defined as devices that radiate heat using piping suspended at regular spaces from the ceiling ( Table 1). The CRCP consists of insulation, piping, and a heat sink. devices that radiate heat using piping suspended at regular spaces from the ceiling ( Table 1). The CRCP consists of insulation, piping, and a heat sink.

This study
The Cooling Radiant Ceiling Panel hangs on the ceiling and exchanges thermal radiation by radiant heat transfer. Since the piping is biased in the room side direction, radiation occurs in a target space.

Current Evaluation Method of the Cooling Capacity of the CRCP
In the Lawrence Berkeley National Laboratory (LBNL) study, the amount of heat transferred from the CRCP surface to the room wall surface and the room air is referred to as the cooling power of the CRCP. In this case, the amount of cooling of the CRCP is the amount of convection and radiation heat exchange due to the difference in the surface temperature and the air temperature in the installation area of the CRCP. Table 2 presents the formula used in calculating the cooling capacity of the CRCP based on the LBNL study and ASHRAE [7,8]. In both studies, the cooling capacity is divided into a convection component and a radiation component, which is a function of the temperature of each surface. The cooling capacity of the CRCP could be obtained by combining the convection and radiation components. devices that radiate heat using piping suspended at regular spaces from the ceiling ( Table 1). The CRCP consists of insulation, piping, and a heat sink.

This study
The Cooling Radiant Ceiling Panel hangs on the ceiling and exchanges thermal radiation by radiant heat transfer. Since the piping is biased in the room side direction, radiation occurs in a target space.

Current Evaluation Method of the Cooling Capacity of the CRCP
In the Lawrence Berkeley National Laboratory (LBNL) study, the amount of heat transferred from the CRCP surface to the room wall surface and the room air is referred to as the cooling power of the CRCP. In this case, the amount of cooling of the CRCP is the amount of convection and radiation heat exchange due to the difference in the surface temperature and the air temperature in the installation area of the CRCP. Table 2 presents the formula used in calculating the cooling capacity of the CRCP based on the LBNL study and ASHRAE [7,8]. In both studies, the cooling capacity is divided into a convection component and a radiation component, which is a function of the temperature of each surface. The cooling capacity of the CRCP could be obtained by combining the convection and radiation components.

Current Evaluation Method of the Cooling Capacity of the CRCP
In the Lawrence Berkeley National Laboratory (LBNL) study, the amount of heat transferred from the CRCP surface to the room wall surface and the room air is referred to as the cooling power of the CRCP. In this case, the amount of cooling of the CRCP is the amount of convection and radiation heat exchange due to the difference in the surface temperature and the air temperature in the installation area of the CRCP. Table 2 presents the formula used in calculating the cooling capacity of the CRCP based on the LBNL study and ASHRAE [7,8]. In both studies, the cooling capacity is divided into a convection component and a radiation component, which is a function of the temperature of each surface. The cooling capacity of the CRCP could be obtained by combining the convection and radiation components.  Table 3 presents the previous studies focused on evaluating the cooling capacity of CRCP. According to Kochendorfer's study, the problem is that the amount of cooling in the same panel may vary due to irregular surroundings, while the amount of cooling capacity may vary depending on the non-uniform environment [9]. Notably, previous studies have highlighted that the boundary condition could be changed by the envelope. The effect of the envelope is that the indoor boundary condition of the surface exposed to the outside air may be changed based on the thermal load generated by the envelope space where the CRCP is installed. This means that the amount of cooling of the CRCP can be changed [10]. In order to solve this problem, it is necessary to establish the non-uniform boundary condition at the environment of the actual building. Table 3. Previous studies on irregular surrounding conditions.

Author or Organization Issues
C. Kochendorfer [9] Evaluates the cooling capacity in an environment where there is no heat exchange with the outside The amount of cooling capacity may vary, based on the surface temperature of the wall N. Fonseca [3] N. Fonseca et al. [10] Estimates the cooling capacity of the panel based on the external load M. De Carli et al. [11] The cooling capacity of the ceiling radiation panel is estimated using unsteady state analysis in consideration of the external load Z. Tian et al. [12] The cooling capacity of the panel is calculated based on the window to the wall ratio considering the influence of the window-installed envelope EN14240 [13] Since the material enclosing the CRCP could affect the cooling capacity of the panel, the adiabatic boundary condition is maintained to prevent the influence of the envelope NT VVS 078 [14] Same as EN14240

Concept for Establishing Boundary Condition
The amount of cooling capacity may vary depending on the analysis of the plenum and the room. The concept for establishing the boundary condition in this study can be summarized through the analysis of the plenum and the analysis of room. Thus far, it is still unclear how the surface temperature caused by the boundary conditions affects the convection and radiative heat transfer in the room fitted with CRCP. In addition, the effect of the boundary condition influencing factors on the cooling capacity of CRCP has not been clarified. Thus, it is necessary to study the relationship between the cooling capacity and the influencing factors on the wall, ceiling, and floor surface temperature in the analysis of the plenum and room.

Heat Transfer Model of Room Equipped with CRCP
In order to establish the boundary condition for evaluating the cooling capacity of the CRCP, the influencing factors affecting the boundary conditions must be analyzed. In this study, the heat transfer process in the plenum and the room is analyzed by establishing the heat transfer model of the test room where the CRCP is installed based on the heat transfer model of ASHRAE (schematic of heat balance processes in zone) to figure out the influencing factors [15]. The heat transfer model in the space with the CRCP based on the ASHRAE model, consists of three parts: a plenum, a ceiling with a panel, and a room. In the actual building, the heat transfer occurs as shown in Figure 3. In the test room, convective heat transfer occurs between the room air and the bottom face of the ceiling, while radiative heat transfer occurs between the zone inside face and the bottom face of the ceiling. Finally, Figure 4 shows the simplified heat transfer model of the test room equipped with the CRCP.

Influencing Factors of Cooling Capacity
From the heat transfer model of a test room with CRCP, the influencing factors in the space are identified. Depending on the characteristics of these factors in the heat transfer model, the boundary condition influencing factors could be categorized into three groups: Group 1 Site and Building, Group 2 Wall Composition, and Group 3 Internal Heat Gain. By analyzing the factors belonging to each group and removing the redundant factors and simplifying them, the factors in each group can be set as in Table 4: Table 4. Group and influencing factors in each group.

Groups Influencing Factors
Group 1: Site and building Sol-air temperature Number of surfaces exposed to the outdoors Group 1 includes the sol-air temperature and the number of surfaces exposed to outdoor air. Considering the situation of the test room in which the sun cannot be accessed, the outside temperature was a substantial source of the sol-air temperature and the temperature increased with the increase in the amount of heat introduced into the room through the walls. The number of faces exposed to the outdoors is a factor that can influence the temperature difference between the outside environment and the test room. Group 2 consists of the material and composition of the wall, the heat transmission coefficient, and the surface heat absorption rate [15]. The material and composition of the wall can be determined by the structure of the wall, the build of the material, and the presence of

Influencing Factors of Cooling Capacity
From the heat transfer model of a test room with CRCP, the influencing factors in the space are identified. Depending on the characteristics of these factors in the heat transfer model, the boundary condition influencing factors could be categorized into three groups: Group 1 Site and Building, Group 2 Wall Composition, and Group 3 Internal Heat Gain. By analyzing the factors belonging to each group and removing the redundant factors and simplifying them, the factors in each group can be set as in Table 4: Table 4. Group and influencing factors in each group.

Groups Influencing Factors
Group 1: Site and building Sol-air temperature Number of surfaces exposed to the outdoors Group 1 includes the sol-air temperature and the number of surfaces exposed to outdoor air. Considering the situation of the test room in which the sun cannot be accessed, the outside temperature was a substantial source of the sol-air temperature and the temperature increased with the increase in the amount of heat introduced into the room through the walls. The number of faces exposed to the outdoors is a factor that can influence the temperature difference between the outside environment and the test room. Group 2 consists of the material and composition of the wall, the heat transmission coefficient, and the surface heat absorption rate [15]. The material and composition of the wall can be determined by the structure of the wall, the build of the material, and the presence of

Influencing Factors of Cooling Capacity
From the heat transfer model of a test room with CRCP, the influencing factors in the space are identified. Depending on the characteristics of these factors in the heat transfer model, the boundary condition influencing factors could be categorized into three groups: Group 1 Site and Building, Group 2 Wall Composition, and Group 3 Internal Heat Gain. By analyzing the factors belonging to each group and removing the redundant factors and simplifying them, the factors in each group can be set as in Table 4: Table 4. Group and influencing factors in each group.

Groups Influencing Factors
Group 1: Site and building Sol-air temperature Number of surfaces exposed to the outdoors Group 1 includes the sol-air temperature and the number of surfaces exposed to outdoor air. Considering the situation of the test room in which the sun cannot be accessed, the outside temperature was a substantial source of the sol-air temperature and the temperature increased with the increase in the amount of heat introduced into the room through the walls. The number of faces exposed to the outdoors is a factor that can influence the temperature difference between the outside environment and the test room. Group 2 consists of the material and composition of the wall, the heat transmission coefficient, and the surface heat absorption rate [15]. The material and composition of the wall can be Energies 2020, 13, 1684 6 of 21 determined by the structure of the wall, the build of the material, and the presence of either outdoor or indoor insulation. The surface heat absorption rate is a coefficient associated with radiative heat transfer, and can be used to determine how the interior surface finishing of the wall affects the boundary conditions. Group 3 consists of a boundary condition influencing factor associated with the amount and the components of internal heat gain [15][16][17][18].

Values and Ranges of Influencing Factors
The values and ranges of these factors affect every heat transfer process and are derived as shown in Table 5 based on previous research and international standard analysis [3,4,7]. The sol-air temperature of Group 1 is the ambient temperature of the building located at 36 degrees north latitude. The wall heat transmission coefficient of 2.02 W/m 2 K is the minimum design requirement for the wall heat transfer rate and 0.15 W/m 2 K is the insulation standard for passive houses. Table 5 is the reference for setting the values, ranges, and values of each factor. The Group 2 wall material is assumed to vary, whereas the rate of heat transfer is the same (e.g., concrete wall, insulation wall). In addition, the composition of the wall varies from the inside to the outside based on the construction method. By setting the heat transmission coefficient of all the walls to 2.02 W/m 2 K as in Table 6, only the influence of the boundary condition due to the wall structure, which is dependent on the material of the wall and the composition of the insulation, can be evaluated. If the material and composition of the wall are the same, the boundary condition of the wall can vary as a result of the change in the thermal conductivity based on the wall thickness. Table 7 shows that the thermal conductivity of the wall utilized for this study varies from 2.02 to 0.15 W/m 2 K by intervals of 0.07.

Concept of the Boundary Condition Derivation Model
In order to derive the boundary condition based on the influence factor changes, the CRCP installed space was modeled via simulation. Table 8 shows that the boundary condition derivation model is implemented in six steps. In step 1, the simulated input variables were derived from the influencing factor and the range of values or values of each factor. In step 2, the room fitted with CRCP and the adjacent room were modeled as in Figure 5, such that it is possible to set the desired influence factors in the space to be analyzed. In step 3, a CRCP and an air conditioning system were modeled in the main room and the adjacent room to derive the thermal load. Specifically, the air conditioning system was operated while the CRCP remained off, and the thermal load removed by the air conditioning system was estimated. The cooling capacity of the CRCP necessary to maintain the set temperature can be estimated. In steps 4 and 5, the process of finding the correct supply water temperature to satisfy the set room temperature, was reiterated by simulation using a fixed flow rate. In step 6, the boundary condition of the test room where the CRCP is installed was finally derived when the CRCP was in operation. Table 8. Modeling process to evaluate the boundary condition. Step Modeling Process

Simulation Model for Evaluating Boundary Condition
The simulation model for evaluating boundary conditions was implemented with the aim of obtaining the boundary condition when the CRCP is in operation under the various influencing factors. EnergyPlus V8-1-0 and LabVIEW were used for simulation. EnergyPlus is suited for radiation and convection heat transfer on the surface, since a heat balance-based solution is possible [19]. The boundary condition based on the varying factors can be used as the output data. In addition, the iterative calculation processes to find the supply water mentioned in step 5 of Section 4.1 were performed by linking LabVIEW and EnergyPlus. The size of the space is 4.0 m (width) × 4.0 m (depth) × 3.0 m (height), as shown in Figure 5. The CRCP and air conditioning system are installed in each zone. In this study, the CRCP was modeled using a self-developed product. Figure 6 shows the system flow chart of the CRCP installed in the test room.

Boundary Conditions on Wall Composition
• Material and composition of wall Table 9 shows the changes in surface temperature that were obtained in the analysis to determine the effect of the wall thickness and layer composition on the boundary conditions.

Simulation Model for Evaluating Boundary Condition
The simulation model for evaluating boundary conditions was implemented with the aim of obtaining the boundary condition when the CRCP is in operation under the various influencing factors. EnergyPlus V8-1-0 and LabVIEW were used for simulation. EnergyPlus is suited for radiation and convection heat transfer on the surface, since a heat balance-based solution is possible [19]. The boundary condition based on the varying factors can be used as the output data. In addition, the iterative calculation processes to find the supply water mentioned in step 5 of Section 4.1 were performed by linking LabVIEW and EnergyPlus. The size of the space is 4.0 m (width) × 4.0 m (depth) × 3.0 m (height), as shown in Figure 5. The CRCP and air conditioning system are installed in each zone. In this study, the CRCP was modeled using a self-developed product. Figure 6 shows the system flow chart of the CRCP installed in the test room.

Simulation Model for Evaluating Boundary Condition
The simulation model for evaluating boundary conditions was implemented with the aim of obtaining the boundary condition when the CRCP is in operation under the various influencing factors. EnergyPlus V8-1-0 and LabVIEW were used for simulation. EnergyPlus is suited for radiation and convection heat transfer on the surface, since a heat balance-based solution is possible [19]. The boundary condition based on the varying factors can be used as the output data. In addition, the iterative calculation processes to find the supply water mentioned in step 5 of Section 4.1 were performed by linking LabVIEW and EnergyPlus. The size of the space is 4.0 m (width) × 4.0 m (depth) × 3.0 m (height), as shown in Figure 5. The CRCP and air conditioning system are installed in each zone. In this study, the CRCP was modeled using a self-developed product. Figure 6 shows the system flow chart of the CRCP installed in the test room.

Boundary Conditions on Wall Composition
• Material and composition of wall Table 9 shows the changes in surface temperature that were obtained in the analysis to determine the effect of the wall thickness and layer composition on the boundary conditions.  Figure 6. Flow chart representing the system configuration of the simulation model.

Boundary Conditions on Wall Composition
• Material and composition of wall Table 9 shows the changes in surface temperature that were obtained in the analysis to determine the effect of the wall thickness and layer composition on the boundary conditions. The walls with the same heat transmission coefficient transmit the same amount of heat regardless of the wall composition. Therefore, in the steady-state under the same sol-air temperature, the wall internal surface temperatures having the same heat transmission coefficients have the same surface temperature.
However, in the unsteady-state where the sol-air temperature changes in the daytime, the indoor surface temperature of the surface exposed to the outside air also changes due to the heat storage. Figure 7 shows that the surface temperature of the lightweight wall is similar to the profile of the sol-air temperature, unlike in heavy wall. Since the outside insulated wall is influenced more by the heat storage, the change in the surface temperature of the indoor surface is smaller than that of the inside insulated wall. The walls with the same heat transmission coefficient transmit the same amount of heat regardless of the wall composition. Therefore, in the steady-state under the same sol-air temperature, the wall internal surface temperatures having the same heat transmission coefficients have the same surface temperature.
However, in the unsteady-state where the sol-air temperature changes in the daytime, the indoor surface temperature of the surface exposed to the outside air also changes due to the heat storage. Figure 7 shows that the surface temperature of the lightweight wall is similar to the profile of the solair temperature, unlike in heavy wall. Since the outside insulated wall is influenced more by the heat storage, the change in the surface temperature of the indoor surface is smaller than that of the inside insulated wall.

•
Heat transmission coefficient of wall Table 10 shows the evaluation of the surface temperature of the wall with the same material properties composed of a single material. The change in the heat transmission coefficient was made by changing the thickness of the insulation from 0.017 to 0.3168 mm.

•
Heat transmission coefficient of wall Table 10 shows the evaluation of the surface temperature of the wall with the same material properties composed of a single material. The change in the heat transmission coefficient was made by changing the thickness of the insulation from 0.017 to 0.3168 mm. The surface temperature of the wall exposed to the outside air is decreased from 31.9 to 26.1 • C when the heat transmission coefficient of the wall changes from 2.02 to 0.15 W/m 2 K when the sol-air temperature is 60 • C (Figure 8). This tendency is as a result of heat transmission coefficient being the determinant factor in the surface temperature of the wall in steady state conditions. When the sol-air temperature is 30 • C the ambient temperature is not affected by the outside air, regardless of the heat transmission coefficient. Figure 9 shows that for the wall that is not exposed to the outside air, with exception of the front wall, the sol-air temperature has almost no influence. The surface temperature of the wall exposed to the outside air is decreased from 31.9 to 26.1 ℃ when the heat transmission coefficient of the wall changes from 2.02 to 0.15 W/m 2 K when the sol-air temperature is 60 ℃ (Figure 8). This tendency is as a result of heat transmission coefficient being the determinant factor in the surface temperature of the wall in steady state conditions. When the sol-air temperature is 30 ℃ the ambient temperature is not affected by the outside air, regardless of the heat transmission coefficient. Figure 9 shows that for the wall that is not exposed to the outside air, with exception of the front wall, the sol-air temperature has almost no influence.

•
Surface heat absorption rate of wall Table 11 and Figure 10 show that the wall with high heat absorption rate on the surface exposed to the outside air has a large surface temperature change according to the change of the sol-air temperature.  The surface temperature of the wall exposed to the outside air is decreased from 31.9 to 26.1 ℃ when the heat transmission coefficient of the wall changes from 2.02 to 0.15 W/m 2 K when the sol-air temperature is 60 ℃ (Figure 8). This tendency is as a result of heat transmission coefficient being the determinant factor in the surface temperature of the wall in steady state conditions. When the sol-air temperature is 30 ℃ the ambient temperature is not affected by the outside air, regardless of the heat transmission coefficient. Figure 9 shows that for the wall that is not exposed to the outside air, with exception of the front wall, the sol-air temperature has almost no influence.

•
Surface heat absorption rate of wall Table 11 and Figure 10 show that the wall with high heat absorption rate on the surface exposed to the outside air has a large surface temperature change according to the change of the sol-air temperature.  Surface heat absorption rate of wall Table 11 and Figure 10 show that the wall with high heat absorption rate on the surface exposed to the outside air has a large surface temperature change according to the change of the sol-air temperature. The higher the surface absorption rate of the wall, the more the boundary conditions are affected by the outside air. When the sol-air temperature is less than 30 • C there is hardly any influence by the outside air, regardless of the surface heat absorption rate of the wall. However, as the ambient temperature increases, the influence of the surface heat absorption rate also increases. The higher the surface absorption rate of the wall, the more the boundary conditions are affected by the outside air. When the sol-air temperature is less than 30 ℃ there is hardly any influence by the outside air, regardless of the surface heat absorption rate of the wall. However, as the ambient temperature increases, the influence of the surface heat absorption rate also increases.

Boundary Conditions on Site and Building
• Sol-air temperature In the wall with the same thermal conductivity, the inner surface temperature of the exposed surface is higher, since the sol-air temperature is higher (Table 12). As the sol-air temperature rises from 30 to 60 ℃ in the wall with 2.02 W/m 2 K, the indoor surface temperature of the exposed surface also rises from 25.7 to 31.9 ℃ Conversely, when the heat transmission coefficient of the wall is as low as the level of the passive house, the sol-air temperature has little effect on the surface temperature ( Figure 11).

Boundary Conditions on Site and Building
• Sol-air temperature In the wall with the same thermal conductivity, the inner surface temperature of the exposed surface is higher, since the sol-air temperature is higher (Table 12). As the sol-air temperature rises from 30 to 60 • C in the wall with 2.02 W/m 2 K, the indoor surface temperature of the exposed surface also rises from 25.7 to 31.9 • C Conversely, when the heat transmission coefficient of the wall is as low as the level of the passive house, the sol-air temperature has little effect on the surface temperature ( Figure 11).  Figure 11. Wall heat transmission coefficient and surface temperature for sol-air temperature.
• Number of surfaces exposed to outdoor air The larger the number of exposed surfaces to the outside environment under the same sol-air temperature, the greater the increase in surface temperature. However, the temperature changes of each surface are very small (approximately less than 1 ℃). The surface temperature increases linearly with the ambient temperature, regardless of the number of surfaces exposed to the outside environment.

Boundary Conditions on Internal Heat Gain
In order to understand the effect of internal heat gain on the boundary condition, the evaluation cases are presented as in Tables 13-15, and the surface temperature and the temperature change trend in each case were derived.

Amount of internal heat gain
As the internal heat gain increases, there is a corresponding surface temperature change ( Figure  12). In particular, when the internal heat gain is 411 W, the surface temperature change is 1.8 ℃ as the radiation component of the internal heat gain changes from 0 to 1. However, when the internal heat gain is as small as 59 W, the surface temperature change due to the radiation component is within a maximum of 0.5. When the radiative component of the internal heat gain is between 0.6 and 0.7, it shows a constant surface temperature, regardless of the internal heat gain.  Figure 11. Wall heat transmission coefficient and surface temperature for sol-air temperature.
• Number of surfaces exposed to outdoor air The larger the number of exposed surfaces to the outside environment under the same sol-air temperature, the greater the increase in surface temperature. However, the temperature changes of each surface are very small (approximately less than 1 • C). The surface temperature increases linearly with the ambient temperature, regardless of the number of surfaces exposed to the outside environment.

Boundary Conditions on Internal Heat Gain
In order to understand the effect of internal heat gain on the boundary condition, the evaluation cases are presented as in Tables 13-15, and the surface temperature and the temperature change trend in each case were derived.

Amount of internal heat gain
As the internal heat gain increases, there is a corresponding surface temperature change ( Figure 12). In particular, when the internal heat gain is 411 W, the surface temperature change is 1.8 • C as the radiation component of the internal heat gain changes from 0 to 1. However, when the internal heat gain is as small as 59 W, the surface temperature change due to the radiation component is within a maximum of 0.5. When the radiative component of the internal heat gain is between 0.6 and 0.7, it shows a constant surface temperature, regardless of the internal heat gain. Table 13. Cases for amount of internal heat gain.

Group 1 Group 2 Group 3
Sol-Air Temp. If there is no surface exposed to the outdoor environment, then the surface temperature will differ only by the effect of internal heat gain. In the case that the radiation component of internal heat gain is 1, the surface temperature tends to increase with increasing internal heat generation. Alternatively, when the radiation component of the internal heat is less or equal to 0.6, the amount of internal heat and the surface temperature tend to be inversely proportional to each other. The surface temperature difference according to the radiative component of the internal heat was maximum at 3.2 ℃ which is more than that of the outdoor exposed surface. This is because the amounts of long wave absorption and re-radiation of the surface due to the internal heat gain is increased.

Component of internal heat gain
The surface temperature based on the change in radiation and convection component through various amounts of internal heat gain of between 59 and 411 W was obtained. When the radiation component of the internal heat generation was 1, the surface temperature increased with the size of the internal heat generation. On the other hand, when the radiation component of the internal heat generation was 0, the magnitude was inversely proportional to the surface temperature of the internal heat generation. Figure 13 presents the surface temperature derivation results based on the size of internal heat generation.  If there is no surface exposed to the outdoor environment, then the surface temperature will differ only by the effect of internal heat gain. In the case that the radiation component of internal heat gain is 1, the surface temperature tends to increase with increasing internal heat generation. Alternatively, when the radiation component of the internal heat is less or equal to 0.6, the amount of internal heat and the surface temperature tend to be inversely proportional to each other. The surface temperature difference according to the radiative component of the internal heat was maximum at 3.2 • C which is more than that of the outdoor exposed surface. This is because the amounts of long wave absorption and re-radiation of the surface due to the internal heat gain is increased.

Component of internal heat gain
The surface temperature based on the change in radiation and convection component through various amounts of internal heat gain of between 59 and 411 W was obtained. When the radiation component of the internal heat generation was 1, the surface temperature increased with the size of the internal heat generation. On the other hand, when the radiation component of the internal heat generation was 0, the magnitude was inversely proportional to the surface temperature of the internal heat generation. Figure 13 presents the surface temperature derivation results based on the size of internal heat generation.  Figure 13. Internal heat gain and surface temperature for radiant fraction (exposed to outdoor air).
In this case, since the change in surface temperature due to significant outdoor temperature can be excluded, it is possible to more accurately determine whether the change in the surface temperature estimated above is as a result of internal heat generation. In the case where there is no surface exposed to the outside environment, the amount of the internal heating increases, as shown in Figure 14. This tendency was least when the radiative component of the internal heat gain was 0.8.  Figure 13. Internal heat gain and surface temperature for radiant fraction (exposed to outdoor air).
In this case, since the change in surface temperature due to significant outdoor temperature can be excluded, it is possible to more accurately determine whether the change in the surface temperature estimated above is as a result of internal heat generation. In the case where there is no surface exposed to the outside environment, the amount of the internal heating increases, as shown in Figure 14. This tendency was least when the radiative component of the internal heat gain was 0.8.

Impact Analysis of Factors
When evaluating the cooling capacity, the boundary conditions have different values and ranges for each influencing factor. Therefore, the impact of influencing factors on the boundary condition should be simulated. The effect of the boundary condition influencing factors on the indoor surface temperature is derived when the scale of the factors is equally set. Figure 15 shows that the values of influence factors belonging to Groups 1-3 are normalized from 0 to 1. When the value of each influence factor changes from 0 to 1, the change in surface temperature is also normalized from 0 to 1.

Impact Analysis of Factors
When evaluating the cooling capacity, the boundary conditions have different values and ranges for each influencing factor. Therefore, the impact of influencing factors on the boundary condition should be simulated. The effect of the boundary condition influencing factors on the indoor surface temperature is derived when the scale of the factors is equally set. Figure 15 shows that the values of influence factors belonging to Groups 1-3 are normalized from 0 to 1. When the value of each influence factor changes from 0 to 1, the change in surface temperature is also normalized from 0 to 1.
When evaluating the cooling capacity, the boundary conditions have different values and ranges for each influencing factor. Therefore, the impact of influencing factors on the boundary condition should be simulated. The effect of the boundary condition influencing factors on the indoor surface temperature is derived when the scale of the factors is equally set. Figure 15 shows that the values of influence factors belonging to Groups 1-3 are normalized from 0 to 1. When the value of each influence factor changes from 0 to 1, the change in surface temperature is also normalized from 0 to 1.  Figure 16 shows that the factors that have the greatest influence on the boundary condition of the CRCP installed space is the number of the surfaces exposed to the outdoor environment. When the number of surfaces exposed to the outside air rises from 0 to 2, there is significant increase in the surface temperature. Thus, the factor 'number of surfaces exposed to the outdoors' is assigned (1.0) as the reference of impact. Subsequently, the influences are significant in the following order: 'sol-air temperature' (impact 0.91), 'heat transmission coefficient' (impact 0.85), 'surface heat absorption rate' (impact 0.63), and 'component of internal heat gain' (impact 0.28). Since the influence of the surface temperature on the factors related to the site and building belonging to Group 1 is high, when evaluating the cooling capacity of the CRCP, the surface temperature should be set based on the adjacent space and the outside air temperature. The influencing factors, which are in Group 2, are less influential compared to Group 1. The internal heat gain factors of Group 3 have relatively little effects on the surface temperature. In a room with more than one surface exposed to the outside air, the thermal load caused by the external load is dominant, compared to the other factors.  Figure 16 shows that the factors that have the greatest influence on the boundary condition of the CRCP installed space is the number of the surfaces exposed to the outdoor environment. When the number of surfaces exposed to the outside air rises from 0 to 2, there is significant increase in the surface temperature. Thus, the factor 'number of surfaces exposed to the outdoors' is assigned (1.0) as the reference of impact. Subsequently, the influences are significant in the following order: 'sol-air temperature' (impact 0.91), 'heat transmission coefficient' (impact 0.85), 'surface heat absorption rate' (impact 0.63), and 'component of internal heat gain' (impact 0.28). Since the influence of the surface temperature on the factors related to the site and building belonging to Group 1 is high, when evaluating the cooling capacity of the CRCP, the surface temperature should be set based on the adjacent space and the outside air temperature. The influencing factors, which are in Group 2, are less influential compared to Group 1. The internal heat gain factors of Group 3 have relatively little effects on the surface temperature. In a room with more than one surface exposed to the outside air, the thermal load caused by the external load is dominant, compared to the other factors.

Evaluation of Cooling Capacity According to Boundary Condition
In order to verify the boundary condition and the influences of the cooling capacity evaluation derived from this study, the cooling capacity of the panel was estimated using various boundary conditions. The simulation of the cooling capacity evaluation is based on the boundary condition of the panel modeled using the EnergyPlus simulation highlighted in Section 4.1. For the calculation of the cooling capacity, the evaluation case consisted of at least three different boundary conditions in the simulation model. The cases consist of the top two factors with the most influence. Finally, the cooling capacity of the panel could be calculated at each boundary condition. Table 16 shows the configuration of the simulation case.  Figure 16. Impact of factors on boundary condition.

Evaluation of Cooling Capacity According to Boundary Condition
In order to verify the boundary condition and the influences of the cooling capacity evaluation derived from this study, the cooling capacity of the panel was estimated using various boundary conditions. The simulation of the cooling capacity evaluation is based on the boundary condition of the panel modeled using the EnergyPlus simulation highlighted in Section 4.1. For the calculation of the cooling capacity, the evaluation case consisted of at least three different boundary conditions in the simulation model. The cases consist of the top two factors with the most influence. Finally, the cooling capacity of the panel could be calculated at each boundary condition. Table 16 shows the configuration of the simulation case.

Evaluation Results of Cooling Capacity
When the sol-air temperature increases from 30 to 60 • C the nominal cooling capacity of the panel increases from 10.75 W/m 2 to between 79.37 and 90.12 W/m 2 ( Figure 17). This tendency is more pronounced when the sol-air temperature is higher. When the number of exposed surfaces increases from 0 to 2, the nominal cooling capacity of the panel increases from 76.77 to 95.32 W/m 2 . As a result, the nominal cooling capacity of the panel is highly evaluated with the same CRCPs ( Figure 18). When the number of exposed surfaces is three or more, the supply water temperature of the CRCP becomes excessively low, and the supply water-return temperature difference becomes larger than the design standard. In this case, it is impossible to eliminate all the thermal loads generated in the room by the panel only, and a secondary system is required. When the sol-air temperature increases from 30 to 60 ℃ the nominal cooling capacity of the panel increases from 10.75 W/m 2 to between 79.37 and 90.12 W/m 2 ( Figure 17). This tendency is more pronounced when the sol-air temperature is higher. When the number of exposed surfaces increases from 0 to 2, the nominal cooling capacity of the panel increases from 76.77 to 95.32 W/m 2 . As a result, the nominal cooling capacity of the panel is highly evaluated with the same CRCPs ( Figure 18). When the number of exposed surfaces is three or more, the supply water temperature of the CRCP becomes excessively low, and the supply water-return temperature difference becomes larger than the design standard. In this case, it is impossible to eliminate all the thermal loads generated in the room by the panel only, and a secondary system is required.

Verification Experiment of Boundary Condition Derivation Model
In order to verify the simulation, the boundary condition for cooling capacity analysis was experimentally implemented using the same method, and the results were compared with the simulation results ( Figure 19) [20][21][22][23]. Table 17 shows that the verified boundary conditions in the experiment includes the flow rate, supply water temperature, the wall surface temperature, and the internal heat gain [24,25].

Verification Experiment of Boundary Condition Derivation Model
In order to verify the simulation, the boundary condition for cooling capacity analysis was experimentally implemented using the same method, and the results were compared with the simulation results ( Figure 19) [20][21][22][23]. Table 17 shows that the verified boundary conditions in the experiment includes the flow rate, supply water temperature, the wall surface temperature, and the internal heat gain [24,25].     Figure 20 shows that the cooling capacity of the CRCP is 100.1-153.0 W/m 2 at a constant flow rate of 0.153 kg/s and a supply water temperature of 11-17 • C Table 18 shows that the nominal cooling capacity of CRCP is 93.2 W/m 2 . As a result of the experiment, the difference in cooling amount from the simulation model is approximately 3.08 W/m 2 . In the simulation, the cooling capacity of the CRCP may vary based on the modeling technique, in addition to the detail of adhesion between various members and joints constituting the panel [26]. The cooling capacity estimated in the simulation and experiment has an agreement of 96.7%. As a result of the experiment, the validity of the simulation model was verified.
Energies 2020, 13, x FOR PEER REVIEW 19 of 21 Figure 20 shows that the cooling capacity of the CRCP is 100.1-153.0 W/m 2 at a constant flow rate of 0.153 kg/s and a supply water temperature of 11-17 ℃ Table 18 shows that the nominal cooling capacity of CRCP is 93.2 W/m 2 . As a result of the experiment, the difference in cooling amount from the simulation model is approximately 3.08 W/m 2 . In the simulation, the cooling capacity of the CRCP may vary based on the modeling technique, in addition to the detail of adhesion between various members and joints constituting the panel [26]. The cooling capacity estimated in the simulation and experiment has an agreement of 96.7%. As a result of the experiment, the validity of the simulation model was verified.

Conclusion
In this study, the boundary condition for the evacuation of cooling capacity of the CRCP according to various boundary condition influencing factors was derived. According to the results of this study, it can be seen that the cooling capacity of the CRCP varies depending on the boundary conditions. In particular, when the space in which the CRCP is installed is exposed to a large area of outside air, and when the outdoor temperature is high, the nominal cooling capacity of the panel could be overestimated. The reason is that excessive radiation and convective heat transfer between the test room and the CRCP may occur. Particular attention should be paid to the calculation of the cooling capacity, especially when installing the ceiling radiant cooling panels adjacent to the walls exposed to the outside air. In addition, the same results are obtained even when the surface heat absorption rate and thermal conductivity of the wall are high. The results of this study, in which the cooling capacity of the CRCPs

Conclusions
In this study, the boundary condition for the evacuation of cooling capacity of the CRCP according to various boundary condition influencing factors was derived. According to the results of this study, it can be seen that the cooling capacity of the CRCP varies depending on the boundary conditions. In particular, when the space in which the CRCP is installed is exposed to a large area of outside air, and when the outdoor temperature is high, the nominal cooling capacity of the panel could be overestimated. The reason is that excessive radiation and convective heat transfer between the test room and the CRCP may occur. Particular attention should be paid to the calculation of the cooling capacity, especially when installing the ceiling radiant cooling panels adjacent to the walls exposed to the outside air. In addition, the same results are obtained even when the surface heat absorption rate and thermal conductivity of the wall are high. The results of this study, in which the cooling capacity of the CRCPs can be estimated in various boundary conditions, are expected to be a significant design basis for applying the panels. In the future, the airflow analysis of the space with CRCPs using CFD will be carried out according to the boundary conditions derived from this study [27].