Confusion about the source code stiff_chem_dpm.f90

MFiX
1

Dear developers, hello.
I currently have many questions and confusion regarding the source code stiff_chem_dpm.f90. I am currently conducting simulations using the CFD-DEM method. I have successfully output the reaction heat calculated in stiff_chem_dpm.f90, and I found that the computed reaction heat is not assigned to the global variable q_source, nor is it assigned to any other variable used in the energy equation calculation. Moreover, I noticed that q_source does not include the reaction heat calculated in stiff_chem_dpm.f90.

Therefore, I would like to ask all the developers:
Is the reaction heat calculated in stiff_chem_dpm.f90 actually used in the energy equation?
Does it contribute to the simulation’s temperature field?

2

Q_source is used in des_reaction_model.f, which is called only when you are not using stiff solver.
If you are using stiff solver, everything is solved in stiff_chem_dmp.f90. a_YDOT are the sources for the ODEs we solved in stiff solver, including gas density, temperature, species, and solid mass, temperature and species.
Please refer to MFiXFractionalStepMethodImplementation_110623.pdf for detailed information of stiff solver in MFiX.

1 Like
3

Thank you for your reply. Your explanation has given me a further understanding of MFIX, and it is very helpful to me. Thank you again!