1. Introduction
Passive heating systems such as the Trombe wall stand out as essential building blocks on the road to zero energy buildings. Quasi-steady models [
1,
2], which allow relatively fast calculations for the initial design, are forced into simplifications by adopting general values for the heat transfer coefficients. Additionally, they cannot offer information about design details that could improve the passive element’s energy performance. Energy balance models [
3,
4,
5] and building energy simulation (BES) models [
6,
7] consider one-dimensional heat transfer, make general assumptions about the flow pattern and the heat transfer within the air gap while using coefficients for taking into account the temperature gradient along the air gap [
2].
Although several researchers are considering several modified forms of the Trombe wall [
8], issues regarding the optimal value of key geometric features of conventional Trombe walls as presented in an extensive literature review [
9] have not yet been resolved. For example, the optimal width of the air gap is the 1/10 of its height according to [
10]. Thus, the optimal width is between 0.2 and 0.3 m [
10]. However, other researchers propose that the width of the gap must be limited from 2 to 5 cm in order to create a small space for air [
6], or even 6 cm [
11]. In terms of the size of the ventilation slots, the relevant research focuses mainly on the estimation of their sizes as a percentage of the Trombe wall surface [
10], while it has not been examined whether their configuration can affect the performance of the wall. The optimal thickness of the Trombe wall depends on the climatic conditions and according to [
12] the optimum thickness is between 15 to 40 cm. Finally, the effect that the distance between the ventilation slots has on the operation of the Trombe wall does not appear to have been investigated in detail. Therefore, there are still gaps in research on the basic geometric features of the conventional Trombe wall.
Computational Fluid Dynamic (CFD), although a time consuming and computationally demanding method for designing and evaluating the operation of a Trombe wall, can guide the design of details (shape of Trombe wall absorption surface, opening geometry, etc.). Furthermore, CFD can offer insight into the flow inside the air gap and around ventilation slots. Consequently, the results can feed quasi-steady models with tested heat transfer coefficients and predict the flow pattern, allowing energy balance models to calculate these coefficients correctly. In addition, CFD can be used to assess and/or feed with info and/or boundary conditions energy balance and BES models [
13].
The majority of CFD studies in recent years use the Finite Volume method for solving the partial differentiate RANS equations [
13,
14,
15,
16,
17,
18], while in older works, researchers also use the Finite Difference Method [
4,
19]. The majority of the published research papers use steady-state calculations [
13,
15,
17,
19], and only a few of them, due to high computational requirements, use unsteady calculations [
14,
20,
21]. The high computational cost, as well as the particular characteristics of the Trombe wall operation to which each researcher wishes to focus his attention, lead several researchers to 2D calculations [
16,
17,
18,
20,
22]. Nevertheless, in recent years, research with 3D CFD simulations of the Trombe wall are also available [
13,
14,
15,
19]. Some CFD studies simulate the operation of solar air heaters in conjunction with the space they serve [
13,
14,
17,
21]. Other CFD studies focus on the operation and optimization of specific parts of the passive system, very often focusing on the air gap without simulating the solid parts [
15,
16,
18,
19].
What is of a major importance is the way researchers treat the radiation (solar and thermal). Several published research works do not simulate all the incident solar radiation, considering known values of radiation or heat fluxes on the storage wall [
15,
16,
18,
19,
21,
22]. Among them, some focus on the heat transfer with convection and conduction while others simulate the heat transfer of thermal radiation with models like S2S [
16,
18]. Nevertheless, there are works that include in the computational domain the transparent cover and the storage wall and they simulate their operation considering as boundary conditions directly the solar incident radiation, using for the simulation the Monte Carlo [
20] method or more often the DO radiation model [
11,
14,
17]. The basic common characteristic of these works is that they consider the transparent cover as gray material with average optical properties along the whole spectrum. This approach cannot simulate accurately the operation of the transparent cover. Another drawback is the fact that they simulate the storage wall as an ‘opaque’ element, which in the DO model means that the participation of this component to heat exchange through the specular component of reflected radiation is limited. In addition, there is no possibility to set different values for absorption and emission coefficients. According to the extended review presented in [
8] the majority of CFD studies use the k-ε model in order to simulate the turbulence effect in the transport phenomena, using grids that vary from 4500 to 1,500,000 computational cells according to the specific needs, the computational domain size and the characteristics of each study.
Solar air heaters are a practical system for producing warm air using solar energy and can be considered as aerothermal systems. The principle of their operation is similar to that of solar thermal panels. From numerous studies concerning solar collectors, information can be collected and assessed in relation to the increase of their efficiency. This can be realized by increasing heat transfer through the proper treatment of the absorption surface (roughness, fins, etc.) [
23]. In [
24] the heat transfer enhancement is studied using fins and rods on the collector surface. Kabeel et al. [
25] examined the use of fins, the use of multiple air paths, and the use of wires perpendicular to the air flow in order to increase the turbulence and thus the transfer coefficient as well as the heat exchange surface.
Studies of solar chimneys could be useful in understanding the simulation of the Trombe wall operation. In [
26] a solar chimney system with earth tubes is studied where the CFD model is validated with the equivalent experimental configuration. A similar study of a solar chimney is made in [
27].
In the present work, a 2D CFD model is developed for simulation of Trombe wall operation serving a test room. The transport phenomena are modeled with the finite volume method while the radiation transport is modeled with the DO model. Special attention has been given on radiation heat transfer considering non-gray cover material with detailed calculation of optical properties of the composite cover in four wavelength bands (Ultra Violet, Visual, Near Infrared and Infrared). This, simulates accurately the transportation of the radiation. Additionally, we treat the storage wall as a ‘semi-transparent’ material improving its participation in the heat transfer calculation. To our knowledge this treatment of radiation heat transfer has not been used in the past by other researchers for the simulation of Trombe wall operation. After the validation of a default geometry, other 9 Trombe wall configurations are examined investigating the influence of air gap width, storage wall width, ventilation slots size, ventilation slots distance, and of the ventilation slots geometry on the Trombe wall basic characteristics (flow pattern/flow rate inside the air gap and heat transfer coefficients). Thus, the present study contributes, through the development of the CFD model, to the current understanding of how the basic geometric features of the conventional Trombe can affect its performance.
2. Materials and Methods
2.1. Mathematical Model
The transport phenomena developed in the solar air heater system are described with the Reynolds Averaged Navier-Stokes equations which are solved with the finite volume method [
28,
29]. Specifically, the equations of continuity, conservation of momentum in two directions (x, y), energy, transport radiation, and the equations corresponding to the turbulence model are solved in a steady-state regime. The simplified approach of a 2D model is considered since in a cross-section including the ventilation slots with buoyancy driven flow, the, buoyancy forces in the horizontal direction may be neglected. This is a common approach adopted by many researchers [
22]. The flow is considered incompressible since, in case of natural convection, the expected value of velocities is not expected to reach a Mach number greater than 0.1. For the modeling of thermal buoyancy, we used the Boussinesq approximation since it allows us to consider the density constant in all the solved equations except for the buoyancy term in the momentum equations. Thus, the computational effort is reduced while a faster convergence is achieved. For the simulation of radiative transfer, the Discrete Ordinate (DO) model was used [
30,
31] with angular discretization 4 × 4 and pixelation 3 × 3. Radiation transport equations are solved in two directions for four wave bands: (a) Ultra-violet radiation–UV (λ = 0.1–0.39 μm), (b) Visual radiation -VIS (λ = 0.38–0.74 μm), (c) Near Infrared radiation–NIR (λ = 0.74–1.2 μm), and (d) Infrared radiation–IR (λ = 1.2–100 μm). Energy and radiation equations are also solved for the solid parts of the computational domain. The ability of the transparent cover to allow the propagation of shortwave solar radiation and not the longwave thermal radiation due to the correct modeling of its spectral optical properties is crucial for the analysis of Trombe wall performance.
The flow was considered turbulent because the expected velocities within the air gap are expected to lead to values greater than 10
9, as evidenced by the presented results. The turbulence effect is modeled through the high Reynolds k-ε model [
32]. The standard high Re k-ε turbulence model was selected because it offers numerical stability, quick convergence, and, using wall functions, does not require very dense space discretization close to the solid boundaries (as in the case of low Reynolds models). At the same time, it offers satisfactory accuracy for calculations in large spaces [
33]. This is the reason that it is used extensively in calculations that include a Trombe wall coupled with the room they serve [
8]. The limited requirements of k-ε model in the number of computational cells were crucial for its selection since a very detailed discretization for the radiation model can be applied. Although this increases the computational cost, it offers an accurate solution of the radiation transmission which is the most important factor for the specific calculations.
2.2. Numerical Model
A grid with 22,000 rectangular computational cells was used. A growth ratio of 1.025 was used for the grid thickening, securing that the dimensionless distance y+ of the first computational cell from the wall boundaries lies between 10 and 40 since wall functions are used. This decision was based on a grid independence test which was carried out using three grids of 14,640, 22,000, and 34,200 computational cells, ensuring that in all grids the requirement regarding the value of y+ will be met. The value of air mass flow inside the air gap had a 3% difference between the grid of 14,640 computational cells and the grid of 22,000. Refining further the grid to 34,200 computational cells the change of air mass flow was lower than 1%. Thus, a grid of 22,000 computational cells with adequate density in critical areas was selected for our calculations.
For the pressure/momentum coupling the SIMPLEC algorithm was used. The discretization for convection terms of momentum and turbulence quantities (k and ε) is implemented using a 2nd Order Upwind (SOU) scheme, which takes into account the values of upstream cells, providing second order accuracy. It is considered adequate for parabolic flows, as is a flow driven by thermal buoyancy. For the convection terms of energy equation, the 3rd order Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL) is used, combining central difference and upwind schemes characteristics. Finally, for Radiation Transport Equations (RTEs) a 1st Order Upwind scheme (FOU) is selected since it provides adequate accuracy along with fast convergence. All the diffusive terms are discretized with a central difference scheme which is a second order of accuracy considering the values contribution of all neighbor cells as required by the nature of diffusive terms. The calculations were performed by the ANSYS Fluent software V14.5. The convergence criterion for all the parameters except energy and radiation was set to 10−4, while for energy was set to 10−6 and for radiation to 10−5. In the 1st geometry, the appropriateness of the convergence criterion was tested for the continuity equation, by setting its convergence criterion to 10−6 and monitoring as a control parameter the mass flow rate inside the air gap. The convergence criterion refinement modified the calculated mass flow rate by only 0.44%, while the execution elongation roughly doubled, without any significant improvement to the accuracy of the calculations. For that reason, the continuity criterion convergence was preserved in the initial value for all the other simulations.
2.3. Case Study and Parametric Study
The examined Trombe wall consists of a heat storage wall with a high absorbing solar radiation coefficient which is covered by a transparent cover. In the gap between them, air can circulate. This air enters the gap through ventilation slots in the lower part of the storage wall and leaves the gap through the upper ventilation slots returning to the room.
2.3.1. Studied Geometry
The 2D model simulates the operation of a Trombe wall cross-section which includes the computational domains of the transparent cover, the 7.6 cm width air gap, the storage/collection wall having a 10 cm thickness, the lower and upper ventilation slots which have a 12.4 cm diameter and the room coupled with the Trombe wall. The room has dimensions of 3 m (length) by 2.8 m (height) as presented in
Figure 1a, while in
Figure 1b the used computational grid is presented. The examined cross section is in the North-South axes, with the Trombe wall covering the south wall.
2.3.2. Boundary Conditions
The transparent cover is modeled as a ‘sandwich’ consisting of an outer surface, an inner solid semi-transparent element, and an inner surface. The outer surface (computational domain boundary) is a semi-transparent wall where a mixed thermal boundary condition is applied (combined convection and radiation heat transfer with environment and conduction with the solid material). In the cover outer surface are defined: (a) the external air temperature, (b) the convection heat transfer coefficient with the external air and (c) the equivalent sky temperature for the calculation of the thermal radiation exchange with the external environment. Additionally, on the outer surface, the solar incident radiation (beam and diffuse) is defined in four wavelength bands. The inner surface of the cover is a semi-transparent wall thermically and optically coupled with the solid material of the cover and the air in the gap between the cover and the storage wall restoring a conjugated heat transfer treatment. The storage wall is also modeled as a ‘sandwich’ consisting of an inner (gap-facing surface) surface, an interior solid semi-transparent material, and an outer (room-facing surface). Both storage wall surfaces are thermically and optically coupled with the interior solid material and with the air in the airgap and the room. The solid material of the storage wall has an absorption coefficient high enough to ensure the extinguish of incident radiation in the first computational cell. The energy and radiation equations are solved both in the fluid and solid zones. All the other room internal surfaces are considered adiabatic.
2.3.3. Optical Properties
The transparent cover consists of two panes of glass with an internal air gap. However, in the simulation, it is treated as a single material with the same thickness and equivalent thermal and optical properties.
Given the transmittance,
τti and reflectance,
ρti of the individual glass panes [
34,
35,
36,
37,
38,
39,
40], the transmittance,
τt, absorptance,
αt, and reflectance,
ρt, of the entire cover (two glazing system and air gap) for vertical incident radiation [
34] are calculated as follows:
Then, with the following relations, the equivalent refractive index,
n, and the equivalent extinction coefficient,
αs, for the whole structure are calculated.
where,
r, the cover reflection, and
τα the transmittance due to absorption losses.
where
α, the cover absorptance, and d, the cover thickness
These calculations are made for the optical properties of the cover system in each of the examined wavelength bands.
According to the above relationships, the simulation of the double-glazed cover 4-15-5 with air gap and low-e surface is materialized with an optically equivalent cover with optical properties given in
Table 1.
For the storage wall the extinction coefficient αs is taken equal to 1000 and the refractive index, n, equal to 1.418 in the whole examined wavelength range.
2.3.4. Thermal Properties
The following
Table 2 gives the thermal properties of the materials used.
2.3.5. Climatic Conditions
The external temperature was taken 14.97 °C. The wind velocity 2.3 m/s. The incident radiation normal to the transparent cover, beam, and diffusive, are given in the following
Table 3 for each considered wavelength band.
2.3.6. Parametric Study
Apart from the default geometry, presented in the previous paragraphs, henceforth it will be referred to as grid 1, nine (9) more geometries were examined in order to study the Trombe wall geometry effect on the wall performance. The nine geometries’ characteristics are summarized in the following
Table 4. In detail, the cases that have been examined are the following: (a) One case with smaller air gap width and 8 cases with bigger width, (b) Three cases with different storage wall thicknesses, (c) Two cases with different ventilation slots’ distances, (d) Two cases with different ventilation slots’ diameter and (e) A case with different upper ventilation slot geometry.
4. Discussion
The purpose of the Trombe wall is to heat up a space utilizing incident sunlight. The heat transfer to the coupled room is achieved either due to convection through the warm air that enters the room from the upper ventilation slot (therefore high flow and high temperature of the exhaust air are required) and through the room-facing storage wall surface due to convection and radiation (therefore the high temperature of this surface is required). However, it is important that this heat penetrates as well as possible into the room and only the warm air stream from the upper ventilation slot can contribute to this. Therefore, in addition to high flow and temperature, equally important is the flow field that develops and the degree to which the room is actually ventilated by the air-gap without leaving large areas in which the existing air simply recirculates. The operation of the Trombe wall associated with heat storage and its time lag performance cannot be studied by steady-state simulation. The findings of the CFD simulations will be assessed with respect to the aforementioned criteria.
It is observed that a small gap width leads to a relatively small flow. Since this reduction is not proportional to the available cross-section it leads to higher air velocities inside the air gap which in turn enhance the heat transfer by convection from the surfaces of the cover and the wall. This along with the fact that the same amount of radiation is used for the heating of a smaller volume of air leads to higher exit air temperature. Inversely, the increase of the airgap width increases the mass air flow while decreasing the gap exit temperature and increasing the temperature of the room facing storage wall surface. This is due to the fact that lower velocities observed inside the gap reduce the rate of heat transfer by convection from the wall surface. A further increase in width leads to an increase in flow but not to an increase in the room facing storage wall surface temperature since larger amounts of air in contact with the warm surfaces gaining larger amounts of heat by convection. It seems that the 7.6 cm air gap can balance the heat transfer mechanisms ensuring high temperatures both at the air exit and at the surface that is in contact with the room with a relatively small flow.
Increasing the thickness of the Trombe decreases the pressure drop at the exit and inlet ventilation slots and thus for the same thermal buoyancy higher mass flow rates are achieved. At the same time increasing the thickness of the storage wall leads to lower air temperature at the exit ventilation slot because a bigger amount of air absorbs the same available solar radiation (through radiation and convection from the gap’s solid surfaces). At the same time the temperature of the gap facing storage wall surface increases because a thicker storage wall leads to the increased thermal resistance of the wall preventing heat transferred by conduction to the room facing wall surface. However, the increase of storage wall thickness allows the air from the gap to penetrate deep into the room improving thus the Trombe wall effectiveness, a conclusion consistent with the results of previous research [
19].
The only significant change resulting from the variation of the diameter of the ventilation slots is related to the velocity distribution at the air exit from the gap. The small observed reduction in temperature in the exit cross-section equals the corresponding small increase in mass flow rate. The result is that there is no significant effect on the efficiency with which the Trombe wall ventilates the coupled room. Finally, the surface temperature of the storage wall does not appear to be significantly affected.
Increasing the distance between the ventilation slots both the air mass flow and the storage wall temperature increase since a larger solar collection area is available, while the storage wall temperature increase is small. Inversely, the temperature at which the air leaves the air gap decreases by a small percentage while its distribution in the exit slot cross section is not significantly affected. What is important is that the ability of Trombe wall to ventilate the coupled room is improved.
The sloping exit slot configuration leads to higher flow rates for specific choices of air gap width and storage wall thickness since it reduces the local pressure drop due to flow direction change. This way improves the operation of Trombe wall and allows the warm air stream to penetrate deep into the coupled room without leaving large areas with recirculation loops.
5. Conclusions
A 2D CFD model was developed for the simulation of the Trombe wall cross-section operation. The developed model is able to model adequately the transparent cover behavior in different radiation wavelengths (four wavelength radiation bands were considered) which is crucial for the simulation of Trombe wall operation as well as the adequate participation of the storage wall in the heat transfer phenomena. The CFD model presentation was accompanied by the description of the calculation of equivalent cover optical properties since the latter is modeled as a single material. The model was applied to 10 geometrical configurations of a Trombe wall and the results were presented and discussed in terms of achieved air mass flow rate, temperature, velocity distributions, and flow fields including in the computational domain the room served by the examined Trombe wall.
According to the presented research findings when high air exit temperatures are required, a small gap width of 5 cm is suggested which offers increased temperature by 7% with respect to a gap width of 7.6 cm and by 11% with respect to a gap of 10 cm. When high room-facing storage wall surface temperatures are required a gap width of 7–8 cm is suggested which gives surface temperature increase by 12% with respect to the 5 cm gap.
However, these options leave a significant part of the serviced room without heating through the warm air stream and in addition lead to low mass air flow, reduction of the gap width by 50% leads to a reduction of the supply by 25%. This can be corrected by increasing the distance between the slots, which presents an increase in the air mass flow rate by 0.34%/cm, or by increasing the storage wall thickness. A storage wall with a thickness of 30 cm gives an increase of air mass flow rate by 73% with respect to a wall with a thickness of 15 cm. Finally, a much better choice is the proper shaping of the exit ventilation slot with an inclination of 30 degrees to achieve an increase of flow rate by 35% with respect to the horizontal slots and gives optimal penetration of warm air flow into the room maintaining relatively high air temperatures.