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

در تحقیق پیش رو جریان یک بعدی در کانال روباز، با استفاده از فرم متوازن شده الگوی عددی شار متوسط وزن‌دار، شبیه سازی شده است که معادلات آب کم عمق را حل می نماید. به منظور پیشگیری از به وجود آمدن نوسانات با منشا غیر فیزیکی، از minmod به عنوان یک تابع محدود کننده شار استفاده شده است. ﺲ از آن، ﺑﺮای دﺳﺘﯿﺎﺑﯽ ﺑﻪ ﯾﮏ ﺷﮑﻞ ﻣﺘﻮازن از اﻟﮕﻮی ﺷﺎر ﻣﺘﻮﺳﻂ وزندار ﺑﺎ ﻫﺪف ﺷﺒﯿﻪ ﺳﺎزی ﺟﺮﯾﺎن آب ﮐﻢ ﻋﻤﻖ، ﯾﮏ ﻓﺮآﯾﻨﺪ ﺑﺎزﺳﺎزی ﺧﻄﯽ روی ﻋﺒﺎرت ﺷﺎر اﻋﻤﺎل ﺷﺪه و ﺑﻪ وﺳﯿﻠﻪ اﻟﮕﻮی ﺑﻪ دﺳﺖ آﻣﺪه در ﺣﺎﻟﺖ ﻫﺎی ﻣﺨﺘﻠﻒ ﺷﺒﯿﻪ ﺳﺎزی ﻫﺎﯾﯽ ﺻﻮرت ﮔﺮﻓﺘﻪ اﺳﺖ. اﯾﻦ ﺷﺒﯿﻪﺳﺎزی ﻫﺎ ﻋﺒﺎرتاﻧﺪ از ﺷﺒﯿﻪﺳﺎزی آﺑﺮاﻫﻪ در ﺣﺎﻟﺖ ﺳﮑﻮن ﺑﻪ ﻣﻨﻈﻮر ﺻﺤﺖﺳﻨﺠﯽ ﺗﻮازن اﻟﮕﻮی ﺷﺎر ﻣﺘﻮﺳﻂ وزن دار و ﻋﺪم واﺑﺴﺘﮕﯽ آن ﺑﻪ ﺷﮑﻞ ﺗﻐﯿﯿﺮﭘﺬﯾﺮیﻫﺎی ﻧﺎﮔﻬﺎﻧﯽ ﮐﻒ، ﺷﺒﯿﻪﺳﺎزی ﺷﮑﺴﺖ ﺳﺪ در دو ﺣﺎﻟﺖ ﺑﺴﺘﺮ ﺧﺸﮏ و ﺑﺴﺘﺮ ﺗﺮ ﮐﻪ ﺑﺎ ﺑﯿﺸﯿﻨﻪ ﺧﻄﺎی ﻧﺴﺒﯽ 3.42 ﻧﺸﺎن دﻫﻨﺪه ﺗﻮاﻧﺎﯾﯽ ﺷﺒﯿﻪﺳﺎزی ﺟﺮﯾﺎن روی ﺑﺴﺘﺮﻫﺎی ﺗﺮ و ﺧﺸﮏ و ﻗﺎﺑﻠﯿﺖ ﺗﺴﺨﯿﺮ ﺷﻮک اﺳﺖ. ﻫﻤﭽﻨﯿﻦ ﺷﺒﯿﻪ ﺳﺎزی ﺟﺮﯾﺎن در آﺑﺮاﻫﻪ ﺑﺎ ﺗﻐﯿﯿﺮﭘﺬﯾﺮیﻫﺎی ﻧﺎﮔﻬﺎﻧﯽ ﮐﻒ ﻧﺸﺎن دﻫﻨﺪه ﻗﺎﺑﻠﯿﺖ اﯾﻦ اﻟﮕﻮ ﺑﺮای ﻣﺪﻟﺴﺎزی ﺟﺮﯾﺎن روی ﺗﻐﯿﯿﺮﭘﺬﯾﺮیﻫﺎی ﻧﺎﮔﻬﺎﻧﯽ ﮐﻒ ﺑﻮده و ﺑﯿﺸﺘﺮﯾﻦ ﻣﻘﺎدﯾﺮ ﺧﻄﺎی ﻧﺴﺒﯽ در ﺷﺒﯿﻪ ﺳﺎزی ﻫﺎی ﺻﻮرت ﮔﺮﻓﺘﻪ ﺑﺮاﺑﺮ 3.47 اﺳﺖ. ﺷﺒﯿﻪﺳﺎزی ﻫﺎی ﺻﻮرت ﮔﺮﻓﺘﻪ و ﻣﻘﺎدﯾﺮ ﺧﻄﺎ ﮔﻮﯾﺎی ﮐﺎراﯾﯽ، ﭘﺎﯾﺪاری و دﻗﺖ اﻟﮕﻮی ﻣﺘﻮازن ﺷﺎر ﻣﺘﻮﺳﻂ وزندار ﺑﻪدﺳﺖ آﻣﺪه در اﯾﻦ ﺗﺤﻘﯿﻖ اﺳﺖ.

Evaluation of well-balanced form of Weighted Average Flux scheme for simulation of flow in open channels

Water is the most strategic liquid in the world. The life of all humankind and animals and plants are relying on the water. The water should supply to the location of demands. One of the most common water transmission ways is open channels. If a sudden change occurs in the channel section, it can affect the whole water flow in the channel. These changes can happen naturally, like aggregation of sediments in a section of the channel. Moreover, the changes may build by humans, like sharp and broadcrested weirs. Thus, it is necessary to simulate openchannel flows to predict possible changes in water surface profile and velocity. Basically, researchers follow three approaches to simulate water flows: the analytic, the experimental, and the numerical approaches. Analytical approach for solving the flow equations is not sufficient due to the complexity and nonlinearity of the equations so there are several restrictions in the modeling. On the other hand, experimental approach is time consuming and expensive. Since the highperformance computers have been developed, researchers attracted to numerical approaches. There are different numerical solutions which are used to solve the water flow equations such as Finite Difference Method, Finite Element Method, Finite Volume Method, etc. The Finite Volume Method is one of the most applicable methods in several computational aspects of engineering, such as computational fluid dynamics and heat transfer problems. In this method, it is necessary to have a strong approximation of numerical flux term for solving flow equations. The Riemann solver provides a reliable approximation for the numerical flux term. The Riemann problem for a set of PDEs is an initial value problem for such PDEs in which the initial condition has a special form. In order to apply numerical solutions, one can use the exact Riemann solver or approximate Riemann solver. The exact Riemann solver uses NewtonRaphson method that takes noticeable cost in time and money and the results rely on the first guess of NewtonRaphson. Therefore, researchers prefer the approximate Riemann solvers such as Harten Lax van Leer (HLL) scheme, Harten Lax van Leer Constant (HLLC) scheme and Weighted Average Flux (WAF) scheme that have acceptable results and running time. WAF scheme can be categorized as a branch of Finite volume method. The scheme was first applied to the Euler equation. This scheme is one of the approximation solution (besides HLL and HLLC methods) of the Riemann problem. Then, Toro used the WAF scheme to simulate twodimensional shallow water equations. Subsequently, WAF has been utilized to simulate flow over different kinds of open channels. Although the scheme shows reasonable results, it is noticeable that the numerical scheme is not wellbalanced essentially. Thus, a wellbalanced WAF scheme should be developed to simulate flow in open channels accurately without nonphysical fluctuations in flow surface. The aim of this research is using the ability of the WAF scheme to simulate shallow water and applying some consideration on the scheme to prevent nonphysical fluctuations in water surfaces. In this paper, a wellbalanced form of WAF which is combined with HLL for estimating flux has been employed to simulate onedimensional flow open channels. MINMOD as an effective slope limiter has been used in order to prevent nonphysical oscillations. Moreover, RungeKutta has been employed as the time integration method to renew depths and velocities. Several different cases have been used to show that the scheme has an excellent shockcapturing ability and can handle the wet and dry condition of channel bed. Importantly, the linear reconstruction for the scheme has been applied to have secondorder accuracy and to prevent the negative depth effect on computations. The scheme is shown to be wellbalanced by evaluating stationary solutions at steady state conditions. Besides, the capability and accuracy of the scheme are verified by the comparison of scheme numerical results with literature analytical and experimental results. The numerical results have shown that the scheme can satisfy the continuity equation and prevent negative depth. For real applications of the scheme, the simulations of flow over sharp changes and dambreak show that RMSEs are in acceptable ranges and there is no nonphysical fluctuations on the water surface profile. Simulating dam break on wet and dry beds, show that the scheme is capable in shock capturing as well as solving wettingdrying problems. In addition, flow with wide range of Froude number over different forms of broad crested weirs, have been employed to verify the robustness, accuracy and stability of the scheme. Hence all of these results prove that the presented wellbalanced scheme is able to simulate different cases of shallow water equation examples accurately.