Update: I tried using:
eps = C_D(cell, tp);
with Thread *tp as primary liquid phase.
No error occurs.(because "dispersed" option only uses primary phase dissipation rate)
Theory Guide for Prince&Blanch Aggregation Kernel says it is using dissipation rate for liquid phase as well.
Is this approach appropriate?