Problem in drag

MFiX
#1

I am using the default drag model “Hill Koch” but some error is coming, I could not find to fix it. please suggest to me how can i fix it.
I also reset the Gas Pressure linear solver 100, even it does not work.
Message from set_bc0.f:564
Outflow BC defined: IJK_P_g remaining undefined.
t= 0.000000 Wrote SPx: B; .RES;
t= 0.0000 Dt= 0.5000E-03 NIT= 25 Sm= 21.481 Hl= 0.0000 CPU= 1.s
Mass imbalance%(0,MMAX): 0.1549E-08 -0.9096E-12
Timestep walltime, fluid solver: 3.129 s
MFiX process has crashed
Previous MFiX run is resumable.
Error from drag_gs.f:1277
Error: Solver crash!
The MFiX solver has terminated unexpectedly
Error information:
float overflow in __drag_gs_mod_MOD_drag_koch_hill
at drag_gs.f:1277
MFiX process has stopped

image
image1725×909 102 KB

lhillkoch_2023-10-22T090326.938997.zip (57.7 MB)

#2

I think this is a bug when the solids fraction is above 0.4. Please try the attached. Copy in your project directory and build the custom solver. This will be fixed in the next point release. Thanks for the report.

drag_gs.f (81.9 KB)