Factor Acéntrico ω · Fisicoquímica · Ingeniería Química

Fisicoquímica · Ingeniería Química · 3er Semestre

Factor Acéntrico de Pitzer

Un parámetro molecular que transforma la forma en que describimos el comportamiento termodinámico de sustancias reales.

⚗️ Termodinámica de Fluidos 📐 Ecuaciones de Estado 🧬 Geometría Molecular 🏭 Diseño de Procesos Kenneth S. Pitzer · 1955

Desplázate

Sección 01

Contexto Histórico y Motivación

📜 El problema con los fluidos reales

Las ecuaciones de estado clásicas —desde van der Waals (1873) hasta Redlich-Kwong (1949)— describían bien los gases monoatómicos y moléculas simétricas, pero fallaban sistemáticamente con hidrocarburos ramificados, alcoholes y compuestos polares de interés industrial.

El problema era estructural: esas ecuaciones asumían que la geometría molecular no importa, cuando en realidad la forma y la polaridad de la molécula alteran profundamente las interacciones intermoleculares y, por ende, la presión de vapor saturada.

👨‍🔬 Kenneth S. Pitzer (1955)

El químico-físico estadounidense Kenneth Sanborn Pitzer, de la Universidad de California Berkeley, analizó datos experimentales de presión de vapor saturada para docenas de compuestos y observó que las curvas de log(P^sat/P_c) en función de la temperatura reducida no coincidían entre sí, excepto para gases nobles y metano.

Pitzer cuantificó esa desviación en un único número adimensional al que llamó factor acéntrico, ω (omega). Publicó su trabajo en el Journal of the American Chemical Society en 1955 y transformó para siempre el diseño de procesos termodinámicos.

💡 Principio de estados correspondientes simple (Leiden, 1881): Establece que todos los fluidos, si se expresan en variables reducidas (T_r = T/T_c, P_r = P/P_c), se comportan de la misma manera. Pitzer demostró que esto solo es cierto para fluidos esféricos (gases nobles). Para el resto, se necesita un tercer parámetro: ω.

Sección 02

Definición Formal

◆ Definición Operacional de Pitzer (1955)

ω ≡ − log₁₀ ( P^sat / P_c ) _{| T_r = 0.7} − 1.000

Donde P^sat es la presión de vapor saturada del fluido puro, P_c es la presión crítica, y la evaluación se realiza estrictamente a temperatura reducida T_r = T/T_c = 0.7.

🔍 ¿Por qué T_r = 0.7?

Pitzer eligió T_r = 0.7 de forma pragmática: es una temperatura dentro del rango líquido para la mayoría de sustancias de interés industrial, lejos de la región crítica (donde los datos son imprecisos) y donde la presión de vapor es suficientemente medible con buena exactitud experimental.

Matemáticamente, a T_r = 0.7 los gases nobles y el metano dan log₁₀(P^sat/P_c) ≈ −1.000, de modo que ω ≈ 0, estableciendo la línea base del principio de estados correspondientes.

📏 Condiciones de referencia

ω = 0 Fluido simple de referencia: Argón (Ar), Kriptón (Kr), Xenón (Xe). Molécula esférica sin momento dipolar.

ω > 0 Todos los fluidos reales de interés ingenieril. Cuanto mayor ω, mayor desviación del comportamiento esférico.

ω < 0 Extremadamente raro. El helio (ω ≈ −0.390) y el hidrógeno (ω ≈ −0.216) lo presentan por efectos cuánticos.

📌 Interpretación física El factor acéntrico mide la no-esfericidad y la no-centralidad de las fuerzas intermoleculares. Un valor de ω más alto indica que la molécula es más alargada, asimétrica o posee mayor polaridad permanente (dipolo-dipolo, puentes de hidrógeno), lo que se refleja en una presión de vapor saturada menor a T_r = 0.7 de lo que predice el principio de estados correspondientes simple.

Sección 03

Significado Molecular

⚛️

Gases nobles

Esfera perfecta
sin dipolo

⛽

0.0–0.2

Hidrocarburos simples

Ligera asimetría

🧪

0.2–0.5

Polares moderados

Cetonas, ésteres

💧

>0.5

Alta polaridad

Alcoholes, ácidos

🔬 Escala del Factor Acéntrico — Compuestos representativos

Gradiente de 0 (fluido simple) a >1 (molécula compleja y asociada)

0.00.10.20.3 0.40.50.60.7 0.80.91.0+

Argón

ω = 0.001

N₂

Nitrógeno

ω = 0.037

CH₄

Metano

ω = 0.011

C₂H₆

Etano

ω = 0.099

C₃H₈

Propano

ω = 0.152

CO₂

Dióx. Carbono

ω = 0.225

C₆H₆

Benceno

ω = 0.210

nC₇H₁₆

n-Heptano

ω = 0.349

Acetona

Propanona

ω = 0.307

CH₃OH

Metanol

ω = 0.559

C₂H₅OH

Etanol

ω = 0.645

H₂O

Agua

ω = 0.345

📊 Comparativa visual — Factor ω normalizado (máx = n-C₁₀H₂₂)

Argón (Ar)

0.001

Nitrógeno (N₂)

0.037

Metano (CH₄)

0.011

Etano (C₂H₆)

0.099

Propano (C₃H₈)

0.152

n-Heptano

0.349

Benceno (C₆H₆)

0.210

Metanol

0.559

Etanol

0.645

n-Decano

0.490

Sección 04

Tabla de Valores ω y Propiedades Críticas

Datos tomados de: Smith, Van Ness & Abbott (8ª ed.) y Perry's Chemical Engineers' Handbook (9ª ed.). T_c en K, P_c en bar, V_c en cm³/mol.

Sustancia	Fórmula	T_c (K)	P_c (bar)	V_c (cm³/mol)	Z_c	ω	Categoría
Argón	Ar	150.8	48.72	74.9	0.291	0.001	Noble
Kriptón	Kr	209.4	55.02	91.2	0.288	0.005	Noble
Xenón	Xe	289.7	58.40	118.8	0.286	0.013	Noble
Metano	CH₄	190.6	46.10	98.6	0.286	0.011	Simple
Nitrógeno	N₂	126.2	33.98	89.5	0.290	0.037	Simple
Oxígeno	O₂	154.6	50.46	73.4	0.288	0.022	Simple
Dióx. de carbono	CO₂	304.2	73.83	94.0	0.274	0.225	Simple
Etano	C₂H₆	305.3	48.72	145.5	0.279	0.099	Simple
Propano	C₃H₈	369.8	42.48	200.0	0.276	0.152	Simple
n-Butano	n-C₄H₁₀	425.1	37.96	255.0	0.274	0.200	Simple
n-Pentano	n-C₅H₁₂	469.7	33.70	313.0	0.270	0.252	Cadena
n-Hexano	n-C₆H₁₄	507.6	30.25	368.0	0.266	0.301	Cadena
n-Heptano	n-C₇H₁₆	540.2	27.40	428.0	0.261	0.349	Cadena
n-Decano	n-C₁₀H₂₂	617.7	21.03	617.0	0.247	0.490	Cadena
Benceno	C₆H₆	562.2	48.98	259.0	0.271	0.210	Simple
Tolueno	C₇H₈	591.8	41.06	316.0	0.264	0.262	Simple
Cloroformo	CHCl₃	536.4	53.70	239.0	0.287	0.222	Polar
Acetona	C₃H₆O	508.2	47.01	209.0	0.233	0.307	Polar
Acetato de etilo	C₄H₈O₂	523.3	38.80	286.0	0.255	0.363	Polar
Metanol	CH₃OH	512.6	80.97	117.0	0.220	0.559	Asociado
Etanol	C₂H₅OH	513.9	61.48	167.0	0.241	0.645	Asociado
1-Propanol	C₃H₇OH	536.8	51.75	218.5	0.253	0.629	Asociado
Agua	H₂O	647.1	220.64	55.9	0.229	0.345	Asociado
Ácido acético	CH₃COOH	592.7	57.86	171.0	0.200	0.467	Asociado

⚠️ Nota importante: Los valores de ω dependen de los datos experimentales de presión de vapor usados. Distintas fuentes (Reid, Prausnitz & Poling; Smith, Van Ness & Abbott; NIST) pueden reportar valores ligeramente diferentes. Siempre cita la referencia bibliográfica de tus datos.

Sección 05

Correlaciones de Pitzer

El factor acéntrico por sí solo no tiene utilidad: su valor radica en su uso como tercer parámetro en las correlaciones que extienden el principio de estados correspondientes a fluidos reales.

📐 Correlación generalizada para el Factor de Compresibilidad Z

El principio de estados correspondientes extendido de Pitzer expresa Z como una combinación lineal en ω del comportamiento del fluido simple de referencia y una función de corrección:

◆ Correlación de Pitzer — Ec. de estado generalizada Z = Z⁽⁰⁾(T_r, P_r) + ω · Z⁽¹⁾(T_r, P_r)

Donde Z⁽⁰⁾ es la función para fluidos simples (Ar, Kr, Xe) y Z⁽¹⁾ es la función de desviación. Ambas son funciones tabuladas de T_r y P_r (tablas de Lee-Kesler).

🌡️ Entalpía residual

◆ Entalpía residual generalizada (H^R / RT_c) = (H^R)⁽⁰⁾ + ω · (H^R)⁽¹⁾

Permite calcular calores de vaporización, capacidades caloríficas residuales y diferencias entálpicas entre estados reales.

🔄 Entropía residual

◆ Entropía residual generalizada (S^R / R) = (S^R)⁽⁰⁾ + ω · (S^R)⁽¹⁾

Fundamental en el cálculo de trabajo de compresor, eficiencias isoentrópicas y análisis de ciclos termodinámicos.

💧 Correlación de Pitzer para Presión de Vapor Saturada

El mismo ω se usa para correlacionar la propia presión de vapor en todo el rango de temperatura (no solo en T_r = 0.7):

◆ Ecuación de Pitzer para P^sat log₁₀(P^sat / P_c) = f⁽⁰⁾(T_r) + ω · f⁽¹⁾(T_r)

f⁽⁰⁾ (Lee-Kesler) f⁽⁰⁾ = 5.92714 − 6.09648/T_r − 1.28862 ln(T_r) + 0.169347 T_r⁶

f⁽¹⁾ (Lee-Kesler) f⁽¹⁾ = 15.2518 − 15.6875/T_r − 13.4721 ln(T_r) + 0.43577 T_r⁶

⚙️ Uso en Ecuaciones de Estado Cúbicas (EOS)

Las modernas ecuaciones de estado cúbicas (Soave-Redlich-Kwong, Peng-Robinson) incorporan ω en el parámetro α(T) que corrige la atracción intermolecular:

◆ SRK — Soave (1972) α(T_r) = [1 + m(1 − T_r^0.5)]²
m = 0.480 + 1.574 ω − 0.176 ω²

◆ PR — Peng-Robinson (1976) α(T_r) = [1 + κ(1 − T_r^0.5)]²
κ = 0.37464 + 1.54226 ω − 0.26992 ω²

✅ La ecuación de Peng-Robinson es la más usada industrialmente (refinería, gas natural, GLP) por su excelente desempeño en densidades de fase líquida y su capacidad de predecir equilibrio líquido-vapor con mayor precisión que SRK para hidrocarburos pesados.

Sección 06

Calculadora Interactiva

🧮

Herramientas de Cálculo — Factor Acéntrico

Introduce la presión de vapor saturada a T_r = 0.7 y la presión crítica.

P^sat a T_r = 0.7 Presión de vapor saturada en bar

P_c — Presión crítica Presión crítica en bar

ω = —

Calcula P^sat/P_c usando la correlación de Lee-Kesler a cualquier T_r.

ω — Factor Acéntrico Adimensional

T_r — Temperatura reducida T/T_c (0.4 – 0.99)

P_c — Presión crítica (bar)

Estima el factor de compresibilidad Z usando la correlación de virial truncada de Pitzer: Z ≈ 1 + B·P/(RT) con B = B⁽⁰⁾ + ω·B⁽¹⁾ (correlación de Pitzer-Curl).

ω — Factor Acéntrico

T_r = T / T_c

P_r = P / P_c

Calcula el parámetro α(T_r) para las ecuaciones de estado SRK y Peng-Robinson.

ω — Factor Acéntrico

T_r — Temperatura reducida

Sección 07

Aplicaciones en Ingeniería Química

🏭 Diseño de Procesos

Los simuladores industriales (Aspen Plus, HYSYS, ProII) usan ω como entrada obligatoria para seleccionar el modelo termodinámico adecuado. Sin ω, no se pueden calcular: entalpías de mezcla, calores de condensación, presiones de burbuja y rocío.

🔧 Compresores y Turbinas

El cálculo de trabajo isoentrópico de compresores de gas natural requiere Z y S^R que dependen de ω. Un error de 0.01 en ω puede traducirse en un error de ~0.5% en la potencia estimada del compresor.

🌡️ Destilación

La presión de vapor calculada con ω determina las volatilidades relativas en columnas de destilación. Es crítico en la separación de fracciones petroleras, bioetanol y refinación de aromáticos.

❄️ Refrigeración

Los refrigerantes modernos (R-134a, R-410A, R-32) son seleccionados y diseñados considerando ω para optimizar el COP (coeficiente de operación) de ciclos de refrigeración y bomba de calor.

⚗️ Supercríticos

La extracción con CO₂ supercrítico (ω = 0.225) requiere modelado preciso en la región cercana al punto crítico. La ecuación PR es preferida por su fiabilidad a T_r próximo a 1.0.

🧬 Biocombustibles

El biodiesel y el bioetanol tienen valores altos de ω (0.6–1.2) por la presencia de grupos éster e hidroxilo. Esto exige modelos más avanzados como UNIFAC o PC-SAFT para mezclas reactivas.

🚨 Limitaciones del factor acéntrico: Para fluidos con fuerte asociación (ácido acético dimeriza en fase vapor) o mezclas con puentes de hidrógeno cruzados (agua-alcoholes), las correlaciones basadas en ω pueden dar errores >15% en propiedades de equilibrio. En esos casos se recomienda CPA EOS, PC-SAFT o SAFT-VR.

Sección 08

Problemas Propuestos

1 Nivel 1 · Básico

Determinación directa del factor acéntrico

Para el propano (C₃H₈), se reporta que la presión de vapor saturada a T = 258.9 K (que corresponde exactamente a T_r = 0.7) es de 6.468 bar. Determina el factor acéntrico ω del propano usando directamente la definición de Pitzer. Compara con el valor tabulado (ω = 0.152) y calcula el error porcentual.

T_c = 369.8 K P_c = 42.48 bar P^sat(T_r=0.7) = 6.468 bar

2 Nivel 2 · Intermedio

Presión de vapor con correlación Lee-Kesler

El n-butano tiene T_c = 425.1 K, P_c = 37.96 bar y ω = 0.200. Calcula la presión de vapor saturada a T = 340 K usando la correlación de Pitzer-Lee-Kesler. Compara el resultado con el valor experimental de 6.90 bar reportado en la literatura.

T_c = 425.1 K P_c = 37.96 bar ω = 0.200 T = 340 K

3 Nivel 2 · Intermedio

Factor de compresibilidad por correlación virial de Pitzer

El dióxido de carbono (CO₂, ω = 0.225, T_c = 304.2 K, P_c = 73.83 bar) se encuentra a 350 K y 50 bar. Usando la correlación del segundo coeficiente del virial de Pitzer-Curl: B·P_c/(RT_c) = B⁽⁰⁾ + ω·B⁽¹⁾, donde B⁽⁰⁾ = 0.083 − 0.422/T_r^1.6 y B⁽¹⁾ = 0.139 − 0.172/T_r^4.2, calcula Z y el volumen molar del gas. ¿Cuánto difiere del gas ideal?

T = 350 K P = 50 bar ω = 0.225 T_c = 304.2 K P_c = 73.83 bar

4 Nivel 3 · Avanzado

Cálculo de α(T) y parámetros EOS para Peng-Robinson

Para el tolueno (ω = 0.262, T_c = 591.8 K, P_c = 41.06 bar), calcula: (a) el parámetro κ de PR, (b) la función α(T_r) a T = 400 K, (c) el parámetro a(T) de la ecuación PR [a = 0.45724 R²T_c²/P_c · α(T)], y (d) el parámetro b [b = 0.07780 RT_c/P_c]. Expresa tus resultados con unidades consistentes (R = 83.14 cm³·bar/(mol·K)).

T_c = 591.8 K P_c = 41.06 bar ω = 0.262 T = 400 K R = 83.14 cm³·bar/(mol·K)

5 Nivel 4 · Desafío

Análisis comparativo de EOS — n-Heptano

El n-heptano (ω = 0.349, T_c = 540.2 K, P_c = 27.40 bar) se comprime desde 1 bar hasta 30 bar a T = 500 K (isotérmico). Calcula: (a) Z usando correlación virial de Pitzer-Curl, (b) α(T) de SRK y de PR, (c) estima el trabajo técnico de compresión isotérmica (por mol) asumiendo comportamiento del virial truncado con Z ≈ cte (promedio de Z inicial y final). Discute qué modelo esperarías que sea más preciso para un hidrocarburo C7 a estas condiciones y por qué.

P₁ = 1 bar P₂ = 30 bar T = 500 K ω = 0.349

Sección 09

Preguntas Frecuentes

¿Por qué el agua tiene ω = 0.345 si es una molécula muy polar?

El factor acéntrico del agua (0.345) parece bajo para una molécula tan polar, pero hay que recordar que ω se mide empíricamente a través de la presión de vapor a T_r = 0.7. El agua tiene una T_c = 647.1 K muy alta, así que T_r = 0.7 corresponde a T = 453 K (180°C), donde la presión de vapor del agua es ~10 bar. La fuerte red de puentes de hidrógeno eleva tanto T_c que el punto de referencia queda en una región donde el comportamiento no se desvía tanto como cabría esperar. Por eso las EOS cúbicas con solo ω fallan para agua en aplicaciones de proceso y se requieren modelos más avanzados.

¿Cómo se determina ω experimentalmente si no tenemos datos a T_r = 0.7 exacto?

En la práctica, rara vez se tiene un punto experimental exactamente a T_r = 0.7. Se usa regresión con la ecuación de Antoine o Clausius-Clapeyron sobre los datos disponibles para interpolar P^sat al punto requerido. Alternativamente, se usa la ecuación de Lee-Kesler en modo inverso: con datos de P^sat a cualquier T_r, se busca el ω que minimiza el error cuadrado medio sobre toda la curva de vapor saturado. Esta segunda estrategia da ω más robusto y confiable.

¿El factor acéntrico cambia con la temperatura o la presión?

No. El factor acéntrico ω es una propiedad del componente puro, definida a condiciones específicas (T_r = 0.7, sobre la curva de vapor saturado). Es una constante de la sustancia, al igual que T_c, P_c y V_c. No cambia con el estado del sistema. Lo que cambia con T y P son las propiedades residuales y el factor Z, que se calculan usando ω como parámetro fijo.

¿Se puede definir un ω para mezclas?

No directamente, pero en la práctica de simulación se calcula un factor acéntrico aparente de mezcla usando reglas de mezcla lineales: ω_mix = Σ(x_i · ω_i). Esta regla solo es válida en el contexto de las correlaciones de Pitzer y para estimaciones rápidas. Las EOS modernas prefieren calcular Z para la mezcla directamente aplicando reglas de mezcla cuadráticas (vdW-1f) a los parámetros a y b de cada componente.

¿Cuándo debo usar PR en lugar de SRK?

La ecuación de Peng-Robinson predice densidades de líquido con menor error que SRK (≈3–5% vs ≈8–15%), especialmente para hidrocarburos C5+. PR es preferida para: gas natural y LNG, refinería (fracciones C5-C20), y cálculos de balance de energía con cambios de fase. SRK sigue siendo útil para gas natural liviano (C1-C3) a presiones bajas y para sistemas muy supercríticos donde la diferencia de densidades es pequeña. Para sistemas polares o asociados, ni SRK ni PR bastan: usa CPA, PC-SAFT o el modelo NRTL-RK.

¿Cómo se relaciona ω con el momento dipolar μ de la molécula?

No existe una relación algebraica directa y universal entre ω y el momento dipolar μ, porque ω captura tanto la no-esfericidad geométrica como la polaridad electrostática. Una molécula lineal no polar (n-decano, ω = 0.490) tiene ω mayor que el agua (ω = 0.345), aunque μ del agua es mucho más grande. El n-decano tiene ω alto por su cadena larga, no por polaridad. Sin embargo, dentro de una familia química homóloga (alcoholes, cetonas), existe correlación positiva entre μ y ω. En la práctica ingenieril se usan correlaciones como la de Halm-Stiel que incorporan μ como corrección a ω para fluidos muy polares.

Sección 10

Glosario de Términos

Factor Acéntrico (ω)Parámetro adimensional propuesto por Pitzer (1955) que cuantifica la desviación del comportamiento de un fluido respecto al de un fluido esférico de referencia (Ar, Kr, Xe), a través de la presión de vapor saturada a T_r = 0.7.

Temperatura crítica (T_c)Temperatura por encima de la cual no puede existir una fase líquida, independientemente de la presión. Punto en que desaparecen las diferencias entre fase líquida y gaseosa.

Presión crítica (P_c)Presión necesaria para licuar un gas a su temperatura crítica. Define, junto con T_c, el punto crítico de la sustancia.

Temperatura reducida (T_r)Razón adimensional T/T_c. Permite comparar sustancias a "estados equivalentes" dentro del principio de estados correspondientes.

Factor de compresibilidad (Z)Relación Z = PV/(nRT) que mide la desviación del comportamiento de un gas real respecto al gas ideal. Z=1 para gas ideal; Z<1 indica atracción dominante; Z>1, repulsión.

Principio de estados correspondientesHipótesis (Leiden, 1881) que establece que todas las sustancias en estados reducidos idénticos (T_r, P_r, V_r) exhiben el mismo comportamiento. Pitzer lo extendió añadiendo ω como tercer parámetro.

Propiedad residual (H^R, S^R)Diferencia entre la propiedad del fluido real y la del gas ideal a las mismas T y P. Cuantifica el efecto de las interacciones intermoleculares sobre las propiedades termodinámicas.

Ecuación SRKEcuación de estado cúbica de Soave-Redlich-Kwong (1972): P = RT/(V−b) − a·α/(V(V+b)). Incorpora ω a través del parámetro m de la función α(T).

Ecuación Peng-Robinson (PR)Ecuación de estado cúbica (1976): P = RT/(V−b) − a·α/[V(V+b)+b(V−b)]. Más precisa que SRK para densidades de líquido. Estándar industrial en simuladores de procesos.

Correlación Lee-KeslerCorrelación modificada de Benedict-Webb-Rubin (1975) que utiliza ω para extender el principio de estados correspondientes a regiones de alta presión y densidad. Proporciona tablas de Z⁽⁰⁾ y Z⁽¹⁾.

Factor acentrico IQ