1. Introduction
The global climate is warming, making the summer months begin to experience frequent extreme heat. On average, the global population experienced three times as many extreme heat days in 2020 than it did in 2008. Extreme heat has been relentless across the Northern Hemisphere this summer and threatens people’s lives. In addition, the current repeated prevention and control measures related to the novel coronavirus pneumonia epidemic tend to be normal. A large number of medical workers wear protective clothing to participate in the work. The high–temperature environment in clothing endangers the health of medical staff. At present, there have been studies on the detection and monitoring of human physiological signals by wearable electronic devices [
1,
2,
3,
4], but even more important is cooling the body when it becomes too hot. Therefore, workers operating in high–temperature environments are in urgent need of a portable cooling device combined with special clothing (such as welding suits, medical protective clothing, etc.). Taking epidemic prevention front–line medical staff as an example, there are traditional cooling methods such as turning on fans, turning on air conditioners in isolation rooms, and touching ice cubes. There are also new clothes with positive pressure protective features for cooling down. However, the above cooling methods have problems such as low refrigeration efficiency, high energy consumption, and poor portability. In addition, for the large temperature difference between day and night in the field, there are still few studies on portable refrigeration and heating two–way–control equipment. Semiconductor refrigeration produces heat while cooling. Taking advantage of this feature, the portable device can be used in hot conditions during the day and can be used at night for heating inside a field tent. Therefore, this paper uses semiconductor refrigeration technology to design and build a small and low–power consumption portable refrigeration device, and study the performance of the device based on CFD simulation.
Current research on portable cooling equipment is focused on the “microclimate zone” within clothing [
5], which is a thin layer of air between the interior of the clothing and the surface of the body skin. The main factors affecting physical comfort in this zone include the flow rate, temperature, and humidity of the fluid in the zone [
6]. The adjustability of the cooling performance, such as the adjustment of air velocity, flow rate, or temperature, should be considered in the device design. Gas and liquid heat transfer are the most common cooling methods, as shown in
Figure 1. The principle of gas heat transfer is to use a fan to blow the gas into the clothing and discharge it from the opening of the clothing so as to accelerate the evaporation of sweat on the surface of the body and achieve cooling effect, while the principle of liquid heat transfer is to use low–temperature liquid to take away the heat on the contact surface to achieve the cooling effect. They both convert the external electrical energy into the kinetic energy of the fluid medium [
7,
8]. The rapid flow of the fluid accelerates the heat transfer process, which is used to control the temperature distribution in the “microclimate zone”. The limitation of the above two methods is that they can only achieve separate refrigeration, and the size and weight of the liquid cooling equipment is greater.
In the study of gas cooling, Zhao et al. [
9] produced an air–cooled suit. A miniature fan and a battery were integrated into the garment, and the miniature fan enhanced forced convection to increase the evaporation rate of the sweat. It can accelerate the temperature reduction of the body surface and microclimate zone. Zhao et al. [
9] made five jackets with the same fan airflow size and only the installation position was located in different positions of the garment with closable openings in order to study the cooling effect of different fan positions, as shown in
Figure 2. The same battery and micro–fan are integrated in all jackets. Each ventilation unit was a small fan with a diameter of 10 cm. They can circulate ambient air around the trunk at a flow rate of approximately 0.012 m/s. The test results show that the fan at different locations has a more pronounced local impact, but only with the cooling function, and the cooling rate is slow and takes about 15 min to reach a steady state.
For the study of traditional water–cooled clothing, Guo et al. [
10] conducted relevant research on liquid cooling clothing. The theoretical model of heat transfer from the human skin surface to the environment shows that the maximum working duration of liquid–cooled clothing is 3.36 h at an ambient temperature of 45 °C, a flow rate of 224.5 mL/min, and a maximum cooling rate of 243.2 W/m
2, and the experimental results were within 10% error from the calculated results [
10]. Although Guo et al. have balanced the cooling effect and weight of the liquid cooling suit, its weight still reached 2 kg. At the same time, the cooling effect decreases with the increase of liquid temperature. Rahman et al. [
11] creatively used a non–toxic gallium–based liquid metal [
12] as the cooling medium for the liquid cooling suit system, which makes it possible to carry more heat in a smaller volume, as shown in
Figure 3. This active cooling system is one–third lighter than that of conventional liquid cooling systems, while the cooling effect is four times longer than that of conventional systems [
6]. However, it still requires liquid storage devices, and there is a risk of liquid leakage. In addition, Lu et al. [
13] and Kang et al. [
14] combined the function of ventilation and heat dissipation to develop a personal lightweight cooling system with better results, as shown in
Figure 3. The system uses phase change material for heat dissipation. Its power consumption is greater than that of air cooling and liquid cooling. Meanwhile, the continuity of phase change refrigeration is not good. It can only provide a cooling time of 90 to 110 min in hot and humid environments and 20 to 30 min in cold and dry environments. Due to the degradation of thermal physical properties during the use of phase change materials, its durability needs to be further studied.
In terms of heating in microclimates, flexible heaters based on stretchable polymers with conductive compositions have seen a number of advances. Current studies in this direction include a metal nano–network [
15], copper–plated fibers [
16], and graphite [
17]. However, these flexible heaters tend to increase local resistance during deformation, resulting in local high temperatures. The core component of the portable refrigeration device designed in this paper is a commercial semiconductor cooling plate, which can greatly reduce the cost of making the refrigeration device. The semiconductor refrigeration method involved can be classified as a heat pump in the usual sense, and it is also commonly referred to as a thermoelectric refrigeration chip, as shown in
Figure 4. The main principle is the use of some semiconductor materials that have the Peltier effect to achieve the purpose of cooling, and it is a reversible process. It can convert heat into electricity. Zhang et al. [
18] reviewed the research progress and proposed the research direction. For refrigeration–related applications, Alwan et al. [
19] installed a cooling chamber with thermoelectric elements on top of the solar still to enhance condensation. Miao et al. [
20] also designed a new type of plate–type ground source radiation system for space heating and cooling in buildings by utilizing the characteristics of semiconductor refrigeration technology. Compared to other refrigeration components, it has the advantage of high reliability and small size and does not require high space size. It needs timely and effective heat dissipation of its heating surface, and this heat can also be used to heat microclimates when heating functions are needed. The semiconductor refrigeration technology has been applied to target cooling devices, refrigeration storage cabinets, lithium battery temperature control devices, and other fields. Semiconductor refrigeration technology in which a semiconductor cooling plate directly wraps around the heat dissipation surface is used for thermal conductivity cooling. It can easily reach sub–zero temperatures in a short period, and the choice of the power of the semiconductor cooling plate makes the time to reach the maximum temperature difference different. The device behaves similarly to phase change refrigeration in terms of overall power consumption, but the lowest temperature it can reach with the same power is difficult to achieve with other refrigeration principles. It meets the conditions of use in the application of a portable semiconductor refrigeration device and facilitates miniaturized design.
The purpose of this paper is to develop a portable semiconductor refrigeration device. It can also realize the functions of refrigeration and heating at the same time and has a very wide application prospect. The comparative analysis of the five aspects of large temperature difference, high efficiency, low energy consumption, miniaturization, and low noise is shown in
Figure 5. Based on the above situation, semiconductor refrigeration technology can overcome the limitations of existing technologies applied to portable refrigeration devices after considering miniaturization and an excellent refrigeration effect under the same power consumption. In this paper, we use computer modeling and CFD simulation to assist in the design and build and test of the actual device based on the simulation results. Finally, we provide a small portable refrigeration device based on semiconductor refrigeration technology which is lightweight and can actively cool down the temperature, making it possible to regulate the temperature inside the garment more flexibly and actively.
In this paper, we use computer modeling and CFD simulation to assist in the design and build and test of the actual device based on the simulation results. Finally, we provide a small portable refrigeration device based on semiconductor refrigeration technology which is lightweight and can actively cool down the temperature, making it possible to regulate the temperature inside the garment more flexibly and actively. The workflow is shown in
Figure 6.
3. Performance Analysis Based on CFD Simulation
3.1. Assumptions of Physical Model
Simulating the performance of semiconductor refrigeration devices easily leads to a large amount of calculations and even non–convergence of results when more conditions are considered. Therefore, before establishing the control equations, the temperature field distribution of portable refrigeration devices based on semiconductor refrigeration technology needs to be decomposed into several basic models, including heat conduction, convective heat transfer, and radiation heat transfer. In order to simulate the solution process as close as possible to the actual situation while reducing the computational effort and facilitating convergence, the model with a small amount of influence is simplified or ignored. The model assumptions in CFD simulation in this paper are shown in
Figure 9.
3.2. Control Equations
The finite element analysis technique is an important means of modern engineering research [
25] and improving the quality of human life [
26]. Since the development of CFD technology, it has been widely used in the analysis of fluid flow, heat conduction, and other related physical phenomena. It is widely used in the comparative analysis of engineering field [
27]. Nguyen, T. [
28] selected the k–e turbulence model to study the effects of the straight fins and V–fins on the numerical simulation results of heat dissipation. The research showed that the effect of the fin shape on the heat transfer was not different when other conditions are equal. In terms of use for cooling, Jing–Ming Dong et al. [
29] also used CFD to study a miniature steam ejector refrigeration system, which is embedded with a capillary pump loop that can result in a compact design.
The Fluent software version used in this paper is Ansys 2021 R2, which is a software application from ANSYS, Inc in the United States. Due to the fast air velocity inside the device, the state of the gas at the outlet of the micro–fan is not smooth, and the flow is in a non–stationary state. After checking the relevant cases, the turbulence model in Fluent software and the standard
k−
ε model with good generality can be used to solve the simulation of the gas flow in the semiconductor cooling device in this design, where the transport equations for the standard
k−
ε turbulence model are as follows in Equations (8) and (9).
where
Gk is the turbulent energy term due to the laminar velocity gradient;
Gb is the turbulent energy term due to buoyancy;
YM is the contribution term to the dissipation rate due to the expansion of the turbulent pulsations in the compressible flow to the global flow;
C1ε,
C2ε,
C3ε are three constants;
σk is the
k–equation turbulence Prandtl number;
σε is the
ε–equation turbulence Prandtl number;
Sk is a user–defined turbulent energy term;
Sε is a user–defined turbulence dissipation source term;
ui is the velocity of
i–phase flow at moment
t;
xi is the spatial position of the flow mass point in phase
i at time
t;
xj is the spatial position of the flow mass point in phase
j at time
t;
μ is the hydrodynamic viscosity; and
μt is the turbulent viscosity coefficient.
In this simulation calculation,
C1ε is 1.44,
C2ε is 1.92,
σk is 1.0,
σε is 1.3, and the turbulent viscosity coefficient is calculated by Equation (10):
where
Cμ is the constant 0.09. The default values of the above parameters are those set in the FLUENT software.
3.3. CFD Models
The scientificity of model establishment determines the degree to which the simulation results reflect the real situation, and at the same time, the deviation of the simulation results can be controlled within a suitable range. Further, numerical calculations and performance simulations can be performed for the internal temperature field of portable refrigeration devices based on semiconductor refrigeration technology. The fluid and solid domains of the device are first modeled using Solid works software. Among them, the fluid domain includes cold air fluid and hot air fluid, and the solid domain includes a semiconductor cooling plate, copper fins, and aluminum fins and shell. Redundant parts that are not useful for the simulation process (circuit parts such as potentiometers) to reduce the meshing effort were removed. The total length of the model is 123.5 mm, the total width is 89 mm, the total height is 29.7 mm, and the volume is 125,887.34 mm3.
The above modeled solid domain and fluid domain were meshed. In this meshing, adaptive sizing is used with a cell size of 0.001 m, a span angle center set to fine, the default grid average surface area of 9.9993 × 10
−5 m
2, and the minimum edge length of 1.1514 × 10
−6 m. The fluid domain is divided into heat dissipating air fluid as well as cooling air fluid, and the mesh division results are shown in
Figure 10a (left is heat dissipating air fluid, right is cooling air fluid). The mesh division results of aluminum fins, copper fins and semiconductor cooling plate are shown in
Figure 10b, and the overall mesh situation and cross–section are shown in
Figure 10c,d. The number of grid nodes is 236,325, the number of cells is 639,385, and the quality of the grid is good to meet the simulation requirements.
In order to verify the independence of different mesh densities, three mesh sizes are selected: 0.0005 m, 0.001 m, and 0.002 m. Because the meshing settings of hot surface simulation and cold surface simulation are the same, one simulation can be selected for verification. The hot surface simulation is taken as the verification object. The verification results are shown in
Table 2. The results show that the variation of mesh density does not affect the accuracy of the final calculation results but also can predict the performance of the device. It can be concluded that the mesh density above is reasonable.
Each material was created in Fluent, including plastic shells, copper, and bismuth telluride (Bi2Te3), the substrate material for the semiconductor cooling plate, and material parameters were entered to match the material properties to each solid and fluid domain. In particular, this simulation involves energy exchange, so it is necessary to open the energy equation. The heat source is set for the semiconductor cooling plate, and the value of the energy source term is calculated based on the semiconductor power and refrigeration workload and the semiconductor refrigeration wafer size.
For the setting of the entrance boundary conditions, the velocity is set to a uniformly distributed constant. The wind speed in the x–axis direction is set to v = 7 m/s at the outlet of the micro–fan based on the calculation results. The y–axis and z–axis directions are ideal flow rates, both set to zero (where the x–axis direction is perpendicular to the end face of the micro–fan outlet and parallel to the cold/hot air flow path axis). The initial temperature is 20 °C and the exit boundary condition is kept as default for the free port.
3.4. Simulation Results
The results of the solution are shown in
Figure 11 and
Figure 12 of the following series, including the overall temperature, section temperature, and outlet flow rate. It can be seen that the temperature of the hot air surface can reach 54.7 °C near the semiconductor cooling plate and about 40 °C at the exit, and the wind speed can also reach 8.6 m/s. The temperature of the cold air surface can reach −9.3 °C near the semiconductor cooling plate and about 8 °C at the exit, and the wind speed can reach 8.3 m/s.
The simulation results of the semiconductor cooling device can meet the performance requirements of air temperature and air speed compared to the design performance requirements. To a certain extent, it can show that the device is feasible and provides some guidance for building the device entity and realizing the designed performance requirements.
5. Conclusions
At present, the principles of air cooling, liquid cooling, and phase–change cooling have been applied to portable refrigeration equipment. The research shows that they can achieve certain cooling effect. However, through comparative analysis, the above technology has problems such as poor refrigeration capacity, large volume, loud noise, and so on. The principle of semiconductor refrigeration can overcome the above shortcomings, so this paper proposes and designs a new idea of semiconductor refrigeration technology applied to portable refrigeration equipment. The design scheme is constructed by three–dimensional modeling software and imported into the CFD module of ANSYS software for simulation analysis. In the simulation, the working conditions of the portable refrigeration device are set up, the model is simplified, and the performance of the device is simulated by computer. The CFD simulation results show that the device has good performance. According to the design scheme, the actual equipment is then made. The device is encapsulated in a 3D–printed shell through circuit control to realize adjustable safety protection and other related functions. The overall power consumption of the device is 82 W. An infrared thermal imaging camera and anemometer are used to measure the wind temperature and speed of the equipment. The temperature of cold air surface and hot air surface reached a stable state of 13.9 °C and 47.2 °C, respectively, within about 4 min. The experimental results are basically consistent with the simulation results, which meet the design performance requirements and verify the feasibility of miniaturization design of the device. The limitation of the current study is that the unit is a stand–alone unit and does not integrate the unit with equipment requiring cooling or heating (e.g., clothing, tents, etc.). This is also one of our future research directions. The expected idea is to attach the device to the garment so that the device can be removed and replaced easily.