A passive containment cooling system (PCCS) has been developed as advanced safety feature for innovative power reactor (iPOWER). Passive systems are inherently less stable than active systems and the PCCS encountered the flashing-induced instability previously identified. The objective of this study is to develop stability maps for flashing-induced instability using MARS (Multi-dimensional Analysis of Reactor Safety) code. Firstly, we conducted a series of sensitivity analysis to see the effects of time step size, nodalization, and alternative MARS user options on the onset of flashing-induced instability. The riser nodalization strongly affects the prediction of flashing in a long riser of the PCCS, while time step size and alternative user options do not. Based on the sensitivity analysis, a standard input and an analysis methodology were set up to develop the stability maps of PCCS. We found out that the calculated equilibrium quality at the exit of the riser as a stability boundary above 5 kW/m(2) was approximately 1.2%, which was in good agreement with Furuya's results. However, in case of a very low heat flux condition, the onset of instability occurred at the lower equilibrium quality. In addition, it was confirmed that inlet throttling reduces the unstable region. (C) 2019 Korean Nuclear Society, Published by Elsevier Korea LLC.