Célula a Combustível

import numpy as np # Definição das constantes L = 1.0 # Comprimento da célula a combustível n = 10 # Número de pontos discretizados na direção x dx = L / n # Espaçamento entre os pontos discretizados T_inf = 25 # Temperatura do fluido de entrada T0 = 800 # Temperatura inicial da célula A = 0.1 # Área da superfície da célula Q_dot = 1000 # Taxa de calor gerada pela reação de oxidação do etanol # Cálculo da taxa de transferência de calor convectiva h = 50 # Coeficiente de transferência de calor convectiva T_avg = (T0 + T_inf) / 2 Q_conv = h * A * (T0 - T_inf) # Cálculo da taxa de transferência de calor gerada pela reação química Q_chem = Q_dot * A # Cálculo das temperaturas dos pontos discretizados T = np.zeros(n) T[0] = T_inf T[-1] = T0 for i in range(1, n-1): T[i] = (Q_chem + Q_conv) / (2 * h * A) + (T_inf + T0) / 2 # Exibição dos resultados print("Temperaturas na célula a combustível:") for i in range(n): print(f"Ponto {i+1}: {T[i]:.2f} °C")