Abstract: This paper studies the modeling and numerical simulation of the impact jet and impact wall in aircraft deicing spray operations, and analyzes the velocity and pressure distribution characteristics of the impact wall. The phase field function φ is employed to explain the jet deflection phenomenon simulated by the phase field and level set models. The reason for the uneven pressure distribution gradient around the impact jet is explored, and the simulation results of different multiphase flow numerical models are compared based on the Realizable k-ε turbulence model. Our research results show the cross-section of the impact jet develops conically along the central axis, with the flow velocity rapidly decreasing to zero at the stagnation point and the pressure reaching its maximum value. The wall jet gradually gains maximum flow velocity when it deviates from the stagnation point by 0.15～0.2 m. The simulation results of different multiphase flow models are compared with the commonly used turbulence coefficient in experiments. Considering the velocity, pressure distribution characteristics and jet turbulence coefficients simulated by the models, it is concluded that the mixture model and the phase-transfer mixing model are more suitable for the study of aircraft deicing jets and other similar large-flow, multi-fluid microcluster-like flows.