Development of a Bio-Solar House Model for Egyptian Conditions

: The need for heating and cooling in traditional housing is becoming increasingly disadvantageous regarding high energy costs. But what is more concerning is the impact on our environment. The main goal of this paper is studying the prospects of using renewable energy for heating and cooling houses through an integrated bio-solar system in order to solve the energy scarcity problem. For this purpose, a simulation model for a bio-solar house made from di ﬀ erent materials (walls made of bricks with straw bales and a roof made of concrete with straw bales) was developed successively in accordance with the energy balance and renewable energies such as biogas and solar energy were applied. This approach enabled an enhancement of the main factors a ﬀ ecting the performance of a building in terms of saving energy. The model was able to predict the energy requirements for heating and cooling of houses, the energy gained by a solar collector and by a biogas digester as well as the energy requirement for heating the biogas digester. Also, the purpose of this paper is to validate this developed simulation model by measuring energy requirements for heating of houses and solar radiation for solar collectors. The model is a simulation model for the bio-solar house with its three main parts—a straw house, a solar collector and a biogas digester. This paper demonstrates the values of the performed measurements and compares them to the theoretical, predicted values. The comparison indicates that the predicted energy requirements for the heating of buildings were a close approximation to the measured values. Another relevant deduction of the validation was the fact that the solar collector delivered the highest heat gain on 21st of June. The results also show that the energy requirements for heating the biogas digester ranged from 5.37 to 9.66, to to 13.11 and 11.01 to 15.69 kW for and of respectively, during the winter season (January) in Toukh city.


Introduction
Energy scarcity is becoming a problem that cannot be ignored. Egypt's energy balance for 2007 indicates that the largest share of final energy consumption occurs in the industrial sector (34.2%), followed by transportation (24.2%), residential buildings (18.8%) and agriculture and mining (4.7%), which altogether makes 81.9% of the total consumption. According to the International Energy Agency [1], the energy consumption measured in the residential buildings of different countries in 2004 was divided into 4% for cooking, 5% for lighting, 17% for water heating, 20% for appliances and 54% for space heating. Energy Information Administration pointed out that the consumed amounts of petroleum started to exceed the produced amounts in Egypt from the last of 2011. This contrast became especially conspicuous in 2014 where on a daily basis only 710,000 barrels were produced problem. To achieve that, a bio-solar system model of a building which can predict the energy gains and losses from buildings under Egyptian climate conditions was developed and alternative energy types of the bio-solar system were studied.

Model Development
Heat balance was performed on bio-solar houses to determine system energy gains, losses and the demand for air conditioning of the system. To maintain the required conditions inside of the building space, the conditioning system must overcome the energy loads which are imposed by the climatic conditions outside of the building and also the factors inside of the building itself.

System Description
The bio-solar house system consists of three main parts-the straw house, a solar collector and a biogas digester as shown in Figure 1. Straw bales are used as a building material or as insulation material for the original construction material. Using biogas and solar energy permits three possible air conditioning systems of the building:

•
Buildings heated by solar energy.

•
Buildings heated by both solar and biogas energy (biogas digester heated by solar energy).

•
Buildings heated by biogas energy that is heated by solar energy.
Energies 2020, 13, 817 3 of 30 energy scarcity problem. To achieve that, a bio-solar system model of a building which can predict the energy gains and losses from buildings under Egyptian climate conditions was developed and alternative energy types of the bio-solar system were studied.

Model Development
Heat balance was performed on bio-solar houses to determine system energy gains, losses and the demand for air conditioning of the system. To maintain the required conditions inside of the building space, the conditioning system must overcome the energy loads which are imposed by the climatic conditions outside of the building and also the factors inside of the building itself.

System Description
The bio-solar house system consists of three main parts-the straw house, a solar collector and a biogas digester as shown in Figure 1. Straw bales are used as a building material or as insulation material for the original construction material. Using biogas and solar energy permits three possible air conditioning systems of the building: • Buildings heated by solar energy.
• Buildings heated by both solar and biogas energy (biogas digester heated by solar energy).
• Buildings heated by biogas energy that is heated by solar energy.

Mathematical Model
The following assumptions for the development of the present model are made: • The model is based on steady-state conditions.
• The temperature of air is uniform across the whole building.
• The density and heat capacity of the air are constant.
• The temperature gradients within the constructions are negligible.
• The building walls have negligible heat capacity.
• The Egyptian time zone is GMT/UTC ± 2, so LSM for Egypt is 30° according to Reference [17]. • The distance between absorbent plate and glass cover is assumed to be 3 mm.

Mathematical Model
The following assumptions for the development of the present model are made: • The model is based on steady-state conditions. • The temperature of air is uniform across the whole building.

•
The density and heat capacity of the air are constant.

•
The temperature gradients within the constructions are negligible. The Egyptian time zone is GMT/UTC±2, so LSM for Egypt is 30 • according to Reference [17].

•
The distance between absorbent plate and glass cover is assumed to be 3 mm. • Clearness number is assumed to be 0.7 along the lines of [18].

•
The absorber plate temperature is assumed to be 90 • C or 363.15 • K and the glass cover temperature is assumed to be equal to the ambient temperature.

•
The internal gains in the digester are to be neglected.

•
The collector tilt angle is 45 • .

Houses Energy Requirements
Energy balance was performed on the house. The energy exchange due to both external and internal energy loads is shown in Figure 2. The external energy loads are transmission, fenestration, perimeter, ventilation and infiltration. The internal loads are occupants, lighting and appliances. The internal energy loads and the fenestration always count as sources of heat gain. Ventilation, infiltration, perimeter and heat transmission loads cause heat loss in winter months and heat gain in summer months. Equation (1) presents the energy balance components of the house. • Clearness number is assumed to be 0.7 along the lines of [18].
• The absorber plate temperature is assumed to be 90 °C or 363.15 °K and the glass cover temperature is assumed to be equal to the ambient temperature. • The internal gains in the digester are to be neglected.
• The collector tilt angle is 45°.

Houses Energy Requirements
Energy balance was performed on the house. The energy exchange due to both external and internal energy loads is shown in Figure 2. The external energy loads are transmission, fenestration, perimeter, ventilation and infiltration. The internal loads are occupants, lighting and appliances. The internal energy loads and the fenestration always count as sources of heat gain. Ventilation, infiltration, perimeter and heat transmission loads cause heat loss in winter months and heat gain in summer months. Equation (1) presents the energy balance components of the house. with with Q-total energy at any time in the house in J Q s -heat gains by people in W Q fens. -heat gains by fenestration in W Q light -heat gain by light in W Q appliances -heat gain by appliances in W Q conv. -heat exchange with the air by convection in W Q cond. -heat exchange with the wall by conduction in W Q p -heat exchange by perimeter of floor in W Total fenestration heat gain by the sun radiation splitting into three parts can be calculated as followed Equations (2)-(5) [17]: SHGC-solar heat gain coefficient IAC-inside shading attenuation coefficient (= 1.0 without an inside shading device) U-overall U-factor in w m −2 k −1 T out -temperature of the outer surface in K T in -temperature of the inner surface in K The radiation emitted by the sun travels through the vacuum of space unaltered. For the determination of the surface direct irradiance (ED) the following Equation (6) is needed [17]: The direct normal irradiance can be determined with Equation (7) [13]: where β-solar altitude angle, degrees C n -clearness number A-apparent solar constant in W m −2 B-atmospheric extinction coefficient (dimensionless) The angle that a ray of light falling on the irradiated surface creates with the line normal to the irradiated surface is called the angle of incidence θ. Its value is determined by [17] demonstrated in Equation (8): cos θ = cos β cos γ sin Σ + sin β cos Σ with γ-surface-solar azimuth angle in degrees -surface tilt angle from the horizontal plane (also called slope) in degrees The solar altitude angle β is calculated as follows with Equation (9) [17]: where L-local latitude angle (negative in the southern hemisphere) in degrees δ-solar declination angle in degrees H-hour angle in degrees Energies 2020, 13, 817 6 of 27 ASHRAE [19] mentioned that the earth's equatorial plane is tilted at an angle of 23.45 • to the orbital plane, so the solar declination δ is determined according to the following Equation (10) provides sufficient accuracy [19]: with n-day of the year The hour angle (H) is determined according to Equation (11) as follows: where AST-apparent solar time in hours The apparent solar time (AST), as determined by a solar time sundial, varies somewhat from the mean time of a clock kept running at a uniform rate. This variation is called the equation of time (ET) and is approximated by the following Equations (12) and (13) [19]: With ET expressed in minutes and The apparent solar time (AST) is calculated as follows (Equation (14)) [17]: where LST-local standard time in hours LON-longitude of site in degrees of Greenwich LSM-longitude of local standard time meridian in degrees of Greenwich (negative in the western hemisphere) The standard meridian longitude is related to the time zone as follows in Equation (15) [17]: TZ-time zone, expressed in hours ahead or behind coordinated universal time UTC. The surface-solar azimuth angle (γ) is defined as the angular difference between the solar azimuth φ and the surface azimuth ψ as shown in Equation (16) [17]: The solar azimuth (Φ) is calculated as follows: The apparent solar radiation (A) and the atmospheric extinction coefficient (B) are calculated by the following equations Equations (18) and (19) For vertical surfaces the diffuse irradiance E d coming from the sky is calculated following Equation (20) [18]: However, for surfaces other than vertical the calculation is carried out as shown in Equation (21): with C-sky diffuse factor (dimensionless) Y-ratio of sky diffuse on a vertical surface to the sky diffuse on a horizontal surface (dimensionless) Sky diffuse factor C is given by Equation (22) [17]: The ratio of the sky diffuse irradiation on vertical surfaces to the sky diffuse irradiation on horizontal surfaces can be calculated with the following Equation (23) [20]: The ground-reflected irradiance E r is computed as [19] and demonstrated in Equation (24): The conduction of heat between the inner surface and the outer surface of the building was calculated by Equation (25): with k-thermal conductivity coefficient in W m −1 K −1 A-the surface area in m 2 z-thickness of wall surface in m The heat transferred through convection can be calculated using Newton's Law of cooling, which is shown in Equation (26).
with h-heat transfer coefficient in W m −2 K −1 Nusselt number (Nu) correlations are traditionally used to predict a heat transfer coefficient within the Equation (27): with L c -characteristic length of the surface in m The characteristic length for a rectangular block is given as follows in Equation (28) [21]: When it comes to free convective surfaces, the Nusselt number is related to another dimensionless number, the Rayleigh number (R a ), through empirical correlations. The Rayleigh number is: with g-gravitational acceleration, 9.81 m s −2 β-coefficient of thermal expansion in K −1 v-kinematic viscosity of the fluid in m 2 s −1 µ-dynamitic viscosity of the air in m −1 s −1 C p -specific heat of air in J kg −1 K −1 The extinction coefficient is given by Equations (30) and (31): For free convection, the Nusselt number is given for horizontal or vertical planes, pipes, rectangular blocks and spheres as follows in Equations (32) and (33) [21]: If Ra between 10 4 and 10 8 Nu = 0.13Ra 1 3 (33) If Ra between 10 8 and 10 12 In cases where wind is present (i.e., forced convection), different flat plate correlations could be used but they run the risk of being unsuitable. The following Nusselt number correlation for mixed laminar and turbulent flow regions (for 5 × 10 5 < Re < 10 8 ) can be used [22].
Re-Reynold's number The previous equation is valid for Prandtl numbers between 0.6 to 60. The Reynold's number, Re, is a dimensionless number representing the ratio of inertial to viscous forces in the boundary layer of the fluid. It can be calculated according to Equation (35): Even though occupant density is low, occupancy loads should be estimated nevertheless. Sensible heat gain per sedentary occupant is assumed to be 67 W. To prevent gross oversizing, the number of occupants should not be overestimated. Recent census studies recommend that the total number of occupants should be based on two persons for the first bedroom, plus one person for each additional bedroom. The occupancy load should then be distributed equally among the living areas because the maximum load occurs when most of the residents occupy these areas [17].
The heat transfer according to perimeter can be estimated by Equation (37) for both unheated and heated slab floors [17]: with F-heat loss coefficient per foot of perimeter in W m −1 K −1 P-perimeter of floor in m The heat gain from electric lighting may be calculated by Equation (38) [20]: with W-total light wattage in W F ul -lighting use factor F sa -lighting special allowance factor The lighting use factor is the ratio of the wattage in use, for the conditions under which the load estimate is being made, to the total installed wattage.

Digester Energy Analysis
The digester has a cylinder shape. According to Fourier's law and Newton's law of cooling for convection, the thermal resistances are calculated as mentioned in Equation (26). Furthermore, according to Newton's law of cooling, the convection heat coefficient is given by Equation (27). From the air out of the digester the heat transfer is from free convection so, the Nusselt number is given by Equations (32) and (33). But, from the digestate inside the building the heat transfer is from forced convection. The Nusselt number is given in Equations (40) and (41) according to [15]: • Laminar flow: • Turbulent flow: The energy required of the digester depends on energy needed for the conditioning system and the bio-solar houses system used as shown below: • If the used system includes a building heated with solar energy and unheated biogas energy or a building heated with solar energy and biogas energy which is heated with solar energy, then the total energy needed from the digester on daily basis should be: = (Avg. yearly rate required for conditioning) × (24-Avg. yearly sunshine hours)/1000. • If the used system involves a building heated with biogas energy which is heated with solar energy, then the total energy needed from the digester on a daily basis should be: = (Avg. yearly rate required for conditioning) (24)/1000.
The volume of the biogas digester depends on this needed energy. The heating value of non-purified biogas is 6.5 kWh m −3 [23]. Necessary biogas quantity will be:
Therefore: added water rate = required feedstock rate × (total solids − 8)/8 Whereas 8 is the optimum total solids ratio for the AD process. By considering the density of digesting mix equal to water density (1000 kg m −3 ), the volume of digester will be: digester volume = (water + feedstock rate daily) * retention time/1000

Solar Collector Performance Analysis
The basic equation for calculating the performance of collectors is given in Equation (42) [24]: with Q u -useful heat gained by collector per unit of aperture area in W m −2 I t -total irradiation of collector in W m −2 τ-transmittance of cover times α-absorptance of plate U L -upward heat loss coefficient in W m −2 K −1 t p -Temperature of the absorber plate in K t out -temperature of the atmosphere in K ASHRAE [25] suggested that an additional term, the collector heat removal factor FR, should be introduced in order to permit the use of the fluid inlet temperature. This additional term is added as shown in Equation (43).
with FR-ratio of the heat delivered by the collector to the heat that would be delivered if the absorber were at t fi .
Assuming that the collector efficiency is (q u /I tθ ) = 80%, the FR ratio can be calculated in the following manner [18]: The total solar irradiation E t Equation (45) is the sum of the direct component from Equation (6) and the diffuse component in Equation (20) and Equation (21) coming from the sky plus whatever amount of reflected shortwave radiation E r may reach the surface from the earth or from adjacent surfaces [24]: The overall heat loss coefficient is given by the following expression [24]: with The coefficient of heat loss from the upper surface of the solar collector with a single cover (U t ) can be determined with Equation (47).
with h pc -heat transfer coefficient of natural convection h rpc -radiation heat transfer coefficient from the absorber to the glass cover h cw -convective heat transfer coefficient between the glass and the environment h rcw -radiation heat transfer coefficient between glass cover and sky For the heat transfer coefficient of natural convection (h pc ) the following Equation (48) applies: with L-distance between the absorber plate and glass cover The convective heat transfer coefficient between the absorber plate and the glass is calculated from the Nusselt number according to the following Equation (49) [26]: The radiation heat transfer coefficient from the absorber to the glass cover (h rpc ) can be stated as in Equation (50) [27]: The convective heat transfer coefficient between the glass and the environment (h cw ) is calculated according to Equation (51) [28]: Radiation heat transfer coefficient between glass cover and sky (h rcw ) is determined by Equation (52): A simple relationship, which ignores vapor pressure of the atmosphere, may be used to estimate the apparent sky temperature demonstrated in Equation (53) [29]: The empirical relations for U b and U e coefficients are mentioned by [29] as follows in Equations (54) and (55):

Model Approach
All computational procedures of the model were carried out using Microsoft Visual Basic 2013. The program's graphical user interface with the data input window below can be seen in Figure 3. The computer program focused on heat balance for predicting the temperature and energy gain or loss of the house. Figure 4 shows the user approach in the model. The parameters used in the model which were obtained from literature are listed in Table 1. Table 1. Parameters used in the energy balance.

Parameter
Units Value Source Degree 90 Energies 2020, 13, 817 12 of 30   The model is designed to consist of three parts: The first part was devoted to houses' performance form as shown in Figure 5. This form was used to determine the difference between building with straw and other materials by calculating the energy quantity needed in a specific building (specific area) for conditioning.  The model is designed to consist of three parts: The first part was devoted to houses' performance form as shown in Figure 5. This form was used to determine the difference between building with straw and other materials by calculating the energy quantity needed in a specific building (specific area) for conditioning. The model is designed to consist of three parts: The first part was devoted to houses' performance form as shown in Figure 5. This form was used to determine the difference between building with straw and other materials by calculating the energy quantity needed in a specific building (specific area) for conditioning.  The second part was devoted to solar collectors' form as shown in Figure 6. This form calculates the thermal energy rate that can be gained and stored when using solar collectors in specific areas.
Energies 2020, 13, 817 14 of 28 The second part was devoted to solar collectors' form as shown in Figure 6. This form calculates the thermal energy rate that can be gained and stored when using solar collectors in specific areas. The third part was devoted to the biogas digester form as shown in Figure 7. This form calculates the thermal energy that can be gained from the biogas digester by calculating the amount of biogas produced from a certain type of feedstock and its heating value (kWh) and then equivalent LPG (liquefied petroleum gas) cylinders. The third part was devoted to the biogas digester form as shown in Figure 7. This form calculates the thermal energy that can be gained from the biogas digester by calculating the amount of biogas produced from a certain type of feedstock and its heating value (kWh) and then equivalent LPG (liquefied petroleum gas) cylinders. The third part was devoted to the biogas digester form as shown in Figure 7. This form calculates the thermal energy that can be gained from the biogas digester by calculating the amount of biogas produced from a certain type of feedstock and its heating value (kWh) and then equivalent LPG (liquefied petroleum gas) cylinders.

Procedures
The objective of this study is the requirement for the surface of the collector and the volume of the digester in order to provide the bio-solar house with the energy required for heating.

Description of the Experimental Bio-Solar House
The experiment was carried out on a single-family detached house (Figure 9 The energy used for heating the house was recorded daily in order to calculate the energy requirement for heating. The recorded energy consumption was used to validate the predicted

Procedures
The objective of this study is the requirement for the surface of the collector and the volume of the digester in order to provide the bio-solar house with the energy required for heating.

Description of the Experimental Bio-Solar House
The experiment was carried out on a single-family detached house (Figure 9 respectively. The surface of windows and doors is 6.0 and 12.0 m 2 , respectively. The bio-solar house is provided with 8 LED lamps (12 Watt each) for lighting. The number of persons for each house is four.
The energy used for heating the house was recorded daily in order to calculate the energy requirement for heating. The recorded energy consumption was used to validate the predicted energy from the energy balance model.
This system provides a solar collector and a biogas digester. Table 2 shows the specifications of the bio-solar house.  This study was done using collector performance formed in the model. Table 3 shows the data of the used collector. The solar radiation was recorded from the weather station in Moshtohor. The recorded solar radiation was used to validate the predicted solar radiation from the model. This study was done using collector performance formed in the model. Table 3 shows the data of the used collector. The solar radiation was recorded from the weather station in Moshtohor. The recorded solar radiation was used to validate the predicted solar radiation from the model.

Case Studies on the Bio-Solar House
The cases, which have been studied on the bio-solar house are the following: 1.
Heating the house with solar energy.

2.
Heating the house and the biogas digester with solar energy during daylight and with biogas at night (temperature of the biogas digester is 309.15 K and retention time is 35 days).

3.
Heating the house and biogas digester with solar energy during daylight and with biogas at night (temperature of biogas digester is 322.15 K and retention time is 17.5 days).

4.
Heating the house with biogas while the biogas digester is heated by solar energy (temperature of biogas digester is 309.15 K and retention time is 35 days).

5.
Heating the house with biogas while the biogas digester is heated by solar energy (temperature of biogas digester is 322.15 K and retention time is 17.5 days). Table 4 shows the required area of solar collectors and the necessary volume of the biogas digester to provide the required energy for different case studies for the bio-solar house. It could be seen that the solar collector areas required for the first, second, third, fourth and fifth case study were 177.17, 179.06, 180.31, 4.25 and 8.50 m 2 . The volume of the biogas digester in order to provide the required energy for the bio-solar house were 9.76, 4.56, 2.28, 9.15 and 4.56 m 3 , respectively.

Solar Collector Area and Biogas Digester Volume Requirements
The results indicate that only the fourth and the fifth case can be applied. In the first, second and third case the collector area is larger than the roof area of the bio-solar house. The results also show that the energy requirements for heating the biogas digester ranged from 5.37 to 9.66, 8.02 to 11.56, 7.58 to 13.11 and 11.01 to 15.69 kW for second, third, fourth and fifth case of heating, respectively, during the winter season (January) in Toukh city. The results also show that the energy requirements for heating the biogas digester ranged from 5.37 to 9.66, 8.02 to 11.56, 7.58 to 13.11 and 11.01 to 15.69 kW for second, third, fourth and fifth case of heating, respectively, during the winter season (January) in Toukh city. Figure 11 shows the energy requirement for cooling and the energy gained by the solar collector and biogas digester for a different bio-solar house during the summer season (July) in Toukh city. The results demonstrate that the energy requirements for cooling the house ranged from 466.08 to 1031.76 kW for all case studies of the bio-solar house. The energy gained by the solar collector ranged from 396.61 to 430.52, 401.84 to 436.19, 405.30 to 439.95, 7.09 to 7.70 and 11.79 to 12.80 kW for first, second, third, fourth and fifth bio-solar house, respectively. The average energy gained by the biogas digester was 148.87, 148.87, 148.87, 297.75 and 297.75 kW for first, second, third, fourth and fifth case, respectively.

Results of Case Studies during Summer Season
The results also show that the energy requirement for heating the biogas digester ranged from 1.14 to 2.27, 3.52 to 5.18, 1.54 to 3.06 and 4.75 to 7.00 kW for second, third, fourth and fifth bio-solar house, respectively.
Energies 2020, 13, 817 20 of 28 Figure 11 shows the energy requirement for cooling and the energy gained by the solar collector and biogas digester for a different bio-solar house during the summer season (July) in Toukh city. The results demonstrate that the energy requirements for cooling the house ranged from 466.08 to 1031.76 kW for all case studies of the bio-solar house. The energy gained by the solar collector ranged from 396.61 to 430.52, 401.84 to 436.19, 405.30 to 439.95, 7.09 to 7.70 and 11.79 to 12.80 kW for first, second, third, fourth and fifth bio-solar house, respectively. The average energy gained by the biogas digester was 148.87, 148.87, 148.87, 297.75 and 297.75 kW for first, second, third, fourth and fifth case, respectively.

Results of Case Studies during Summer Season
The results also show that the energy requirement for heating the biogas digester ranged from 1.14 to 2.27, 3.52 to 5.18, 1.54 to 3.06 and 4.75 to 7.00 kW for second, third, fourth and fifth bio-solar house, respectively. 1000 kg for pig slurry, maize silage, maize grain, grass silage, fresh grass, cattle-dung, cattle-slurry, poultry, whole wheat crop, wheat grain, food waste, sewage sludge, old bread, chicken litter/dung, horse manure, clover, chaff and residuals from vegetables, respectively. The highest value of biogas production (686.0 m 3 per 1000 kg) was found by the use of clover waste, while the lowest value of biogas production (12 m 3 per 1000 kg) was found for pig slurry. 1000 kg for pig slurry, maize silage, maize grain, grass silage, fresh grass, cattle-dung, cattle-slurry, poultry, whole wheat crop, wheat grain, food waste, sewage sludge, old bread, chicken litter/dung, horse manure, clover, chaff and residuals from vegetables, respectively. The highest value of biogas production (686.0 m 3 per 1000 kg) was found by the use of clover waste, while the lowest value of biogas production (12 m 3 per 1000 kg) was found for pig slurry.      Figure 14 shows the predicted and the measured energy requirement for heating the house during the winter season (January). It can be seen that the predicted energy requirements were in a reasonable agreement with those measured, seeing as they ranged from 89.65 to 390.12 kW experimentally, while ranging 85.07 to 373.63 kW theoretically. Figure 14 shows the predicted and the measured energy requirement for heating the house during the winter season (January). It can be seen that the predicted energy requirements were in a reasonable agreement with those measured, seeing as they ranged from 89.65 to 390.12 kW experimentally, while ranging 85.07 to 373.63 kW theoretically. The variations between the predicted and measured energy requirement for heating the house during the winter season (January) are shown in Figure 15. The best fit for the relationship between the predicted and the measured value could be described as follows in Equation 56:

Energy Requirement for Heating the House
where Ep is the predicted energy requirement for heating the house in KW EM is the measured energy requirement for heating the house in KW The variations between the predicted and measured energy requirement for heating the house during the winter season (January) are shown in Figure 15. The best fit for the relationship between the predicted and the measured value could be described as follows in Equation (56): where E p is the predicted energy requirement for heating the house in KW E M is the measured energy requirement for heating the house in KW Energies 2020, 13, 817 24 of 28 Figure 15. Comparison between predicted and measured energy requirement for heating the house during winter season. Figure 16 shows the predicted and the measured solar radiation from solar collector. It can be noticed that the average solar radiation predicted by the model was in good agreement with the radiation measured by the system, seeing as it ranged from 161.00 to 561.85 W m −2 experimentally while ranging from 264.88 to 470.68 W m −2 theoretically. The stratification of solar radiation has been predicted well for the system.   Figure 16 shows the predicted and the measured solar radiation from solar collector. It can be noticed that the average solar radiation predicted by the model was in good agreement with the radiation measured by the system, seeing as it ranged from 161.00 to 561.85 W m −2 experimentally while ranging from 264.88 to 470.68 W m −2 theoretically. The stratification of solar radiation has been predicted well for the system.  Figure 16 shows the predicted and the measured solar radiation from solar collector. It can be noticed that the average solar radiation predicted by the model was in good agreement with the radiation measured by the system, seeing as it ranged from 161.00 to 561.85 W m −2 experimentally while ranging from 264.88 to 470.68 W m −2 theoretically. The stratification of solar radiation has been predicted well for the system. The variations between the predicted and measured solar radiation from the solar collector are shown in Figure 17. The best fit for the relationship between the predicted and the measured value can be determined with the following equation: The variations between the predicted and measured solar radiation from the solar collector are shown in Figure 17. The best fit for the relationship between the predicted and the measured value can be determined with the following equation:  Figure 17. Comparison between predicted and measured solar radiation. Figure 18 shows the predicted and the measured hourly average solar radiation from the solar collector. The average hourly solar radiation predicted by the model was in good agreement with the radiation measured by the system. It ranged from 39 to 634, 22 Figure 18 shows the predicted and the measured hourly average solar radiation from the solar collector. The average hourly solar radiation predicted by the model was in good agreement with the radiation measured by the system. It ranged from 39 to 634, 22

Conclusions
A simulation model for the bio-solar house was developed successively according to energy balance to optimize the main factors affecting the performance of buildings built with different materials (the walls are made of bricks with straw bales and the roof is made of concrete with straw bales). Applying renewable energies such as biogas and solar energy to this system results in energy savings and the reduction of energy related costs. The model was able to predict the energy requirements for heating and cooling of the houses, the energy gained by the solar collector, the energy gained by the biogas digester and the energy requirements for heating the biogas digester. Also, carrying out an experiment to validate the model results through measuring energy requirement for heating house and the solar radiation for the solar collectors. The most important results obtained within the frames of validation can be summarized as follows: • The required solar collector area and the biogas digester volume for providing the bio-

Conclusions
A simulation model for the bio-solar house was developed successively according to energy balance to optimize the main factors affecting the performance of buildings built with different materials (the walls are made of bricks with straw bales and the roof is made of concrete with straw bales). Applying renewable energies such as biogas and solar energy to this system results in energy savings and the reduction of energy related costs. The model was able to predict the energy requirements for heating and cooling of the houses, the energy gained by the solar collector, the energy gained by the biogas digester and the energy requirements for heating the biogas digester. Also, carrying out an experiment to validate the model results through measuring energy requirement for heating house and the solar radiation for the solar collectors. The most important results obtained within the frames of validation can be summarized as follows: • The required solar collector area and the biogas digester volume for providing the bio-solar house with the required energy were 177. 17 The highest value of biogas production (686.0 m 3 per 1000 kg) was reached using the clover waste, whereas the lowest value of biogas production (12 m 3 per 1000 kg) was reached using pig slurry. Funding: Open Access Funding by TU Wien.