DEM with reaction : ERROR des\\calc_epg_des.f:205

Hello, everyone,
I want to add DES chemical reaction to a fluidized bed.
I started from copy the “fuild_bed_dem_2d” tutorial case, and imitated the “usr_rates.f” and “usr_rates_des.f” in “silane_pyrolysis_3d_pic” case.

First, I added two gas phase reactions, which works well.
However, when I began to add one DES reaction(C+CO2–>2CO), this erro kept prompting out.
%E5%9B%BE%E7%89%87

When removed reaction, everthing worked well.

I checked the source file following the erro message, and I think it may have something to do with volume fraction.(Maybe not added to one? But shouldn’t be this erro prompted out)

Here is the my project files:
fb_2019-09-27T162730.249344.zip (23.2 KB)

What should I do?

By the way, below is the message showing in the text command line windows before the erro message:

INFO des\calc_epg_des.f:165
Error 1100: Unphysical gas phase volume fraction calculated. A .vtp
file will be written and the code will exit. Fluid cell details:
Fluid Cell IJK: 2851 I/J/K: (13,2,2)
EP_G = 0.000 CUT_CELL_AT = F
PINC: 14 VOL = 0.3000E-06
Global ID Position
70 (0.1102, 0.0008, 0.0008)
71 (0.1118, 0.0008, 0.0008)
72 (0.1134, 0.0008, 0.0008)
73 (0.1150, 0.0008, 0.0008)
74 (0.1166, 0.0008, 0.0008)
75 (0.1182, 0.0008, 0.0008)
76 (0.1198, 0.0008, 0.0008)
215 (0.1102, 0.0025, 0.0008)
216 (0.1118, 0.0025, 0.0008)
217 (0.1134, 0.0025, 0.0008)
218 (0.1150, 0.0025, 0.0008)
219 (0.1166, 0.0025, 0.0008)
220 (0.1182, 0.0025, 0.0008)
221 (0.1198, 0.0025, 0.0008)
This is a fatal error. A particle output file (vtp) will be written
to aid debugging.
ERROR des\calc_epg_des.f:205

This case does not have interpolation enabled by default. Enabling it on the DEM pane could fix your error.

1 Like

Thank you, I will try it as soon as possible.

Still not working.

I tried “Field to particle and particle to field” interpolation and set scheme as “Square DPVM” with width 3e-3.

Hi I got the same problem with you.
I am using the MFIX 21.1.3.
I happened to fix this problem by defining a value for minimum fluid volume fraction in the advanced part of dem panel of solid column.

Regards
Yuanhe

1 Like

thanks, I will try it.