بررسی روش‌های انتگرال‌گیری زمانی در حل عددی معادله‌های دوبعدی آب‌های کم‌عمق

حل عددی معادله‌های دوبعدی جریان نیازمند گسسته‌سازی ترم‌های مکانی و زمانی است و با انتگرال‌گیری عددی از ترم‌های زمانی، به‌روزرسانی معادله‌ها انجام می‌شود. در راستای افزایش دقت و پایداری حل، بیش‌تر از روش‌های با دقت مرتبه دو برای انتگرال‌گیری زمانی ترم‌های معادله‌ها استفاده می‌شود. در این مطالعه دو روش متداول در انتگرال‌گیری زمانی شامل روش rungekuttaمرتبه 3و روش تجزیه عملگر strang، که دارای دقت مرتبه دو هستند موردبررسی قرارگرفته‌اند. برای داشتن قضاوتی صحیح از عملکرد این دو روش نسبت به هم، شیوه‌های عددی کاملاً یکسانی در برخورد با ترم‌های مکانی و ترم منبع‌های معادله به‌کار گرفته شد. سپس با استفاده از دو مدل به‌دست‌آمده که تنها روش برخورد با ترم‌های زمانی در آن‌ها متفاوت است، مسائل یک‌بعدی و دوبعدی دارای نتایج آزمایشگاهی یا حل تحلیلی اجرا شد. نتایج نشان می‌دهد که در مسائل یک‌بعدی با نوسان زیاد، روش رانگ کوتا مرتبه3 در ابتدا خطای کم‌تری داشته، اما با کاهش میزان نوسان، هر دو روش به دقت مشابهی خواهند رسید. در مسائل دوبعدی، پارامتر خطا در روش strangبرای تکرارهای مختلف روند یکنواختی خواهد داشت. اما مقدار خطای روش رانگ کوتا مرتبه3 کم است که با تکرار محاسباتاندکی افزایش می‌یابد.  بااین‌وجود هر دو مدل، مسائل جریان یک و دوبعدی را به‌درستی مدلسازی کرده که به‌خوبی بیانگر چینش صحیح روش‌های به‌کاررفته در آن‌ها در برخورد با ترم‌های مکانی، زمانی و ترم منبع‌‌ها است.

Assessment of Time Integration Methods in the Numerical Solution of Two-Dimensional Shallow Water Equations

The 2D shallow water equations are used in flow simulation of rivers, floodplains, coastal currents, etc. In the research, updating or the socalled numerical integration of temporal terms of twodimensional equations using firstorder methods is more stable but less accurate. In contrast, highorder accuracy methods have numerical stability problems and cause divergence (Brouwer et al., 2014). For this reason, secondorder accurate methods with median properties are widely used. Despite much research on how to deal with spatial terms, according to a review by the authors, there is less research on how to deal with the temporal terms of equations. In addition, studies on time integration methods are limited to solving 1D problems. In this research, two different time integration methods of RungeKutta third order (RK3 method) and Strang splitting operator method (Strang method), which have a secondorder of accuracy and are commonly used in various research (Huang et al., 2013), have been investigated. Therefore, two models have been obtained in which the applied time integration methods are different, but the ways adopted to deal with spatial and sources terms of equations are same. Then, 1D and 2D reference problems are implemented using these two models and their results are presented in order to recognize the appropriate time integration method for solving 2D shallow water equations.