0
Research Papers

# Nonisothermal Transient Flow in Natural Gas Pipeline

[+] Author and Article Information
M. Abbaspour

Gregg Engineering, Inc., 403 Julie Rivers Drive, Sugar Land, TX 77479

K. S. Chapman

National Gas Machinery Laboratory, Kansas State University, 245 Levee Drive, Manhattan, KS 66502

J. Appl. Mech 75(3), 031018 (May 02, 2008) (8 pages) doi:10.1115/1.2840046 History: Received December 06, 2006; Revised August 31, 2007; Published May 02, 2008

## Abstract

The fully implicit finite-difference method is used to solve the continuity, momentum, and energy equations for flow within a gas pipeline. This methodology (1) incorporates the convective inertia term in the conservation of momentum equation, (2) treats the compressibility factor as a function of temperature and pressure, and (3) considers the friction factor as a function of the Reynolds number and pipe roughness. The fully implicit method representation of the equations offers the advantage of guaranteed stability for a large time step, which is very useful for gas pipeline industry. The results show that the effect of treating the gas in a nonisothermal manner is extremely necessary for pipeline flow calculation accuracies, especially for rapid transient process. It also indicates that the convective inertia term plays an important role in the gas flow analysis and cannot be neglected from the calculation.

<>

## Figures

Figure 3

Comparison of present work (a) with Kiuchi model (b) for Δt=1.0min

Figure 4

Comparison of present work (a) with Kiuchi model (b) for Δt=0.1min

Figure 5

Comparison of present work (a) with Kiuchi model (b) for Δt=0.01min

Figure 6

Oscillation on flow rate at closing valve process for Δt=0.01min

Figure 7

Variation of inlet flow rate without considering temporal inertia term in momentum equation Δt=0.01min

Figure 8

Solution for isothermal model by considering convective inertia term for (a) Δt=1.0min, (b) Δt=0.1min, and (c) Δt=0.01min

Figure 9

Comparison of the result for impact of convective inertia term and without convective inertia for Δt=0.1min

Figure 10

Comparison of the result for impact of convective inertia term and without convective inertia for different flow rates and Δt=0.1min

Figure 1

Mesh of the solution

Figure 2

Pipe information and boundary condition for flow through the valve

Figure 11

Variation of temperature (a) and Joule–Thompson coefficient (b) along the pipe with respect to time

Figure 12

Temperature comparison for current study and Coulter–Bardon equation

Figure 13

Temperature boundary condition at Node 1

Figure 14

Solution for nonisothermal condition for (a) Δt=1.0min, (b) Δt=0.1min, and (c) Δt=0.01min

Figure 15

Comparison of the result for impact of convective inertia term (isothermal) and nonisothermal condition for Δt=0.1min

Figure 16

Temperature (a) and compressibility factor (b) for nonisothermal condition and Δt=0.1min

## Errata

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections