Dracocephalum moldavica is widely used as an ornamental, medicine, and perfume in industry. Real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) is widely and accurately utilized for gene expression evaluations. Selecting optimal reference genes is essential for normalizing RT-qPCR results. However, the identification of suitable reference genes in D. moldavica has not been documented. A total of 12 reference genes in D. moldavica were identified by PEG6000 (15%) treatment under hypertonia conditions in different tissues (roots, stem, leaves, flower, seeds and sepal) and during three stages of flower development, then used to validate the expression stability. There were four algorithms (delta Ct, geNorm, NormFinder, and BestKeeper) used to analyze the stability. Finally, the RefFinder program was employed to evaluate the candidate reference genes’ stability. The results showed that ACTIN, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and EF1α (elongation factor-1α) were stable reference genes under the PEG6000 treatment. Heat shock protein 70 (HSP70) was the most stable gene across different flower development stages. ADP-ribosylation factor (ARF) was the most stable gene in different tissues and total samples. This study provides reliable gene expression studies for future research in D. moldavica.