Opened 16 years ago

Closed 16 years ago

#19 closed defect (fixed)

Incompatibilidade de resultados a partir de v0.9.46

Reported by: Rodolfo Rodrigues Owned by: Rafael de Pelegrini Soares
Priority: major Milestone:
Component: EMSO Version:
Keywords: Cc:

Description

É 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.

Attachments (5)

series_reactors.mso (12.4 KB) - added by Rodolfo Rodrigues 16 years ago.
cstr_sample_0.9.45.rlt (242 bytes) - added by Rodolfo Rodrigues 16 years ago.
cstr_sample_0.9.46.rlt (278 bytes) - added by Rodolfo Rodrigues 16 years ago.
cstr_sample_0.9.47.rlt (270 bytes) - added by Rodolfo Rodrigues 16 years ago.
cstr_sample_0.9.48.rlt (276 bytes) - added by Rodolfo Rodrigues 16 years ago.

Download all attachments as: .zip

Change History (9)

Changed 16 years ago by Rodolfo Rodrigues

Attachment: series_reactors.mso added

Changed 16 years ago by Rodolfo Rodrigues

Attachment: cstr_sample_0.9.45.rlt added

Changed 16 years ago by Rodolfo Rodrigues

Attachment: cstr_sample_0.9.46.rlt added

Changed 16 years ago by Rodolfo Rodrigues

Attachment: cstr_sample_0.9.47.rlt added

Changed 16 years ago by Rodolfo Rodrigues

Attachment: cstr_sample_0.9.48.rlt added

comment:1 Changed 16 years ago by Rodolfo Rodrigues

Summary: Compatibilidade de resultados a partir de v0.9.46Incompatibilidade de resultados a partir de v0.9.46

comment:2 Changed 16 years ago by Rodolfo Rodrigues

Component: EMLEMSO
Owner: changed from Paula Bettio Staudt to Rafael de Pelegrini Soares

comment:3 Changed 16 years ago by Rafael de Pelegrini Soares

Status: newassigned

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 m3, 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 m3.

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

comment:4 Changed 16 years ago by Rafael de Pelegrini Soares

Resolution: fixed
Status: assignedclosed

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

Note: See TracTickets for help on using tickets.