به منظور بررسی پارامترهای جریان محترق دوفازی گاز- جامد در چرخه بالستیک داخلی، مدل سازی عددی بر اساس یک مدل تئوری ارائه می شود. این مدل شامل معادلات پیوستگی جرم، اندازه حرکت و انرژی برای هر فاز به همراه معادلات تنش بین دانه ای بین ذرات جامد در بستر منخلخل می باشد. از فرآیند شیمیایی احتراق بستر جامد چشم پوشی می شود اما دبی جرمی گاز حاصل از احتراق بستر به کمک روابط تجربی آهنگ سوزش محاسبه می گردد. هنگامی که احتراق بستر جامد آغاز می گردد، فشار، دما و چگالی گاز به سرعت افزایش می یابد و فشار مخلوط دوفاز را افزایش می دهد. این افزایش فشار در نهایت سبب حرکت مرز متحرک می شود. به منظور بررسی و مطالعه این تغییرات یک روش عددی صریح بر مبنای الگوریتم مک کورمک برای حل دستگاه معادلات حاکم بر مسئله ارائه می شود. از آنجایی که مسئله دارای مرز متحرک است، از شبکه بندی خاصی در این روش عددی استفاده می شود. تغییرات فشار در طول زمان، سرعت مرز متحرک و سایر مشخصه های جریان دو فاز بدست می آیند و در نهایت با سایر مراجع موجود مقایسه می شوند.