ESVAZIAMENTO DE RESERVATÓRIOS ESFÉRICO, CILÍNDRICO E CÔNICO/ DETERMINAÇÃO DE COEFICIENTE DE DESCARGA скачать в хорошем качестве

ESVAZIAMENTO DE RESERVATÓRIOS ESFÉRICO, CILÍNDRICO E CÔNICO/ DETERMINAÇÃO DE COEFICIENTE DE DESCARGA

LIVROS CITADOS NO VÍDEO: https://www.e-papers.com.br/produto/m... https://www.e-papers.com.br/produto/b... 1) DEDUÇÃO DA EQUAÇÃO BASE - 0:47 2) RESERVATÓRIO ESFÉRICO - 10:22 3) RESERVATÓRIO CILÍNDRICO - 24:04 4) RESERVATÓRIO CÔNICO - 36:24 5) DETERMINAÇÃO DO COEFICIENTE DE DESCARGA NO SCILAB A PARTIR DE DADOS EXPERIMENTAIS - 44:17 Acompanhe também conteúdos em:   / professormarcosfpmoreira     / cristianismo_na_industria   https://professormarcosfpm.wixsite.co...   / professormarcosfpmoreira      / @cristianismo_na_indústria   O Prof. Marcos Moreira é professor Associado do Departamento de Engenharia Química da Universidade Estadual do Oeste do Paraná (Unioeste, Toledo-PR, 24°43'29,77''S/53°45'11,84''O). http://buscatextual.cnpq.br/buscatext... É autor dos livros Scilab aplicado à Engenharia Química; " Operações Unitárias da Engenharia Química utilizando o Excel/VBA; Balanços Globais de Massa na Engenharia Química e Matemática Aplicada à Engenharia Química. Para mais informações acesse: https://www.e-papers.com.br/product_a... Para entrar em contato escreva para [email protected] OBS: Onde você encontrar "(MENOR QUE)" nos códigos altere para o símbolo matemático adequado!!! // PROGRAMA PARA RESERVATÓRIO ESFÉRICO clc clear global x_exp y_exp alfa y0 cte // dados de entrada t_exp=[0 0.5341823 2.0862089 4.5770556... 7.9228045 12.033766 16.81333... 22.156425 27.947376 34.056811... 40.336944 46.61389 52.673862... 58.234548 62.868741 65.577883]; t_exp=t_exp'; h_exp=[0.15:-0.01:0]; h_exp=h_exp'; t0=t_exp(1); h0=h_exp(1); d=0.008; D=0.15; R=D/2; A_sai=%pi*d^2/4; g=9.81; cte=-(2*g)^0.5*A_sai/(%pi); // modelo function t_teorico=t_cal(h,C) t_teorico=((2/3)*D*h^1.5-(2/5)*h^2.5-(2/3)*D*h0^1.5+(2/5)*h0^2.5)/(C*cte); endfunction // Realizando a busca valor=1e1000; for C=0.001:0.001:1 F=(t_cal(h_exp,C)-t_exp)'*(t_cal(h_exp,C)-t_exp) if abs(F) (MENOR QUE) valor then coeficiente_de_descarga=C valor=F end end C_otimo=coeficiente_de_descarga; // imprimindo valores correlacao=correl(t_cal(h_exp,C_otimo),t_exp); printf('Coeficiente de descarga = %.2f
',C_otimo); printf('R^2 = %.6f
',correlacao^2); // Gráfico plot(t_cal(h_exp,C_otimo),h_exp,'b-',t_exp,h_exp,'ro'); legend(["calculado", "experimental"],1); xtitle("AJUSTE MODELO","tempo (s)","Nível (m)"); // PROGRAMA PARA RESERVATÓRIO CILÍNDRICO clc clear global t_exp h_exp cte h0 // dados de entrada t_exp=[0 17.395074 34.790147 52.185221 69.580295 86.975368... 104.37044 121.76552 139.16059 156.55566]; t_exp=t_exp'; h_exp=[0.15 0.1176835 0.0987007 0.0827789 0.0685675 0.0555059... 0.0432932 0.031744 0.0207339 0.0101745]; h_exp=h_exp'; h0=h_exp(1); t0=t_exp(1); d=0.008; D=0.15; L=0.3; A_sai=%pi*d^2/4; g=9.81; cte=-(2*g)^0.5*1.5*A_sai/(2*L); // modelo function h_teorico=h_cal(t,C) global h0 cte h_teorico=D-((D-h0)^1.5-cte*C*t)^(2/3); endfunction // Realizando a busca valor=1e100; for C=0.001:0.001:1 F=(h_cal(t_exp,C)-h_exp)'*(h_cal(t_exp,C)-h_exp) if abs(F) (MENOR QUE) valor then coeficiente_de_descarga=C valor=F end end C_otimo=coeficiente_de_descarga; // imprimindo valores correlacao=correl(h_cal(t_exp,C_otimo),h_exp); printf('Coeficiente de descarga = %.2f
',C_otimo); printf('R^2 = %.6f
',correlacao^2); // Gráfico plot(t_exp,h_cal(t_exp,C_otimo),'b-',t_exp,h_exp,'ro'); legend(["calculado", "experimental"],1); xtitle("AJUSTE MODELO","tempo (s)","Nível (m)"); // PROGRAMA PARA RESERVATÓRIO CÔNICO clc clear global x_exp y_exp alfa y0 cte // dados de entrada t_exp=[0 2.4591706 4.9183412 7.3775118... 9.8366824 12.295853 14.755024 17.214194... 19.673365 22.132535 24.591706]; t_exp=t_exp'; h_exp=[0.15 0.1438097 0.1371915... 0.130056 0.122279 0.1136787 0.1039717... 0.0926701 0.0787958 0.0597161 0]; h_exp=h_exp'; t0=t_exp(1); h0=h_exp(1); d=0.008; D=0.15; R=D/2; H=0.15; A_sai=%pi*d^2/4; g=9.81; cte=-(2*g)^0.5*2.5*A_sai/(%pi*(R/H)^2); // modelo function h_teorico=h_cal(t,C) h_teorico=(h0^(5/2)+C*cte*t)^(2/5); endfunction // Realizando a busca valor=1e1000; for C=0.001:0.001:1 F=(h_cal(t_exp,C)-h_exp)'*(h_cal(t_exp,C)-h_exp) if abs(F) (MENOR QUE) valor then coeficiente_de_descarga=C valor=F end end C_otimo=coeficiente_de_descarga; // imprimindo valores correlacao=correl(h_cal(t_exp,C_otimo),h_exp); printf('Coeficiente de descarga = %.2f
',C_otimo); printf('R^2 = %.6f
',correlacao^2); // Gráfico plot(t_exp,h_cal(t_exp,C_otimo),'b-',t_exp,h_exp,'ro'); legend(["calculado", "experimental"],1); xtitle("AJUSTE MODELO","tempo (s)","Nível (m)");