Cálculo do Número Ideal de Placas para Isolamento Térmico de Tubulação - Estudo em MATLAB

clear all clc Determinação de parâmetros pertinentes para os cálculos Variável simbólica do comprimento do tubo syms L Pi varpi31416 Raios r1006 r2007 Coeficientes k112 h150 Temperaturas Tinf1300 celsius Tinf225 celsius Temperaturas Kelvin Tinf1kelvinTinf127315 celsius Tinf2kelvinTinf227315 celsius Material A kA008 epsilonA090 Prandtl Pr07296 Grashof g981 beta322103 T3ideal50 celsius T3idealkelvinT3ideal27315 nu1678105 Radiação sigma567108 Tfilme375 celsius TfilmekelvinTfilme27315 Ar kar002624 Determinação das resistências térmicas já conhecidas R11h12varpir1 R2logr2r12varpik1 Determinação do loop for que vai fazer a iteração para encontrar o número de placas Variáveis auxiliares C0 n0 nplacas0 varloop0 Determinação do número de placas while varloop0 for i111000000 r3r2i0002 rayleigh8PrgbetaT3idealkelvinTinf2kelvinr33 nu2 if rayleigh 1010 rayleigh 102 C0675 n0058 elseif rayleigh 102 rayleigh 102 C102 n0148 elseif rayleigh 102 rayleigh 104 C0850 n0188 elseif rayleigh 104 rayleigh 107 C0480 n0250 elseif rayleigh 107 rayleigh 1012 C0125 n0333 end hpeliculaCrayleighnkar2r3 R3logr3r22varpikA R41hpelicula2varpir3 R514sigmaepsilonATfilmekelvin32varpir3 ReqR1R2R3R4R5R4R5 dotQTinf1Tinf2Req T3realTinf1dotQR1R2R3 if T3real50 nplacasi break end end varloop1 end fprintfO número ideal de placas do material A é3d nplacas fprintfCom essa idealização de placas encontramos T332f ºC T3real T3T3real T1Tinf1dotQR1 T2Tinf1dotQR1R2 Plots r1plotlinspace0r1 r2plotlinspacer1r2 r3plotlinspacer2r3 T1plotlinspaceTinf1T1 T2plotlinspaceT1T2 T3plotlinspaceT2T3 figure plotr1plot103T1plotyLineWidth15 hold on grid on plotr2plot103T2plotrLineWidth15 plotr3plot103T3plotbLineWidth15 legendRegião Interna ao TuboParede do DutoMaterial A locationbest xlabelRaio mm ylabelTemperatura ºC graphgca graphFontSize12