This study employed a non-equivalent quasi-experimental pre-test/post-test control-group design to study the effect of the PhET simulation intervention on students' engagement, satisfaction, and academic achievement in the learning of direct current electric circuit concepts among Bhutanese students. We analysed the pre- and post-test scores and perceptions of 57 ninth-grade students, divided into experimental group (EG, n = 29) and control group (CG, n = 28), from one high school in Paro District, Bhutan. The EG students were taught with the PhET simulation intervention, while the CG students were taught with the traditional chalk-talk method. The pre- and post-test scores were collected with the Electric Circuits Conceptual Evaluation (ECEE) inventory. Mean, standard deviation, a two-sample t-test, and multiple linear regression (MLR) were computed using R and RStudio. The t-test revealed a statistically significant difference in the mean post-test scores of CG and EG students. MLR analysis further confirmed that this difference was due to the PhET simulation intervention, ruling out the influence of other confounding variables. Additionally, an instrument called the PhET Engagement-Satisfaction Questionnaire was developed to assess EG students' engagement level and satisfaction with the PhET simulation intervention. Confirmatory factor analysis and Cronbach's alpha calculation confirmed its validity and reliability. Data from the PhET Engagement-Satisfaction Questionnaire unveiled significant impact of the PhET simulation intervention on students' engagement level and their overall satisfaction, reinforcing prior research. However, further research with a larger sample size, incorporating lesson observations, interviews, and our measurement tool, is necessary to ascertain whether the findings it yields align with the present study's findings.