Dear developers, hello. I noticed that in the solve_energy_eq.f source code, the heat source term for particle reactions is represented by HOR_S , and the reaction heat calculated through the subroutine is assigned to RXNS_Qs . However, I have not found any part in the source code where RXNS_Qs is assigned to HOR_S , which has left me quite confused. Therefore, I would like to ask whether the reaction heat calculated by the subroutine has been included in the heat source term.
HOR_S is the energy source term for solid in TFM, which is calculated directly in rrates0.f.
However, RXNS_Qs is used for DPM (like DEM, CGDEM, et al), which is calculated in des/calc_rrates_des.f. and added to the total source term in des/calc_thermo_des.f.
@HangZhou Thank you very much for taking the time to reply to me amidst your busy schedule. Your answer has helped me gain a better understanding of the MFIX source code. According to your explanation, HOR_S is the source term in the energy equation calculated using the TFM model. However, I would like to ask where in the DPM model the fluid temperature is calculated. I haven’t found any other source files for solving the energy equation besides solve_energy_eq.f. Therefore, I’m unsure how the homogeneous reaction heat affects the fluid temperature.
I sincerely hope you can once again offer your guidance and insights.
If you are using the stiff solver, everything is updated in chem/stiff_chem_dpm.f90 for DPM models.
If you are not using the stiff solver, the gas energy equation is solved in solve_energy_eq.f. That is, solve_energy_eq.f solves energy equations for both gas and solid in TFM and only for gas in DPM.
@HangZhou Thank you for your reply. I will try it as soon as possible.