Based on the double shockwave approximation procedure,combining with RGFM and level-set method,local Riemann problem for strong nonlinear equations of state such as JWL equation of state was constructed and then solved to suppress successfully the numerical oscillation caused by high-density ratio and high-pressure ratio across the explosion product-water interface.A fifth order WENO finite difference scheme and the third-order TVD Runger-Kutta method were utilized for spatial discretization and time advance,respectively.A novel enclosed-type MPI-based parallel methodology for RGFM procedure on uniform structured meshes was presented to realize the parallelization of three dimensional RGFM code for air blast and underwater explosion,which had dramatically improved the practical scale of computing model.For air blast cases,the peak overpressure at same location and same proportional distance of gauge points for both TNT and aluminized explosives was monitored and analyzed for revealing the rule of peak pressure propagation and the influence of aluminum powder combustion on peak overpressure of explosion wave.The overall process of three dimensional bubble pulsation generated by underwater explosion of both explosives mentioned above was successfully simulated with high order numerical scheme.The numerical results obtained indicated that the attenuation of explosion wave formed by aluminized explosives is slower than that caused by TNT.The influence of aluminum powder combustion on bubble pulsation was also investigated by comparing TNT with aluminized explosives.