Examen Parcial - II Unidad FENOMENOS DE TRANSPORTE
Examen Parcial - II Unidad FENOMENOS DE TRANSPORTE
Examen Parcial - II Unidad FENOMENOS DE TRANSPORTE
SECCIÓN: A
CICLO: IV
UNIDAD: 1
INTEGRENTES:
SEMESTRE: 2023 - II
En la siguiente tabla se presentan algunas propiedades del Helio en fase gas y el correspondiente
número de Prandtl [Cryogenic system, 2nd. Edition, R.F. Barron, Oxford University Press, 1985].
Calcule la conductividad térmica por medio de la ecuación de Chapman-Enskog y el
correspondiente número Prandtl para seis temperaturas en el intervalo de 5 a 300 K. También,
reporte la diferencia porcentual entre los valores que obtuvo y los que se presentan en esta tabla:
Solución:
El problema nos pide calcular la conductividad térmica y el número de Prandtl, esto con la
finalidad de compararlo con los que nos brindan en la tabla (D.1) y calcular la diferencia porcentual
para cada uno de ellos, para esto tendremos en cuenta las siguientes formulas:
√𝑇/𝑀
𝑘 = 1.9891 × 10−4
𝜎 2 Ω𝑘
En esta primera formula se utiliza para calcular la conductividad térmica, tenemos los valores
contantes de la masa molar (4.003 𝑔/𝑚𝑜𝑙) y el valor de sigma que viene a ser de 2.576 𝐴̇, la
temperatura sale de la tabla a la que vamos a trabajar, y la integral de colisión que se viene a
calcular de la siguiente formula:
𝑇
𝑇∗ =
𝜀/𝐾
En este punto para poder calcular la temperatura adimensional tenemos que el valor de 𝜀/𝐾 es de
10.2 𝐾, con todos estos datos podemos calcular la conductividad térmica, de esta forma lo
calculamos de una manera experimental y hacemos la comparación con la que tenemos en la tabla
para poder calcular la diferencia porcentual. Luego de la misma tabla vamos a tomar los valores
del calor especifico y la viscosidad para calcular el número de Prandtl, esto viene dado por la
siguiente formula:
𝐶𝑝 × 𝜇
𝑃𝑟 =
𝑘
Para calcular este valor utilizaremos el valor de la conductividad térmica que calculamos
anteriormente, también debemos tener en cuenta el factor de conversión ya que al resolver la
ecuación tenemos por unidades a 𝑐𝑎𝑙/𝑐𝑚 × 𝑠 × 𝐾 esto debemos convertirlo a 𝑚𝑊/𝑚 × 𝐾, para
esto tenemos el siguiente factor de conversión:
1𝑐𝑎𝑙 𝑚𝑊
= 418 680
𝑐𝑚 × 𝑠 × 𝐾 𝑚×𝐾
Para finalizar usamos el software Matlab para insertar todos los datos y las formulas para que nos
de los resultados y nos estime la diferencia porcentual:
% Constantes
kB = 1.38064852e-23; % Constante de Boltzmann en J/K
NA = 6.02214076e23; % Número de Avogadro
sigma = 2.576; % Diámetro de colisión en Amstrongs
m = 4.003/NA/1000; % Masa molar del helio en Kg
R = 8.31446261815324; % Constante universal de los gases ideales en J/mol*K
Tc = 10.2; % Temperatura crítica del helio en K
% Cálculos
T_star = T/Tc; % Temperatura reducida
Omega_k = 1.16145*T_star.^(-0.14874) + 0.52487*exp(-0.77320*T_star) + 2.16178*exp(-
2.43787*T_star); % Integral de colisión
k = (0.00019891*sqrt(T/4.003))./((sigma^2).*Omega_k); % Conductividad térmica en
W/m*K
K = k*418680; % Conductividad térmica en mW/m*K
Mu = mu*1e-6; % Viscosidad en Pa*s
cp = Cp*1e3; % Cp en J/Kg*K
Z = rho.*R.*T./m/1e5; % Factor de compresibilidad
Npr = mu.*Cp./K; % Número de Prandtl para un gas real
% Diferencia porcentual
diff_k = (K - 25.9)./25.9*100;
diff_Npr = (Npr - 0.684)./0.684*100;
% Resultados
disp('Temperatura reducida, T* = '); disp(T_star);
disp('Integral de Colisión, Omega_k = '); disp(Omega_k);
disp('Conductividad térmica, k (cal/m*s*K) = '); disp(k);
disp('Conductividad térmica, k (mW/mK) = '); disp(K);
disp('Diferencia porcentual en k (%) = '); disp(diff_k);
disp('Número de Prandtl, Npr = '); disp(Npr);
disp('Diferencia porcentual en Npr (%) = '); disp(diff_Npr);
Cryogenic system, 2nd. Edition, R.F. Barron, Oxford University Press, 1985