Error when running ErgunTestMPI with HerschelBulkley

Hi all,

I have a problem running the tutorial case ErgunTestMPI with the following transport property.

transportModel HerschelBulkley;
nu nu [ 0 2 -1 0 0 0 0 ] 1e-04;
rho rho [ 1 -3 0 0 0 0 0 ] 1000;
HerschelBulkleyCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0] 1;
tau0 tau0 [ 0 2 -2 0 0 0 0 ] 1e-03;
k k [ 0 2 -1 0 0 0 0 ] 1e-06;
n n [ 0 0 0 0 0 0 0 ] 1;
}

At Time = 0.006, the process stucked and then was killed, as shown below.

Time = 0.006

Courant Number mean: 2.25314e+19 max: 3.27342e+21
Starting up LIGGGHTS
Executing command: 'run 100 pre no post no'
Step Atoms KinEng 1 Volume centerOf dragtota dragtota dragtota
10501 10000 0.52049761 0.0069654779 8.8303085e-05 0.011480873 -516.59689 288.15303 -143.29341
10510 10000 141.51224 0.20657257 0.0041844205 0.012145453 -516.59689 288.15303 -143.29341
10520 10000 719.6193 0.22200545 0.19090371 0.012751432 -516.59689 288.15303 -143.29341
10530 10000 1752.9829 0.22805005 2.6168253 0.013214352 -516.59689 288.15303 -143.29341
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 4823 on node Geotech521 exited on signal 9 (Killed).

However, if the yield stress decreased by a factor of 10, as shown below, there was no error.transportModel HerschelBulkley;
nu nu [ 0 2 -1 0 0 0 0 ] 1e-04;
rho rho [ 1 -3 0 0 0 0 0 ] 1000;
HerschelBulkleyCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0] 1e-1;
tau0 tau0 [ 0 2 -2 0 0 0 0 ] 1e-04;
k k [ 0 2 -1 0 0 0 0 ] 1e-07;
n n [ 0 0 0 0 0 0 0 ] 1;
}

Thanks in advance!