Computational Heat and Mass Transfer and Fluid Flow through Complex Passages for Emerging Applications
View this Special IssueResearch Article  Open Access
Modeling Validation and Control Analysis for Controlled Temperature and Humidity of Air Conditioning System
Abstract
This study constructs an energy based model of thermal system for controlled temperature and humidity air conditioning system, and introduces the influence of the mass flow rate, heater and humidifier for proposed control criteria to achieve the controlled temperature and humidity of air conditioning system. Then, the reliability of proposed thermal system model is established by both MATLAB dynamic simulation and the literature validation. Finally, the PID control strategy is applied for controlling the air mass flow rate, humidifying capacity, and heating, capacity. The simulation results show that the temperature and humidity are stable at 541 sec, the disturbance of temperature is only 0.14°C, 0006 in steadystate error of humidity ratio, and the error rate is only 7.5%. The results prove that the proposed system is an effective controlled temperature and humidity of an air conditioning system.
1. Introduction
The air conditioning system can regulate indoor air quality to provide a more comfortable and safer environment. It is also important in the hightech industry and medical industry. Liu et al. [1] use computational fluid dynamics (CFD) to create a thermally comfortable and healthy environment by air distributed information from numerical simulation. Safaei and Goshayshi [2] finite volume method, which is based on the difference of temperature and momentum, has effect on mechanical ventilation or pressure different and finally results in indoor and outdoor flow. Safaei et al. [3] investigate the energy transferring through two different heated side walls, especially transferring in rooms and buildings. However, as the controller of airconditioning system has long time constant, multivariable, and high control complexity by interaction between parameters, it is thus difficult to obtain accurate mathematical model. Due to the high energy consuming, the research on high efficiency energysaving air conditioning system has become a topical subject. An effective control law is very important for developing high energy efficiency airconditioning system. The discussion on the influence of airconditioned zone energy input and output on indoor temperature and humidity, as well as the analysis of the dynamic characteristics of air conditioning system, is important in analyzing air conditioning system and designing control law. Therefore, it is necessary to develop an approximate thermodynamic model and validate the feasibility of model at the design stage of air conditioning system. Many recent studies have used the theoretical or experimental method to explore the dynamic model of air conditioning system as the basis of air conditioning system control. Tsao et al. [4] developed high efficiency makeup air units (MAUs) for technology factories in subtropical zone by modeling thermodynamic system and proposed eight schemes for analysis, among which was the dry coil cooling (DCC) based energysaving scheme. However, the thermal model was built for airconditioning equipment only, and the effect of energy on the state of airconditioning space was not discussed. Soyguder et al. [5] built a thermal model of air conditioning system and used PID expert fuzzy control to adjust the airconditioned zone temperature change. However, that system lacked humidity control for airconditioned zone. Tashtoush et al. [6] developed a thermal model of air conditioning system and airconditioned zone and used PID to control the outflow air temperature and humidity for constant temperature and humidity control. However, only the outflow air temperature and humidity were controlled, while the operation of the equipment was not discussed. The outflow air temperature and humidity matching problem of the chilled water was not considered. Rehrl and Horn [7] proposed a simple thermodynamic model of airconditioning equipment and combined it with feedback method for predictive control of air conditioning system. However, in terms of humidifier, the control mode was designed according to the characteristic curve of valve and humidifying capacity, and the test room was built instead of thermal model of airconditioned zone for control. Platt et al. [8] built a model of air conditioning system and used genetic algorithm to optimize air conditioning control. However, only the temperature change of airconditioned zone was controlled and the thermal model of airconditioning equipment was not built, and the genetic algorithm was used as control base. Yao et al. [9] predicted the variance in chilled water by building the thermal model of chilled water coil and validated it by simulation. However, they did not apply the control law to the airconditioned zone. Soyguder and Alli [10] used the air conditioning system thermal model and applied inputoutput data to PID ANFIS to control the fan speed. Barbosa and Mendes [11] built the thermal model of various equipment types of central air conditioning system and used weather data to simulate the onoff control of chiller. Kasahara et al. [12] modeled the airconditioned zone and controlled temperature and humidity air conditioning system and used the thermal model and PI control to change the air output for constant temperature and humidity. Riederer et al. [13] built the thermal model of room system but did not regard the room as uniformly mixed constant volume. The air temperature change was deduced from the positions of sensors, and it was validated by experimental simulation. Kang et al. [14] proposed the HVAC model which introduces CO_{2} concentration to consider the comfort of room, and linearquadratic regulator (LQR) is applied for optimizing and stabilizing the system. AmosAbanyie et al. [15] investigated the thermal mass effect, window size, and night ventilation on peak indoor air temperature (PIAT) via E plus (E+) simulation. Homod [16] analyzed the advantage and disadvantage from types of HVAC model.
Therefore, in order to implement constant temperature and humidity control, the controlled temperature and humidity air conditioning system was used in this study. The air feeder, cooling coil heater, and humidifier constituted the airconditioned zone at constant temperature and humidity. The thermal model of controlled temperature and humidity air conditioning system and airconditioning equipment was built to discuss the influence of various types of equipment in the air conditioning cabinet on the temperature and humidity of airconditioned zone. Finally, considering the matched cases, the PID control law was used to control the air mass flow rate (), heater input percentage (), and humidifier input percentage (). This study investigates the controlling based on thermal and humidity conditions only without considering the energy transferring such as CFD, thus to discuss the stabilization time, overshoot and stability of controlled temperature, and humidity air conditioning system under PID control.
2. Modelling
2.1. Description
The main function of controlled temperature and humidity air conditioning system is to make the airconditioned zone at constant temperature and humidity by applying different control strategies, so as to meet different industrial requirements, such as chemical laboratory, measuring and testing laboratory, process of hightech industry, and electronic machine room. Figure 1 shows a typical controlled temperature and humidity air conditioning system.
First the difference between the practical status and target value is obtained by detection of indoor air state. Then the flow of return air (RA) is controlled by exhaust air (EA) and the outdoor air (OA) is introduced to generate mixing air (MA) in airhandling unit. The mixing air is processed by cooling and dehumidifying by cooling coil (C.C.), heating by heating coil (H.C.), and humidifying by humidifier (H.D.) to mix the desired state for airconditioned zone gradually.
As shown in Figure 2, the algorithm first compares the difference between given , and room , , respectively, and delivery to PID controller. The control parameters , , and of C.C, H.C., and H.D. are changed and begin to affect the room temperature and humidity. The controlled , are compared again with room , and begin the next loop.
2.2. Mathematical Model of Controlled Temperature and Humidity of Air Conditioning System
Figure 1 shows that, under the control of chilled water coil, heater, and humidifier of the controlled temperature and humidity air conditioning system, the airconditioned zone is provided with appropriate energy to keep its temperature and humidity corresponding to the set conditions. The modeling of the system of the exhaust air, mixing air, chilled water coil, heater, and humidifier components is discussed. The design of comprehensive mathematical model of a controlled temperature and humidity air conditioning system is unpractical and unreasonable due to the complexity from its multivariable. Therefore the mathematical model will be simplified by the following assumptions.(1)The specific heat and density of air are constants.(2)There is no air leakage in the process.(3)The boundary is insulated.(4)The air is regarded as the ideal gas.
2.2.1. Exhaust Air and Mixing Air Model
Generally, introduction of the outside air conditioning will raise cooling load on the airconditioning zone. However, an appropriate amount of outside air can dilute the hazardous material in the airconditioned zone, thus providing necessary fresh air for breathing, removing pollutant and controlling the temperature and humidity of work site, and upgrading the indoor air quality. The exhaust air damper is linked to outside air damper and controls the air discharge and outside air volume. Figure 3 shows the mixing process of indoor return air and outside air in the air mixing zone of the air conditioning cabinet, the outside air and return air volumes are adjusted by controlling the outside air damper opening, and the air is mixed in the mixing zone of air conditioning cabinet.
If the mixing process is a steady flow process, the energy conservation for the mixing process is expressed as and the balance of the mass of water is expressed as
2.2.2. Cooling Coil Model
The energy in the air is transferred into the chilled water cycle, and the chilled water cycle delivers the energy to the evaporator. The evaporation of refrigerant in the evaporator takes the energy away, reducing the zone temperature and humidity. According to the air conditioning system, the chilled water temperature is generally 7~12°C. When the return air passes through the chilled water coil, if the chilled water temperature is lower than the return air dewpoint temperature, the vapor in humid air is condensed for dehumidification. The mixed air enters the chilled water coil, taking the sensible heat () and latent heat () of air away. The cooled and dehumidified air is fed to the heater. If the heat exchange rate is 100% when the air passes through the chilled water coil, it is steady flow process, and the moisture condensed on the coil is removed completely and immediately by gravity as shown in Figure 4.
Disregarding the variance in kinetic energy and potential energy, the energy conservation of chilled water is expressed as The humid air energy variation of (3) is divided into the energy of dry air and water vapor, expressed as The sensible heat of air () only influences the air temperature and the latent heat of air () only influences the air vapor, so (4) can be expressed as (5) and (6), respectively, as follows: Equation (6) can be reduced to vapor mass conservation equation, expressed as
2.2.3. Heating Coil Model
It is switched on when the room temperature is lower than the set point. By heating the temperature of supplied air with heater, the heat energy is transferred to the airconditioned zone, and let the temperature reach the set point. The heater is divided into hotwater heating, electric heater, and steam heating. In this study, the hot water coil is used as the heater of controlled temperature and humidity air conditioning system. The air through chilled water passes through the heating coil to take the heat into the air, and the heated air is fed to the humidifier as shown in Figure 5. The heating capacity of heater is pure sensible heat change, the heater output power () heats the air, and the humidity ratio of airconditioned zone will not be influenced.
If the heating process is steady flow process, the energy conservation relationship of heater as shown in Figure 5 is expressed as
2.2.4. Humidifier Model
The humidifier output power () converts electric energy into latent heat, the moisture is added in the system by steam humidification so that the humidity of the air supply changes, and the system processed air is fed into the airconditioned zone. The humidifier humidification is pure latent heat modification; therefore, only humidity changes, and the temperature of airconditioned zone is not influenced. If the humidifier is steady flow process, it can be described by a simple vapor mass conservation equation as
2.2.5. Integration to a Zone Hygrothermal Model
In order to make the airconditioned zone reach the set temperature and humidity, the controlled temperature and humidity air conditioning system must be used. The controlled temperature and humidity air conditioning system consists of chilled water coil, heater, and humidifier. The mathematical model of airconditioned zone is simplified by the following assumptions.(1)The airconditioned zone is an open and unsteady flow system in constant volume.(2)The air temperature and humidity in the airconditioned zone are distributed uniformly.(3)The specific heat and density of air are constants.(4)There is no air leakage in the process.(5)The wall surface of airconditioned zone is insulated.(6)The air is regarded as ideal gas.
According to thermodynamic energy conservation equation, the thermal model [17] can be expressed as If the kinetic energy and potential energy variations in the process are ignored and the system has not applied work, (10) is reduced to
Figure 6 shows that the airconditioned zone state is significantly influenced by the persons, office equipment in operation, and various machine tools in the space and the heat transfer to the space is shown in Figure 2 so that the air state in the space changes. The temperature and humidity increase continuously, and the state keeps deviating from comfort zone. Different air conditioning loads influencing indoor state can be reduced to sensible heat load () and latent heat load (). The energy conservation equation (12) for airconditioned zone can be deduced from (11): According to the law of conservation of mass of vapor, the airconditioned zone humidity change can be expressed as where and are the energy and vapor mass, respectively, taken in or out of the system by various equipment types in the air conditioning cabinet, so (2) are combined with the thermal model of various equipment types of the air conditioning system in Section 2.2. The energy conservation and vapor mass conservation equations represent the indoor temperature and humidity variations, respectively, expressed as
3. Results and Discussion
In order to validate the controlled temperature and humidity air conditioning system model deduced in Section 2, MATLAB is used to create the dynamic simulation program of the system, and set the simulation condition working procedure of [6]. The changes in the temperature and humidity ratios of airconditioned zone are compared with the literature, so as to validate the accuracy of the thermal model and the dynamic simulation program. When the dynamic simulation program is determined, a wellknown control law is still required to direct the action of controlled temperature and humidity air conditioning system, so that the system function can be perfected. Therefore, the PID control is used and the cooling and dehumidification capability matching of chilled water is considered. ZieglerNichols rule is applied to obtain PID control parameters, so as to analyze the effectiveness of PID control of controlled temperature and humidity air conditioning system.
3.1. Modeling Validation
The dynamic simulation of the proposed model, based on the hypothesis in chapter 2, is built and processed by Matrix Laboratory (MATLAB) and is shown in Figure 7. The program simulates eight input parameters, which are zone sensible heat (), zone latent heat (), initial temperature of airconditioned zone (), zone initial humidity ratio (), air mass flow rate (), heater capacity (), outside air (), and humidifier capacity (), respectively. When the above parameters are imported into the simulation, the program calculates the energy variation in the space resulting from various equipment types.
This study developed an open loop control test and set different air mass flow rates. The effectiveness of dynamic model of the system is validated by comparing the results of this study with the simulation results of the literature. This test neglected the external wall and heat transfer. The set conditions are shown in Table 1.

The simulation results are shown in Figures 8 and 9. If there is only zone sensible heat, the stable temperature decreases as the air mass flow rate increases, and the humidity ratio converges to 0.00763 . According to (14), when there is sensible heat load, the difference between zone temperature and chilled water temperature decreases with the indoor temperature. The temperature is stabilized until the heat taken away by the chilled water equals the sensible heat load. Therefore, with a larger air mass flow rate, the stable temperature is lower. When the system lacks air mass flow rate, the sensible heat load cannot be taken away adequately, so that the cooling is failed, and the return air temperature difference increases continuously. However, as the temperature difference changes, the heat load taken away by the chilled water increases and it reaches balance with sensible heat load at last. In terms of humidity ratio without latent heat, according to (5), regardless of , the outflow air humidity ratio is eventually converged. Relevant parameters are shown in Table 2. As seen, when is large, the stabilization time is short. The simulation results are consistent with the trend in the literature. This study did not consider the wall heat transfer, and the simulation results were compared with the literature. At 4000 sec when stabilization was most approximately reached, the maximum error was 1°C, and the error rate was 4%. The humidity ratio converged to 0.763 with maximum error was 2.4%.
3.2. Simulation Results
PID control is a closed loop control system, commonly used in industry. The PID algorithm can provide a considerable control effect with good adjustment. Sensors are applied for detecting the present state of airconditioned zone and returning them back to comparator. The state of airconditioned zone reaches the set state through proportional control, integral control, and differential control. As it is simple and understandable and its adjustment does not need precise system model, it is widely used in various domains. The PID control inputoutput relationship is expressed as the following equation: where , , and are the proportionality, integral, and differential coefficients of PID controller. The proportional control regulates and controls the system according to the difference between input and output signals. However, the steadystate error cannot be removed by only proportional control. In order to eliminate steadystate error, the integral control is adopted. The integral term takes the integral of time and the integral term increases with time, so that the control output increases until the steadystate error is eliminated. In the process of regulation, the system may have oscillation even run out of control due to inertial or lag component, and the integral term magnifies the oscillation on the contrary. Therefore, the differential term is used to predict the variation trend to improve the interference effect. However, the three control parameters influence each other. Hence, how to obtain control parameters becomes an important but difficult problem in PID control.
The PID controller setup parameters in this paper are designed according to ZieglerNichols empirical rule, which is regarded as a good parameter optimization in PID control [6]. The parameters are adjusted continuously to make system oscillation, so as to know the system response. The control parameter settings are derived from empirical equation, as shown in Table 3. In this paper, considering the matching of cooling and dehumidification effects of chilled water coil, the air mass flow rate () is used as control parameter, and the heater and humidifier control is added. With the assistance of heater output power percentage () and humidifier output power percentage (), the loss of chilled water due to excessive cooling or dehumidification is compensated.

According to the system parameters of desired controlled temperature and humidity of air conditioning system, a PID controlled real time of the proposed air conditioning system can be built by Simulink as shown in Figure 10. Simulation begins with given initial conditions, and the controller adjusts P, I, and D, respectively, to control room temperature and humidity.
As shown in Table 3, the chilled water is equipped with two PID controllers, which feedback indoor temperature and humidity ratio data. Two air mass flow rate data are exported under control. As the dehumidification and cooling only depend on the chilled water coil, considering the system matching problem, this paper uses large value as the actual output. Therefore, if sensible heat or latent heat is removed excessively, the heater and humidifier are used for compensation.
This paper observes the convergence and stability of PID control based on constant temperature and humidity system thermal model in simulated summer environment. Moreover, it validates the feasibility of PID control in controlled temperature and humidity air conditioning system. The set conditions are described below.(1)°C; .(2)°C; .(3)Zone volume is 36 m^{3}.(4)Sensible heat load is 1.3 kW.(5)Air mass flow rate is 0~0.3 kg/s.(6)Heater is 0~0.6 kW.(7)Humidifier is 0~1 kW.
At the beginning, the overheat room temperature is detected by PID controller of chilled water coil, triggering the increase of mass flow rate to cool down the room temperature till reaching 17.81°C. The overcool temperature stimulates the PID controller of heater and the heater is working till setting temperature. The crossfunction of these two PID controllers results in a sharp control curve as shown in Figure 11.
In humidifier, due to the fact the humidity ratio controlling needs only adjusting the chilled water coil, only one PID controller, which belong to chilled water coil, is working and results in a smooth curve as shown in Figure 12.
According to the indoor temperature and humidity ratio variations in Figures 6 and 7, under the PID control, the humidity ratio can be controlled well. However, in terms of temperature stability, the maximum disturbance is 0.14°C. Because there are two PID controllers, when the temperature and humidity approach the set values, and the output air between the two controllers, resulting in oscillation. In addition, as the airconditioned zone cooling response speed is higher than dehumidification, when the temperature is reached, the chilled water still needs to continue dehumidification, so that the indoor temperature decreases to 17.81°C before it rises to the set value. The maximum overshoot of temperature is 4.19°C. The temperature and humidity are stable at 541 sec, but the steady state of humidity ratio is 0.0074 , the steadystate error is 0.0006 , and the error rate is 7.5%.
The simulation data are displayed in the psychrometric chart, as shown in Figure 13. The air state changes in the airconditioned zone in the PID control mode. First, it is high temperature and high humidity state and only the chilled water is used for cooling and dehumidification. Therefore, the air state almost develops towards the outflow air temperature and humidity point straightly. Soon after the temperature is lower than the set point, the heater is started up to lower the temperature to the set point. The matching of various equipment types stabilizes the air state inside the airconditioned zone at the set point.
4. Conclusion
This study developed a thermal model for controlled temperature and humidity air conditioning system, including airconditioned zone, air mixing stage, chilled water coil, heater, and humidifier. Based on the thermal model, MATLAB was used to simulate open loop control system, and it is compared with references to validate the feasibility of the system. As the wall heat transfer was disregarded, the maximum error between room temperature simulation and literature was 1°C, and the maximum error rate between humidity ratio simulation and the literature was only 2.4%. Finally, the PID controller was used for control. When the system matching was considered, the ZieglerNichols empirical rule was applied to control the air mass flow rate, humidifier output power, and heater output power. The results showed that the PID control could reach stabilization at 541 sec at high temperature and high humidity. The maximum overshoot was 4.19°C, the humidity ratio error rate was only 7.5%, and the maximum temperature disturbance was 0.14°C. The goal for constant temperature and humidity could be attained effectively.
Nomenclature
:  Mass flow rate of the air stream (kg/s) 
:  Specific heat (kJ/kg°C) 
:  Temperature (°C) 
:  Humidity ratio 
:  Enthalpy (kJ/kg) 
:  Enthalpy of air through chilled water (kJ/kg) 
:  Latent heat of vaporization of vapor (kJ/kg) 
:  Maximum input power of heater (kJ) 
:  Maximum output power of humidifier (kJ). 
:  Air 
da:  Dry air 
:  Air through chilled water 
:  Heater 
:  Latent heat 
:  Sensible heat 
:  Water 
MA:  Mixing air 
OA:  Outside air 
RA:  Return air 
WO:  Return water 
Wi:  Supply water. 
:  Heater input power percentage 
:  Outside air damper opening 
:  Humidifier output power percentage. 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by the Department of Refrigeration, Air Conditioning and Energy Engineering, National ChinTi University of Technology, as well as National Science Foundation of Taiwan and under Grant no. NSC1002622E167016CC3.
References
 W. Liu, S. Mazumdar, Z. Zhang et al., “Stateoftheart methods for studying air distributions in commercial airliner cabins,” Building and Environment, vol. 47, no. 1, pp. 5–12, 2012. View at: Publisher Site  Google Scholar
 M. R. Safaei and H. R. Goshayshi, “Investigation of turbulence mixed convection in air filled enclosures,” Journal of Chemical Engineering and Materials Science, vol. 2, no. 6, pp. 87–95, 2011. View at: Google Scholar
 M. R. Safaei, M. Goodarzi, and M. Mohammadi, “Numerical modeling of turbulence mixed convection heat transfer in air filled enclosures by finite volume method,” International Journal of Multiphysics, vol. 5, no. 4, pp. 307–324, 2011. View at: Publisher Site  Google Scholar
 J. Tsao, S. Hu, T. Xu, and D. Y. L. Chan, “Capturing energysaving opportunities in makeup air systems for cleanrooms of hightechnology fabrication plant in subtropical climate,” Energy and Buildings, vol. 42, no. 11, pp. 2005–2013, 2010. View at: Publisher Site  Google Scholar
 S. Soyguder, M. Karakose, and H. Alli, “Design and simulation of selftuning PIDtype fuzzy adaptive control for an expert HVAC system,” Expert Systems with Applications, vol. 36, no. 3, pp. 4566–4573, 2009. View at: Publisher Site  Google Scholar
 B. Tashtoush, M. Molhim, and M. AlRousan, “Dynamic model of an HVAC system for control analysis,” Energy, vol. 30, no. 10, pp. 1729–1745, 2005. View at: Publisher Site  Google Scholar
 J. Rehrl and M. Horn, “Temperature control for HVAC systems based on exact linearization and model predictive control,” in Proceedings of the 20th IEEE International Conference on Control Applications (CCA '11), pp. 1119–1124, September 2011. View at: Publisher Site  Google Scholar
 G. Platt, J. Li, R. Li, G. Poulton, G. James, and J. Wall, “Adaptive HVAC zone modeling for sustainable buildings,” Energy and Buildings, vol. 42, no. 4, pp. 412–421, 2010. View at: Publisher Site  Google Scholar
 Y. Yao, Z. Lian, and Z. Hou, “Thermal analysis of cooling coils based on a dynamic model,” Applied Thermal Engineering, vol. 24, no. 7, pp. 1037–1050, 2004. View at: Publisher Site  Google Scholar
 S. Soyguder and H. Alli, “Predicting of fan speed for energy saving in HVAC system based on adaptive network based fuzzy inference system,” Expert Systems with Applications, vol. 36, no. 4, pp. 8631–8638, 2009. View at: Publisher Site  Google Scholar
 R. M. Barbosa and N. Mendes, “Combined simulation of central HVAC systems with a wholebuilding hygrothermal model,” Energy and Buildings, vol. 40, no. 3, pp. 276–288, 2008. View at: Publisher Site  Google Scholar
 M. Kasahara, Y. Kuzuu, T. Matsuba, Y. Hashimoto, and K. S. Kamimura, “Physical model of an air conditioned space for control analysis,” ASHRAE Transactions, vol. 106, part 2, pp. 307–317, 2000. View at: Google Scholar
 P. Riederer, D. Marchio, J. C. Visier, A. Husaunndee, and R. Lahrech, “Room thermal modelling adapted to the test of HVAC control systems,” Building and Environment, vol. 37, no. 89, pp. 777–790, 2002. View at: Publisher Site  Google Scholar
 C. S. Kang, J. I. Park, M. Park, and J. Baek, “Novel modeling and control strategies for a HVAC system including carbon dioxide control,” Energies, vol. 7, pp. 3599–3617, 2014. View at: Google Scholar
 S. AmosAbanyie, F. O. Akuffo, and V. KutinSanwu, “Effects of thermal mass, window size, and nighttime ventilation on peak indoor air temperature in the warmhumid climate of ghana,” The Scientific World Journal, vol. 2013, Article ID 621095, 9 pages, 2013. View at: Publisher Site  Google Scholar
 R. Z. Homod, “Review on the HVAC system modeling types and the shortcomings of their application,” Journal of Energy, vol. 2013, Article ID 768632, 10 pages, 2013. View at: Publisher Site  Google Scholar
 Y. A. Cengel and M. A. Boles, Thermodynamics an Engineering Approach, McGraw Hill, New York, NY, USA, 6th edition, 2006.
Copyright
Copyright © 2014 JingNang Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.