[en] Highlights: • A LES-VOF model is conducted to simulate atomization of coaxial swirling jets. • Structure and flow field of coaxial swirling jets are investigated. • Merging process occurs at the nozzle exit and generates additional perturbation. • The Rayleigh mode instability dominates the breakup of ligaments. - Abstract: Spray atomization process of a liquid-liquid coaxial swirl injector in bipropellant thruster has been investigated using volume of fluid (VOF) method coupled with large eddy simulation methodology. With fine grid resolution, detailed flow field of interacted liquid sheet has been captured and analyzed. For coaxial swirling jet, static pressure drop in the region between the liquid sheets makes two liquid sheets to approach each other and merge. A strong pressure, velocity and turbulent fluctuations are calculated near the contact position of two coaxial jets. Simulation results indicate that additional perturbations are generated due to strong radial and axial shear effects between coaxial jets. Observation of droplet formation process reveals that the Rayleigh mode instability dominates the breakup of the ligament. Droplet diameter and distribution have been investigated quantitatively. The mean diameter of the coaxial jets is between that of the inner and the outer jets. Compared with the individual swirling jets, wider size distributions of droplets are produced in the coaxial jets.