CONSIDERACIONES SOBRE ECUACIONES

PARA FLUJO EN ESTADO ESTACIONARIO EN

GAS NATURAL TUBERÍAS
En este trabajo, una discusión sobre las particularidades de las ecuaciones de caída
de presión que se utilizan en el diseño de las tuberías de gas natural se llevará a
cabo. Se presentan varias versiones de acuerdo con los diferentes regímenes de
flujo bajo consideración y a través de la presentación de estas ecuaciones se discute
el soporte físico básico para cada una, así como su factibilidad.
Palabras clave: flujo de gas natural, caída de presión, gasoductos, Renouard.
INTRODUCCIÓN
El diseño de gasoductos y redes es comúnmente presentado a través de una serie
de procedimientos numéricos y recomendaciones, y generalmente las ecuaciones
de flujo son recomendadas por los diversos autores de acuerdo con el diseño y
cálculos comunes practicos, sin un análisis profundo del razonamiento físico básico
que está detrás de cada una de esas ecuaciones. En este trabajo una discusión
sobre las particularidades de las ecuaciones de caída de presión que se utilizan en
él se realizará el diseño de gasoductos y varias versiones presente. El desarrollo de
la ecuación de flujo es comúnmente encontrado en varios libros y publicaciones en
Mecánica de fluidos o conectado a tecnologías de utilización de gas industrial
(Pritchard et al., 1978; Katz y Lee, 1990), en consecuencia, es inútil presentar una
vez más tal derivación. El lector puede consultar el trabajo de Mohitpour y col. (2000)
donde dicho análisis se presenta para constantes y flujo de fluido compresible en
estado inestable.
Nomenclatura
C, C '= constante genérica
d = densidad relativa del gas, adimensional
D = diámetro interno de la tubería, m
E = término de energía potencial, ecuación (34), Pa2
f = coeficiente de fricción de Darcy, adimensional
g = aceleración gravitacional, m / s2
H = altura de los puntos 1 y 2, m
K = constante, sin dimensiones
L = longitud del tubo, m
M = masa molecular, kg / kmol
n = exponente para el caudal de gas (rango de valores entre 1.74y 2)
P, P '= presión absoluta, Pa
P1 = presión absoluta en la entrada de la tubería, Pa
P2 = presión absoluta a la salida de la tubería, Pa
Pavg = presión promedio de flujo, Pa
Pst = presión estándar, 1.01325 × 105 Pa
st Q & = caudal de gas volumétrico en condiciones estándar, m3 / s
R = resistencia al flujo por unidad de longitud de tubería
R = constante de gas universal, 8314.41 J / (kmol K)
Re = número de Reynolds del flujo de gas, adimensional
T = temperatura absoluta, K
Tavg = temperatura promedio de flujo, K
Tst = temperatura estándar, 288.15 K
zavg = factor de compresibilidad del gas, adimensional
zst = factor de compresibilidad en condiciones estándar, zst ≈ 1
Símbolos griegos
α = coeficiente, adimensional
β = coeficiente, adimensional
max = variación máxima del factor de fricción, adimensional
P = caída de presión, Pa
ε = rugosidad de la pared, m
η = factor de eficiencia, adimensional
μ = viscosidad dinámica del gas, Pa s
ν = viscosidad cinemática del gas, m2 / s
ρ = densidad de gas, kg / m3
τ = esfuerzo cortante, Pa
Subíndices
Aire: aire
Aplicación: aparente
Promedio: promedio
Cg: gas ciudad
Cr: crítico
Ent: entrada
St: estándar
Sta: estándar para aire
W: pared
1: relativo al punto genérico 1
2: relativo al punto genérico 2
ECUACIÓN GENERAL PARA FLUJO EN ESTADO ESTACIONARIO
Considerando la ecuación de momento aplicada a una porción de tubo de longitud
dx, dentro del cual fluye un fluido compresible con una velocidad media u, por
ejemplo, gas natural, suponiendo un estado estable condiciones donde ρ es la
densidad del gas, p es el gas estático absoluto presión, A es el área de la sección
transversal de la tubería (πD2/4) y dH representa una variación en alto, la ecuación
diferencial resultante es,

En la ecuación anterior f es el coeficiente de fricción de Darcy que es relacionado

con el esfuerzo cortante de la pared mediante,

Antes de proceder a la integración de la ecuación. (1) entre dos puntos genéricos 1

y 2 de la tubería, es más conveniente simplificar el término de disipación viscosa
para facilitar su integración, de lo contrario, debería saberse que la función u2 = f (x)
continúa con la integración. Pero como ρu = ṁ/A = C es una constante, de acuerdo
con el condiciones de estado estable, entonces,
Para integrar cada uno de los términos de la ecuación anterior, seguirá una
discusión más detallada.
Término de energía cinética (ρ2u du)
Como ρ = C u entonces,

Término de trabajo de fuerza de presión (ρdP)

Como

La integración de este término se convierte,

Se debe presentar alguna explicación sobre los parámetros, Tavg y Zavg. Para seguir
un enfoque más fácil para la integración, se usaron valores promedio para el factor
de compresibilidad y la temperatura del gas. La temperatura promedio del gas se
puede obtener a través de Mohitpour et al. (2000)
Mientras que la presión promedio se determinará a través de

Finalmente, el Zavg se determinará en base a los valores promedio obtenidos

anteriormente para P y T a través de tablas o fórmulas adecuadas, Smith (1990).

Término de energía potencial (ρ2 gdH)

Se asume nuevamente que (P M / (z R T)) 2 puede determinarse mediante valores

promedio que simplifican enormemente la integración, Mohitpour et al. (2000), y el
resultado final es

Este enfoque se justifica por el hecho de que no existe una relación matemática
simple entre altura, presión y temperatura, y el error introducido no tiene sentido.

Disipación de energía por fricción viscosa

Donde L es la longitud de la tubería entre los puntos 1 y 2. La suma de los términos

determinados previamente rinde,
Y, como de costumbre, el término de energía cinética es insignificante en
comparación con los otros términos,

Normalmente, se usa el caudal volumétrico en condiciones estándar Qst (288.15 K

y 1.1325 × 105 Pa), en lugar de la constante C, y combina la definición de esta
constante y del caudal másico

Y porque para un gas perfecto,

Donde d es la densidad del gas en relación con el aire y Mair = 28.9625 ≈ 29 kg /

kmol, la combinación de Eqs. (12) a (14) da en unidades SI que,
Según Osiadacz (1987), el flujo de gas real en una tubería es inferior al calculado
mediante la ecuación de flujo, es decir, la ecuación. (15) cuando η = 1, debido a la
fricción adicional impuesta por accesorios como curvas, tes, válvulas y también por
otros efectos como corrosión, ensuciamiento y deposición de polvo / óxido. Para
tener en cuenta tales reducciones de flujo adicionales de una manera simple y
efectiva, es una práctica común usar un factor multiplicador correctivo, el factor de
eficiencia η, que generalmente toma valores entre 0.8 y 1. Mohitpour et al. (2000)
sugieren valores de η entre 0.92 y 0.97, aunque la experiencia recomienda que para
tuberías viejas puede ser tan bajo como 0.7, Osiadacz (1987).
A veces, la influencia de η se introduce como una corrección de la longitud de la
tubería, es decir, L/η2 reemplaza a L y, en consecuencia, se utiliza una longitud
equivalente que puede tomar valores entre 1.56 L y L, correspondientes al rango
0.8 < η <1.

Teniendo en cuenta que , y que

, excluyendo por simplicidad el factor η, se puede

dar una forma más simple aunque menos común a la última ecuación,

Reemplazando todas las constantes de la ecuación anterior que están fuera de los
corchetes por C ’,

Para una tubería horizontal y siempre la

energía potencial o el término de elevación se pueden descuidar y, en consecuencia
 1
Con esta última ecuación, la importancia de la relación , conocida como el factor
√𝑓
de transmisión, así como del diámetro interno de la tubería D, en la tasa de flujo de
gas, se evalúa fácilmente. El factor de transmisión es un parámetro importante
porque es una representación simple de la transmisividad del gas dentro del
conducto, Smith et al. (1956) En el otro extremo, el diámetro interno de la tubería es
otro parámetro importante para el diseño de sistemas de tuberías, ya que, por
ejemplo, en el caso de una duplicación de su valor, el flujo aumentará 22.5 ≈ 5.66
veces.
Al reelaborar la última ecuación, el flujo volumétrico estándar de gas se puede
escribir como,

O de hecho,

Y de manera genérica se puede decir que

Donde n e R dependen de la ecuación utilizada para el cálculo del factor de

transmisión, 1 f.
Regímenes de flujo
Antes de discutir las ecuaciones disponibles en la literatura para el cálculo del factor
de fricción o, que es lo mismo, para el cálculo del factor de transmisión, es más
conveniente considerar, incluso en un sentido amplio, los diferentes regímenes de
flujo que generalmente se encuentran en Transporte de gasoducto. Para las líneas
de transmisión típicas con gas a alta presión y caudales de gas moderados a altos,
generalmente se observa una de las dos situaciones siguientes: Totalmente flujo
turbulento (flujo de tubería irregular) o flujo parcialmente turbulento (flujo de tubería
suave).
Considerando la definición del número de Reynolds, la relación entre la velocidad

promedio del gas y su caudal volumétrico, y sabiendo que

para condiciones de estado estacionario

Como , donde

Pero R = 8314.41 J / (kmol K) y recordando que las condiciones estándar se definen

Como Pst = 1.01325 × 105 Pa y Tst = 288.15 K entonces,

Como se sabe, para Re menor que 2100 el flujo es laminar, mientras que para Re
por encima de 2100 el flujo se considera turbulento. Entre los flujos laminar y
turbulento hay una región de transición, para la cual no hay correlaciones de caída
de presión disponibles.
Régimen de flujo laminar
Aunque generalmente las situaciones de flujo de gas en las tuberías son
turbulentas, en aras de la exhaustividad del análisis, la situación del flujo laminar
también se trata en el presente texto. Cuando el flujo laminar se desarrolla
completamente a través de un conducto de sección transversal circular, White
(1999) da el factor de fricción de Darcy, independiente de la rugosidad de la tubería.
Según Shah y London (1978), el flujo laminar dentro de un conducto está
completamente desarrollado más allá de una determinada longitud de la Cuaresma,
determinada por,

Para una región cuando el flujo aún está en desarrollo, es decir, para longitudes
más cortas que la longitud de entrada, el factor de fricción Darcy promedio o
aparente se puede calcular a través de Shah y London (1978)

Donde la longitud adimensional viene dada por x+

(Ecuación 29)
X es la longitud de la tubería y D es el diámetro interno de la tubería.
Cuando x+→ ∞, entonces fapp Re → 64.

REGIONES DE FLUJO PARCIAL Y TOTALMENTE TURBULENTO

En el flujo parcialmente turbulento, el grosor de la subcapa laminar es mayor que la
rugosidad absoluta de la pared del tubo, y luego hay una región laminar que cubre
la pared interna del tubo y una región turbulenta fuera de ella. Es como si hubiera
un flujo turbulento dentro de una tubería de paredes lisas y esa es la razón por la
cual la pared de la tubería se designa como hidráulicamente lisa, la caída de presión
es independiente de la rugosidad de la tubería, Munson et al. (1998) El factor de
fricción de Darcy en esta situación puede calcularse mediante la ecuación de
Prandtl-Von Karman, Mohitpour et al. (2000).
 Figura 1. Influencia de la rugosidad relativa en la transición del régimen de flujo
turbulento de tuberías hidráulicamente lisas a rugosas.
Con el aumento del número de Reynolds, el grosor de la subcapa laminar
disminuye, la rugosidad de la tubería se vuelve cada vez más importante, lo que
altera la subcapa laminar y, después de una breve región de transición, el
coeficiente de fricción se vuelve independiente del número de Reynolds, Munson et
al. (1998), es decir, estas son las condiciones del régimen totalmente turbulento. En
tales circunstancias, este factor de Darcy está determinado por la ecuación de
Nikuradse, Mohitpour et al. (2000)
(Ecuación 31)
La ecuación (30) para tuberías lisas se puede utilizar hasta que la influencia de la
subcapa viscosa se reemplace por la influencia de la rugosidad de la tubería Smith
(1990) y, a partir de entonces, la ecuación. (31) debe ser utilizado. Como se puede
ver en la Fig. 1, la transición de Reynolds depende de la rugosidad relativa y cuanto
más pequeña es, más tarde es la región de transición.
Designando por número crítico de Reynolds, el Reynolds a cuyo valor hay un
cambio brusco del flujo turbulento en tuberías lisas hacia el flujo turbulento en
tuberías rugosas, o Recr de flujo totalmente turbulento, se puede trazar un gráfico
muy conveniente, Fig. 2, que se presenta a través de una línea continua la relación
entre Recr y ε/Do, en otras palabras, muestra el límite entre estos dos regímenes de
flujo.
Dicha línea puede ajustarse mediante la siguiente ecuación,
(Ecuación 32)
Muchos investigadores, según lo mencionado por Smith (1990), adoptan una
modificación de la ecuación de Colebrook-White, utilizando la constante 2.825 en
lugar de 2.51, válida para las regiones turbulentas parciales, de transición y
totalmente turbulentas,

Figura 2. El número crítico de Reynolds Recr y la rugosidad relativa.

Sin embargo, este mismo autor (Smith, 1990), se refiere a que los datos
experimentales para tuberías comerciales no siguen esta ecuación de Colebrook-
White modificada, sino que tienen un cambio brusco entre los regímenes parciales
y totalmente turbulentos y, en consecuencia, las ecuaciones. (30) y (31) deben
usarse en lugar de la ecuación. (33) En tales circunstancias, la ecuación. (32) es
crucial para definir la elección correcta. La figura 1 es una guía óptima para una
comprensión clara de la situación.
REGÍMENES DE FLUJO MÁS FRECUENTES
La velocidad del gas dentro de las tuberías está limitada para reducir la erosión o el
ruido y, en consecuencia, el Reynolds también se limita. Entonces es interesante
comparar este valor límite práctico con el Reynolds crítico previamente definido para
la transición del régimen inter turbulento. La Tabla 1 muestra una comparación entre
estos dos números de Reynolds. Las condiciones de flujo utilizadas para los
cálculos fueron, presión absoluta de 5 atm para dar un límite máximo de Re para
las redes de distribución y utilización; densidad de gas natural 0.65 también por la
misma razón que ρ = P29 d / (ṜT); temperatura 288,15 K; velocidad media del gas
10 m/s, que es para la gran mayoría de las situaciones un valor límite de diseño;
viscosidad absoluta del gas 1.0758 × 10-5 Pa s, rugosidad absoluta de la tubería
0.0191 mm Smith et al. (1956), que es un límite inferior para algunos tubos de acero
y un buen valor máximo absoluto de referencia para todos los tubos de cobre y
polietileno, valor que impone un Reynolds crítico máximo para tubos de acero y un
mínimo para tubos de cobre y polietileno que probablemente aparecen en
condiciones normales de trabajo, ec. (32)
Tabla 1. Comparación entre Recr y Relim

Como se puede observar en la observación de la Tabla 1, el límite de Reynolds

siempre es inferior al Reynolds crítico, lo que significa que en las redes de
distribución y utilización, el subregimo de flujo es parcialmente turbulento, es decir,
tuberías hidráulicamente lisas. Se producirá un flujo totalmente turbulento
principalmente en las redes de primer y segundo nivel donde las presiones son más
altas, lo que impone la misma velocidad máxima, un aumento del volumen
específico y, en consecuencia, del número de Reynolds.
El caso anterior cambia sustancialmente si se utiliza una rugosidad absoluta de
0.046 mm, que es el valor sugerido típico para tubos de acero comerciales (White,
1999) (con una incertidumbre en su definición de ± 30%). En estas circunstancias,
también puede ocurrir un flujo totalmente turbulento para todos los diámetros de
tubería de acero disponibles, consulte la Tabla 1, en términos prácticos, solo las
tuberías de cobre o polietileno pueden garantizar condiciones exclusivas para un
flujo parcialmente turbulento.
ECUACIONES MÁS UTILIZADAS PARA FLUJO CONSTANTE
La ecuación (15) puede tomar una forma más simple cuando las constantes se
reemplazan por sus valores correspondientes y el término de energía potencial se
reemplaza por
(Ecuación 34)
Así que eso,
(Ecuación 35)
Para flujos de baja presión, es decir, cuando la presión manométrica es inferior a 50
mbar, la diferencia de presión P2 1 – P2 2 puede simplificarse, Osiadacz (1987),
(Ecuación 36)
Y teniendo en cuenta que para bajas presiones,
(Ecuación 37)
Donde ′ avg P es la presión promedio dentro del conducto de transporte,
(Ecuación 38)
El resultado es,
(Ecuación 39)
En la ecuación anterior, P1 y P2 pueden estar en presiones absolutas o
manométricas. Como se muestra, las ecuaciones de baja presión no son más que
simplificaciones de las ecuaciones de alta presión, que han sido modificadas
matemáticamente para facilitar los cálculos. Las ecuaciones de alta presión tienen
una aplicación universal en lugar de ecuaciones de baja presión, que tienen un
campo de aplicación más estrecho.
Como la ecuación para el factor de transmisión es genérica,
(Ecuación 40)
Mientras que el número de Reynolds se puede escribir como,
(Ecuación 23)
La ecuación (15) puede reelaborarse para tomar una forma más general
(Ecuación 41)
En las siguientes secciones, las ecuaciones más comunes utilizadas para el cálculo
de las caídas de presión en las tuberías de gas se analizan a la luz de lo que se ha
dicho anteriormente. Todas las ecuaciones se han convertido a unidades SI, aunque
varias de ellas se conocen en sistemas de unidades bastante diferentes.
ECUACIONES PARA FLUJOS TURBULENTOS EN TUBERÍAS
HIDRÁULICAMENTE LISAS
Panhandle A Ecuación
En este caso, el factor de transmisión viene dado por,
(Ecuación 42)
Lo que significa que α = 3.436 y β = 0.07305. Reelaborando la ecuación. (41) con
estos valores, se obtiene la ecuación Panhandle A en unidades S.I.
(Ecuación 43)
Según la sugerencia que se encuentra con frecuencia en la literatura técnica,
Mohitpour et al. (2000), que para esta ecuación, μ = 1.0415 × 10-5 Pa s, la ecuación
anterior se convierte en,
(Ecuación 44)
AGA Ecuación parcialmente turbulenta La ecuación de la American Gas Association
para parcialmente turbulento utiliza el factor de transmisión dado por la expresión
de Prandtl-von Karman,
(Ecuación 30)
Y reemplazándolo en la ecuación. (41) da,
(Ecuación 45)
BLANCO
Para simplificar el cálculo del factor de transmisión y asegurar al mismo tiempo que
los factores de fricción resultantes son cercanos a los obtenidos de una ecuación
de referencia como la de Prandtl-von Karman, un buen enfoque es utilizar la
ecuación de White (White, 1979) para f. Esta ecuación además de tener una
pequeña banda de error entre + 0.24% y -3.23%, para los resultados calculados,
cuando el Reynolds pasa de 2000 a 2 × 106, también tiene la ventaja de lo explícito.
Cabe destacar que, entre las diversas ecuaciones presentadas por Branco et al.
(2001) la ecuación blanca, que se muestra a continuación, es la que presenta
desviaciones más pequeñas en comparación con la ecuación de Prandtl-von
Karman (o AGA parcialmente turbulenta),
(Ecuación 46)
La ecuación de flujo correspondiente se obtiene mediante la introducción de la
ecuación blanca. (46) en la ecuación. (41) dando
 (Ecuación 47)
La figura 3 es una comparación de los coeficientes de fricción utilizados en las
ecuaciones de Panhandle A, AGA parcialmente turbulentas y blancas. Como se
puede ver, la f correspondiente a la última ecuación tiene valores más altos que los
de Panhandle A en números bajos de Reynolds. Este hecho es una demostración
de cuán limitado es el rango de aplicación de varias de las ecuaciones disponibles,
ya que muchas de ellas se desarrollaron para situaciones particulares, Schroeder
(2001). Como se mencionó anteriormente, es evidente la buena superposición de
los valores del factor de fricción obtenidos con el AGA parcialmente turbulento (o
Prandtl-von Karman) y las ecuaciones de White.

Figura 3 Comparación de factores de desgaste para Panhandle A, AGA p/t y

ecuaciones blancas.
FLUJO TURBULENTO EN TUBERÍAS ÁSPERAS
Panhandle B Ecuación
El factor de transmisión para esta ecuación viene dado por,
(Ecuación 48)
Y entonces la ecuación Panhandle B se obtiene introduciendo en la ecuación (41),
α = 8.245 y β = 0.01961,
(Ecuación 49)
Usando la viscosidad típica de μ = 1.0758 × 10-5 Pa s, la ecuación anteriorv se
convierte en,
(Ecuación 50)
Ecuación de Weymouth
De acuerdo con Mohitpour et al. (2002) esta ecuación sobreestima el cálculo de la
caída de presión y, por eso, se usa con mayor frecuencia en el diseño de redes de
distribución a pesar de ser menos exacta que otras ecuaciones. En este caso, el
factor de transmisión viene dado por,
(Ecuación 51)
Con el diámetro interno de la tubería en metros.
Reemplazando este factor de transmisión en la ecuación. (41) el resultado es,
(Ecuación 52)
Ecuación totalmente turbulenta de AGA
Esta es una de las ecuaciones más recomendadas y utilizadas para este tipo de
flujo, pudiendo estimar con alta precisión los valores de flujo y caída de presión, si
la rugosidad de la tubería se conoce con exactitud. El factor de transmisión viene
dado por la expresión Nikuradse, Eq. (31)
(Ecuación 31)
Que después de ser reemplazado en la ecuación. (41) da,
(Ecuación 53)
Ecuación de Colebrook-White modificada
Esta ecuación combina los tres regímenes de flujo, parcialmente turbulento,
transición y totalmente turbulento y se recomienda (Mohitpour et al., 2002) cuando
el sistema está operando en la región de transición entre ambos regímenes, aunque
otros autores tienen opiniones diferentes sobre este tema (Smith, 1990; Gersten et
al., 2000). El factor de transmisión viene dado por una modificación de la ecuación
de Colebrook-White donde la constante 2.51 fue reemplazada por 2.825 (18) para
lograr una mejor concordancia con los datos experimentales en números de
Reynolds más altos (Gersten et al., 2000).
(Ecuación 33)
Y, en consecuencia, la ecuación de flujo de Colebrook-White modificado será,
(Ecuación 54)
Gersten y col. Ecuación
Esta ecuación es el resultado reciente del trabajo de GERG (Groupe Europeen de
Recherches Gazières) fundado en 1961 y compuesto por miembros de ocho países
europeos. Es una ecuación válida para ambos regímenes, de manera similar a lo
que sucedió con la ecuación de Colebrook-White modificada previamente
considerada, pero según sus defensores, esta ecuación reproduce con alta fidelidad
la transición entre flujo parcial y totalmente turbulento como se muestra en los
resultados experimentales, es decir Un cambio brusco en el factor de transmisión
con el número de Reynolds. El factor de transmisión es dado por Gersten et al.
(2000)
(Ecuación 55)
Mientras que la ecuación de flujo correspondiente es
(Ecuación 56)
El parámetro t es equivalente al factor de eficiencia η presente en la ecuación. (22)
y de manera similar, es igual a la unidad en ausencia de caídas de presión
localizadas. El parámetro n no tiene importancia para los flujos de tubería
hidráulicamente suaves, así como para el flujo totalmente irregular, pero por el
contrario es fundamental para la región de transición entre estos dos flujos.
Para n = 1 hay una transición suave, como lo sugiere la ecuación de Colebrook-
White modificada, pero en el otro extremo, la transición abrupta como se muestra
en los datos experimentales de gas natural, está bien descrita para n = 10.
Comparación de factores de fricción
La Figura 4 tiene una comparación entre varios factores de fricción utilizados en las
ecuaciones mencionadas anteriormente para un flujo totalmente turbulento. Para
llevar a cabo tal comparación, se usó una rugosidad típica de 0.046 mm, así como
tres diámetros de tubería de 1", 15" y 40". Como se puede ver en la ecuación
Panhandle B, aunque solo depende del número de Reynolds sigue de cerca el AGA
para mediano y grande diámetros. También se ve que la ecuación de Weymouth
sobreestima el coeficiente de fricción en comparación con los resultados AGA y esta
diferencia aumenta con la reducción del diámetro de la tubería. En lo que respecta
a la evolución del factor de fricción Colebrook-White modificado, esto se calculó
para dos diámetros de tubería 1" y 40" y, como se esperaba, los resultados de esta
ecuación son similares a los obtenidos con la ecuación AGA para Reynolds alto.

Figura 4 Coeficientes de fricción utilizados en Panhandle B, AGA completamente

turbulento, Colebrook-White y Weymouth.
En la mayoría de la literatura y para ambas ecuaciones, Panhandle A y B, el
exponente de la relación Tst / Pst es diferente de la unidad contraria a las versiones
presentadas aquí. Tal diferencia es comprensible si se comprende que en el
presente texto se ha utilizado una forma general de la ecuación numérica de
Reynolds, la ecuación. (23), mientras que en los desarrollos que aparecen en la
literatura, la ecuación de Reynolds. (24) fue utilizado. En esta última ecuación, la
temperatura y la presión estándar se reemplazaron por 288,15 K y 1 atm, mientras
que esto no se hizo con la misma relación, Tst / Pst, en las ecuaciones de flujo. (15)
o (41). Existe una contradicción en la forma en que se desarrolló la ecuación y esto
resulta en el hecho de que las ecuaciones de Panhandle A y B que se encuentran
en la literatura solo dan resultados correctos si la temperatura y presión estándar
que se utilizarán con ellas son las mismas utilizado previamente en la definición del
número de Reynolds, es decir, 288.15 K y 1 atm.
OTRAS ECUACIONES PARA LA ECUACIÓN DE DISTRIBUCIÓN
IGT de flujo de estado estable
En este caso, el factor de transmisión viene dado por,
(Ecuación 57)
Y usando ahora la ecuación. (41) donde α = 2.3095 y β = 0.100, se obtiene la
siguiente ecuación,
(Ecuación 58)
Donde el factor de compresibilidad promedio Zavg., se supone como 1.
Ecuación de Mueller
Para este caso, el factor de transmisión viene dado por,
(Ecuación 59)
Y siguiendo la misma metodología que en la sección anterior,
(Ecuación 60)
Una vez más, el factor de compresibilidad promedio Zavg., Se supone como 1.
Ecuación de Fritzsche
Ahora el factor de transmisión viene dado por,
(Ecuación 61)
Y reelaborando la ecuación. (41)
(Ecuación 62)
Usando el enfoque de μ = 1.0415 × 10-5 Pa s y suponiendo que zavg es igual a la
unidad, la última ecuación se convierte en
(Ecuación 63)
Ecuación de polos
En esta ecuación, así como en las siguientes, y de manera similar a la Ecuación de
Weymouth, el factor de transmisión solo depende del diámetro interior de la tubería.
Para la ecuación de Polo, la transmisión factor tomará los valores mostrados en la
Tabla 2 (Mohitpour et al.,2000).
1
Tabla 2. Para la ecuación de Polo.
√2

Como la ecuación de Polo es válida para flujos de baja presión, para obtener dicha
ecuación, el factor de transmisión de acuerdo con lo que se presenta en la Tabla 2
se introduce en la ecuación (39). En la literatura técnica, esta ecuación de Polo se
presenta descuidando el término de energía potencial E, el factor de compresibilidad
se supone uno, Tst = 273.15 K, Tavg = 277.8 K y Pst = 1.015598 × 105 Pa. En tales
condiciones, la ecuación. (39) se convierte,
(Ecuación 64)
Donde C se da en la Tabla 3.
Para la media de P'avg, véase la ecuación (39), se adoptó un valor de media de P'avg
= Pst + 390 [Pa] para que los valores de C fueran idénticos a los encontrados en la
literatura técnica de Mohitpour et al. (2000) para la ecuación de Polo.
Tabla 3.C constante en la ecuación de polo
Ecuación de Spitzglass (presión media)
En esta ecuación, el factor de transmisión viene dado por,
(Ecuación 65)
Donde D [m] es el diámetro interior de la tubería. Reemplazando f en la ecuación.
(35) da,
(Ecuación 66)
Sin embargo, la expresión que generalmente se encuentra en la literatura se
simplifica más ya que la energía potencial no se tiene en cuenta, el factor de
compresibilidad se supone igual a uno, Tst = Tavg = 288.8 K y Pst = 1.01325 × 105
Pa, dando como resultado la siguiente ecuación para transporte a media presión,
(Ecuación 67)
Ecuación de Spitzglass (baja presión)
En este caso, la ecuación. (65) del factor de transmisión se reemplaza en la
ecuación. (39) para baja presión en lugar de la ecuación. (35) para mayor presión.
Entonces,
(Ecuación 68)
Una vez más y de manera similar a lo que sucedió con la ecuación de Pole, la
ecuación de Spitzglass de baja presión generalmente se presenta en la literatura
técnica en una forma más simple porque la energía potencial no se tiene en cuenta,
el factor de compresibilidad se supone nuevamente igual a uno, también T st = Tavg
= 288.15 K y Pst = 1.01325 × 105 Pa.
Finalmente, P′avg = Pst + 1210 [Pa] y, al usar todos estos enfoques, la ecuación de
Spizglass para flujos de baja presión se convierte en,
(Ecuación 69)
Ecuación Renouard (presión media)
Las ecuaciones de Renouard se usan con frecuencia en Portugal y España (Brucart,
1987; Andrés et al., 1989; Becco, 1989) para el dimensionamiento de líneas de gas.
En la situación actual, el factor de transmisión según Brucart (1987) viene dado por,
(Ecuación 70)
Y reemplazándolo en la ecuación. (41) con α = 2.4112 y β = 0.09 da la siguiente
ecuación,
(Ecuación 71)
Esta última ecuación solo es válida en los casos en que Q [m3/h]/D[mm]<150
(Brucart, 1987), es decir, que Re <4 × 106, de acuerdo con la ecuación. (25) Para
números de Reynolds más altos, el factor de transmisión debe calcularse a través
de,
(Ecuación 72)
1
Que da valores de 9% más altos que los obtenidos a través de la ecuación (70),
√𝑓
en otras palabras, usar el factor de transmisión tal como se define en la ecuación
(70) da como resultado una sobreestimación de las pérdidas por fricción.
Renouard (1952) también sugiere para tuberías interiores, ramas y columnas
ascendentes (Re <104), el uso de la siguiente ecuación,
(Ecuación 73)
Resultando en factores de fricción superiores a los obtenidos por medio de la
ecuación (70).
Usando una viscosidad absoluta de μ = 1.0757 × 10-5 Pa s, descuidando el término
de energía potencial y suponiendo Tst = Tavg = 288.15 K, Pst = 1.01325 × 105 Pa y
que Zavg = 1, Eq. (71) se vuelve a trabajar para dar una ecuación aplicable a flujos
de presión media,
(Ecuación 74)
En la literatura técnica, la ecuación de Renouard viene escrita de forma explícita en
términos de la diferencia de presión y, en tales circunstancias, la ecuación. (74) se
convierte,
(Ecuación 75)
Mientras que en la literatura técnica (Brucart, 1987) aparece como,
(Ecuación 76)
Donde la constante 4810 toma el valor 48600 si se utilizan las siguientes unidades,
P - kg/cm2, L - km, Q - m3/h y D - mm. Cabe destacar que en el desarrollo que
condujo a la ecuación (76), Renouard (1952) usó un valor constante para la
velocidad cinemática del gas de la ciudad, Vcg = 2.2 × 10-5 m2/s, y una temperatura
de flujo de 15 ºC. Esto significa que al fijar el valor de la viscosidad cinemática del
gas, la densidad también se mantiene fija, lo cual es físicamente inexacto cuando
se consideran flujos de gas compresibles a presión media o alta, porque la
viscosidad cinemática de los gases depende en gran medida de la presión. El
número de Reynolds de flujo de gas ahora se define como
(Ecuación 77)
Y cada vez que se usa un gas con una viscosidad cinemática V diferente del gas
urbano, un factor de corrección multiplicador (V/Vcg)0.18 debe aplicarse en la
ecuación. (76), Renouard (1952). El producto entre dicho factor de corrección y la
densidad relativa del gas (V/Vcg)0.18 generalmente se denomina densidad corregida.
Para obtener una comprensión más profunda del razonamiento y las
simplificaciones matemáticas llevadas a cabo por Renouard (1952), es necesario
regresar y reelaborar la ecuación. (22) para el número de Reynolds de flujo de gas.
Tal ecuación se puede escribir así como
(Ecuación 78)
Lo que significa que esta ecuación es idéntica a la ecuación. (77) siempre que la
viscosidad utilizada sea una viscosidad aparente calculada a través de
(Ecuación 79)
El factor de corrección que se aplicará a la ecuación. (76) cuando se usa un gas
diferente del gas de la ciudad utilizado originalmente por Renouard (1952), o porque
se está considerando una temperatura diferente, debe calcularse a través de una
de las siguientes proporciones
(Ecuación 80)
Comparando Eqs. (75) y (76) y teniendo en cuenta lo expuesto anteriormente, el
error obtenido en el cálculo del diferencial de presión mediante el uso de la
ecuación. (76) en lugar de la ecuación. (75) viene dado por
(Ecuación 81)
Y si en la ecuación (75) se supone que la constante, 4088, cambia de acuerdo con
la viscosidad del gas y en la ecuación. (76) se utiliza una densidad de gas corregida
de acuerdo con la ecuación. (80), la ecuación. (81) se convierte,
(Ecuación 82)
Llevando después de la simplificación al siguiente resultado,
(Ecuación 83)
Lo que significa que ambas ecuaciones. (75) y (76) son equivalentes. Sin embargo,
este resultado sugiere que se debe usar la ecuación (75) en lugar de la ecuación.
(76), porque esa ecuación no solo se dedujo de una manera físicamente correcta,
sino que no requiere el uso de la densidad corregida, proporcione en la ecuación.
(71) se aplica la densidad típica del gas que se utiliza. Para una corrección aún
mayor, sería necesario usar una viscosidad calculada para condiciones de flujo
promedio.
El uso de la ecuación (76) sin la corrección de la viscosidad, aunque es bastante
común en la práctica diaria, lleva a una sobreestimación en el cálculo de la caída de
presión, ya que la diferencia de presión P21 – P22 es aproximadamente 6 a 9% mayor
que el valor obtenido a través de la ecuación. (75)
Ecuación Renouard (baja presión)
Para baja presión, la ecuación. (70) todavía se usa para el factor de transmisión,
pero la ecuación de caída de presión se basa en la ecuación. (39) Entonces se
obtiene que,
(Ecuación 84)
Como en la situación anterior, suponiendo que μ = 1.0757 × 10 -5 Pa s, ignorando el
término de energía potencial y considerando que Tst = Tavg = 288.15 K, Pst = 1.01325
× 105 Pa y Zavg = 1 se puede obtener una ecuación más simple
(Ecuación 85)
P′avg se define como P′avg = Pst + 1390 [Pa] (este valor de presión se utiliza para
convertir las ecuaciones de alta a baja presión que se encuentran en la literatura).
Tal tratamiento matemático es equivalente a reemplazar P21 – P22 por (P1 – P2)2P′avg
en las ecuaciones. (74) a (76). El factor de eficiencia utilizado con las ecuaciones
de Renouard es generalmente 0.91, que es equivalente a usar, en la ecuación de
flujo, una longitud equivalente de 1.2 L.
Régimen de flujo laminar
El flujo laminar puede ocurrir en tuberías de gas doméstico, lo que justifica la
deducción de la ecuación de flujo para esta situación. El factor de transmisión viene
dado por,
(Ecuación 86)
1
Donde K es igual a donde el valor de fapp Re viene dado por la ecuación.
√𝑓𝑎𝑝𝑝 𝑅𝑒
(28) Reemplazando la definición del factor de transmisión para flujo laminar en la
ecuación. (41) y resolverlo para
(Ecuación 87)
Despreciando el término de energía potencial E, suponiendo que T avg = Tst = 288.15
K, Zavg = 1 y μ = 1.0757 × 10-5 Pa y adaptando la última ecuación exclusivamente
para cálculos de flujo de baja presión haciendo P′avg = Pst + 390 [Pa], el resultado
es,
(Ecuación 88)
Que es la ecuación de flujo simplificada para el régimen laminar de baja presión.
La Figura 5 presenta una comparación entre varios factores de fricción para
regímenes de flujo parcial y totalmente turbulentos. En las situaciones en que el
diámetro y la rugosidad de la tubería son esenciales para el cálculo de f, se supuso
para fines de cálculo que D = 4" y una tubería de acero con ε = 0.046 mm de
rugosidad (White, 1999).
A partir del análisis de la figura 5, está claro que la ecuación de Muller proporciona
valores f aproximados a los de la ecuación parcialmente turbulenta AGA (AGA p/t),
para números de Reynolds de hasta 4×104, mientras que la ecuación IGT
proporciona valores bastante cercanos a los de AGA p/t precisamente más allá de
Re = 4×104. La combinación de las ecuaciones de Muller e IGT es una buena
alternativa al uso de la ecuación AGA p/t, ya que esas dos son ecuaciones
explícitas, mientras que la ecuación AGA p/t es implícita, por lo que es más difícil
de tratar. Sin perder de vista la figura 5, se observa que los coeficientes de fricción
de las ecuaciones de Renouard son inferiores a los obtenidos de AGA p/t para los
números de Reynold menores de 7.4 ×103, lo que significa que el valor de f se
subestima cuando se utiliza ecuación Renouard. Sin embargo, para los números de
Reynolds superiores a 7,4 × 103, el factor de fricción para la ecuación de Renouard
es superior al valor obtenido con la ecuación de AGA p/t; en este caso hay una
sobreestimación del factor de fricción que coloca al diseñador en el lado seguro
durante el diseño de equipos y tuberías. La ecuación de Fristzsche tiene un
comportamiento general bastante similar a la ecuación de Renouard, aunque da
mayores factores de fricción cuando Re> 7.4 × 103.

Figura 5 Evolución de los factores de fricción utilizados en varias ecuaciones de

flujo discutidas previamente
Las ecuaciones de Polo y Spitzglass, debido a que solo son función del diámetro de
la tubería, deberían, en teoría, ser las más adecuadas para un flujo totalmente
turbulento en el que parecen sobreestimar el valor de f en comparación con la
ecuación AGA totalmente turbulenta. En el flujo parcialmente turbulento, para
números bajos de Reynolds y diámetros de tubería pequeños, las ecuaciones de
Pole y Spitzglass aún pueden sobreestimar el factor de fricción, ya que este factor
aumenta en gran medida con la disminución del diámetro interno de la tubería. De
hecho, para D = 3/4" f es igual a 0.0438 o 0.06579 de acuerdo con las ecuaciones
de Pole o Spitzglass, sin embargo, para diámetros de tubería grandes y Reynolds
pequeños, estas ecuaciones pueden subestimar el valor de f como se ve claramente
en la figura 5. Una vez más, está claro que el soporte físico de las ecuaciones es
débil y, en consecuencia, su uso debe evitarse siempre que sea posible. El uso de
ecuaciones de flujo turbulento para situaciones de flujo laminar impondrá una
sobreestimación del coeficiente de fricción, aunque en la mayoría de las situaciones
prácticas el número mínimo de Reynolds es 1400. Como se puede ver en la Fig. 5
para un conducto largo, cuando el efecto de longitud de entrada se puede descuidar
y, en consecuencia, fapp Re = 64, el factor de fricción obtenido a través de las
ecuaciones de Renouard, Fritzsche e IGT es inferior al uno obtenido a través de
correlaciones adecuadas para el flujo laminar, pero solo para Reynolds menor de
1400. Lo mismo ocurre con las ecuaciones AGA p/t y Muller cuando el número de
Reynolds es inferior a 1000. Sin embargo, para tuberías más cortas, una situación
común en edificios y pequeñas redes de gas, el coeficiente de fricción laminar es
más alto que los obtenidos por las ecuaciones turbulentas. Para un conducto con
una longitud de 80 D, el producto fapp Re está cerca de 94 y de acuerdo con lo que
se observa en la figura 5, en esta situación, el factor de fricción laminar siempre es
superior al valor obtenido por medio de cualquier otra ecuación turbulenta. Sin
embargo, los diámetros mínimos prácticos en instalaciones domésticas están muy
por encima de tamaños sugeridos por estos cálculos laminares y, en consecuencia,
este análisis es puramente académico. La Figura 6 presenta la variación máxima
detectada para el factor de fricción en relación con el valor promedio calculado por
medio de las ecuaciones de flujo laminar parcialmente turbulentas, Renouard,
Muller, IGT y AGA en función del número de Reynolds. Hay una dispersión mínima
de alrededor del 5% para Reynolds más cercana a 1.6 × 104, pero dicha dispersión
aumenta para el número extremo de Reynolds, ya sean muy pequeños o muy
grandes.

Figura 6 Variación máxima del factor de fricción en relación con su valor promedio
cuando se calcula a través de las siguientes ecuaciones: AGA p/t, Renouard,
Muller, IGT y para el régimen de flujo laminar.
De este análisis se puede decir que en el cálculo y diseño de redes de distribución
y / o utilización de gas en régimen parcialmente turbulento, que actualmente
depende en gran medida de los procedimientos numéricos, dicho procedimiento de
cálculo podría llevarse a cabo a través de las siguientes ecuaciones explícitas (para
cálculos de presión): para Reynolds por encima de 2100, White Eq. (47), aunque
esta ecuación seguirá estando implícita en términos de D y Q &, mientras que en el
régimen laminar Re <2100, Eq. (88) es la opción correcta recomendada, aunque
por razones prácticas no es frecuente el uso de diámetros de tubería muy pequeños
(<18 mm), lo que hace que el uso de esta última ecuación no sea muy común. En
todas estas ecuaciones recomendadas, la temperatura de flujo promedio depende
de la elección del diseñador y la viscosidad absoluta del gas se calcula de acuerdo
con dicha temperatura.

Figura 7 Comparación entre las ecuaciones AGA, Colebrook-White y Gersten et al

(2000).
La Figura 7 presenta el rendimiento de varias ecuaciones de flujo en la transición
entre regímenes de flujo parcial y totalmente turbulento para una tubería de 4 "con
0.046 mm de rugosidad absoluta. La mejor ecuación para traducir el
comportamiento experimental es el de Gersten et al. (2000) con el transición abrupta
entre los dos regímenes de flujo. Ambas ecuaciones AGA, usando la ecuación (32)
como criterio de selección, siguen de cerca la ecuación de Gersten et al. (2000) (y
en consecuencia los datos experimentales) y al mismo tiempo son formalmente más
simples de usar En esta Fig. 7 también es evidente la insuficiencia de la ecuación
de Colebrook-White modificada para traducir correctamente el comportamiento
experimental en la región de transición.
Conclusiones
A pesar de la gran cantidad de trabajo realizado sobre el flujo de gas natural, es
bastante útil unir en un solo texto la multiplicidad de ecuaciones que se encuentran
dispersas en la gran cantidad de literatura técnica disponible. Luego se muestra que
la diferencia básica entre las varias ecuaciones presentadas y discutidas, no está
en su aspecto formal, sino en las ecuaciones de coeficientes de transmisión
correspondientes. Sin embargo, hay una ecuación única en este universo que difiere
de todas las demás, no solo por medio de la ecuación del coeficiente de transmisión,
sino también por la forma en que se deduce la ecuación de flujo, la ecuación de
Renouard.
La comparación entre las diferentes ecuaciones de flujo utilizando las ecuaciones
AGA como base de referencia fue fundamental para la definición, para cada
ecuación bajo análisis, del rango de aplicación correspondiente y los errores que
provienen de su uso. Es, en lo que respecta a los autores, la primera vez que todas
estas ecuaciones se presentan juntas en unidades S.I. Las principales conclusiones
de este trabajo son:
1. Es posible, a través de una ecuación muy simple, correlacionar la rugosidad
relativa de la tubería ε/D con el Reynolds crítico correspondiente a la
transición del flujo de gas natural parcialmente turbulento;
2. Según las limitaciones prácticas para el flujo de la tubería de gas natural
(Kennedy, 1993) y teniendo en cuenta la aspereza típica de la tubería, es
posible concluir que para presiones inferiores a 5 bar y para tuberías de cobre
y polietileno, el único régimen de flujo disponible es el parcialmente flujo
turbulento y, en tales circunstancias, la rugosidad de la tubería no es
necesaria para fines de cálculo. Sin embargo, para las tuberías de acero, la
existencia de un flujo totalmente turbulento es común y no es correcto usarla
en los cálculos de ecuaciones de flujo que ignoran la influencia de la
rugosidad interna de la tubería;
3. Las ecuaciones de flujo de baja presión son simplificaciones de la ecuación
de flujo que se obtuvo a través de la integración de la ecuación de la
conservación del momento lineal para fluidos compresibles. El resultado es
una ecuación de alcance limitado, ya que solo se puede aplicar a un rango
de presión limitado, es decir, < 50 mbar;
4. El uso de una fórmula para el cálculo del número de Reynolds donde Tst e Pst
se reemplazan previamente por los valores numéricos correspondientes,
impone que en la deducción de las ecuaciones de flujo el exponente de la
relación (Tst/Pst) sea diferente de la unidad. Por otro lado, si durante los
cálculos realizados con esta ecuación de flujo los valores utilizados para Tst
y Pst son diferentes de los utilizados en el cálculo de Re, habrá un error en
los resultados numéricos;
5. La utilización de las ecuaciones AGA, junto con la ecuación para la
determinación de Recr recomendada en el presente texto, es una buena
 alternativa a la de Gersten et al. (2000), porque las ecuaciones AGA son
formalmente más simples y se traducen igualmente bien el cambio abrupto
encontrado en la transición del flujo de gas natural parcialmente turbulento.
6. Se debe usar la ecuación (75) en lugar de la ecuación. (76) al calcular
tuberías de gas natural a través de la ecuación de Renouard.