Neste trabalho estudamos as transições de fases cinéticas e o comportamento crítico da produção de entropia em modelos de spins com interação entre primeiros vizinhos e sujeitos a duas dinâmicas de Glauber, as quais simulam dois banhos térmicos a diferentes temperaturas. Para tanto, é admitido que o sistema corresponde a um processo markoviano contínuo no tempo o qual obedece a uma equação mestra. Dessa forma, o sistema atinge naturalmente estados estacionários, que podem ser de equilíbrio ou de não-equilíbrio. O primeiro corresponde exatamente ao modelo de Ising, que ocorre quando o sistema se encontra em contato com apenas um dos reservatórios. Dessa forma, há uma transição de fase na temperatura de Curie e o balanceamento detalhado é seguramente satisfeito. No segundo caso, os dois banhos térmicos são responsáveis por uma corrente de probabilidade que só existe visto que a reversibilidade microscópica não é mais verificada. Como conseqüência, nesse regime de não-equilíbrio o sistema apresenta uma produção de entropia não nula. Para avaliarmos os diagramas de fase e a produção de entropia utilizamos as aproximações de pares e as simulações de Monte Carlo. Além disso, admitimos que a teoria de escala finita pode ser aplicada no modelo. Esses métodos foram capazes de preverem as transições de fases sofridas pelo sistema. Os expoentes e os pontos críticos foram estimados através dos resultados numéricos. Para a magnetização e a susceptibilidade obtemos = 0,124(1) e = 1,76(1), o que nos permite concluir que o nosso modelo pertence à mesma classe de Ising. Esse resultado refere-se ao princípio da universalidade do ponto crítico, que é verificado devido o nosso modelo apresentar a mesma simetria de inversão que a do modelo de Ising. Além disso, as aproximações de pares também mostraram uma singularidade na derivada da produção de entropia no ponto crítico. E as simulações de Monte Carlo nos permitem sugerir que tal comportamento é uma divergência logarítmica cujo expoente crítico associado vale 1.