جزییات یک مدل عددی برای شبیه سازی جریان آشفته در درون بنادر رودخانه ای مربعی شکل در این مقاله ارایه می شود. برای مدلسازی جریان تور بلان از معادلات هیدرودینامیکی انتگرال گیری شده در عمق رینولدز استفاده شده است. برای گسسته سازی معادلات حاکم از روش تفاضلات محدود و برای حل آنها از روش نیمه ضمنی ADI و الگوریتم توماس استفاده شده است. در این پژوهش سعی بر این بوده است که با به کارگیری مدل های صفر و دو معادل های توربلانس (طول اختلاط و k-ε) و نیز شرایط مرزی بسته مختلف بتوان تا حد امکان به شکل دقیق الگوی جریان در درون بنادر مزبور دست یافت.