key: cord-1035090-irhbuuq5 authors: Sun, Jian; Zhang, Mingkan; Gehl, Anthony; Fricke, Brian; Nawaz, Kashif; Gluesenkamp, Kyle; Shen, Bo; Munk, Jeff; Hagerman, Joe; Lapsa, Melissa title: COVID 19 vaccine distribution solution to the last mile challenge: experimental and simulation studies of ultra-low temperature refrigeration system date: 2021-11-08 journal: Int J Refrig DOI: 10.1016/j.ijrefrig.2021.11.005 sha: 8bef9014bfea73cac95d36ea81150e740ea8374c doc_id: 1035090 cord_uid: irhbuuq5 Most COVID-19 vaccines require ambient temperature control for transportation and storage. Both Pfizer and Moderna vaccines are based on mRNA and lipid nanoparticles requiring low temperature storage. The Pfizer vaccine requires ultra-low temperature storage (between −80°C and −60°C), while the Moderna vaccine requires −30°C storage. Pfizer has designed a reusable package for transportation and storage that can keep the vaccine at the target temperature for 10 days. However, the last stage of distribution is quite challenging, especially for rural or suburban areas, where local towns, pharmacy chains and hospitals may not have the infrastructure required to store the vaccine. Also, the need for a large amount of ultra-low temperature refrigeration equipment in a short time period creates tremendous pressure on the equipment suppliers. In addition, there is limited data available to address ancillary challenges of the distribution framework for both transportation and storage stages. As such, there is a need for a quick, effective, secure, and safe solution to mitigate the challenges faced by vaccine distribution logistics. The study proposes an effective, secure, and safe ultra-low temperature refrigeration solution to resolve the vaccine distribution last mile challenge. The approach is to utilize commercially available products, such as refrigeration container units, and retrofit them to meet the vaccine storage temperature requirement. Both experimental and simulation studies are conducted to evaluate the technical merits of this solution with the ability to control temperature at −30°C or −70°C as part of the last mile supply chain for vaccine candidates. Along with development and manufacturing, vaccine distribution also presents great challenges for supply chain since most COVID-19 vaccines require ambient temperature control for transportation and storage. Two types of vaccine have been developed by manufacturers (Pfizer and Moderna). Both vaccines are based on mRNA and lipid nanoparticles requiring low temperature storage. The Pfizer vaccine requires ultra-low temperature storage (−80 °C to −60 °C), while the Moderna vaccine requires −30 °C storage. Pfizer has designed a reusable package for transportation and storage that can keep the vaccine at the target temperature for 10 days and can accommodate between 1,000 and 5,000 doses. However, the last stage of distribution is quite challenging (Owens 2020) , especially for rural or suburban areas, where local towns, pharmacy chains and hospitals may not have the infrastructure required to store the vaccine at the required temperature ( Heilweil 2020) . Also, the need for a large amount of ultralow temperature refrigeration equipment in a short time period creates tremendous pressure on equipment suppliers (Sari, 2018) . In addition, there is limited data available to address ancillary challenges of the distribution framework for both transportation and storage stages, including safety concerns due to human exposure to large amounts of CO 2 from dry-ice sublimation, issues due to the pressure increase caused by dry-ice sublimation, and the potential damage caused by ultra-low temperatures. As such, there is a need for a quick, effective, secure and safe product/solution to mitigate the challenges faced vaccine distribution logistics. To address this gap, ORNL is collaborating with industry partner to develop a local vaccine storage solution using refrigerated container technology to ensure a significant increase in dry ice life to provide reliable temperature control for vaccine storage. The overall objective of this study is to assess the technical merits of utilizing container refrigeration units with the ability to control container temperature at -30 °C as part of the last mile supply chain for vaccine candidates. The study will be applicable to over the road trailer refrigeration units as well. As shown in Figure 1 , the test platform was set up at Oak Ridge TN. It is composed of a refrigeration storage container unit with size of 12.19 m x 2.44 m x 2.62 m (L x W x H), a refrigeration system, and vaccine packages. The refrigeration storage container unit is fully instrumented with,  Thermocouples, locations of thermocouples depend on different tests  CO 2 & O 2 sensors inside container box, both near front at reefer. The O 2 and CO 2 sensors were used in testing with a 5-psi diaphragm compressor to pump air from inside the container to outside. This will be done by drilling a hole and using clear tubing to move air to the sensors.  Scales to measure dry ice weights during sublimation  Omega 5V differential pressure monitor to measure the pressure difference from inside and outside the container. Figure 1 : Test platforms at Oak Ridge Tennessee The test unit ( Figure 1 ) is a refrigeration unit of lightweight aluminum frame construction, designed to be bolted on to the front of a container and serve as the container's front wall. The refrigeration system is self-contained all electric unit which consists of compressor, evaporator, condenser, fans and controls to provide precise temperature control from -35 °C to 30 °C. The unit operates on nominal 380/460 volt, 3-phase, 50/60 hertz (Hz) power. Control system power is provided by a transformer which steps the supply power down to 12 and 24 volts, single phase. The controller is a microprocessor operating automatically to select cooling, holding or heating as required to maintain the desired set point temperature within very close limits. The unit may also be equipped with an electronic temperature recorder. The controller has a keypad and display for viewing or changing operating parameters. The display is also equipped with lights to indicate various modes of operation. As shown in Figure 2 , two types of packages are used in this study: generic insulated box and pharmaceutical box. The details of these two packages are listed in table1. Polystyrene foam Polystyrene foam The tested refrigeration container units and vaccine packages were fully instrumented to monitor the operational performance. The measurement includes temperature distribution inside refrigeration container unit, temperature inside vaccine packages, supply and return air temperature of refrigeration system, pressure difference between inside and outside of refrigeration container unit, CO 2 and O 2 concentration level inside refrigeration container unit. The detailed list of the measurement points and specification of the instrumentation are given in table 2. The placement of the boxes is shown in Figure 3 . Twenty of the styrofoam boxes (B1 to B20) were placed inside the refrigeration container, and one (B0) was placed inside lab at air conditioning space for comparison purposes. When the dry ice was delivered, 21 boxes (B0 to B20) were filled with dry ice. Two of the boxes (B20 and B9) on the inside of the container were placed on a scale. For this test, thermocouples were placed strategically throughout the refrigeration storage container to monitor the temperature distribution inside. 50 thermocouples were placed of the container's vertical surface A, B, C, D, and E. A scan rate of 1 minute was used to collect data. The placement of the sensors is plotted in Figure 3 . A local dry ice manufacturer was contacted to deliver dry ice to fill each of packages as shown in Figure 4 . An empty small payload box was set inside each storage package with dry ice surrounding it. A thermocouple was placed in the center of each payload box to represent the vaccine temperature. Four thermocouples were placed inside two storage packages ( B20 and B9) at different heights to measure the temperature change of storage package along with the dry ice sublimation process. The general test procedure included the following steps:  Plug refrigeration unit into 460 V 60 Hz power source, turn the unit on.  Select a unit set point and defrost set.  Receive dry ice delivery and load each styrofoam container with dry ice around 22.5kg.  Place each box in the proper location as shown in Figure 4 with a thermocouple "buried" within the dry ice. Using shipping tape, seal off each outer carboard box. After boxes are placed appropriately, close the container rear doors.  Turn on 5 V power supply connected to the diaphragm pump. Start the data acquisition and ensure a scan rate of 1 minute has been selected. Check that the O 2 and CO 2 readings are ~20.8% and 0%, respectively. Ensure that the fresh air ventilation (FAV) is 100% closed and record the starting weight of dry ice inside the container and ambient. Lab tests were conducted to quantitively evaluate the performance of the proposed refrigeration storage container method, identify its potential risks, and improve and optimize the strategy. Four aspects of this proposed solution were studied in detail: temperature distribution of refrigerated storage container, temperature distribution inside the vaccine dry-ice package, dry-ice sublimation, and CO 2 concentration inside the refrigerated storage container Compared with regular storage, the refrigerated storage container can create −30 °C surrounding temperatures of the vaccine packages and slow the dry-ice sublimation process to extend the useful lifetime of maintaining the vaccines in dry-ice packages. Maintaining the temperature surrounding vaccine dry-ice packages within the required temperature tolerance is essential for this refrigeration storage container solution, which makes thermal design optimization of a cargo layout for vaccine cartons a necessity. Before optimizing the cargo layout, lab tests were performed to better understand the temperature distribution of the refrigerated storage container. Figure 5 shows the temperature profiles of four temperature measurements (C2, C5, C6, C9 in Figure 3 ). Figure 6 plots the temperature profile of one defrost period. In this test case, a typical defrost strategy of container refrigeration system: time-based strategy (once every 24 hours), was applied. The defrost was conducted through shutting down the refrigeration system and turning on the electrical heater installed under the evaporator coil. The test results indicate the following:  In general, the refrigerated storage container maintained the set point temperature during testing except for the defrost period.  The container temperatures increased to the range of 15 °C to 25 °C in most of the defrost period. During door open events, the temperature rose above 35 °C.  Because of defrost operation, the refrigerated container air temperature control excursion lasted approximately 2 h in this test which depends on the cargo load and container insulation. During this defrost period, the air temperature rose. The defrost operation ( e.g., termination temperature limits) can be adjusted to minimize the container air temperature rising. Because the air was supplied by the refrigeration system from the floor and returned to the top, the temperatures close to the unit top (C2), as shown in Figure 6 , were around 0.8 °C to 2.2 °C higher that the temperatures close to the unit floor (C9). During the defrost period, this temperature difference was around 19 °C. Other test data also indicate the air temperature inside refrigerated container unit increased with the height. The temperature distribution of the refrigerated storage container based on the test data presents several suggestions and concerns as follows:  Placing the vaccine boxes below a certain height is recommended to take advantage of the relative lower temperature near the supply air.  A height limit can be found through simulation based on the temperature gradient.  As noted earlier, the air temperature will rise during the defrost mode of the refrigeration unit, and this increased the surrounding air temperature inside the vaccine box. The potential risk of the defrost operation on the dry-ice package temperature and vaccine temperature will be discussed in the next section. As described in previous section, a small payload box was placed inside each storage package with dry ice surround the payload box ( Figure 4) . A thermocouple was placed in the center of each payload box to represent the vaccine temperature. Figure 7 shows the temperature profiles of air inside the payload box in some dry-ice packages. The initial mass of dry ice in each storage packages are listed in table 3. Because the temperatures of each payload box follow similar temperature profiles, a generic temperature profile of a vaccine payload box can be defined as shown in Figure 8 . The temperature profile can be divided into two regions: a temperature-tolerable zone and a temperature non-tolerable zone. If the amount of dry ice left in the vaccine package can't maintain the payload box temperature within the required temperature range of between −80 °C and −60 °C, it is in the temperature non-tolerable zone which occurs at approximately loss of 86% dry-ice in this test. Otherwise, it is in the temperature-tolerable zone, which can be further divided into two stages: a stable period and an unstable period. In the stable period, the dry ice can keep the payload box at the dry-ice boiling point (sublimation point) temperature of −78.5 °C. When the payload box temperature starts to rise but is still within the required temperature range, it is in the unstable period of the temperature-tolerable zone. The switch time from stable period to unstable period occurs at approximately loss of 31% dry-ice in this test. As shown in Figure 9 , four thermocouples were installed inside vaccine storage package to monitor their temperature variation during dry ice sublimation. The effects of the refrigeration system defrost operation on the payload box temperature are shown in Figures 7, 8 and 9 . In general, the defrost operation presents little influence on the payload box temperature at the stable stage of the temperature-tolerable zone, relatively small impact (<1 °C) on the payload box temperature at the unstable stage of the temperature-tolerable zone, and relatively large impact (>1 °C) on the payload box temperature in the temperature non-tolerable zone. Maintaining the payload box temperature within the temperature-tolerable zone is a key performance index to evaluate the ultralow-temperature refrigeration solution. According to our observations, to maintain the vaccine temperature within the required temperature range, the vaccines should be surrounded by dry ice during the transportation and storage to avoid them being exposed to the CO 2 gas released from dry-ice sublimation. The dry-ice sublimation rate was monitored through dry-ice weight variation during the tests. The purpose of monitoring the dry-ice sublimation rate was to quantify the extended dry-ice lifetime in dry ice vaccine package stored in refrigeration containers compared with in normal indoor conditions (21.1 o C). Test data for the dry-ice weight over time are shown in Figure 10 . Table 4 lists several dry-ice sublimation performance parameters, including sublimation rate. In both tests, the dry-ice lifetime was extended significantly when vaccine dry-ice packages were stored in refrigerated conditions vs. room ambient.  In test 1, the dry-ice lifetime was extended from 120 h to 237 h under the conditions of −30 °C refrigeration temperature and generic insulated box.  In test 2, the dry-ice lifetime was extended from 123 h to 306 h under the conditions of −35 °C refrigeration temperature and pharmaceutical box. The test results show that the surrounding temperature of the vaccine dry-ice package significantly affects the dry-ice sublimation rate. Besides the surrounding temperature of the vaccine dry-ice box, other factors influence the dry-ice lifetime include package materials, insulation, leakage situation. Therefore, storage package design is also very critical for the vaccine storage. According to National Institute for Occupational Safety and Health (NIOSH) (FSIS , 2020), the recommended CO 2 exposure limit for 15 min is 3%. A level of 4% is designated by NIOSH as "immediately dangerous to life or health." The NIOSH-recommended CO 2 exposure limits in the workplace is 0.5%, or 5,000 parts per million, for a 40 h work week. Because of the large amount of dry ice being used in the vaccine storage, maintaining a safe level of CO 2 concentration during loading/unloading is critical for this proposed vaccine distribution last-mile solution. Two ventilation strategies have been tested: fresh air ventilation (FAV) and door open strategies. The goal was to reduce the CO 2 concentration build-up inside the container by introducing fresh air or opening the door. Figure 11 shows the CO 2 concentrations in the refrigerated container during the entire dry-ice sublimation. A FAV strategy was applied during the test by opening the outtake of FAV. FAV is to provide ventilation for container. The air exchanger depends on static pressure differential, which will vary depending on the container and how the container is loaded. A vent position sensor is installed to determine the position of the fresh air vent ( upper or lower) and sends data to the controller display. The test results indicate that the FAV strategy is impractical for keeping the CO 2 concentration at a safe level for vaccine package loading/unloading. FAV took around 22 h to drop the CO 2 concentration from 20.6% to 6.4%; further reducing the CO 2 concentration to 0.6% can take around 6 days (147 h), which makes the strategy impractical. The door open strategy, as shown in Figure 12 , is a simple and effective solution to manage the CO 2 concentration within safe levels. The test results show that opening one of the rear doors (there are two rear doors) can reduce the CO 2 concentration close to a safe level (5%) within 20 min. In addition to the experimental approach, modeling and simulation, for example computational fluid dynamic (CFD) modeling, has been also used to address the challenges faced by the proposed vaccine distribution solution: refrigeration storage container with vaccine dry-ice package. The modeling and simulation can not only provide a visualization of the temperature distribution in the container, but also offer an optimal distribution of the boxes with an economic cost. The modeling results can provide guidance of experimental setups to save both time and cost. A 3D CFD model of a refrigerated container was developed using a commercial software, ANSYS/FLUENT (version 17.2) ( ANSYS, 2017) . The model was used to investigate the airflow, temperature distribution, and CO 2 distribution in the container The refrigerated container model used in this study has dimension of L = 11.65 m, H = 2.56 m, and W = 2.28 m. A Cartesian coordinates system was employed in which x, y, and z are the directions for width, height, and length of the container, respectively. A refrigeration unit was located at the back of the container with the return air grille located near the ceiling of the container. Insulated sandwich panels formed the ceiling and walls of the container. The interior surface of the wall and ceiling panels consisted of a stainless steel sheet, and the exterior side was covered with an aluminum sheet. The insulating material between the steel and aluminum sheets is polyurethane (PUR). The PUR panel was 136 mm thick for the ceiling and 66 mm thick for the walls. The stainless steel and aluminum sheets were both 2 mm thick. The air inlets (supply air) were located at the floor of the container, where a series of T-shaped floor beams were used as diffusers to disperse the airflow. Twenty boxes filled with dry ice, with L box = 0.54 m, H box = 0.5 m, and W box = 0.5 m, were deployed on the floor of the container. Because CO 2 leaks from the vaccine boxes to the container, in the present work, two-phase flow, air and CO 2 , was modeled. A mixture model was employed to describe the multi-phase flow and heat transfer. The governing equations of the mixture model include continuity, momentum, and energy equations. Where, ⃗ is the mass averaged (mixture) velocity and is the mixture density, p is the pressure and μ t is the turbulence dynamic viscosity. μ m is the viscosity of the mixture, λ eff is the effective conductivity and h l is the enthalpy for phase l. The detail calculation of these variables is elaborated by Sun et al, 2021. Air and CO 2 were treated as incompressible ideal gases in this model. The flow was turbulent in all conditions evaluated, so a standard k-ε model was adopted in the model. Therefore, the turbulence dynamic viscosity can be written as (4) The k-ε model can be found in the reference ( ANSYS 2017). The validated CFD model was used to simulate a refrigerated container exposed to a sunny day from morning (07:00) to evening (16:00). Several assumptions have been made for this simulation as following,  The temperature changes of three external surfaces of a refrigerated container caused by solar radiation referred to existing studies at Hakata Port, Japan (Shinoda, 2016 , Budiyanto, 2017 .  Because of the lack of temperature data on the container door surfaces and the relatively small surface area of the doors compared with the other surfaces, the door surfaces were assumed to be perfectly insulated in the present case study.  The simulation started at 07:00, when the air temperature in the container was −33 °C. Initially, the CO 2 concentration was 409.8 ppm, which is the average atmospheric CO 2 concentration (Lindsey, 2020) .  After 07:00, CO 2 started to leak from the vaccine packages to the container at a rate of 0.0832 kg/h per package based on the experimental measurement.  The supply air was blown from the T-shaped floor to the container starting at 07:00 at a rate of 0.38 m 3 /s (800 CFM) and a temperature of −37.3 °C.  In this simulation, other operations, such as defrost and opening/closing of the door, were ignored to simplify the model. To investigate the temperature distribution in the container and the insulating panels, five vertical lines were defined, which are line 1 to line 4 (green lines in Figure 13 ) and the center line (red line). The center line is the cross-line of the planes x = 1.14 m and z = 5.825 m, located in the center of the container. The temperature distribution at that moment is confirmed by Figure 14a , which depicts the temperature distribution on the two cross-sections at 08:00. Subsequently, because of the solar radiation and higher ambient temperature, the container temperature started to rise, so the container temperature at 12:00 was much higher than at 08:00. Also, a constant temperature region below y = 0.25 m existed only at 16:00, whereas temperature kept increasing after leaving the T-shaped floor at 08:00 and 12:00 because the circulating flows had not been developed at 08:00 and 12:00. Therefore, at 08:00 and 12:00, supply air flowed directly upward and could not generate a constant temperature region. To visualize the temperature distribution near the vaccine packages, Figures 14d , e, and f show the temperature profiles at the cross-section y = 0.4 m at 08:00, 12:00, and 16:00, respectively. These profiles indicate that at 08:00, near all the packages, the temperature was lower than −36 °C. Because of the temperature of the external surfaces, the temperature increased to about −35 °C at 12:00. At 16:00, the temperature near the packages was even lower than at 12:00 because of the circulating flows. In addition to temperature distributions within the refrigerated container, the CO 2 concentration was also studied because of the use of dry ice in the vaccine packages. Because of CO 2 accumulation, the CO 2 concentration kept increasing from 08:00 to 16:00, from about 10,000 to 85,000 ppm. Because of the lack of circulating flows, at 08:00, the CO 2 concentration was increasing from the door to the refrigeration unit (Figure 15a and b) following the airflow direction. Previous study shows that the CO 2 concentration above 1,000 ppm could potentially be harmful to human health [9] , so even at 08:00, the container was not safe to human health. Therefore, extra fresh air ventilation is essential to ensure the safety of people entering the container. Some key conclusions from this study are as follows:  At 16:00, in the container, a circulating flow pattern was observed. Such a pattern did not exist at 08:00.  The circulating flow acted as a limiting factor to resist heat transfer in the container when y < 0.25 m. Therefore, a constant temperature region existed when y < 0.25 m at 16:00. As a result, the temperature around the box at 16:00 was lower than at 12:00. However, at 08:00 and 12:00, such a constant temperature region did not exist in the container.  The PUR panels played a critical role in keeping the temperature in the container low even on a summer day.  The circulating flow also acted as a cap to resist CO 2 spreading in the container, leading to a low CO 2 concentration region when y < 0.25 m.  The CO 2 concentration was about 85,000 ppm in the container at 16:00 because of the accumulation of CO 2 . Even at 08:00, the CO 2 concentration was still above 10,000 ppm, which is much higher than what is safe for human health. Therefore, extra fresh air ventilation is essential to ensure the safety of people entering the container for off-loading the product. The purpose of this study was to mitigate the last-mile challenge faced by COVID-19 vaccine distribution, especially for rural or suburban areas, where local towns, pharmacy chains, and hospitals may not have the infrastructure required to store the vaccine at the required temperature. This study proposes an effective, secure, and safe ultralow-temperature refrigeration solution that uses commercially available refrigeration container units that are operated in the global cold-food chain and can be modified quickly to support the vaccine distribution supply chain according to the more demanding requirements of vaccine storage. To evaluate the technical merits and potential risk of this proposed solution, tests and computational fluid dynamics simulation studies were conducted to evaluate temperature distributions of refrigeration storage containers and vaccine packages, as well as the dry-ice sublimation process. According to the test and simulation results, several conclusions and observation are summarized as follows.  For the tested refrigeration storage container with supply air at the floor and return air at the top, refrigeration storage container temperature gradient analysis though testing and simulations suggest placing the vaccine boxes at the floor to minimize the refrigeration storage container air temperature variation, especially during the defrost period, and maximize the dry-ice lifetime.  The refrigeration storage container temperature increased significantly during the defrost operation.  Along with dry-ice sublimation, temperatures inside the vaccine package were not uniformly distributed. The dry-ice section maintained its sublimation temperature of −78.5 °C, whereas the temperature of the CO 2 gas section rose because of the heat transfer or leakage with air outside the vaccine box.  A generic temperature profile of vaccine payload box was defined to include two regions: temperature tolerable zone and temperature non-tolerable zone. If the amount of dry ice left in the vaccine package cannot maintain the payload box temperature within the required temperature range, it is in the temperature non-tolerable zone. Otherwise, it is in the temperature-tolerable zone, which can be further divided into two stages: a stable period and an unstable period. In the stable period, the dry ice can keep the payload box at the dry-ice boiling point (sublimation point) temperature of −78.5 °C. When the payload box temperature starts to rise but is still within the required temperature range, it is in the unstable period of the temperature-tolerable zone.  In general, the defrost operation presents little influence on the payload box temperature at the stable stage of the temperature tolerable zone, relatively small impact (<1 °C) on the payload box temperature at the unstable stage of the temperature tolerable zone, and relatively large impact (>1 °C) on the payload box temperature in the temperature non-tolerable zone. Therefore, ensuring enough dry ice surrounds the vaccine inside the dry-ice package is critical to maintain the vaccine temperature within the required temperature range.  The surrounding temperature of vaccine packages was shown to significantly affect the dry-ice sublimation rate. Under the lab test setup, compared with typical 21.1 °C ambient storage, using a −35 °C storage temperature extended the dry-ice lifetime by approximately 2.49 times. Other factors influencing dry-ice lifetime include vaccine box design, insulation, leakage rate, and dry-ice fill density.  The existing FAV strategy is not a practical solution to maintain a safe CO 2 concentration level for vaccine package loading/unloading procedures because of the long time to reduce the CO 2 concentration to near the safe level of 0.5%.  We recommend opening the rear door for at least 30 min before accessing the refrigerated container to load/unload vaccine package. In addition, we recommend having portable CO 2 monitoring to ensure that the CO 2 concentration is below the safe limit before accessing the refrigeration storage container. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Funding for this research was provided by the US Department of Energy, Office of Energy Efficiency and Renewable Energy. The authors would like to thank Erika Gupta and Samuel Petty, Program Manager of Building Technologies Office, for their support of this work. The authors would like to thank the support from the Oak Ridge National Laboratory Operations and Facilities team (Gerald Barth, Margret Smith, Brandy Milun). The authors would like to thank Doug Auyer, David Brisson from Carrier Global Corporation to provide some test data. The authors would like to thank Yarom Polsky) and Carrier Global Corporation Safety and Health Group, Carbon dioxide health hazard informaiton sheet How quick can US distribute a Covid-19 vaccine? Here are the four biggest logistical challenges, Recode Atmospheric Carbon Dioxide, Climate.gov Vaccinating rural America won't be easy Cooling load calculation of cold storage container for vegetables case study C Campus-UISI Energy Saving Effect of Roof Shade for Reefer Container in Marine Container Terminal COVID-19 vaccine distribution solution to the last-mile challenge: Cryogenic refrigeration to facilitate the maintenance of ultralow temperature