Assessment of sensitivity of depth integrated solutions to longitudinal discontinuities on the channel bed

Download
2019
Mohammadi, Ramez
Depth integrated equations can be easily solved over large domains to provide flood inundation maps. In urban and rural areas however, there may be numerous natural or artificial bottom boundary discontinuities in the form of rapid variations in bed elevation. Such discontinuities cause abrupt changes in the source terms of the governing equations and can significantly affect stability and accuracy of the numerical solution. This study concentrates on the assessment of sensitivity of the governing equations to longitudinal discontinuities and proposes solutions to alleviate associated numerical complications. The presence of dry regions in the domain is also considered which requires additional modifications in the code to deal with the wet/dry fronts. In this thesis, Godunov’s type Finite Volume Method (FVM) is used for the numerical solution of the shallow water equations. Weighted Averaged Flux (WAF) method, that is 2nd order extension of the first order Godunov’s scheme, based on HLL Riemann solver, is used to compute the fluxes. The source term treatment is based on the “hydrostatic reconstruction” and first and second order well-balanced schemes are obtained. Creation of probable negative water depths during linear piecewise reconstruction of the water surface near dry areas has been prevented in the 2nd order well-balanced scheme. Eight test cases with their available analytical solutions that are widely used in the literature are selected to validate the developed codes. The results of the developed codes are also compared with sets of experimental data if available. Test cases are solved by pure WAF method without well-balancing property, WAF method with first order well-balanced scheme, and WAF method with second order well-balanced scheme. Two test cases are specifically selected to show the capability of the developed codes in solving flows in regions with wet and dry areas. In the end, a new test case is introduced to observe behavior of the numerical solutions in a bed consisting of a series of sharp corners due to positive and negative steps. It is observed that well-balanced schemes can produce water surface profiles without spurious oscillations and steady-state horizontal hydrostatic water surface is recovered without noisy fluctuations.

Suggestions

Depth integrated equations applied to longitudinal discontinuities on the channel bed
Razmand, Maral; Aydın, İsmail; Department of Civil Engineering (2015)
Depth integrated equations can be solved over large domains to provide flood inundation maps. In urban and rural areas however, there may be numerous natural or artificial bottom boundary discontinuities in the form of rapid variations in the bed elevation. Such discontinuities cause abrupt changes in the source terms of the governing equations and can significantly affect stability and accuracy of the numerical solution. A 1D code is developed for shallow water equations using HLL approximate Riemann solve...
Water surface profile computation around square blocks by imaginary representation of the obstructions in the shallow water equations
Hofioni, Sayed Omar; Aydın, İsmail; Department of Civil Engineering (2021-7)
Depth integrated shallow flow equations can be solved easily over large domains to provide flood inundation maps. However, there may be numerous natural or artificial obstructions that increase resistance to flow due to drag forces induced on the submerged volumes. It is possible to numerically solve the governing equations on a suitable grid system for the flow field around such obstructions and determine the flow depths. But, when the number of obstructions increases with arbitrary size and distribution, ...
Inverse problems for parabolic equations
Baysal, Arzu; Çelebi, Okay; Department of Mathematics (2004)
In this thesis, we study inverse problems of restoration of the unknown function in a boundary condition, where on the boundary of the domain there is a convective heat exchange with the environment. Besides the temperature of the domain, we seek either the temperature of the environment in Problem I and II, or the coefficient of external boundary heat emission in Problem III and IV. An additional information is given, which is the overdetermination condition, either on the boundary of the domain (in Proble...
Development of a two-dimensional Navier-Stokes solver for laminar flows using cartesian grids
Şahin, Mehmet Serkan; Aksel, Mehmet Haluk; Department of Mechanical Engineering (2011)
A fully automated Cartesian/Quad grid generator and laminar flow solver have been developed for external flows by using C++. After defining the input geometry by nodal points, adaptively refined Cartesian grids are generated automatically. Quadtree data structure is used in order to connect the Cartesian cells to each other. In order to simulate viscous flows, body-fitted quad cells can be generated optionally. Connectivity is provided by cut and split cells such that the intersection points of Cartesian ce...
Impact of physical and chemical heterogeneities on optimal remediation design and costs
Aksoy, Ayşegül (null; 2001-01-24)
The impact of physical and chemical aquifer heterogeneities on optimal remediation design and costs is investigated by linking a genetic algorithm optimization library with a contaminant transport simulation model. Various levels of physical and chemical (sorption) aquifer heterogeneities are examined. In the first level, heterogeneity is limited to the hydraulic conductivity (K) field. Then systems with heterogeneity in both K and the distribution coefficient (Kd) are considered. The final level of heterog...
Citation Formats
R. Mohammadi, “Assessment of sensitivity of depth integrated solutions to longitudinal discontinuities on the channel bed,” Thesis (M.S.) -- Graduate School of Natural and Applied Sciences. Civil Engineering., Middle East Technical University, 2019.