codigomatlabmecflui

Páginas: 6 (1470 palabras) Publicado: 23 de septiembre de 2015
%FACULTAD DE INGENIERIA UNAM
%MECANICA DE FLUIDOS
%PRESENTAN: CORRELFULLC7+ 1.O

%El propósito de este programa es la caracterización de la fracción pesada
%del componente C7+, como temperatura crítica, presión crítica, factor
%acéntrico, peso molecular, volumen crítico, temperatura de ebullición y
%factor de compresibilidad Z. mediante el uso de diferentes correlaciones,
%se proponen 6:Riazi-daubert; Standing; Rowe; Mogulas; Kesler-Lee;Winn
%Sim-daubert. y finalmente se presenta una gráfica comparándolas.
%El funcionamiento del programa se basa en un menú gráfico en el que se
%encuentran las 7 opciones, se elige alguna y se siguen las instrucciones.

%Limpia operaciones previas
clear all
clc
%%Comienza el programa con la opción 1
mypick = menu('selecciona unaopción','Riazi-Daubert', ' Rowe', ' Standing', 'Magoulas-Tassios','Kesler-Lee','Winn y Sim-Daubert', 'Gráfica Comparativa')
if mypick==1
disp(' Ecuación de Riazi-Daubert ')
fprintf('\n\n');
disp('este metodo se aplica si se conoce el peso molecular del C7+ y la gravedad especifica.')
fprintf('\n\n');
peso = input('Ingresa el peso molecular del C7+: ')gravedad = input('Ingresa la gravedad específica: ')

format long

disp('Calculando la Temperatura Crítica')
Tc = 544.4*((peso)^0.2998)*((gravedad)^1.0555)*exp((-1.3478*10^-4)*(peso)-0.61641*(gravedad) + 0)
disp('Calculando la Presión Crítica')
Pc = ((4.5203)*10^4)*((peso)^-0.8063)*((gravedad)^1.6015)*exp((-1.8078)*(10^-3)*(peso)+(-0.3084*gravedad))
disp('Calculando elVolumen Crítico')
Vc = (((1.206)*10^-2)*(peso^0.20378)*(gravedad^-1.3036)*exp(((-2.657*10^-3)*peso+(0.5287*gravedad))+ (2.60120*( 10^-3))*(peso)*(gravedad)))
disp('Calculando la Temperatura de Ebullición')
Tb = 6.77857*(peso^0.401673)*(gravedad^-1.58262)*exp((3.77409*10^-3)*(peso)+(2.984036*gravedad)+((-4.25288*10^-3)*(peso)*(gravedad)))


%Muestra los Resultados enPantalla
fprintf('\n\n El valor de la Temperatura crítica (Tc) es: %6.6f en ºR \n ', Tc)
fprintf('\n\n El valor de la Presión Crítica (Pc) es: %6.6f en (lb/in^2) abs \n ', Pc)
fprintf('\n\n El valor del Volumen crítico (Vc) es: %6.6f en (ft^3/lb) \n ', Vc)
fprintf('\n\n El valor de la Temperatura de Ebullición (Tb) es: %6.6f en ºR \n ', Tb)fprintf('\n\n');
disp(' Haciendo uso de la ecuación de Edmister para estimar el factor acéntrico')

w = ((3*(log10(Pc/14.70)))/(7*((Tc/Tb)-1)))-1
fprintf('\n El valor del factor acéntrico es: %6.6f \n, ', w)
fprintf('\n gracias por utilizar este programa ')
end %fin de opción 1

%%Comienza la opción número 2.
if mypick==2
disp('ecuación de Rowe ')
fprintf('\n\n')
disp('Este método calcula la Tc, Tb y Pc dado el peso molecular M y gravedad especifica')
fprintf('\n\n');

disp('Utilizando las correlaciones de Rowe')


Rpeso = input('Ingresa el peso molecular del C7+: ')
Rgravedad=input('Ingresa la gravedad específica: ')

format long
disp('Calculando el número de átomosdel C7+');
n = ((Rpeso-2.0)/14)

disp('Calculando el coeficiente a')
a = 2.95597 -( 0.090597*((n)^(2/3)))

disp('Obteniendo la temperatura crítica')
RTc = 1.8*(961-(10^(a)))

disp('Calculando el coeficiente y')
Y = -0.0137726826*a+0.6801481651


disp('Obteniendo la presión crítica')
RPc = ((10)^(4.89165+Y))/RTc

disp('Calculando latemperatura de ebullición Tb')
RTb = 0.0004347*(RTc^2)+ 265

%Muestra los resultados en pantalla
fprintf('\n\n El valor de la Temperatura crítica (Tc) es: %6.6f en ºR \n ', RTc)
fprintf('\n\n El valor de la Presión Crítica (Pc) es: %6.6f en (lb/in^2) abs \n ', RPc)
fprintf('\n\n El valor de la Temperatura de ebullición (Tb) es: %6.6f en ºR \n ', RTb)...
Leer documento completo

Regístrate para leer el documento completo.

Conviértase en miembro formal de Buenas Tareas

INSCRÍBETE - ES GRATIS