本发明公开了一种改善激波稳定性的数值模拟方法、装置、设备和存储介质。其中，方法包括将基于欧拉坐标的磁流体方程模型转换为拉格朗日坐标的目标磁流体方程模型；在拉格朗日坐标中，建立物理量穿过激波时需要满足的激波跳跃条件，获得目标磁流体方程模型在接触界面处的速度和磁张力，其中，物理量包括速度和磁张力；将速度和磁张力逆变换为欧拉坐标，在欧拉坐标中，针对激波前后不同的磁场，建立只有对流项磁场的激波跳跃关系，获得相等的对流磁场表达式；采用对流磁场表达式模拟多维磁流体数值。本发明提供的方案避免了数值模拟时在激波波头附近出现非物理凸起的缺陷，提高了数值模拟的可靠性。