MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/210727 210727 Research Article Analytical Model of Fluid Flow through Closed Structures for Vacuum Tube Systems Park Chan-Hee 1 Cheon Dae-Sung 1 Park Joonam 2 Shateyi Stanford 1 Korea Institute of Geoscience and Mineral Resources (KIGAM) Daejeon 305-350 Republic of Korea kigam.re.kr 2 Wonkwang University Iksan Jeonbuk 570-749 Republic of Korea wku.ac.kr 2015 1692015 2015 25 02 2015 30 07 2015 1692015 2015 Copyright © 2015 Chan-Hee Park 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.

An analytical model is developed to evaluate airtightness, which is one of the most important requirements of vacuum tube transportation systems. The main objective of the model is to anticipate the pressure inside a closed structure, which initially decreases and then rises with time owing to the inflow of the air outside. The model is formulated by using Darcy’s law and by solving the differential equations that consider the air permeability of the material and physical configuration of the tube. Equations are derived for a tube section in two cases: one with constant thickness and another with variable thickness. Although the developed model must be verified experimentally, results simulated by using the model with several assumed sections were consistent with the findings of a previous study. The mathematical model developed here to predict the behavior of the pressure inside a vacuum tube structure could be effectively used to evaluate the airtightness of vacuum tube transportation systems. Such data could provide background technical information for the practical design of the system.

1. Introduction

The air resistance that acts against a transportation vehicle increases with speed. Hence, minimizing air resistance could enhance vehicle speed without excessive consumption of propulsion energy. Vacuum tube systems, which constitute an innovative type of transportation mode where high speeds can be achieved for a vehicle by reducing the air pressure inside a tube structure to minimize air resistance (Figure 1), are increasingly attracting attention for their effectiveness and low emission levels. During the past several decades, transportation systems that use reduced air resistance such as the Vactrain system, Cargo Cap, and PCP (Pneumatic Capsule Pipeline) have been proposed . Although these systems differ in their detailed characteristics, they share the basic idea of maintaining the pressure inside a tube structure at a low level. Although an actual vacuum tube system has not yet been realized, according to a recent study , the speed of the vehicle could reach up to 700 km/h provided that the pressure inside the tube is maintained at 0.1 atm. The system is therefore considered as a strong candidate that could provide an alternative to existing aviation systems.

Concept of vacuum tube transportation system.

One of the most important requirements that must be met by the infrastructure of a vacuum tube system is a sufficient level of airtightness. Otherwise, the capacity of the vacuum pumps would need to be inevitably increased, resulting in an increase in construction costs. In addition, the operation interval of the pumps would need to be reduced, increasing maintenance costs as well.

In a vacuum tube, the inside pressure is much lower than the atmospheric pressure outside. Air will then inflow into the vacuum tube structure from the outside unless the tube is completely airtight. The flow rate of air going through a porous medium can be calculated using Darcy’s law , given the pressure values on both sides. The pressure inside a vacuum tube, however, does not remain constant but increases with time because the tube is a closed hollow structure with finite volume; in contrast, the pressure outside the tube, that is, the atmospheric pressure, remains constant (Figure 2). The pressure inside the tube is therefore treated as a function of time formulated through a differential equation. Mathematical expressions that describe the pressure change inside a vacuum tube with variable thickness are formulated in this study. The differential equation is set up by applying Darcy’s law and considering the physical configuration of the tube structure and initial pressure condition.

Schematic view of a circular vacuum tube with constant thickness.

2. Flow Model for Tube Structure with Constant Thickness

Park et al.  have carried out both analytical and experimental studies on evaluating the airtightness of concrete tube structures with different joint configurations. Let us consider a circular tube section with constant thickness as shown in Figure 2, where the initial pressure inside the tube is greatly reduced. As noted above, unlike the pressure outside the tube, which remains constant, the pressure inside the tube increases as air flows into it. Given that the temperature remains constant, the inside pressure increases in proportion to the amount of air inflow. That is,(1)P1·Q1t=dP2tdt·V2,where Q1 is the flow rate (m3/s) and P1 and P2 are the pressures (Pa) outside and inside the tube, respectively.

Meanwhile, the flow rate of compressible fluid that flows into a tube structure can be expressed by applying Darcy’s law as follows:(2)Q1t=kAP12-P2t22P1μh,where k is the intrinsic permeability of the tube material (m2), μ is the viscosity of the fluid (kg/m·s), h is the thickness of the tube (m), and A is the area across which the air flows per unit length of the tube (m2).

Substituting (2) into (1), the gradient of the pressure inside the tube with respect to time is given as follows:(3)dP2tdt=kAP12-P2t22V2μh=kAP122V2μh-kA2V2μh·P2t2.

By solving the differential equation of form y=b-a·y2, the pressure inside the tube is expressed as(4)P2t=P1·1+C1·exp-kAP1/μhV2·t1-C1·exp-kAP1/μhV2·t,where C1 is a constant. The pressure inside the tube when t=0, P2(0), will be the target pressure level for operating the vacuum tube transportation system. C1 is determined as C1=(β-1)/(β+1) when β  (=P2/P1) is defined as the initial pressure ratio. If the initial condition is set so that the pressure inside the tube is 10 percent of the outside pressure; that is, β=0.1; the pressure inside the tube is expressed as (5)P2t=P1·1-0.8182·exp-kAP1/μhV2·t1+0.8182·exp-kAP1/μhV2·t.

Using (5), the pressure change inside a vacuum tube structure with time can be anticipated. For example, Figure 3 shows the relation between time and the pressure inside a tube with a physical configuration of r0=1.0 m and h=0.15 m for several different values of intrinsic permeability using (5).

Pressure change with time at several levels of air permeability.

3. Flow Model for Tube Structure with Variable Thickness

The mathematical model described in the previous section was established for the special case where the thickness of the tube section was constant. The actual thickness of a vacuum structure, however, may not always be constant; instead, it may vary along the radial direction because a tube section of a structural member is designed in consideration of various loads (e.g., dead, live, and special loads), in addition to cost and constructability (see Figure 4 for examples).

Examples of tube sections with variable thickness.

Let us consider a porous medium with a thickness that varies along one direction (the z-direction) as shown in Figure 5. That is, the width b and height l of the medium are constant but the thickness h varies along the z-direction. Assume that a fluid flow of Q1 occurs in the y-direction through the medium owing to a pressure drop between P1 and P2. The flow rate dQ1 corresponding to infinitesimal area dA=bdz satisfies (6)dQ1=kμhz·P12-P222P1·dA.

Flow of fluid through a porous medium with variable thickness.

The total flow rate Q1 can then be formulated as(7)Q1=kμ·P12-P222P1·0lbhzdz.

Figure 6 shows a schematic view of a tube structure with a length of b, of which the thickness is variable along the radial direction. Figure 7 shows a more detailed view of the section where the distances from the centroid to the inner wall and the inner wall to the outer wall are a(θ) and h(θ), respectively. The gradient of the pressure change at time t inside the tube is (8)dP2tdt=P1V2Q1t=P1V2·k·bμ·02πrθhθdθ·P12-P2t22P1,where rθ=aθ+(1/2)h(θ).

Tube structure with variable thickness.

Physical definitions of variables for cross section with variable thickness.

Upon solving the differential equation, we obtain(9)P2t=P1·1+C1·exp-g·t1-C1·exp-g·t,where(10)g=P1V2·k·bμ·02πrθhθdθ.

4. Illustrative Example

This section provides illustrative applications of the model developed in this study. Let us first consider a tube section, denoted as Section B, which has a square outer section with a length of d on each side and a circular inner section with a radius r0 (Figure 8). Assuming that the angle from the vertical line passing through the centroid is θ and the thickness at that location is h(θ),(11)hθ=d2cosθ-r00θπ4,rθ=r0+hθ20θπ4.

Definition of Section B.

The pressure inside the tube at time t is then expressed as(12)P2t=P1·1+C1·exp-g1·t1-C1·exp-g1·t,where (13)g1=P1V2·k·bμ·02πrθhθdθ,02πrθhθdθ=8·0π/4d+2r0cosθ2d-2r0cosθdθ=16d·ln2d·cosθ-4r+sinθ·16r2-4d2/2d-4r·cosθ16r2-4d2-4θ.

Let us consider one more tube section, denoted as Section C, whose outer section is rectangular and inner section is elliptical (Figure 9). In this case, r(θ) and h(θ) are defined as follows: (14)r0θ=a·bcosθ·a2+b2·tan2θ,hθ=dv-2r0θ·cosθ2cosθ0θarctandvdhdv-2r0θ·sinθ2sinθarctandvdhθπ2,rθ=r0θ+hθ2.

Definition of Section C.

Figure 10 shows the pressure change as a function of time for tube structures with three different sections when the intrinsic air permeability is k=1.0×10-16 m2. Section A was a circular section with constant thickness as in Figure 2, with a diameter of 2.0 m and a thickness of 0.15 m. Section B had a square outer section with a length of d=2.3 m on each side and a circular inner section with a radius of r0=1.0 m. Section C had a rectangular outer section and an elliptical inner section with parameters a=1.2 m, b=1.0 m, dv = 2.3 m, and dh = 2.7 m (Figure 9). The pressure change plots shown in Figure 10 are for the case where the initial pressure inside the tubes was 10% of the atmospheric pressure. A comparison of the plots shows that the pressure inside the tube rose most rapidly for the tube with Section A and most slowly for the tube with Section C. This indicates that Sections C and A are the most and least airtight, respectively, among the three sections. Note that the average thickness of Section A is smaller than that of the other two sections and that the volume of the space inside the tube is the greatest for Section C. The trend in Figure 10 is therefore consistent with the results of a previous study  showing that an increase in the volume or thickness of the tube enhances its airtightness.

Comparison of pressure behaviors for tubes with different sections.

Although further experimental verification is needed, the mathematical model developed here to predict the behavior of the pressure inside a vacuum tube structure could be effectively used to evaluate the airtightness of vacuum tube transportation systems. The data thus obtained would provide background technical information for the practical design of the system.

5. Conclusions

An analytical model is developed to evaluate airtightness, which is one of the most important requirements of vacuum tube structures. The model anticipates the pressure inside a closed structure, which initially decreases and then rises with time. This is achieved by using Darcy’s law and solving differential equations that consider the air permeability of the material and physical configuration of the tube. Equations for the pressure change are formulated for tube sections with uniform thickness and variable thickness. The mathematical model developed is applied to several tube structures with assumed sections. The results show that the tube structure with the largest internal volume is the most airtight. In addition, given the same internal volume, the tube with a section that has a greater average thickness is more airtight. Although further experimental verification is needed, the results simulated by using the developed model are consistent with previous results showing that an increase in the volume or thickness of the tube enhances its airtightness. The mathematical model developed here for predicting the behavior of the pressure inside a vacuum tube structure could be effectively used to evaluate the airtightness of vacuum tube transportation systems. The obtained data could provide background technical information for the practical design of the system.

Disclaimer

The opinions, findings, conclusions, and recommendations expressed in this paper are those of the authors and do not necessarily reflect those of the sponsor.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by a grant from the International Cooperation R&D Program of the Korea Institute of Energy Technology Evaluation and Planning (KETEP), funded by the Ministry of Knowledge Economy of the Republic of Korea (2010T100100963).

Cassat A. Bourquin V. Mossi M. Badoux M. Vernez D. Jufer M. Macabrey N. Rossel P. SWISSMETRO—project development status Proceedings of the International Symposium on Speed-up and Service Technology for Railway and Maglev Systems (STECH '03) 2003 Tokyo, Japan 453 410 Becker D. A. Characterization and use of the new NIST rapid pneumatic tube irradiation facility Journal of Radioanalytical and Nuclear Chemistry 1998 233 1-2 155 160 10.1007/bf02389664 2-s2.0-0031877080 Choi S. K. The necessity of development of super speed maglev tube train Korean Society for Railway 2009 11 4 13 18 Korea Railroad Research Institute Development of new infra-structure technology for ultra high speed tube train (PK09001C) Project Report 2009 Mehta P. Monteiro P. Concrete: Microstructure, Properties, and Materials 2005 New York, NY, USA McGraw-Hill Park J. Kim L.-H. Nam S.-W. Yeo I. Performance evaluation of airtightness in concrete tube structures for super-speed train systems Magazine of Concrete Research 2013 65 9 535 545 10.1680/macr.12.00161 2-s2.0-84880724682