A identificação de modelos é determinante no sucesso das modernas técnicas de controle avançado de processos. Um modelo para o sistema pode ser obtido através de modelagem rigorosa, baseada em equações governantes do sistema ou através da modelagem empírica e estimação de parâmetros. Embora mais rápida e fácil, a modelagem empírica necessita de alguns testes de identificação nos quais as variáveis manipuladas são variadas de modo que resultem em variações nas variáveis controladas. Os testes de identificação podem apresentar custos muito elevados tendo em vista que o sistema pode sair de seu ponto normal de operação, gerando produtos com folga de especificação. Este fato ocorre porque usualmente as perturbações aplicadas nas variáveis manipuladas nas indústrias de processos são independentes umas das outras, aumentando a duração do teste de identificação. Desta forma, neste trabalho foi desenvolvida uma nova metodologia de projeto de perturbações simultâneas para a identificação de modelos dinâmicos baseada na direcionalidade do sistema, com o objetivo de fornecer dados mais ricos para se capturar corretamente o comportamento multivariável do sistema e manter o processo no ponto de operação normal. As perturbações são projetadas conforme as características de um modelo simplificado do processo, ou pré-modelo. Este modelo inicial é obtido essencialmente de dados históricos de planta, selecionados através de uma sistemática análise de correlação desenvolvida neste trabalho A metodologia proposta é composta de duas partes: a primeira parte diz respeito à análise dos dados históricos de planta para obtenção de informações prelimirares as quais são utilizadas no planejamento de perturbações, tais como amplitude do ruído de medida, correlação entre as variáveis de processo, constante de tempo do sistema e matriz de ganhos. E a segunda parte consiste no cálculo da amplitude das perturbações baseado nos resultados da primeira etapa do planejamento. Para sistemas mal-condicionados verificou-se que as perturbações planejadas pela metodologia removem menos a planta de seu ponto de operação gerando resultados mais consistentes em relação às perturbações tradicionais. Já para sistemas bem-condicionados, os resultados são semelhantes. A metodologia foi aplicada em uma unidade piloto experimental e numa unidade de destilação da PETROBRAS, cujos resultados apontam pouca remoção dos sistemas do ponto de operação e modelos consistentes. A validação dos modelos também foi contemplada na dissertação, uma vez que foi proposto um novo critério de validação que considera a derivada dos dados de planta e a do modelo e não apenas os dados de planta e os dados da simulação das saídas do modelo.

O principal objetivo deste trabalho é estabelecer caminhos apropriados para a estimativa da eficiência de pratos valvulados com downcomer e pratos dualflow em colunas de destilação industriais. O conhecimento da eficiência tem importância fundamental no projeto e na avaliação do desempenho de colunas de destilação. Pesquisando a literatura, pôde ser identificada uma árvore de alternativas para compor o modelo de eficiência de prato, dependendo dos modelos de transferência de massa, do modelo de distribuição das vazões de líquido e vapor sobre o prato, do modelo de arraste de líquido, do modelo de equilíbrio da mistura multi-componente, dos modelos de propriedades físicas, do modelo da altura da espuma sobre o prato e da definição de eficiência. Neste trabalho, diferentes métodos para a predição da eficiência em pratos valvulados e pratos dualflow foram compostos e comparados contra dados de três colunas de destilação industriais sob diferentes condições operacionais. Os modelos foram inseridos no simulador Aspen Plus 12.1, em linguagem Fortran, junto com dados geométricos dos pratos, propriedades dos fluidos e dados de operação das colunas de destilação industriais. Para cada coluna foi escolhido o melhor pacote termodinâmico pela verificação dos perfis de temperatura e composições de topo e fundo obtidos via simulação contra os valores reais correspondentes das colunas de destilação industriais. Foi feita uma modificação no parâmetro da fração de furos que estão em jetting no modelo hidráulico da dispersão acima do prato de Garcia (1999). Essa modificação produziu resultados melhores na predição da fração de furos em jetting e eficiência dos pratos dualflow e resultados similares aos de Garcia (1999) para a eficiência dos pratos valvulados.

Este trabalho envolve a avaliação da eficiência energética em um sítio industrial, que contém três plantas de produção de polietilenos, utilizando conceitos da análise pinch associados com heurísticas e técnicas de projetos. Na primeira etapa do trabalho, foi realizada uma análise preliminar, na qual foram verificadas as possibilidades de melhorias nas condições de operação das unidades. Com esta avaliação, obteve-se uma economia de 500kg/h de vapor de baixa pressão com a simples alteração na sua pressão de fornecimento para uma das unidades. Na seqüência, com base no conhecimento dos processos e na avaliação das curvas compostas das unidades, verificou-se o potencial de integração energética das plantas individualmente, selecionando-se uma unidade (Planta 3) para ser reprojetada. Para esta unidade foram eliminadas do problema as correntes utilizadas em condições esporádicas, as correntes fora da faixa de integração e, por fim, aquelas cuja integração seria demasiadamente complexa. Após estas simplificações, verificou-se o Ponto de Estrangulamento Energético da unidade, os trocadores que empregam utilidades de uma forma inapropriada, o consumo atual e o mínimo de vapor. Com base nestas informações foram verificadas as modificações estruturais na rede de trocadores de calor da Unidade 3 que resultariam em uma unidade energeticamente mais eficiente. Para isto, utilizou-se o programa Aspen Pinch (da Aspen Tech) que gerou uma lista com vinte possibilidades de adição de novos trocadores. Considerando-se somente as modificações com economia superior a 200kW e as restrições de segurança, layout e controle, apenas duas oportunidades foram detalhadas. A primeira delas, foi o aproveitamento do vapor de flash resultante da expansão do condensado de média pressão da unidade. A segunda oportunidade envolveu a adição de um novo trocador, que foi projetado utilizando-se o programa Bjac conjuntamente com o Aspen Plus. Ambas as oportunidades foram avaliadas economicamente, apresentando Taxas Interna de Retorno muito superiores à Taxa mínima de atratividade de 15% considerada neste trabalho.

Em plantas de polimerização, uma das principais propriedades relacionadas à qualidade do produto é o índice de fluidez, ou MI. Esta propriedade é geralmente medida em laboratório, em intervalos de algumas horas. Por outro lado, as variáveis de processo estão disponíveis continuamente e são controladas com o objetivo de atingir as propriedades desejadas do produto. O desenvolvimento de um analisador virtual, que permita a estimação do MI de forma contínua, com base nestas variáveis medidas, pode levar a um melhor controle da qualidade do produto final. Este trabalho desenvolve um analisador virtual de MI para uma planta de polietileno de alta densidade (PEAD). Existem vários tipos de modelos que podem ser utilizados na implementação de um analisador virtual, assim como várias formas de correção on-line das suas estimativas. Este trabalho direciona especial atenção aos modelos semi-empíricos, e ao uso de um filtro de Kalman estendido (EKF) na correção das estimativas e dos parâmetros do modelo. O objeto de estudo desse trabalho foi um conjunto de produtos, cuja produção é realizada em dois reatores em série. Foram usados 21 conjuntos de dados representando diferentes períodos de operação da planta. Estes conjuntos de dados continham as principais variáveis de planta e os resultados de MI medidos em laboratório. A partir destes dados, foram estimados os parâmetros do modelo para cada reator. Também foi proposta e testada uma metodologia para a sintonia dos parâmetros do EKF. De posse dos modelos, do filtro e de seus respectivos parâmetros, foram gerados os resultados do analisador.

As constantes pressões por redução de margem de lucro, melhoria de qualidade de produtos e segurança operacional que a indústria de processamento vem sofrendo, tem levado as mesmas a utilizarem ferramentas especializadas de otimização de processos. Nas décadas de 1980 e 1990, esta indústria investiu fortemente na utilização de ferramentas de otimização estacionária em tempo real (RTO) nas formas online e offline. Esta tecnologia já atingiu o seu grau de maturidade. Porém está limitada pelas suas características estacionárias, não tendo capacidade de otimizar processos diante das perturbações freqüentes, tais como alterações de qualidade e quantidade de carga, transições decorrentes de alterações de programação de produção ou de receita de uma produção em batelada ou semi-batelada, dentre outros. Para cobrir este espaço, a otimização dinâmica em tempo real (DRTO) é a tecnologia adequada para reduzir a quantidade de produtos fora de especificação e otimizar o lucro operacional diante destas perturbações. Porém, esta tecnologia ainda não atingiu o seu grau de maturidade, tendo ferramentas comerciais apenas na sua versão offline. Com o objetivo de contribuir com a consolidação desta tecnologia, propõe-se estudar e desenvolver uma arquitetura de sistema de DRTO para operar nas plantas de processo, promover melhorias conceituais nesta tecnologia, e desenvolver uma ferramenta de diagnóstico e sintonia de DRTO. Esta estrutura é bem completa e fornece a flexibilidade necessária para uso industrial. A ferramenta de diagnóstico permite resolver os problemas que ocorrerem ao longo do uso do DRTO. Além disso, é apresentada uma nova metodologia de análise e solução de inviabilidades, baseada em otimização multiobjetivos aplicada à otimização dinâmica. Esta técnica pode ser utilizada tanto online quanto offline, para diagnóstico, e na solução de problemas de otimização dinâmica, conferindo mais robustez ao DRTO, evitando insucessos do otimizador. / Considering the constant pressures for profit margins reduction, improvement of products quality and operational safety, that the processing industry has been submitted, has led them to the use of specialized tools of processes optimization. In the 1980 and 1990 decades, this industry has invested strongly in the use of stationary real-time optimization tools (RTO), in the online and offline forms. This technology already reached its degree of maturity, but limited to its stationary characteristics, having no ability to optimize processes due to frequent disturbances, such as quality and feed flows transitions, consequence of frequent changes in the production scheduling or recipes in batch or semibatch operations, and others. To cover this space, the dynamic real-time optimization (DRTO) is the appropriate technology to reduce the off-spec production and optimize the operational profit when the process is submitted to these disturbances. However, this technology has not reached its maturity, having only commercial tools available only in yuor offline version. In order to contribute to the consolidation of this technology, it is proposed to study and develop a DRTO system architecture to operate in process plants, to promote a conceptual improvements in this technology, and to develop a DRTO diagnostic and tuning tool. This structure is quite complete and provides the flexibility required for industrial application. The diagnostic tool allows us to solve problems that occur during the use of the DRTO system. In addition, a new methodology for infeasibility analysis and solution in DAOP is proposed here, and it is based on the solution of the multiobjective dynamic optimization. This technique can be used in online and offline modes, for diagnostics and troubleshooting dynamic optimization problems, giving more robustness to the DRTO systems, avoiding some optimizer failures.

Cada vez mais muitos processos industriais convencionais de separação, tais como destilação, centrifugação e extração com solventes, entre outros, estão sendo substituídos pelos processos de separação por membranas em virtude desses apresentarem vantagens em relação aos processos convencionais, além de serem considerados uma tecnologia limpa. Entre estas vantagens pode-se citar a economia de energia, pois não envolvem mudança de fase; a possibilidade de processamento de substâncias termolábeis, pois podem operar à temperatura ambiente; a maior seletividade; a facilidade de scale – up, pois são modulares e não extensivos em mão de obra, entre outros. Uma característica limitante destes processos é a existência de uma camada de elevada concentração de partículas ou moléculas que se forma próxima à superfície da membrana e que ocasiona uma redução no fluxo permeado, podendo modificar as características de retenção da membrana. Este fenômeno é conhecido como polarização por concentração. O entendimento dos mecanismos de polarização é considerado de grande importância neste campo e o estudo e validação de modelos matemáticos preditivos torna-se necessário para se obter sucesso na aplicação dos processos com membranas no meio industrial O objetivo geral deste trabalho consiste na análise teórico-experimental do processo de microfiltração tangencial com suspensão aquosa de partículas de sílica. O objetivo específico é a utilização de dados experimentais para a estimação de um parâmetro do modelo matemático que descreve o perfil de concentração na superfície da membrana ( membrana). A configuração tangencial implica que a alimentação, contendo as partículas a serem removidas ou concentradas, flui paralelamente à superfície da membrana reduzindo, assim, o fenômeno de polarização por concentração. Os experimentos foram efetuados em uma unidade de bancada de microfiltração tangencial, também conhecido como fluxo cruzado, em um módulo para membrana plana com canal de escoamento retangular sob condições de escoamento laminar. Os seguintes parâmetros foram variados durante os experimentos: pressão transmembrana, velocidade tangencial e concentração de partículas de sílica Os dados experimentais foram utilizados em uma modelagem matemática, baseada nas equações fundamentais da mecânica dos fluidos e da transferência de massa, que considera as propriedades físicas da suspensão, viscosidade dinâmica e difusividade mássica, variáveis com a concentração de sílica para se estimar a concentração máxima de sílica junto à superfície da membrana. Esta modelagem engloba os parâmetros operacionais, bem como as características da membrana nesta única variável. Com a finalidade de se verificar como a concentração máxima sobre a membrana é afetada pelas condições operacionais e, conseqüentemente, como esta afeta o fluxo permeado, a concentração máxima sobre a superfície da membrana é apresentada em função das condições operacionais variadas durante os experimentos de microfiltração. Finalmente, após a estimação desta variável, pôde-se concluir que ela é influenciada pelas condições operacionais típicas do processo de microfiltração tangencial, tais como velocidade tangencial de escoamento, pressão através da membrana e concentração da alimentação.

A consolidação da automação na indústria química trouxe um desafio: transformar o grande volume de dados provenientes dos processos em informação útil. Com o advento dos sistemas digitais, a expansão no desenvolvimento de sensores e métodos de amostragem possibilitou a aquisição de uma elevada quantidade de dados advindos do campo industrial, de forma que, atualmente, engenheiros e operadores dificilmente carecem de informações sobre as diferentes variáveis de um processo. Os algoritmos para compressão de dados de planta surgiram como alternativa para reduzir o espaço de armazenamento demandado por essas informações. Comprimir dados significa gravar apenas uma porção da informação original, devendo-se preservar, no entanto, as características importantes que elas carregam. Ultimamente, não é apenas o espaço em disco que deve ser priorizado quando se fala em compressão. É necessário que os dados sejam fieis às reais informações do processo e, além disso, devem ser recuperados e transmitidos em velocidade razoável. Neste trabalho são apresentadas duas metodologias para a compressão de dados de processos químicos. Para tanto, são utilizadas técnicas para o cálculo da derivada de sinais ruidosos e curvas polinomiais (splines cúbicas) que preservam as características dos dados. Também são propostas sistemáticas para o ajuste automático dos parâmetros nas rotinas de compressão do sistema PI da OSIsoft®. A fim de avaliar essas propostas, foram utilizados alguns estudos de caso, compreendendo sinais reais de uma planta laboratorial e outros artificiais gerados por simulação computacional, abrangendo dinâmicas distintas e características peculiares. Os resultados mostram que os algoritmos são promissores para gravar dados de processos sem que se perca a essência da informação, a qual armazenada pode ser transformada em conhecimento sem prejuízo da qualidade em termos dinâmicos e estáticos. / The consolidation of the automation in the chemical industry has brought a challenge: translating the large volume of data in useful information. With the advent of digital systems, the expansion in the development of sensors and sampling methods allowed the acquisition of a large amount of data from the process field. Thus, nowadays, engineers and operators rarely suffer from lack of information arising from the several process variables. The algorithms for data compression appeared as an alternative to reducing the storage space demanded by these information. Data compressing means to record only a portion of the original information, preserving, however, the relevant features that they hold. Recently, is not only the disk space that must be prioritized when one talks about compression. It is necessary data to be reliable to the actual process information, and, moreover, must be retrived and transmitted at an aceptable speed. In this work, two methodologies for data compression of chemical processes are presented. Therefore, techniques to estimate derivatives of noisy signals and polynomial curves (cubic splines), which preserve the data features, are used. Also, systematic approaches to automate the tuning of the parameters in the OSIsoft® PI System® compression routines are proposed. In order to evaluate these proposals, some case studies are taken, composed by real signals, from a laboratorial plant, and artificial ones generated by computacional simulation, embracing diversified dynamics and peculiar features. The results show that the proposed algorithms are promising to store processes data without significant loss of information, which can be converted in knowledge without jeopardizing the quality in static and dynamic aspects.

O poli(3-hidroxibutirato) é um biopolímero obtido a partir de fontes renováveis e é um poliéster biodegradável. Possui um grande potencial de aplicação principalmente em artefatos e embalagens descartáveis e também na área médica, devido sua biocompatibilidade. Neste trabalho foi estudada a produção de P(3HB) por Bacillus megaterium em meio mineral, buscando o aumento na produção deste biopolímero, bem como representar matematicamente este bioprocesso. Foram avaliadas a influência do pH e da disponibilidade de oxigênio em cultivos utilizando sacarose como fonte de carbono. Também foi avaliada a influência de diferentes fontes de nitrogênio, com glicerol ou sacarose como fonte de carbono. A influência do pH foi estudada mediante experimentos em biorreator de 4 L com controle da concentração de oxigênio dissolvido acima de 40 % em relação a saturação, com ou sem controle de pH. Observou-se que o acúmulo de P(3HB) foi maior quando o pH não foi controlado, 34 % em relação a biomassa em comparação com 28 % quando pH foi controlado. Foi realizada a modelagem da produção de P(3HB), com proposição de quatro modelos. Os dados experimentais dos cultivos com e sem controle de pH foram utilizados na etapa de estimação de parâmetros. Em comparação com modelos relatados na literatura os modelos propostos apresentaram melhor ajuste, com menor ou igual número de parâmetros, e forneceram estimativas dos parâmetros com menor variabilidade. A avaliação da influência do oxigênio foi realizada através de cultivos em biorreator de 4 L, sem controle de pH, com diferentes condições de transferência de oxigênio, obtidas variando-se a rotação do impelidor entre 100 rpm e 600 rpm. O maior acúmulo de biopolímero foi obtido em experimento com rotação do impelidor 200 rpm, o qual apresentou um baixo coeficiente volumétrico de transferência de oxigênio (kLa = 0,006 s-1) e possibilitou acúmulo de 62 % de P(3HB) em relação a biomassa seca. Para avaliar a melhor fonte de nitrogênio, foram realizados experimentos em estufa incubadora rotatória com cinco fontes de nitrogênio diferentes: cloreto de amônio, citrato de amônio, sulfato de amônio, fosfato de amônio e nitrato de amônio. Com relação à fonte de carbono, sacarose possibilitou maior acúmulo de P(3HB) e biomassa que glicerol. Com relação à fonte de nitrogênio, cloreto de amônio possibilitou maior acúmulo de P(3HB) e biomassa quando utilizado glicerol como fonte de carbono, enquanto que sulfato de amônio aumenta a produção com sacarose. / Poly(3-hydroxybutyrate) is a biodegradable polyester produced from renewable resources. This polymer has a great potential to be used as disposable items and packaging, and also in the medical field due to its biocompatibility. In this work, the P(3HB) production by Bacillus megaterium in mineral medium was studied aiming the increase of biopolymer production, as well the development of a mathematical model for this bioprocess. The influence of pH and oxygen transfer rate in a B. megaterium culture using sucrose as carbon source was evaluated, and the influence of different nitrogen sources using sucrose or glycerol as carbon source was analyzed. The influence of pH-medium was studied in bench-scale bioreactor under controlled and uncontrolled pH. The dissolved oxygen was controlled above 40 % of saturation. It was observed that P(3HB) accumulation was higher in the experiment without pH control (34 % of dry biomass) than in the experiment with pH control (28 % of dry biomass). Four models were proposed to perform the mathematical modeling of this bioprocess. The experimental data of the cultures with and without pH control were used for parameter estimation. Compared to earlier models reported in the literature, the proposed models show better fitting, with less or equal number of parameters, and provide parameters estimative with much lower variability. The evaluation of oxygen transfer rate influence was performed in 4 liters bioreactor cultures without pH control, with different conditions of oxygen transfer rates by ranging impeller speed between 100 rpm and 600 rpm. The higher P(3HB) accumulation (62 % of dry biomass) was attained in the experiment with low volumetric oxygen transfer rate (kLa = 0.006 s-1), at 200 rpm impeller speed. To evaluate the best nitrogen source for sucrose or glycerol as carbon sources, experiments were performed in shaker with five different nitrogen sources: ammonium chloride, ammonium citrate, ammonium sulfate, ammonium phosphate and ammonium nitrate. With respect to carbon source, sucrose allowed higher biomass and P(3HB) accumulation than glycerol. With respect to nitrogen source, ammonium chloride allowed highest biomass and P(3HB) accumulation when glycerol was used as carbon source, whereas ammonium sulfate showed the best results for sucrose.

Reações de polimerização em reatores fase gás é o processo mais utilizado atualmente na produção de poliolefinas. Entretanto, existe grande dificuldade de reproduzir um reator em fase gás em escala laboratorial, devido à presença de uma fase sólida e uma fase gasosa em ausência de solvente, a agitação e remoção de calor. O presente trabalho visa uma melhor compreensão dos processos de copolimerização de eteno e buteno, com catalisadores Ziegler-Natta em reatores fase gás e fase líquida a fim de aplicar os resultados obtidos em um reator fase líquida de laboratório em um reator fase gás industrial. A possibilidade de aplicação destes resultados é estudada através da modelagem matemática fenomenológica dos processos. O processo estudado é a copolimerização de eteno e buteno com catalisadores Ziegler-Natta para produção de polietileno linear de baixa densidade em reatores tipo Spherilene. Após ajuste dos modelos para o reator fase líquida e o reator fase gás, efetuaram-se análises de sensibilidade paramétrica para verificar o comportamento das variáveis de saída. A metodologia utilizada para determinação dos critérios de comparação dos dois modelos foi estabelecer um planejamento experimental no processo fase gás, baseado na análise de sensibilidade, utilizando-se as variáveis de saída desse modelo para se obter as condições de operação do reator fase líquida. Essas condições foram relacionadas com as propriedades finais do polímero do reator em fase gás através da temperatura e relações de C4/C2 e H2/C2. A metodologia apresentada permite avaliar as condições adequadas para operar o reator laboratorial em fase líquida e obter dados que podem ser utilizados no estudo do reator fase gás. Os valores obtidos para as constantes cinéticas estão diretamente ligados ao tipo de processo utilizado, utilização ou não de solvente, remoção de calor e outros critérios que diferenciam o processo em fase líquida do processo em fase gás. Logo, aplicando os fatores de correção às constantes cinéticas obtém-se, no reator em fase líquida, um polímero com propriedades finais muito semelhantes ao produto obtido no reator em fase gás de leito fluidizado, considerando as diferenças de processo. / Polymerizations reactions in gas-phase reactors are the utmost process used nowadays for the production of polyolefin. However, it is very difficult to reproduce a gas-phase reactor in a laboratorial scale, due to the presence of solid and gas phases and absence of diluents, the mixing and heat removal. The present work seeks a better understanding of the copolymerization processes of etene/1-butene with Ziegler-Natta catalyst in phase-gas and liquid-phase reactors in order to apply the results obtained in a liquid-phase laboratory reactor in a gas-phase industrial reactor. The possibility of application of these results is studied through the phenomenological mathematical modeling of the processes. The polymerization process studied is the copolymerization of etene/1-butene with Ziegler-Natta catalyst for the production of low density linear polyethylene in a Spherilene reactor. After adjustment in the liquid and gas-phase model, parametric sensitivity analysis were made to evaluate the behavior of the output variables. The methodology used to determinate the comparison criteria was established by experimental design in the gas-phase model, based on the sensitivity analysis results. The output variables of the gas-phase model were used in order to obtain the operating conditions of the liquid-phase reactor. These operating conditions were linked to the final polymer properties through variables such as temperature, C4/C2 and H2/C2 ratios. The presented methodology allows the evaluation of the operating conditions of the laboratory-scale liquid-phase reactor and the collection of data to be used in gas-phase reactors studies. The value obtained to the kinetic constants is directly connected to the kind of process used, presence or absence of diluents, heat removal and other criteria that differentiate both processes. Therefore, the application of correcting factors to the kinetic constants produce, in the laboratory-scale liquid-phase reactor, polymers with final properties very similar to the polymer obtained in the industrial gas-phase reactor, considering all the differences in the process.

O trabalho proposto apresenta um novo método analítico para a resolução da equação de Schrödinger dependente do tempo. Esse método é empregado na simulação de reações químicas envolvendo compostos orgânicos, a fim de estimar o impacto ambiental decorrente do lançamento de substâncias químicas em corpos hídricos. O método é baseado no emprego de transformações de Bäcklund, que fornecem soluções em forma fechada expressas como funções arbitrárias de argumentos conhecidos. Essa formulação proporciona uma redução significativa no tempo de processamento requerido para a obtenção de resultados numéricos, permitindo a simulação de sistemas químicos contendo até 200 átomos. Simulações simbólicas são apresentadas. / This work presents a new analytical method for solving the time dependent Schrödinger equation. This method is employed to simulate chemical reactions involving organic compounds, in order to estimate environmental damages due to the disposal of chemicals in water bodies. The method is based on Bäcklund transformations, which supply closed form solutions as arbitrary functions of known arguments. This formulation provides a significant reduction on the time processing required to obtain numerical results, allowing to carry out the simulation of chemical systems containing up to 200 atoms. Symbolic simulations are reported.