Application of a subdivision algorithm for solving nonlinear algebraic systems

Fernanda De Castilhos Corazza, José Vladimir de Oliveira, Marcos Lúcio Corazza
2008 Acta Scientiarum: Technology  
A subdivision algorithm is presented and applied to solving commonly found chemical engineering problems described by nonlinear algebraic systems. For this purpose, a web-based library available in the literature was used as the main source to select a wide class of one-and multidimensional problems, comprising phase and chemical equilibrium, conversion in tubular and continuous stirred tank reactors, material and energy balances, etc. The problems are classified according to the literature as
more » ... ow, intermediate and of high degree of numerical difficulty based on specific characteristics, like discontinuities in the functions, multiple solutions with the occurrence of false and unfeasible roots, and the presence of null derivative values. It is shown that the algorithm is efficient and robust, even for multidimensional problems of high numerical difficulty, allowing to find simultaneously all the feasible roots of nonlinear algebraic systems, naturally excluding false and unfeasible solutions, with a relatively low CPU time. These features make the algorithm an interesting alternative to deal with chemical engineering problems in contrast to some methods currently in the literature. RESUMO. Aplicação de um algoritmo de subdivisão para solução de sistemas de equações algébricas não-lineares. Um algoritmo de subdivisão é apresentado e aplicado à solução de problemas descritos por sistemas de equações algébricas não-lineares comumente encontrados na engenharia química. Uma biblioteca disponível, na literatura, foi utilizada como fonte principal para a seleção dos problemas a serem resolvidos com uma ou várias dimensões, compreendendo problemas de equilíbrio químico e de fases, conversão em reatores tubulares e contínuos, balanços material e energético, entre outros. Os problemas foram classificados pela literatura com grau de dificuldade numérica baixa, intermediária e alta, com base em características específicas como a existência de descontinuidades nas funções, múltiplas soluções com raízes falsas. O algoritmo mostrou-se eficiente e robusto, mesmo para problemas multidimensionais de alta dificuldade numérica, permitindo encontrar simultaneamente todas as raízes corretas (fisicamente possíveis) dos sistemas algébricos não-lineares, naturalmente excluindo soluções falsas com um tempo de CPU relativamente baixo. Estas características fazem deste algoritmo uma alternativa interessante para solucionar os problemas da engenharia química em contraste com alguns métodos atualmente disponíveis na literatura. Palavras-chave: análise numérica, equações não lineares, múltiplas soluções, algoritmo de subdivisão, simulação, processos de engenharia química.
doi:10.4025/actascitechnol.v30i1.3198 fatcat:dnddottj4ff2layav6ounsb23u