Experimental Investigation and Energy Performance Simulation of Mongolian Ger with ETS Heater and Solar PV in Ulaanbaatar City

: There are approximately 200,000 households living in detached houses and gers (yurts) with small coal stoves that burn raw coal in Ulaanbaatar city. A proper heating system and improvement of the energy e ﬃ ciency of residential dwellings are vitally important for Ulaanbaatar city to reduce air pollution as well as for the operation of the current central energy system. This study shows the experimental results for two gers with two di ﬀ erent heating systems and di ﬀ erent thermal insulation, for investigating the merits of each. The technical feasibility of the system consisting of an electric thermal storage (ETS) heater with a daytime charging schedule and areal photovoltaic (PV) system was also examined by using a simulation with software developed in MATLAB (R2020a, MathWorks, USA). As a result of the experiment, the indoor comfort level and energy e ﬃ ciency of the ger with added insulation and an ETS heater with nighttime charging were shown to be enhanced compared with those of the reference ger. The ger with added insulation and the ETS heater consumed 3169 kWh for electric appliances and 5989 kWh for the heating season. The simulation showed that the PV self-consumption rate is 76% for the Ger 2 with the ETS heater because of the daytime charging schedule of the ETS heater. The PV system supplied 31% of the total energy consumed, with the remaining 69% from the main grid.


Introduction
Energy-efficient buildings with proper heating systems reduce air pollution, fuel consumption, and electricity costs. In addition, with rapidly growing world energy use, climate change is necessitating energy-efficient buildings, proper heating systems, and energy sources for them [1].
Ulaanbaatar city is one of the coldest cities in the world and has a long space-heating season, which continues annually from October to May. About 70% of the total energy consumed is used for space heating. There are, in total, 184,000 households that have 3-7 kW individual coal-burning stoves. Their pollutant emissions have recently been becoming a big problem and impacting the health of citizens. About 40% of these households live in Mongolian gers (yurts) [2,3].
Therefore, improving the energy efficiency of Mongolian gers and the choice of proper heating systems is vitally important for reducing pollutant emissions [4,5]. This paper investigates the energy performance of Mongolian gers with electric thermal storage (ETS) heaters and added envelope insulation through comparative experimentation with a reference ger, which has a traditional heater and envelope materials. Benoit et al. introduced the design and architectural data for gers in his article about the sun clock, while Tsovoodavaa et al. provided more detailed historical, design, and structural information on the yurt [6,7]. Additionally, her extended article presents the results of thermal dynamic simulations in the Mongolian climate [8].
In 1980, the ETS heater was designed, in which the electricity converted by a heating element to thermal energy is accumulated in a magnesium brick for later use. Various developments in terms of design, numerical and experimental studies based on the designed ETS heater, performance evaluation, and feasibility studies of the ETS heater were carried out [9][10][11][12][13]. Additionally, it plays a big role in the integration of renewable energy generation with the main grid. It also brings various benefits to customers such a lowered electricity cost and to the main grid by load shifting and peak shaving [14,15]. In the second section of this study, a 5 kW ETS heater, which is to be installed in the ger with additional insulation and supplied from a grid-connected 3.12 kW p solar photovoltaic (PV) (hereafter referred as the proposed system), was simulated for evaluating its technical feasibility.

Experimentation with the Insulated Ger with ETS Heater
We installed two Mongolian gers for evaluating the enhancements of energy performance (47 • 52 8" N and 107 • 7 22" E). Figure 1 shows a view of the installed gers at the Urgakhnaran test site. The ger shown in Figure 1a has standard envelope materials and a coal-burning stove (hereafter referred to as Ger 1). The ger in Figure 1b has added envelope materials and the ETS heater (hereafter referred to as Ger 2).
Energies 2020, 13, x FOR PEER REVIEW 2 of 14 insulation through comparative experimentation with a reference ger, which has a traditional heater and envelope materials. Benoit et al. introduced the design and architectural data for gers in his article about the sun clock, while Tsovoodavaa et al. provided more detailed historical, design, and structural information on the yurt [6,7]. Additionally, her extended article presents the results of thermal dynamic simulations in the Mongolian climate [8].
In 1980, the ETS heater was designed, in which the electricity converted by a heating element to thermal energy is accumulated in a magnesium brick for later use. Various developments in terms of design, numerical and experimental studies based on the designed ETS heater, performance evaluation, and feasibility studies of the ETS heater were carried out [9][10][11][12][13]. Additionally, it plays a big role in the integration of renewable energy generation with the main grid. It also brings various benefits to customers such a lowered electricity cost and to the main grid by load shifting and peak shaving [14,15]. In the second section of this study, a 5 kW ETS heater, which is to be installed in the ger with additional insulation and supplied from a grid-connected 3.12 kW solar photovoltaic (PV) (hereafter referred as the proposed system), was simulated for evaluating its technical feasibility.

Experimentation with the Insulated Ger with ETS Heater
We installed two Mongolian gers for evaluating the enhancements of energy performance (47° 52′ 8″ N and 107° 7′ 22″ E). Figure 1 shows a view of the installed gers at the Urgakhnaran test site. The ger shown in Figure  1a has standard envelope materials and a coal-burning stove (hereafter referred to as Ger 1). The ger in Figure 1b has added envelope materials and the ETS heater (hereafter referred to as Ger 2). The thermal and technical details of the envelope materials of these gers are explained in Appendix A. Here, the added insulation material on the envelope of Ger 2 is highlighted. The crown wheel of Ger 2 is insulated by a 4 mm polycarbonate sheet, while 3 mm aluminum foils are aligned through the roof, wall, door, door frame, and floor for blocking the radiation of thermal energy in Ger 2. We used 30 mm thicker felt for the wall and 100 mm thicker Styrofoam for the floor of Ger 2 than Ger 1′s. Figure 2 shows an electrical diagram of the experimental setup and electrical wiring of the two gers. As shown in Figure 2a, the light, television, refrigerator, and outlet, etc.-the major electrical appliances-of Ger 1 were supplied from the main grid (1), and the total electricity consumption of both gers was measured using a power meter (2) with a three-phase current sensor (3). Voltage (10), current (9), temperature (12), and humidity (11) sensors attached to module loggers for the voltage (4), current (5), humidity (6), and temperature (7) measured the electricity consumption, indoor air temperature, and humidity. Figure 2b shows an electrical diagram of Ger 2.
The electricity consumption of the ETS heater was measured by voltage (4) and current (5) loggers with a voltage probe (9) and current clamp (10). The communication base (8) collected all the The thermal and technical details of the envelope materials of these gers are explained in Appendix A. Here, the added insulation material on the envelope of Ger 2 is highlighted. The crown wheel of Ger 2 is insulated by a 4 mm polycarbonate sheet, while 3 mm aluminum foils are aligned through the roof, wall, door, door frame, and floor for blocking the radiation of thermal energy in Ger 2. We used 30 mm thicker felt for the wall and 100 mm thicker Styrofoam for the floor of Ger 2 than Ger 1 s. Figure 2 shows an electrical diagram of the experimental setup and electrical wiring of the two gers. As shown in Figure 2a, the light, television, refrigerator, and outlet, etc.-the major electrical appliances-of Ger 1 were supplied from the main grid (1), and the total electricity consumption of both gers was measured using a power meter (2) with a three-phase current sensor (3). Voltage (10), current (9), temperature (12), and humidity (11) sensors attached to module loggers for the voltage (4), current (5), humidity (6), and temperature (7) measured the electricity consumption, indoor air temperature, and humidity. Figure 2b shows an electrical diagram of Ger 2.
Energies 2020, 13, 5840 3 of 13 parameters via radio signals. The measurement resolution was 2 min. The measurement continued between 28 October 2017 and 8 June 2018 (refer to Table A2 of Appendix B).
Ger 1 used a lignite coal for heating and cooking. We measured the coal weight just before use in the coal stove. The measurement continued for 20 days in the annually coldest month, which is January. The lignite coal was extracted from the Baganuur coal reserve, and the coal weight was converted to kWh units according to Kyle's converter [16].  Table 1 shows the electrical and technical details of the ETS heater installed in Ger 2. The ETS heater is factory equipped with an electronically integrated room-temperature control unit and a clock-signal-controlled relay. The charging mode of the ETS heater is controlled through the charging intensity selected by the customer, and the off-peak time for the electricity load in the main grid (between 21:00 and 05:00). For the period of the off-peak tariff, the K1 relay switches on the thermal relay of the heating element. The charging intensity is set by using a 47 kohm potentiometer. If the signal received from the core temperature sensor PT100 exceeds a set value of resistance, the temperature regulator disconnects. The maximum value of the core temperature is 460 °C. If the room temperature sinks below the temperature set by the 47 kohm potentiometer, the temperature regulator switches on the heater fan, dissipates warm air into the room, and turns off at the set temperature. The 450 kohm resistance regulates the fan speed. Immediate heating on demand can provide a day-acting element. The day-acting element is run using an integrated rocker switch and works on the on-peak tariff.  The electricity consumption of the ETS heater was measured by voltage (4) and current (5) loggers with a voltage probe (9) and current clamp (10). The communication base (8) collected all the parameters via radio signals. The measurement resolution was 2 min. The measurement continued between 28 October 2017 and 8 June 2018 (refer to Table A2 of Appendix B).
Ger 1 used a lignite coal for heating and cooking. We measured the coal weight just before use in the coal stove. The measurement continued for 20 days in the annually coldest month, which is January. The lignite coal was extracted from the Baganuur coal reserve, and the coal weight was converted to kWh units according to Kyle's converter [16]. Table 1 shows the electrical and technical details of the ETS heater installed in Ger 2. The ETS heater is factory equipped with an electronically integrated room-temperature control unit and a clock-signal-controlled relay. The charging mode of the ETS heater is controlled through the charging intensity selected by the customer, and the off-peak time for the electricity load in the main grid (between 21:00 and 05:00). For the period of the off-peak tariff, the K1 relay switches on the thermal relay of the heating element. The charging intensity is set by using a 47 potentiometer. If the signal received from the core temperature sensor PT100 exceeds a set value of resistance, the temperature regulator disconnects. The maximum value of the core temperature is 460 • C. If the room temperature sinks below the temperature set by the 47 potentiometer, the temperature regulator switches on the heater fan, dissipates warm air into the room, and turns off at the set temperature. The 450 resistance regulates the fan speed. Immediate heating on demand can provide a day-acting element. The day-acting element is run using an integrated rocker switch and works on the on-peak tariff.

Simulation of the Proposed System
We simulated the operation of the solar PV combined with the ETS heater for evaluating its energy performance for the heating season. The simulation was carried out by using software developed in MATLAB. The system configuration was imagined consisting of a 5 kW electric heater with 40 kWh of thermal storage and a 3.12 kW PV system. We proposed using the ETS heater specified in Table 1 in the ger insulated according to Appendix A.
The electricity-generating solar PV system was connected to the main grid. The rated nominal capacity of the PV array is 3.12 kW; the inclined angle of the PV array is 45 • .
The energy performance simulation procedure for the ger with the ETS heater and solar PV is presented in Figure 3. The simulation data generally consisted of the total electricity consumption and the estimated output of the solar PV. The data resolution was 2 min, and the data set had the same period as the measurements from experimentation. The sum of the electricity consumption of the electric appliances and the generated electricity consumption of the ETS heater made up the total electricity consumption. Here, we used the electricity consumption of the electric appliances in Ger 2 measured from the experimentation. The load profile of the ETS heater scheduled to charge in the daytime was generated from the estimated daily heating demand of the Mongolian ger, which had an envelope material specified in Appendix A.

Simulation of the Proposed System
We simulated the operation of the solar PV combined with the ETS heater for evaluating its energy performance for the heating season. The simulation was carried out by using software developed in MATLAB. The system configuration was imagined consisting of a 5 kW electric heater with 40 kWh of thermal storage and a 3.12 kW PV system. We proposed using the ETS heater specified in Table 1 in the ger insulated according to Appendix A.
The electricity-generating solar PV system was connected to the main grid. The rated nominal capacity of the PV array is 3.12 kW; the inclined angle of the PV array is 45°.
The energy performance simulation procedure for the ger with the ETS heater and solar PV is presented in Figure 3. The simulation data generally consisted of the total electricity consumption and the estimated output of the solar PV. The data resolution was 2 min, and the data set had the same period as the measurements from experimentation. The sum of the electricity consumption of the electric appliances and the generated electricity consumption of the ETS heater made up the total electricity consumption. Here, we used the electricity consumption of the electric appliances in Ger 2 measured from the experimentation. The load profile of the ETS heater scheduled to charge in the daytime was generated from the estimated daily heating demand of the Mongolian ger, which had an envelope material specified in Appendix A. The estimation method for heating demand is defined in Section 2.3, while Section 2.5 details the estimation of solar PV output (refer to Table A2 of Appendix B).
The electricity flow from the PV system to the main grid can be calculated by subtracting the total electricity consumed in the ger from the output power of the PV system. The amount of electricity consumed in the ger from the PV system can be expressed by subtracting the electricity flow from the PV system to the main grid from the output power of the PV system. The amount of electricity consumed in the ger from the main grid can be expressed by subtracting the electricity of the PV system consumed in the ger from the total electricity consumed in the ger. See Equations (1)- (3). The estimation method for heating demand is defined in Section 2.3, while Section 2.5 details the estimation of solar PV output (refer to Table A2 of Appendix B).
The electricity flow from the PV system to the main grid can be calculated by subtracting the total electricity consumed in the ger from the output power of the PV system. The amount of electricity consumed in the ger from the PV system can be expressed by subtracting the electricity flow from the PV system to the main grid from the output power of the PV system. The amount of electricity consumed in the ger from the main grid can be expressed by subtracting the electricity of the PV system consumed in the ger from the total electricity consumed in the ger. See Equations (1)- (3).
P pv ger,t = P pv,t − P pv grid,t (2) Energies 2020, 13, 5840 where P pv grid,t stands for the electricity flow from the solar PV to the main grid at timestep t (kW), P pv,t stands for the electricity generated from the PV system at timestep t (kW), and P ger,t stands for the electricity consumed in the ger at timestep t (kW). P pv ger,t is the electricity consumed in the ger from the PV system. P grid,t stands for the electricity consumed in the ger from the main grid.
The daily amounts of P pv grid,t , P pv ger,t , and P grid,t were calculated from the sums of 24 h values. The monthly average daily amounts can be expressed by dividing the monthly total value by the days per month. See Equation (4), where X t stands for the amounts of P pv grid,t , P pv ger,t , and P grid,t at timestep t (kW). T stands for the number of hours per day (24). D stands for the number of days per month. d stands for the day number of the month.

Estimation of Heating Demand
There are various modelling techniques, which are generally categorized as bottom-up and top-down, available [17]. Building physics-based modelling is used to estimate heating demand. The heating demand of a ger, in another words, the thermal energy supply from the ETS heater (H ETS ) to the ger, can be expressed by subtracting the conduction (H Con ), ventilation (H Ven ), infiltration (H In f ), etc. heat losses through building envelope elements from the heat gains including irradiance (H Irr ), occupants (H Occ ), and electric appliances (H App ) (refer to Equation (5). The heat loss conducted through walls, doors, crown wheels, ceilings, floors, etc. can be calculated as expressed in Equation (6). There is 15% extra heat, which is the radiation loss through the roofs to outside that is kept inside the ger by the aluminum foil materials. Since the ger has natural ventilation, Equation (7) can express the ventilation heat loss, whereas the heat loss caused by infiltration can be calculated as shown in Equation (8).
H ETS,t = H Con,t + H Ven,t + H In f ,t + H Occ,t + H App,t + H Irr,t H Con,t = A·U· T ind,t − T air,t H Ven,t = C p ·ρ·q v · T ind,t − T air,t H Con,t = C p ·ρ·n·V· T ind,t − T air,t H Occ,t = n adult,t ·b adult,t ·P child + n child,t ·b child,t ·P child (9) where A is the area of exposed surface (m 2 ), U is the overall heat transmission coefficient (W/m 2 K), and T ind,t is the indoor air temperature (20 • C). T out,t is the outdoor air temperature at timestep t. C p is the specific heat capacity of the air (1000 W/kgK); ρ is the density of the air (1.2 kg/m 3 ). q v is the air volume flow (m 3 /s). n is the number of air replacements in the room per second. V is the room volume (m 3 ). n adult,t is the number of adults at the timestep t. n child,t is the number of children at the timestep t. b adult,t and b child,t are binary numbers that express occupied (1) or unoccupied (0). P App,i is the energy rate of the ith electric appliance. F i, t is the heat gain rate for the ith electric appliance (0.734). b i, t is a binary number that indicates either an on or off state for the ith electric appliance. G GHI is the global horizontal irradiance (W/m 2 ). A CW is the surface area of the crown wheel.
Refer to Table A2. The measured ambient air temperature is used for estimating the conduction, ventilation, and infiltration heat losses. Heat gain due to irradiance is estimated by using global horizontal irradiance (GHI). Physical characteristic data for the Mongolian ger and thermal and technical data of its envelope materials are used in the calculations of heat gain and losses.
Estimating the conduction heat loss through the floor is not the same as estimating the conduction heat loss above ground. The estimation of the conduction heat loss through the floor involves two significant difficulties, including the fact that soil has a specific heat, so the heat can both flow and be stored as it flows. Second, the soil temperature changes with both the season and depth from the surface. However, the heat conducted from the ger to the ground (H f ,t ) can be calculated with the following simplified Equation (12). The conductance and the earth temperature are modelled according to Equations (13) and (14) [18].
where U f ,t is the effective heat transmission coefficient (W/m 2 K). R f is the thermal resistivity of the ger floor (m 2 K/W). T e is the effective earth temperature ( • C). T oa,yr is the annual average ambient air temperature ( • C). T amp is the amplitude of the surface temperature ( • C). D is the depth below the surface. α is the thermal diffusivity of the ground soil (10 −6 m 2 /s). t year is the current time. t shi f t is the day of the year with the minimum surface temperature.
The modelling of the earth temperature uses typical meteorological year (TMY) data. Detailed information including the rate of power, time schedule of the occupants, and operation time of major electric appliances are presented in Table 2. Equations (5)- (14) are entered into the software developed in MATLAB.

Load Profile Generation of ETS Heater Daytime Charging
The load profile of the ger with daytime charging of the ETS heater was generated for simulating the proposed system. The daily overall heating demand of the ger was calculated from the estimated heating demand with a 2 min resolution (refer to Equation (15)). The charging hours per day can be calculated with Equation (16). The binary number with a 2 min interval was created from the calculated charging hours of the ETS heater. The ETS charging hours were symmetrically aligned before and after midday (12:00 PM. Thus, the ETS heater has the ability to accumulate the electricity from the PV system as much as possible, and the PV self-consumption increases. The load profile with the daytime charging can be expressed as a multiplication between the created binary numbers and nominal rated power of the ETS heater (refer to Equation (17)).
where T is the number of measurements per hour. P ETS,nom is the nominal rated power of the ETS heater. b charge,t is the binary file, which indicates on (1) and off (0) states for the ETS heater.

Estimation of Solar PV Output
Here, we estimated the output power of the 3.12 kW solar PV array with a 45 • inclined angle for simulating the energy performance of the proposed system.
The plane of array (POA) irradiance, wind speed, and ambient air temperature of Ulaanbaatar were used for modelling the array temperature and output power of the PV. The data resolution was 2 min, while the data collection continued from 28 October of 2017 to 8 June of 2018. Equation (18) was referred to for estimating the output power of the solar PV. P pv,t = P pv,stc ·(G POA /1000)·(1 + γ(T A + 25))·K (18) T A = G POA · e a+b * WS ·T air (19) where P pv,stc is the rated capacity under the standard test condition (STC), γ is the temperature correction coefficient of the maximum power, and K is the loss factor including the incident angle, soil, snow, shading, degradation, conversion, and other unknown losses (0.87). G POA is the irradiance received at the 45 • inclined plane of the array. T A is the array temperature. The estimation of the array temperature is expressed by Equation (19). T air is the ambient air temperature, WS is the wind speed, and b and a are parameters that depend on the array construction and materials as well as on the mounting configuration of the module (−3.56 and −0.075) [19]. Figure 4 shows the indoor heat indices of the two experimental gers against the standard heat index. The indoor air temperature of Ger 1 measured across a wider range than that of Ger 2, while the humidity of Ger 2 was increased compared to that of Ger 1, as shown by the heat indices of Ger 1 and 2 illustrated by black triangles and green triangles.

Experimental Results
where is the number of measurements per hour.
, is the nominal rated power of the ETS heater.
, is the binary file, which indicates on (1) and off (0) states for the ETS heater.

Estimation of Solar PV Output
Here, we estimated the output power of the 3.12 kW solar PV array with a 45° inclined angle for simulating the energy performance of the proposed system.
The plane of array (POA) irradiance, wind speed, and ambient air temperature of Ulaanbaatar were used for modelling the array temperature and output power of the PV. The data resolution was 2 min, while the data collection continued from 28 October of 2017 to 8 June of 2018. Equation (18) was referred to for estimating the output power of the solar PV.
where , is the rated capacity under the standard test condition (STC), is the temperature correction coefficient of the maximum power, and is the loss factor including the incident angle, soil, snow, shading, degradation, conversion, and other unknown losses (0.87).
is the irradiance received at the 45° inclined plane of the array.
is the array temperature. The estimation of the array temperature is expressed by Equation (19).
is the ambient air temperature, is the wind speed, and and are parameters that depend on the array construction and materials as well as on the mounting configuration of the module (−3.56 and −0.075) [19]. Figure 4 shows the indoor heat indices of the two experimental gers against the standard heat index. The indoor air temperature of Ger 1 measured across a wider range than that of Ger 2, while the humidity of Ger 2 was increased compared to that of Ger 1, as shown by the heat indices of Ger 1 and 2 illustrated by black triangles and green triangles.

Experimental Results
The wide fluctuation of temperature can be explained through the difference between the electronically integrated room temperature control of the ETS heater in Ger 2 and the uncontrolled combustion of coal in the small coal stove in Ger 1. Since conduction, ventilation, and infiltration heat gains are direct functions of the indoor air temperature, the heat consumption also follows the temperature rise. The aluminum foil material could have increased the indoor air humidity of Ger 2. Furthermore, the radiation barrier materials should be carefully selected for insulating the ger.  The wide fluctuation of temperature can be explained through the difference between the electronically integrated room temperature control of the ETS heater in Ger 2 and the uncontrolled combustion of coal in the small coal stove in Ger 1. Since conduction, ventilation, and infiltration heat gains are direct functions of the indoor air temperature, the heat consumption also follows the temperature rise. The aluminum foil material could have increased the indoor air humidity of Ger 2. Furthermore, the radiation barrier materials should be carefully selected for insulating the ger.
Energies 2020, 13, 5840 8 of 13 Figure 5 shows the daily energy performance of Ger 1 and 2 on the annually coldest days of the year, which are January 1 to January 20. The grey-and blue-colored bars indicates the daily heating consumption and electricity consumption of the electric appliances of Ger 1, and the daily energy consumption of Ger 2 is presented by green-and aqua-colored bars for the ETS heater and electric appliances. the other hand, the total heating consumption of Ger 1 is 30% higher than the electricity consumption of the ETS heater. The added insulation of Ger 2, and the electronically integrated room-temperature control unit of the ETS heater etc., could reduce the heating consumption. Households living in detached houses and gers should use the various heaters equipped with electronically integrated room-temperature control units and clock-signal-controlled relays. This is because they allow energy saving for space heating and effectively interact with the central energy system (CES).
The ETS heater was set to charge the thermal energy at nighttime. Even the electricity tariff from 21:00 PM to 06:00 AM is exempted, in accordance with the energy regulations set by the government; as 0 ₮ for the period of heating, there is a lack of implementation of the three tariff; people still prefer cheap coal with its low operational costs vs. the initial cost of the ETS heater and electricity meter. Therefore, there are economic challenges. Ger 1 uses 3880 kg/yr of lignite-coal and 490 kg/yr of wood for space heating and cooking [20]. Meanwhile, the electric appliances and ETS heater of Ger 2 used 3169 kWh (electric appliances) and 5989 kWh (ETS heater) of electricity for the period of experimentation. The changing of a heating system from a coal-burning stove to an ETS heater that is supplied from a thermal power plant (TPP) reduces the amount of pollutant emission. Direct combustion in coal stoves emits approximately 787.56 kg/yr of pollutant emissions including total suspended particles (TSP), SOx, CO, and particulate matter (PM) 2.5, while the emission of Ger 2 for the period of experimentation was 40.84 kg. The annual amount of Ger 1's emission and emission amount per kWh of the TPP were obtained from [20]. Figure 6 shows a daily operational profile of the ger with the ETS heater and solar PV for one week, which is 5th of February to 12th of February, as a simulation result. The validation results for The cooking and hot water are supplied by the coal combustion in the coal-burning stove. Thus, the electricity consumption of the electric appliances in Ger 2 is~10% higher than that in Ger 1. On the other hand, the total heating consumption of Ger 1 is 30% higher than the electricity consumption of the ETS heater. The added insulation of Ger 2, and the electronically integrated room-temperature control unit of the ETS heater etc., could reduce the heating consumption. Households living in detached houses and gers should use the various heaters equipped with electronically integrated room-temperature control units and clock-signal-controlled relays. This is because they allow energy saving for space heating and effectively interact with the central energy system (CES).

Simulation Results
The ETS heater was set to charge the thermal energy at nighttime. Even the electricity tariff from 21:00 PM to 06:00 AM is exempted, in accordance with the energy regulations set by the government; as 0 Energies 2020, 13, x FOR PEER REVIEW 8 of 14 Figure 5 shows the daily energy performance of Ger 1 and 2 on the annually coldest days of the year, which are January 1 to January 20. The grey-and blue-colored bars indicates the daily heating consumption and electricity consumption of the electric appliances of Ger 1, and the daily energy consumption of Ger 2 is presented by green-and aqua-colored bars for the ETS heater and electric appliances.
The cooking and hot water are supplied by the coal combustion in the coal-burning stove. Thus, the electricity consumption of the electric appliances in Ger 2 is ~10% higher than that in Ger 1. On the other hand, the total heating consumption of Ger 1 is 30% higher than the electricity consumption of the ETS heater. The added insulation of Ger 2, and the electronically integrated room-temperature control unit of the ETS heater etc., could reduce the heating consumption. Households living in detached houses and gers should use the various heaters equipped with electronically integrated room-temperature control units and clock-signal-controlled relays. This is because they allow energy saving for space heating and effectively interact with the central energy system (CES).
The ETS heater was set to charge the thermal energy at nighttime. Even the electricity tariff from 21:00 PM to 06:00 AM is exempted, in accordance with the energy regulations set by the government; as 0 ₮ for the period of heating, there is a lack of implementation of the three tariff; people still prefer cheap coal with its low operational costs vs. the initial cost of the ETS heater and electricity meter. Therefore, there are economic challenges. for the period of heating, there is a lack of implementation of the three tariff; people still prefer cheap coal with its low operational costs vs. the initial cost of the ETS heater and electricity meter. Therefore, there are economic challenges.
Ger 1 uses 3880 kg/yr of lignite-coal and 490 kg/yr of wood for space heating and cooking [20]. Meanwhile, the electric appliances and ETS heater of Ger 2 used 3169 kWh (electric appliances) and 5989 kWh (ETS heater) of electricity for the period of experimentation. The changing of a heating system from a coal-burning stove to an ETS heater that is supplied from a thermal power plant (TPP) reduces the amount of pollutant emission. Direct combustion in coal stoves emits approximately 787.56 kg/yr of pollutant emissions including total suspended particles (TSP), SO x , CO, and particulate matter (PM) 2.5, while the emission of Ger 2 for the period of experimentation was 40.84 kg. The annual amount of Ger 1's emission and emission amount per kWh of the TPP were obtained from [20]. Figure 6 shows a daily operational profile of the ger with the ETS heater and solar PV for one week, which is 5th of February to 12th of February, as a simulation result. The validation results for the estimated PV output and heating demand are presented in Table A3 of Appendix C. The numerical values indicate the monthly average of the mean average percentage error (MAPE), mean absolute error (MAE), and root mean square error (RMSE). numerical values indicate the monthly average of the mean average percentage error (MAPE), mean absolute error (MAE), and root mean square error (RMSE).

Simulation Results
The amount of electricity from the PV system used in the ger is sufficiently high in the first three days because of the availability of solar irradiance, while on the 8th and 9th of February, PV selfconsumption represents a small fraction.  Figure 7 shows the profile of the energy performance of the Ger with the ETS heater and solar PV for the heating season. The energy performance is represented by monthly average daily energy consumption. For the period of the heating season, the fraction of electricity consumed in the ger from the solar PV is 31% of the total energy consumption. The remaining 69% is the amount of electricity imported from the main grid. The ger used 76% of the output of the PV system, while the PV system exported 24% of the total output of the PV to the main grid. The amount of solar PV exported to the main grid was a small fraction when the highest amount of heating energy was demanded in November, December, January, and February, and the self-consumption of PV increases because of the daytime charging of the ETS heater.
Furthermore, the potential for satisfying the demand for hot water through the surplus power of solar PV has been prospected when there is no heating demand in summer. Therefore, various types of thermal storage materials, for instance, water, ceramic, and phase-change material (PCM), need to be considered in the system configuration of this study. The optimal operation of the ETS heater powered by the solar PV and the main grid has the potential to reduce electricity costs for customers and increase renewable energy installation in the main grid. The amount of electricity from the PV system used in the ger is sufficiently high in the first three days because of the availability of solar irradiance, while on the 8th and 9th of February, PV self-consumption represents a small fraction. Figure 7 shows the profile of the energy performance of the Ger with the ETS heater and solar PV for the heating season. The energy performance is represented by monthly average daily energy consumption. For the period of the heating season, the fraction of electricity consumed in the ger from the solar PV is 31% of the total energy consumption. The remaining 69% is the amount of electricity imported from the main grid. The ger used 76% of the output of the PV system, while the PV system exported 24% of the total output of the PV to the main grid. The amount of solar PV exported to the main grid was a small fraction when the highest amount of heating energy was demanded in November, December, January, and February, and the self-consumption of PV increases because of the daytime charging of the ETS heater. The orange color is the amount of electricity consumed in the ger from main grid. The blue color is the amount of electricity consumed in the ger from the PV system. The yellow color is the electricity flow from the PV system to the main grid.

Conclusions
This study experimentally investigated the potential of a ger with an ETS heater and added insulation material in the cold climate of Ulaanbaatar city by comparison with a reference ger from the viewpoint of energy consumption. In addition to this, the technical feasibility of a ger with an ETS heater scheduled to charge during the daytime and solar PV was evaluated by using a simulation based on in-house software developed in MATLAB. Profile of the energy performance of the ger with the ETS heater and solar PV for the heating season. The energy performance is represented by the monthly average daily energy consumption. The orange color is the amount of electricity consumed in the ger from main grid. The blue color is the amount of electricity consumed in the ger from the PV system. The yellow color is the electricity flow from the PV system to the main grid.