Natural Gas - Transient Flow in Gas Transmission Lines

The transient flow of gases in long pipe lines is a problem of industrial interest. The present discussion deals with the application of the conservation of momentum and material to the transient flow of gas in conduits which are long with respect to their diameter. Solutions of the resulting equations include consideration of variations in friction and specific volume of the fluid as functions of time and position along the conduit. The changes in state along the pipe line are established as a function of time. The methods proposed may be extended to other situations. The flow of compressible fluids under steady conditions' has been considered in detail. Primary emphasis was placed on the development of simple means of taking into account the non-uniform flow' realized in long conduits under such conditions. Binder2 discussed the essential differences between the flow of compressible and uoncompressible fluids. Bonilla" proposed direct methods of solution for equations describing isothermal flow in long pipes. Similar considerations limited to steady flow have been discussed by Ruth.' Palsgrove 5 sugrested a graphical means of solving specific problems relating to the flow of gases in conduits of uniform section. Some of the basic concepts of the flow of a compressible viscous fluid in straight tubes have been described by Kuo.6 Problems of this type require knowledge of the equation of state of the fluid. In many cases it liar been convenient to assume that the fluid follows the relationship ascribed to a perfect gas. The treatment of unsteady flow of compressible fluids is less complete. This state of affairs result* from the fact that the differential equations relating pressure and flow rate with time and position for unsteady flow are nonlinear in character.' Muskat and Wyckoff presented an excellent treatment of the situation in laminar Bow where inertia forces were not of controlling importance. Somewhat earlier Hurst9 discussed the uncteady flow of fluids in petroleum reservoirs. Helhering-ton. MarRoberts. and Huntington'" established experimentally the pressure distribution and accumulative flow rate of natural gas through an unconsolidated sand. Recently Joffe11,12 conhitlcretl the specific problem of the flow of natural gas in relatively long.,. pipe lines for unsteady conditions with emphasis upon the economic feusibility of the storage of gas in such a reservoir. These discussions were based upon approximations of the partial differential equations applying to the condititons of flow. The results were descriptive of the physical situation and indirated the effertiveness of large pipe lines as cyclic storage facilities for natural gas. The advent of analog computers" has increased materially the practicality of solving nonlinear partial differential expressions. Several iterative numerical and graphical methods14,15,16 have been developed to solve nonlinear equations. However, it appears that the pair of simultaneous nonlinear partial differential equations involved in problems of this nature are not directly susceptible to solution by presently available analog computer or formal iterative numerical techniques. Nevertheless. graphical analyses afford a means of obtaining data of engineering value without special equipment. MATHEMATICAL CONSIDERATIONS Since the length of the conduit is large in comparison to the cross section in the cases to be considered, it is possible to treat the characteristics of the flow as independent of position in the cross section of the conduit. This possibility follows since the energy- connected with the velocity profile is small compared to that required to overcome the friction associated with its movement throughout the length of the line. Average values of the pressure. velocity. kinetic energy. momentum, and specific volume may he assigned to the fluid at any particular cross section. From a mathematical viewpoint this may be expressed in the following way for velocity: µ = µ(?X)..........(1) The uncertainties introduced by such assumptions usually are negligible when the ratio of length to diameter is greater than the order of 1.000. On this basis a volume element ma!. be considered to occupy the entire cross section of the conduit. Such an element is presented in Fig. 1.