The approach to calculate the aerodynamic drag of maglev train in the evacuated tube

In order to study the relationships between the aerodynamic drag of maglev and other factors in the evacuated tube, the formula of aerodynamic drag was deduced based on the basic equations of aerodynamics and then the calculated result was confirmed at a low speed on an experimental system developed by Superconductivity and New Energy R&D Center of South Jiaotong University. With regard to this system a high temperature superconducting magnetic levitation vehicle was motivated by a linear induction motor (LIM) fixed on the permanent magnetic guideway. When the vehicle reached an expected speed, the LIM was stopped. Then the damped speed was recorded and used to calculate the experimental drag. The two results show the approximately same relationship between the aerodynamic drag on the maglev and the other factors such as the pressure in the tube, the velocity of the maglev and the blockage ratio. Thus, the pressure, the velocity, and the blockage ratio are viewed as the three important factors that contribute to the energy loss in the evacuated tube transportation.


Introduction
The speed of traditional trains is limited because of the dynamic friction between the wheels of the train and the fixed rail on the ground. When the trains are running at a low speed, most of the energy is consumed by friction. The trains can be levitated above the rail to avoid such friction with the technology of magnetic levitation [1]. There are three types of levitation technologies: electromagnetic suspension (EMS), electrodynamic suspension (EDS), hybrid electromagnetic suspension (HEMS) [2]. Even with these three methods, the velocity of trains could not be improved remarkably because of the aerodynamic drag. When the trains run at low speeds, this drag is not evident. At high speeds, however, the aerodynamic drag is too large to be neglected. Whatever the trains are levitated or not, the aerodynamic drag is the dominate part of drag when it runs at a high speed in the atmosphere near the earth's surface. At a speed range between 400 and 500 km/h, the aerodynamic drag accounts for 80 %-90 % of the total drag including the aerodynamic drag, the eddy resistance force, and the braking force [3]. The train speed is much lower than the airplane speed because airplanes flight in a circumstance of rarefied gas in the high altitude. In view of this fact, the evacuated tube transportation (ETT) was proposed to reduce the aerodynamic drag and improve the speed of the maglev train. Shen [4] and Yan [5] discussed the possibility, strategy, and the technical proposal for developing the ETT in China.
Theoretically speaking, when the inner part of the tube is in the condition of absolute vacuum, the aerodynamic drag for the levitation train inside the tube will be zero. However, this is very hard to realize. An alternative is to draw-off the gas partly and optimize the train shape. Therefore, the influence of the air pressure, the velocity, and the blockage ratio on the Maglev train in the evacuated tube system is a very interesting topic to study. Up to now, some research works have been done to explore what conditions are suitable for future ETT. Raghuathan and Kim [6] reviewed the state of the art on the aerodynamic and aeroacoustic problems of high-speed railway train and highlighted proper control strategies to alleviate undesirable aerodynamic problems of high-speed railway system. Various aspects of the dynamic characteristics were reviewed and aerodynamic loads were considered to study the aerodynamic drag [7]. Wu et al. [8] simulated the maglev train numerically with software STARCD based on the N-S equation of compressible viscosity fluid and k-e turbulence model. The flow field, the pressure distribution and the aerodynamic drag coefficient were also analyzed to illustrate the relationship between the aerodynamic drag and the shape of the train in evacuated tube. The pressure distribution in the whole flow field and the relation between the aerodynamic drag and the basic parameters were derived in [9]. Shu et al. [10] simulated the flow field around the train based on the 3D compressible viscous fluid theory and draw the conclusion that its aerodynamic performance is relevant to the length of the streamlined nose. In [11][12][13], the simulated results showed that the speed, the pressure and the blockage ratio significantly affect the aerodynamic drag of the train in an evacuated tube.
In this paper, the experimental system model developed by Superconductivity and New Energy R&D Center of South Jiaotong University was used to study the aerodynamic drag in the tube. The basic mass conversation equation and the momentum conversation equation [14] were used at first to deduce the expression that describes the relationships between the drag and the main parameters such as the tube pressure, velocity, and the blockage ratio. Then the aerodynamic drag is calculated with that expression. Finally, the calculated results are compared with the experimental data to verify the validity of the deduced expression. Considering an infinitesimal part for any arbitrary circumstance, it is well known that every part follows the laws of mass conversation equation and the momentum conversation equation. These two equations [15,16] are listed as: where q ¼ q ðx; y; z; tÞ q is density of the infinitesimal part, m velocity of the infinitesimal part, m i each component of velocity (different when in different coordinates), F i body force in each direction, l dynamic viscosity, and p pressure of the infinitesimal part. Eq. (1) means that the quality flowing into the infinitesimal part is equal to that out of this part. And Eq. (2) means that the rate of change of arbitrary mass's momentum is equal to the sum of the force acting upon it.
All of the following mathematical derivations in this paper are based on Eqs. (1) and (2). To study the evacuated tube transportation, an experimental model was developed.
This model called evacuated tube system for maglev train (ETSMT) includes three components: the evacuate tube, maglev, and propulsion system.
The tube is made of Perspex with the circumference of 10 m and the circular permanent magnetic guideway (PMG) is placed along the bottom of the tube. The positions of the tube, the train, and the PMG and the used rectangular coordinate system are shown in Fig. 1.
In Fig. 1, a 0 and b 0 stand for the width and the height of the tube respectively; a, b, and c stand for the width, the height, and the length of the maglev train, respectively. h 0 is the height of the PMG and h the levitated height.
As shown in Fig. 1, the train is levitated above the rail with height h and can only move in the x-direction due to the selfguiding characters of high-temperature-superconducting (HTS)-PMG system. For convenience, we suppose: 1) The tube is straight and the train runs along it only in the x-direction.
2) The magnetic flux of permanent magnetic rail is constant in the x-direction. 3) The pressure of any part of inner tube is constant in the z-direction with neglecting the atmospheric molecular mass. 4) The train is in the center of the tube whatever the blockage rate is.
In the Cartesian rectangular coordinate system the velocity is expressed as: where q is the atmospheric density, v x ; v y; v z are the components of velocity in x, y, z axis directions respectively, e x ; e y ; e z are the unit vectors of each axis. When considering the assumptions of (1) and (2), v y and v z are zero. Then Eqs. (1) and (2) are modified as: To demonstrate the evident effect of the pressure in tube, the streamlined nose of the train is not adopted. The schematic diagram in moving direction of train is shown in Fig. 2.
As shown in Fig. 2, the maglev train is levitated above the PMG and there is no dynamic friction between the train and the rail. According to assumption (2), there is no vibration in the z-direction, which means that all of the kinetic energy of the free levitated running train is consumed because of the aerodynamic drag after the train gains the initial kinetic energy. In an ideal situation, when the tube pressure is zero, the aerodynamic drag is equal to zero. Since this condition is almost impossible to realize, the actual practice is to pull the air out of the tube to form a suitable pressure. The purpose of this work is to explore the relationship between the drag and the tube pressure.
The aerodynamic drag of a running maglev in this system is composed of three components: (1) F 1 : the force on windward side of the train due to the collision between air and the train, (2) F 2 : the air friction on four side faces of the train and (3) F 3 : the force caused by the different pressures of windward side and the tailstock side of the train Each force will be discussed in following sections based on Eqs. (4) and (5). 2.1 The calculation of F 1 For simplification, a long section of the air ahead of the train is moving at the same velocity as the train because of the character of the air. An infinitesimal part of the air in area 1 is considered in Fig. 3. We suppose that the velocity of the thin layer of air is to vary after a tiny time dt after collision and the displacement of this layer is dx away from the windward side along x-axis within another tiny time dt.
The velocity of the infinitesimal air before collision is After collision, its velocity is equal to that of the train's. So the kinetic energy of this air is 1 2 qdxdydz Á v 2 x and we have the equation: F 1 is equal to zero when the velocity of the train is smaller than sound velocity because the velocity of atmospheric molecule is equal to the sound velocity after the collision with the windward side of the train. When the velocity of train is larger than sound velocity, the air column with the length of v c Á dt is affected within the period of dt and the velocity of that air column is approximatively equal to train velocity. The force can be calculated by combining Eqs. (6) and (7) and the momentum conversation equation: Thus, where p 0 ¼ 101; 325 Pa, and q 0 ¼ 1:293 kg=m 3 .
Because of the Brownian movement of molecules, it is reasonable to neglect F 1 when the train runs at a low speed.

The calculation of F 2
In this EETSMT, there are four side faces where friction force generates, as shown in Fig. 4, where F 2U ; F 2B ; F 2L and F 2R represent the frictions on upper side, lower side, left side, and right side respectively. The left side is toward the inside of the paper and the right side is toward the outside.
On upper side, considering a infinitesimal volume of dxdydz, the area of contact between the infinitesimal and the air is ds ¼ dxdy.
This component of the aerodynamic drag F 2 exits because of gas viscidity. The regularity of the fluid velocity distribution between the tube wall and train body side is described by the function of variable z [16]: where f 2U is velocity function of length n; b and b 0 are shown in Fig. 1. The relationship between the fluid internal friction stress and the velocity gradient according to Newton's proposal is where s friction stress, l viscosity, of 2U ðnÞ on is the change rate of the velocity from train body side to the tube wall. The friction stress is the force on a unit area with the direction perpendicular to the velocity, and the viscosity l is affected by the temperature instead of the air pressure. So the l is constant when the temperature is unchanged. From the analysis above, the friction of infinitesimal dxdydz is Likewise, F 2B ; F 2L and F 2R can be deduced.
2.3 The calculation of F 3 Figure 5 shows that F 3 is generated by pressure difference between the headstock and the tailstock side of the train. This force is larger when the velocity of the train is greater. In Fig. 5, the pressure of the inner tube is p. The train windward side is x-z side with area S. For a small time interval dt, the train moves with distance dx. And there is no interpenetration of air between different areas 1, 2, and 3 within dt. So the velocity of infinitesimal gas is v x when taking the train as a reference. According to the Bernoulli formula, the relationship between the pressure and the velocity at point A is: And the pressure at point B is So the pressure difference is:

Approximation of the total aerodynamic drag
The blockage rate is defined as where a; b; a 0 and b 0 are also illustrated in Fig. 1. When the train is running, the pressure, density, and flow velocity of arbitrary gas are functions of time and space. According to the assumptions and definitions mentioned above, Eqs. (9), (14), and (17) could be modified as: Fig. 6 The total calculated drag when blockage rate is 0.1, 0.2, 0.5, 0.8, respectively The total aerodynamic drag is expressed as According to Eq. (22), the relations between the total drag and the blockage ratio, the velocity and the pressure are shown in Fig. 6.
According to Fig. 6, we can calculate the total aerodynamic drag with Eq. (22) and the known parameters of blockage ratio b r , velocity v, and pressure p, and easily obtain the relationship between them. Figure 7 shows the ETSMT located in a tube made up of Perspex. It is vacuumized with a vacuum pump and the pressure inside the pipe can be detected by an instrument. We designed a control system to gain a fixed pressure ranging from 2,000 to 101,325 Pa. The experimental steps are listed:

The experimental system of the evacuated tube
1) The HTS maglev was fixed above the PMG with nonferromagnetic material at some height such as 0.01 m and then the liquid nitrogen was poured into the train. After the train was levitated, the non-ferromagnetic material must be removed from the PMG. 2) The opening hole of the pipe was covered and then the vacuum pump was started with the control system to reach the design pressures such as 10,000, 8,000, and 5,000 Pa and etc.
3) The liner induction motor was started and then the train could be drove to move when the maglev train runs near the LIM. Thus, the train speed can be accelerated to a certain value such as 3 m/s. 4) The LIM was stopped when the train's speed reached an expected value. Then the time difference between the position check points A and B was recorded to gain the decreasing train velocity. 5) The velocity was calculated with necessary parameters.
All parameters of this experimental system in Fig. 1 are listed in Tables 1 and 2. When T = 288.15 K, l ¼ 1:78 Â 10 À5 kg/(mÁs) and v c ¼ 340 m/s, the effect of the pressure variation could be neglected.

The comparisons of theoretical and experimental results
The total aerodynamic drag of the running train cannot be measured directly because the train is freely levitated above the PMG. The average velocity between check points A and Fig. 7 The experimental system of ETT     When the pressure is constant in the inner tube, the running time is less and the blockage ratio is larger. That is to say, if the blockage ratio is larger, so is the aerodynamic drag. Both the calculated and the experimental results show such a trend.

Conclusions
(1) When the pressure and the blockage ratio are constant, the aerodynamic drag is a quadratic function of the velocity. When the velocity of the train is bigger than the sound velocity, the formula of the aerodynamic drag becomes more complex. (2) If the blockage ratio is smaller, the drag becomes smaller. In practice, the blockage ratio is impossible to be very small because of the limitation of the pipe's section size. So a suitable blockage ratio should be determined in the design of the ETT system. (3) If the pressure in the tube is zero, the aerodynamic drag equals to zero no matter how the velocity and the blockage rate vary. This ideal condition is difficult to realize because of the technological limitation.
In this work, the comparison between the theoretical and experimental results was made when the velocity of the Maglev train is small. When the system runs in a lower pressure, more efforts must be made to solve more sophisticated technical problems. Thus, the speed, the pressure, and the blockage ratio each must have reasonable values for ETT. In such a case, the magnetic drag between the Maglev train and PMG may be negligible. In future study, we will consider the effect of magnetic drag between the Maglev train and PMG at high speeds.  Fig. 11 The relation between the velocity and the time when the pressure is 3,000 Pa The approach to calculate the aerodynamic drag 207