Floating-point exception in des_reaction_model.f (Line 109)

MFiX Bug report
#1

sum_X_s in des_reaction_model.f computed as zero for des phases with SPECIES_EQ == .false. when mmax > 1 and ANY_SOLIDS_SPECIES_EQ == .true. (at least one phase has SPECIES_EQ == .true.)

Fix: Change Line 110 to:

WHERE(PARTICLE_STATE(:MAX_PIP) == NORMAL_PARTICLE .AND. SPECIES_EQ(PIJK(:MAX_PIP,5))) &

1 Like
#2

Hi Femi, can you please attach your setup files so we can take a closer look?

#3

@jeff.dietiker I have attached the case with the issue.

Before: mfix_with_zero_division_error.zip (375.5 KB)

After: mfix_with_fix.zip (375.2 KB)