arxiv:711.839v1 [physics.ins-det] 6 Nov 27 Conduction Effect of Thermal Radiation in a Metal Shield Pipe in a Cryostat for a Cryogenic Interferometric Gravitational Wave Detector Takayuki Tomaru A, Masao Tokunari B, Kazuaki Kuroda B, Takashi Uchiyama B, Akira Okutomi B, Masatake Ohashi B, Hiroyuki Kirihara B, Nobuhiro Kimura A, Yoshio Saito A, Nobuaki Sato A, Takakazu Shintomi C, Toshikazu Suzuki A, Tomiyoshi Haruyama A, Shinji Miyoki B, Kazuhiro Yamamoto B, Akira Yamamoto A A High Energy Accelerator Research Organization (KEK), 1-1 Oho, Tsukuba, Ibaraki 35-81, Japan B Institute for Cosmic Ray Research, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa 277-8585, Japan C Nihon University, 4-2-1-62 Kudankita, Chiyoda 12-73, Japan A large heat load caused by thermal radiation through a metal shield pipe was observed in a cooling test of a cryostat for a prototype of a cryogenic interferometric gravitational wave detector. The heat load was approximately 1 times larger than the value calculated by the Stefan-Boltzmann law. We studied this phenomenon by simulation and experiment and found that it was caused by the conduction of thermal radiation in a metal shield pipe. 1 Introduction A cryogenic interferometric gravitational wave detector is a next-generation detector for developing the field of gravitational wave astronomy. One of the features of the detector is that it uses cryogenic mirrors in the optical interferometer to reduce thermal noises [1], which are one of the principle sensitivity limitations. At present, a Large-scale Cryogenic Gravitational wave Telescope (LCGT) [2] has been proposed in Japan. As a prototype of the LCGT, the Cryogenic Laser Interferometer Observatory (CLIO) [3] has been constructed. It has two 1 m baselines and four cryogenic sapphire mirrors to form Fabry-Perot cavities. Each sapphire mirror is set in an independent cryostat and cooled by ultra-low vibration cryocooler systems [4, 5, 6]. Figure 1 shows the schematic representation of the CLIO cryostat. A feature of this cryostat is that it has a large laser opening that is facing a 3K beam duct since it is difficult to cool all of beam ducts in huge interferometers. In the case of CLIO, the diameter of the laser opening is 3cm. To reduce the 1
Vibration Isolation System Cryostat Cryocooler for Shield Pipe Double Pendulum Vibration Isolater 1K Beam Duct 3K 8K Laser Beam Heat Link Suspension Fiber 1K Mirror (2K) Vacuum Cryocooler 4K Figure 1: Schematic representation of the CLIO cryostat. To reduce the solid angle between the 3K beam duct and the cryostat to reduce the thermal radiation heat load, a thermal radiation shield pipe with a diameter of 3cm and length of 5m was introduced in the beam duct. It was separated from the cryostat and independently cooled by an ultra-low 8K cryocooler system. solid angle between the 3 K area and the cryostat, i.e., to reduce the thermal radiation heat load from the 3K area to the cryogenic area, an aluminum pipe was installed as a thermal radiation shield in the beam duct. The diameter and length of the pipe are 3cm and 5m, respectively. It is separated from the main cryostat and independently cooled to 1K by a vibration-free 8K cryocooler system. In general, while designing cryostats, the thermal radiation heat load from an area with a temperature T 2 to an area with a temperature T 1 was estimated by using the Stefan-Boltzmann law: P = ǫ 1 ǫ 2 ǫ 1 + ǫ 2 ǫ 1 ǫ 2 σ(t 4 2 T 4 1 )A Ω 2π, (1) where ǫ i (i = 1, 2) is the emissivity; σ, the Stefan-Boltzmann constant; A, the surface area; and Ω, the solid angle between both the areas. This equation implies that the heat generated by thermal radiation is transferred from one surface to the other. The thermal radiation heat load in the CLIO cryostat was also estimated using this equation; however, an extremely large heat load over 3W, which was approximately 1 times larger than the design, was found in its cooling test. Since the cryostat reached the design temperature when the laser opening was closed, the large heat load could be generated due to the thermal radiation through the shield pipe. Here, it should be noted that 3K black body radiation consists of infrared (IR) light with a peak wavelength of 2 µm. Moreover, the IR reflectivity of aluminum is sufficiently large. Therefore, thermal radiation could be reflected and conducted through an aluminum shield pipe. We investigated this phenomenon by simulation and experiment. 2
2 Analysis by Simulation Figure 2: Ray-trace model for conduction of thermal radiation in a cylindrical pipe. To investigate the hypothesis of the conduction phenomenon of thermal radiation in a radiation shield pipe in the CLIO, we simulated it by using a ray-trace model. Figure 2 shows the schematic representation of the model. When IR light is isotropically emitted from the 3 K area, the power incident onto the area dx at a distance x from the end of the pipe is given as follows: dω = dω P(x) = P 2π, (2) 2πa2 dx, (3) (x 2 + a 2 ) 3/2 where P is the total emitted power from the 3K area; dω, the solid angle; and a, the pipe radius. This IR light reflects N(x) times inside the pipe wall with a reflectivity of R. When the length of the pipe is L, the number of reflections N(x) is given by ( L x ) N(x) = 1 + n, (4) 2x where n((l x)/2x) is a function that (L x)/2x is converted to integer. Although the reflectivity R depends on the wavelength and incident angle of light, we assumed it to be constant in this model. From the above equations, the IR power P ref reaching the 4K area after reflection in the pipe is calculated as P ref = P L R N(x)dΩ 2π = P L R N(x) a 2 dx. (5) (x 2 + a 2 ) 3/2 On the other hand, the direct power incident on the 4K area without reflection, which is described by the Stefan-Boltzmann law in equation (1), is rewritten as a 2 P SB = P 2L2. (6) Table 1 lists the calculated results. This result shows that the total power reaching the 4 K area through a shield pipe was several hundreds times larger than the direct power incident without reflection. The conduction effect accounted for several dozen percent of the total emitted power from the 3K area. Therefore, the simulation result supported the hypothesis that the large heat load obtained in the CLIO cryostat was possibly caused by the conduction effect of thermal radiation in the radiation shield pipe. 3
Table 1: Simulation result of conduction effect of thermal radiation by a IR ray-trace model in the case of the CLIO. R shows IR reflectivity of aluminum pipe. (P ref + P SB )/(P SB ) implies the ratio of the total heat load through the shield pipe to the direct incident power without reflection, and P ref /P implies the ratio of the heat load with reflection in the pipe to the total emitted heat power from 3K. R =.9 R =.95 R =.97 (P ref + P SB )/(P SB ) 37 622 898 P ref /P 14% 28% 4% 3 Experimental Verification Figure 3: Experimental setup for measurement of heat load through an aluminum shield pipe by using a cryostat of early prototype of cryogenic interferometer, called CLIK. The method for measurement was calorimetry. To verify the phenomenon of thermal radiation conduction in the shield pipe, we performed a model experiment by using a cryostat of the early prototype of a cryogenic interferometer, called CLIK. Figure 3 shows the experimental setup. We used the calorimetric method. A pure aluminum pipe with a diameter of 7cm, length of 1.35m, and temperature of 1K, was installed in a vacuum duct. A thermal radiative source was set in front of the 3K side of the pipe, and a bolometer was placed at the opposite side of the pipe in the cryostat. Paper was used as the thermal radiative source since it is well known that it has a large emissivity of approximately.9. The bolometer consisted of an IR absorber called ultra-black NiP (UB-NiP), a film heater, a thermometer, and 4
a thin aluminum heat link. The UB-NiP had a large IR absorption rate of approximately 8% and small outgassing ratio in a high vacuum system [7]. A silicon thermometer made by Lake Shore Inc. was attached to the UB-NiP in order to measure its temperature. A film heater was used to calibrate the relation between the heat load and the temperature of the bolometer, and the heat link was a heat path between the bolometer and the 4K stage in the cryostat. The measurement was performed twice and good repeatability within 2% was confirmed. From this measurement, we observed a heat load of 39mW on an average. This value was compensated by using the emissivity of paper and absorption rate of the UB-NiP. This observed heat load was 74 times larger than the value estimated by the Stefan-Boltzmann law in equation (6). To confirm that the experimentally observed heat load was certainly caused by the conduction effect of thermal radiation modeled in section 2, we checked the IR reflectivity of the aluminum shield pipe by using an IR free electron laser (FEL) at Tokyo University of Science [8]. In this measurement, the regularly reflected laser power on a sample surface was compared to the injected laser power. The sample was fabricated by cutting a part of the aluminum shield pipe used in the CLIK experiment. Although large errors of several percent due to the instability of the FEL power and inaccuracy of detectors were present in the measurement data, the reflectivity of the aluminum pipe at a wavelength of 8.33 µm was approximately 98%, and its dependence on the laser injection angle was small. The IR reflectivity of the aluminum pipe estimated from the CLIK experiment and equation (5) was 95.%. Therefore, we concluded that these results showed the legitimacy of the conduction model of thermal radiation in a shield pipe. 4 Discussion 3 Arbitrary Unit 25 2 15 1 5 1 2 3 4 5 x m Figure 4: Plot of integrand in equation (5) in the case of the CLIO. The unit of the vertical axis is arbitrary. The origin of the horizontal axis (x) was set as the end of the aluminum pipe of the 3K side. The reduction in the conduction effect of thermal radiation is an important issue for the cryostat design of cryogenic interferometers. In a simple case in 5
which a part of the shield pipe has a small IR reflectivity, we can easily estimate its reduction rate by using the ray-trace model described in section 2. Figure 4 shows a plot of the integrand in equation (5) in the case of the CLIO. This graph shows that the incident power at x < 1 m mostly affects the heat load in the cryostat. Therefore, it is effective to introduce an IR absorption area at the 3K side of the pipe. When a low reflectivity coating of R is used at a region between x = and x = s, equation (5) is modified to [ s P ref = P R N(x,s) R N(x,L) N(x,s) a 2 dx (x 2 + a 2 ) 3/2 + L s R N(x,L) a 2 (x 2 + a 2 ) 3/2 dx ]. (7) Table 2: Calculation results of (P ref + P SB)/P SB in equation (7). (P ref + P SB )/(P SB ) implies the ratio of the total heat load through the shield pipe to the direct incident power without reflection, and (P ref +P SB)/(P ref +P SB ) implies the reduction rate of the total heat load. In this calculation, the reflectivity R of aluminum pipe and reflectivity R of absorption area were assumed to be 95.% and 1.%, respectively. s [m].5 1. 2.5 (P ref + P SB)/P SB 111 38 7.4 (P ref + P SB)/(P ref + P SB ) 18% 6% 1% Table 2 shows the calculated results of (P ref + P SB)/P SB. This result indicated that the IR absorption area was required for more than half of the shield pipe. 5 Conclusion We studied the conduction effect of thermal radiation in a cryogenic shield pipe. Simulation results showed that thermal radiation could be reflected and conducted through the shield pipe and that the heat load caused by this effect could be several hundred times larger than the estimated value by the Stefan-Boltzmann law. A model experiment showed that the observed heat load was 74 times larger than the direct incident power calculated by the Stefan- Boltzmann law. From these results, a 95.% IR reflectivity of the cryogenic shield pipe was derived in the CLIK experiment, and this value was consistent with the result of the IR reflectivity measurement of the pipe by using the IR laser. Therefore, we confirmed the occurrence of the conduction effect of thermal radiation in a long shield pipe. Although the reduction of this effect was briefly discussed, a precise study of its reduction through experiments is a future issue. 6
Acknowledgements We would like to extend our thanks to Dr. T. Imai and Prof. K. Tsukiyama of the IR FEL Research Center, Tokyo University of Science, for the measurement of the IR reflectivity of an aluminum pipe. This work is partly supported by the cooperative program of the Institute for Cosmic Ray Research, University of Tokyo. References [1] T. Uchiyama et al.: Phys. Lett. A 261 (1999) 5. [2] K. Kuroda et al.: Class. Quantum Grav. 23 (26) S215. [3] S. Miyoki et al.: Class. Quantum Grav. 23 (26) S231. [4] T. Tomaru et al.: Cryocoolers 13 (25) 695. [5] R. Li et al.: Cryocoolers 13 (25) 73. [6] T. Tomaru et al.: Cryogenics 44 (24) 39. [7] T. Tomaru et al., J. Phys.: Conf. Ser. 32 (25) 476.; T. Tomaru et al.: J. Vac. Soc. Jpn. 48 (25) 31. (in Japanese) [8] H. Kuroda: Jpn. J. Appl. Phys. 41 (22) Suppl. 41-1, 1. 7