The implementation flow of processing incompressible Navier-Stokes equations based on velocity-correction schemes is introduced, and velocity-correction schemes based on Picard iteration is developed by introducing Picard linearization to process the convection term in the velocity governing equation. Compared with the traditional method, the projection method after Picard linearization can be solved with a larger time step, which improves the stability of the solution method, and the convergence accuracy meets the requirements, which confirms the reliability of the solution method.

Two phase detonations of aluminum dust are simulated in a multi-fluid model to study particle energy calculation methods. In previous studies heat capacities of solid particles are constants, while realistic heat capacities change with temperature. In this simulation, effects of realistic heat capacities are studied. Numerical results show that detonation parameters are influenced significantly. The results with realistic capacities are close to experiments, while the results with fixed capacities overestimate pressure and detonation velocity. In detonation initiation, run-up distance is mainly decide by ignition energy, while realistic effect makes the distance shorter than that in the fixed heat capacity case.