Prticas Modelao Ambiental

Trabalho Computacional Transferncia de Energia e MassaGuillaume Riflet,Dvidas: guillaume.riflet@ist.utl.pt Ext:3430Pgina: https://fenix.ist.utl.pt/homepage/ist1467301Metodologia de resoluoLer o enunciado do problemaConceber uma soluo computacionalPartir o problema e a soluo em problemas mais pequenos (fazer um diagrama)Resolver as peas do problema, uma de cada vez

2Diagrama por blocos genrico dum problema computacionalInputProcessing /Factory/EngineOutputIdentificar o processamento / escrever o(s) algoritmo(s)Identificar as variveis necessrias para o processamentoPensar na inicializao das variveis e na seleco dos algoritmosPensar no input fornecido pelo utilizador para a resoluo do problema e pensar na interface com o utilizadorPensar no output da soluo, na sua apresentao (tabelas e grficos)Diagrama por blocos genrico dum problema computacionalInputProcessing /Factory/EngineOutputProcessing variables 1AlgorithmsVariables initialization 2Algorithm selectionUser given input / User-Interface (UI)Product(s) 3TablesGraphics

Single or continuous input?Single or continuous output?Identificao dos AlgoritmosProcessing /Factory/EngineProcessing variables 1Algorithms

enunciadoexplicitoimplicitoAlgoritmos e Variveis - esquema numrico explcito no tempo

explicito

Processing /Factory/EngineAlgoritmos - esquema numrico implcito no tempo

implcitoABCTiCall Thomas(A,B,C,Ti)Sistema matricial (n equaes a n incgnitas) e tridiagonal. Resoluo por eliminao de Gauss ou algoritmo semelhante. Neste caso usa-se o eficiente algoritmo de Thomas.Processing /Factory/EngineAlgoritmos ciclos no cdigoProcessing /Factory/EngineFor l = 1 to T ... For i = 1 to N ... Next ...NextCiclo no tempoCiclo no espaoVariveis de processamentoProcessing /Factory/EngineN tamanho do array da barra ou do canal i ndice de clula do canal indo de 1 at NCtdt(1 to N) - array da Concentrao ou Temperatura em t + dtCt(1 to N) array Concentrao ou Temperatura em tt ndice de tempo indo de 1 at LT nmero total de iteraes no tempoDif coeficiente de difusoSk coeficiente de decaimentoA(1 to N), B(1 to N), C(1 to N) e Ti(1 to N) arrays de entrada do algoritmo de Thomas. Resultado escrito em Ti.

Os valores em algumas variveis tm que ser deduzidos a partir dos dados de entrada do utilizador e do enunciadoExerccio IProcessing /Factory/EngineRecorrendo s variveis acima definidas, implemente um programa em VBA que resolva o mtodo explcito da equao da difuso duma propriedade num fluido numa barra unidimensional utilizando arrays. Inicialize todas as variveis no cdigo. Resolva at aula seguinte.