hello everyone
I am a beginner with the MFIX code, and currently I am studying its source code. However, I have been unable to locate the part where the fluid temperature is calculated. Additionally, I cannot determine whether the heat released from homogeneous reactions is included in the source term of the fluid energy equation. I would really appreciate it if experienced developers could offer me some guidance.
Fluid temperature is solved in model/solve_energy_eq.f. Yes, the heat of reaction is included in the source term.
Thanks Jeff,Thank you for helping me resolve my confusion and deepening my understanding of the MFIX source code. Currently, I have encountered a problem: I am unable to output the heat generated by the homogeneous reaction. As a result, I cannot confirm whether the heat from the homogeneous reaction has been included in the heat source term. I would like to ask you again: Is it possible to output the heat from the homogeneous reaction? Or is the reason that it cannot be output because the homogeneous reaction heat is calculated implicitly?I hope you can help me clarify this when you have some free time.
We do not support outputting the heat sources from reactions in the current version of MFiX.
If everything is set correctly, the heat of reactions should be included correctly. To confirm it, I suggest turning off all other heat transfers (like convection, conduction et al) and letting the reactions react, then check if the temperatures changes.
@HangZhou I followed the approach you suggested and found that the homogeneous reaction heat indeed has an effect on the fluid temperature! Perhaps my previous statements were not clearly expressed. I do not wish to output the numerical value of the fluid heat source; rather, I want to output the value of the heat released by the homogeneous reactions themselves.
Currently, I have successfully output the values related to the heterogeneous reactions by using des_usr_var in the file calc_thermo_des.f. Therefore, I would like to ask whether it is possible to output the numerical value of the heat from homogeneous reactions by adding some code. If this can be achieved, I would greatly appreciate your guidance on which source file needs to be modified and what specific code should be added.
The source terms for gas phase can be exported in the array ReactionRates(IJK,N) where IJK is the cell index and N is an index between 1 and nrr. nrr is the number of varibales you want to output. It does not need to be reaction rates even with the name of ReactionRates.
@HangZhou
Thank you for taking the time to reply to me. I would like to know if it is exported from the source file solve_energy_eq.f.
Yes. it is in solve_energy_eq.f