Olá Daniel, tudo bem?

Me chamo Tuiza, sou graduanda em Engenharia Química pela UFRN e estou usando o DWSIM para reproduzir uma simulação de um reator de amônia e usar os dados para construir a pesquisa do meu TCC.

O maior problema da simulação se encontra no fato de está saindo uma quantidade menor de amônia, comparado a quatidade de entrada no reator.

Acredito que o problema esteja na cinética da reação. O autor do artigo que estou tomando como base, usou uma cinética desenvolvida no estudo de DYSON e SIMON (1968): DYSON, D. C., SIMON, J. M., 1968, A kinetic expression with diusion correction for ammonia synthesis on industrial catalyst, I & EC Fundamentals, v. 7, n. 4, pp. 605610. ISSN: 0196-4313.

Na imagem abaixo você pode conferir a equação da taxa. O valor da constante de equilíbrio também estava disponível na literatura da seguinte forma:

logKa = (-2.691122*LOG10(T)-5.519265*10^-5*(T)+1.848863*10^-7*(T^2)+2001.6*(1/T)+2.6899)

Conforme já falei, a reação se trata da síntese de amônia e ocorre na presença do catalisador de ferro. Por isso, usei o modelo cinético catalítico heterogêneo, chegando assim, na seguinte expressão:

NUMERADOR = ((1.7698*10^15)*EXP(-40765/(1.987*T))*((((10^(-2.691122*LOG10(T)-5.519265*10^-5*(T)+1.848863*10^-7*(T^2)+2001.6*(1/T)+2.6899))^2)*((R1^3)^(1-0.5))*R2*((R1^3)^(0.5)))-(((10^(-2.691122*LOG10(T)-5.519265*10^-5*(T)+1.848863*10^-7*(T^2)+2001.6*(1/T)+2.6899))^2)*((P1^2)^(0.5))*R2*((P1^2)^(1-0.5))))

DENOMIDADOR = ((P1^2)^0.5)*((R1^3)^(1-0.5))

Considerei: R1=Hidrogênio; R2 = Nitrogênio e P1 = amônia

Neste caso a expressão está usando o alfa = 0.5, contudo já tentei substituir pelo valor de 0.75, porém os erros permanecem.

Gostaria de saber se você poderia me dar um direcionamento, a fim de chegar na solução deste problema.

Agradeço desde já a sua disponibilidade e atenção!