颗粒/流体相互作用

一个标准的 mfix_level::Evolve 步骤如下所示（为了简洁，我们省略了一些细节）：

// 计算颗粒体积分数：
mfix_calc_volume_fraction(lev, sum_vol)

// 演化流体
if ( use_proj_method )
    EvolveFluidProjection(lev, nstep, dt, ...)
else
    EvolveFluidSimple(lev, nstep, dt, ...)

// 将流体阻力施加到颗粒上
mfix_calc_drag_particle(lev)

// 移动颗粒（使用作用在颗粒上的力）
pc -> EvolveParticles(lev, nstep, dt ...)

这里， lev 表示自适应网格细化（AMR）级别， pc 是指向 MFIXParticleContainer 实例的指针。