---

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