A coupled particle-continuum simulation method is presented to compute the multicomponent mixture plume flow from the attitude-control engines of a satellite through the use of the Navier–Stokes solver and direct simulation Monte Carlo method. In this study, the implicit hybrid flux iteration scheme for the axisymmetric compressible Navier–Stokes equations is constructed to solve the inner flowfield of the nozzle. The Navier–Stokes computing technique is implemented considering slip boundary theory for the near-continuum slip flow near the nozzle exit. Direct simulation Monte Carlo methods for the axisymmetric core plume and three-dimensional far-field plume flows of multispecies mixture gas are developed by using a radial weighted factor, by tracking the molecular motion trajectory, and with secondary Cartesian and unstructured cell processing techniques. A multiregion decomposition and hybrid Navier–Stokes/direct simulation Monte Carlo algorithm with two-way and one-way coupling is established to solve the internal and external flow of the thruster, including the near-field, far-field, backflow, and gas-surface contaminated regions. After constructing the coupled Navier–Stokes/DSMC simulation scheme, the present method is applied to solve the plume flowfield and impinging contamination effects from the satellite attitude-control engine and solar array panel; the simulation results correspond well with the experimental measurements from the low-density wind tunnel and the theoretical predictions. To study the contamination effects produced by the multicomponent mixture plume, the current method is employed to simulate a five-component mixture plume flowfield of tens of meters from two representative attitude-control engines installed in different locations of a satellite in orbit. The results of the nonreacting multicomponent flow with equivalent mole fractions can be achieved with a one-component gas simulation if the species is assigned the properties of the mixture. The deposition rate of the plume flows produced by the two teamwork engines can be superimposed on each other. This method can be used to efficiently predict the impingement contamination effects of the gas-fired mixture plume flow on the satellite and solar array panels.