As the core of a hypersonic propulsion system, the effective mixing efficiency of fuel and air in a supersonic combustor is crucial for its performance. This study focuses on a cold supersonic flow and employs computational fluid dynamics (CFD) techniques combined with Euler−Lagrange method's discretephase model (DPM) for multiphase flows, K−H and R−T (Kelvin−Helmholtz and Rayleigh−Taylor) mixing and atomization models, turbulence models, and surface evaporation models to investigate the injection, atomization, and mixing characteristics of kerosene in supersonic airflow. In order to enhance the mixing efficiency between kerosene and air while reducing flow losses, this study examines a staggered dual-jet injection scheme, with the dual jets arranged at the center of the cavity and having a dual-jet spacing of 10 and 20 mm, respectively. Starting from the interaction mechanism between jets, the impact of different staggered dualjet spacings on the kerosene jet penetration height, span expansion area, angle of the shock wave, and Sauter mean diameter distribution was analyzed. The results show that a short dual-jet spacing (10 mm) leads to greater penetration height, wider span expansion, and a larger angle of the shock wave. When the dual-jet spacing is shorter, the interaction between the fuel jet and the cavity shear layer is stronger, resulting in an improved fuel mixing efficiency. The achievements of this study are consistent with previous experimental measurements and the literature, demonstrating a strong theoretical foundation for optimizing the design of hypersonic engines by deepening the understanding of the fundamental atomization mechanisms of kerosene jets in cold-state supersonic flows. Moreover, these results hold practical significance in improving the efficiency of kerosene combustion and enhancing the performance of flame stabilization devices.