Incompatibilidade de resultados a partir de v0.9.46

[Rodolfo Rodrigues]

É um mesmo arquivo (series_reactors.mso) da biblioteca de modelos didáticos do Fogler (chap2) que ao rodar até a versão 0.9.45 reporta resultados corretos. Mas nas versões superiores que a 0.9.45 reporta resultados incorretos (inclusive em ordem de grandeza).

A partir da versão 0.9.46, ele está reportando resultados diferentes para cada nova versão.

Em anexo segue o arquivo (series_reactors.mso) e os resultados gerados para o primeiro flowsheet (sample_cstr) para cada versão a partir da 0.9.45 (foi testado também com a 0.9.44 e não há diferença de resultados com esta):

• sample_cstr_0.9.45.rlt (resultado correto)
• sample_cstr_0.9.46.rlt (resultado um pouco diferente mas da mesma ordem de grandeza)
• sample_cstr_0.9.47.rlt (resultado incorreto inclusive em ordem de grandeza)
• sample_cstr_0.9.48.rlt (resultado incorreto inclusive em ordem de grandeza)

Para comparação usar a variável Vt (volume total) que é a variável principal tratada nos problemas.

[Rafael de Pelegrini Soares]

Na versão 0.9.46 em diante o solver padrão para problemas NLA é o "nlasolver".

O critério de convergência deste solver considera os limites das variáveis: se a variação não é significativa frente os limites e tolerâncias fornecidas então considera que convergiu.

No teu problema é utilizado o tipo volume, este tem limite superior de 10e6 m^{3, certamente inadequado para o problema onde se tem poucos litros. Ou seja, 500 l ou 15 l são iguais com uma tolerância de 1e-3 quando o valor pode ser da ordem de 10e6 m}3.

Soluções possíveis:

• Utilizar o outro solver não-linear: "sundials"
• Arrumar o types.mso para que o limite superior de volume seja menor
• Trocar o solver padrão para o "sundials"
• Solicitar um relativeAccuracy maior neste exemplo

[Rafael de Pelegrini Soares]

Alterado o limite para o tipo "volume", isto resolve o problema deste exemplo.