There seems to be an issue with the MPPIC method. After the program runs, the particles appear to collapse to the bottom, and their motion is also quite strange.
looks like your IC had only 10% solids volume fraction that collapsed under gravity to form a packed layer… seems fine to me.
Okay, I’ll try adjusting the settings.
I double-checked the case settings. The initial particle volume fraction is set to ic.init-bed.fluid.volfrac = 0.48. I was actually running the example in benchmarks/08-pic/patil, but I had only modified the end time to 5 seconds. That’s why I feel the results from this example are a bit strange.