نشریه بین المللی مهندسی صنایع و مدیریت تولید

  در این مقاله با استفاده از روش تفاضل مرکزی مرتبه چهارم قطری شده و معادلات ناویر-استوکس لایه-نازک  شوک ( Shock-Fitting ) حول مخروط سرپخ شبیه‌سازی شده‌است. به­علت استفاده از عبارت­های مرتبه بالای بسط سری تیلور در انفصال عددی مشتقات، این روش نسبت به روش­های مرتبه پایین دارای دقت بیشتر و نیز خطاهای عددی (خطای پراکندگی) کمتر می­باشد. چگونگی انفصال عددی مشتقات در روی مرزها و نقاط مجاور آن در پایداری این روش نقش بسزایی دارد. با استفاده از این روش در یک شبکه نسبتاً درشت، می­توان نتایجی بسیار نزدیک به نتایج شبکه ریزِ روش "بیم-وارمینگ" که دارای دقت مرتبه دوم است، به­دست آورد. با ریز شدن شبکه در این روش، دقت بالای آن نسبت به روش "بیم-وارمینگ" محسوس­تر خواهد شد. لازم به ذکر است که سرعت همگرایی این روش نسبت به روش "بیم-وارمینگ" بیشتر می­باشد. این روش قابلیت همگرایی تا دقت ماشین ( Machine Accuracy ) را نیز دارد.