Questions about the MFIX source code file stiff_chem_dpm.f90

MFiX
1

Currently, I am using the CFD-DEM method for simulation. By learning from other developers, I learned that by using the code snippet ReactionRates(IJK,N), it is possible to export cell data. I copied the source code file stiff_chem_dpm.f90 into my home directory and added the ReactionRates code to print out the heat of homogeneous reactions in my simulation. However, I was not able to successfully output any numerical values. After further examination of the code, I found that in MFIX’s des_time_march.f and main.f source files, there is no call to the subroutine stiff_chem_rrates_dpm defined in stiff_chem_dpm.f90. Therefore, I believe this may be the reason why the heat of homogeneous reactions is not being calculated.

I’m not sure whether my understanding is correct. If this is indeed the cause of the problem, I would like to ask experienced developers for some advice on how I can successfully output the values of my homogeneous reaction heats.

2

In stiff_chem.f90 at line 625:

    if (dem_solids .or. pic_solids) then
      call stiff_chem_rrates_dpm(a_NEq(1:NEq_d), &
             a_Y(1:ODE_d), m_norm(:), a_YDOT(1:ODE_d))
    else
      call stiff_chem_rrates_tfm(a_NEq(1:NEq_d), &
             a_Y(1:ODE_d), m_norm(:), a_YDOT(1:ODE_d))
    endif
3

@cgw Thank you for your reply. I’ve been quite busy with work lately and only saw your message now. Your response was very helpful to me. Thank you again.