Demo entry 3441744

Estado Estacionario

   

Submitted by anonymous on Dec 29, 2015 at 15:14
Language: Matlab. Code size: 2.3 kB.

%-------------------------------------------------------------------------%
%       Código para generar los calores iniciales del sistema             %
%-------------------------------------------------------------------------%

clear all; clc;

e_i     = 0.0074;           %Espesor del tubo
g_AL    = 209.3/54;         %Condictividad térmica aluminio
R_i     = 0.0266;           %Radio del tubo
R_agua  = 0.045;            %Radio del baño de agua de mar
T_in    = 112.15;           %Temperatura entrada GNL
T_in_a  = 287.15;           %Temperatura entrada agua de mar
q_a     = 1.48;             %Caudal de agua de mar
q       = 0.0464;           %Caudal de GNL in
c_GN    = (2.54*10^3);      %Calor específico GN
c_GNL   = (3.66*10^3);      %Calor específico GNL
c_agua  = 3.902*10^3;       %Calor específico del agua de mar
E       = 200000;           %Energía cambio de fase
J       = 7;                %Número de rodajas
dx      = 6/J;              %Alto de cada rodaja
i       = J*2;              %Para llenar matrices diagonales
ro_agua = 1023;             %Densidad del agua [Kg/m^3]
ro_GNL  = 419.7;            %Densidad del GNL  [Kg/m^3]
ro_GN   = 31.19;            %Densidad del GN   [Kg/m^3]

m_GNL   = ro_GNL*pi*((R_i)^2)*dx;
m_GN    = ro_GN*pi*((R_i)^2)*dx;
m_agua  = ro_agua*pi*(R_agua^2-(R_i+e_i)^2)*dx;

R   = (e_i)/(g_AL*2*pi*R_i*dx);
k   = ((e_i)/(g_AL*2*pi*R_i*dx))*(q)*(c_GN);
k_1 = ((e_i)/(g_AL*2*pi*R_i*dx))*(q)*(c_GNL);
k_a = ((e_i)/(g_AL*2*pi*R_i*dx))*(q_a)*(c_agua);
g   = k_a;
w   = -k_1-1;
y   = -k-1;
z   = -k_a-1;

A = [eye(J)*y eye(J);eye(J) eye(J)*z];
A(1,1) = w;

 for i=1:1:J-1
     A(i+1,i)         = k;
     A(i+J,i+(J+1))   = g;
 end

b = zeros(1,J*2)';
b(1)   = E*q*R - k_1*T_in;
b(J*2) = -k_a*T_in_a;

x = (inv(A))*b;
m = zeros;

 for i=1:1:J
     m(i,1) = x((J+1)-i);
     m(i,2) = x((J*2+1)-i);
 end  
 
 Q_GN = zeros;
 Q_agua = zeros;
 Q_GNL = (c_GNL)*(m_GNL)*(m(J,1));
 Q_GN_ini = zeros;
 Q_agua_ini = zeros;
 
 for i=1:1:J-1
     Q_GN(i+1) = (c_GN)*(m_GN)*(m(J-i,1));
 end
 for i=1:1:J
     Q_agua(i) = (c_agua)*(m_agua)*(m((J+1)-i,2));
 end
 
 Q_GN_ini(1) = Q_GNL;
for i=2:1:J
    Q_GN_ini(i) = Q_GN(i);
end
for i=1:1:J
    Q_agua_ini(i) = Q_agua(i);
end
 M = m-273;

This snippet took 0.01 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).