5 AGRADECIMENTOS À minha família. Vera, Zé, Mama, nada teria sentido sem vocês. Ao meu orientador, Prof. Dr. Reinaldo Giudici. Constantemente ocupado e ainda assim, atencioso, bondoso e paciente. Ao Prof. Dr. Pascal Floquet e Prof. Dr. Xavier Julia do ENSIACET. Aos demais professores, alunos e funcionários do Departamento de Engenharia Química da EPUSP. Àqueles que, circunstancialmente, acompanharam parte da jornada, entre eles: Aos amigos do CRUSP, Luciano, Robson, Paulão, frequentadores do 301E, Cátia, Andrès, Betão, Aninha, Tati, Iara e tantos outros. Ao pessoal do Grupo de terapia do Centro de Saúde Escola Butantã, 'ah, se me vissem agora... bando de malucos!'. Aos não-medicados, Carol, Gian, Alê-Lobisomem, porque a oceanografia também é nossa praia. À minha aluna mais querida, Cristiane: uns saem, outros entram...nessa 'fria'. Jordi. Merci camarade. Em diversos momentos, 'nunca nos amamos tanto'. A todos que diretamente ou não, contribuíram para a realização deste trabalho. Ao financiamento da FAPESP, CAPES e CNPq. Agradeço a todos. 6 DEDICATÓRIA Dedico a vocês meus amores Aurélio, Bruno e Lígia 7 RESUMO Materiais poliméricos são amplamente utilizados e, atualmente, especificações mais rígidas para aplicações especiais tem sido impostas a esses materiais. A microestrutura, em geral, exerce intensa influência na determinação das propriedades macroscópicas dos materiais, sendo assim, o controle da microestrutura e seu relacionamento com essas propriedades é de interesse estratégico. Modelos matemáticos dos processos de polimerização são importantes para relacionar as variáveis de processo com a produtividade e a microestrutura do polímero. Neste trabalho um modelo matemático do processo de copolimerização em emulsão do acrilato de butila com o acetato de vinila foi elaborado. O modelo matemático gerado constitui-se num sistema que engloba equações algébricas (p.ex., de equilíbrio de fases) e diferenciais ordinárias de primeira ordem (p.ex., os balanços de massa dos componentes), num conjunto algébrico diferencial não linear de primeira ordem (DAE). Esse sistema foi solucionado numericamente utilizandose, tanto o método BDF (backward differentiation formulas) quanto as NDF (numerical differentiation formulas), implementado em MATLAB 6.5R13; mais especificamente, obtiveram-se bons resultados utilizando-se a ode113. Os resultados mostram que o modelo representa razoavelmente bem os dados experimentais de trabalhos anteriores do grupo de pesquisa, embora um dos três parâmetros de identificação do modelo tenha sido penalizado com valores abaixo daqueles geralmente reportados na literatura. Apesar disso, o modelo mostrou-se funcional e pode ser útil na simulação do processo. Perfis polinomiais de variações temporais de temperatura também foram aplicados na etapa de otimização simples dos parâmetros das equações desses perfis, objetivando melhorias (redução) no tempo de reação. Observou-se, todavia, que o sistema DAE é de implementação complicada e requer cuidados adicionais na etapa de geração da inicialização consi...