Development of a two-dimensional flow model using a time-splitting scheme and regular rectangular mesh

Document Type : Original Article

Authors

1 Department of Water Sciences and Engineering, Faculty of Agriculture, Razi University, Kermanshah, Iran.

2 Department of Civil Engineering, Faculty of Engineering, Razi University, Kermanshah, Iran.

Abstract

Objective: In this study, a two-dimensional numerical model is presented to simulate the flow in open channels. The governing equations of the flow are depth-averaged Navier-Stokes equations, which are separated and solved using the finite volume method and the time-splitting scheme. One of the advantages of the time-splitting method is the ability to solve the different terms of the equations separately using the best numerical scheme. Another advantage of this method is its simplicity and ease of understanding, because in all stages the equations are solved in one dimension.
 
Method: Generating the mesh is one of the most important initial operations for solving the governing equations. in this research, a regular rectangular grid was used. The water level variable was defined at the center of the cell and the velocity and discharge were defined at the boundaries of each cell. In order to solve the governing equations of the flow in this model, each time step is divided into several stages. Initially, the advection terms of the equation of motion were solved based on the Fromm scheme and its output was considered for the diffusion term. The diffusion term in the momentum equation was solved implicitly using the TDMA algorithm. Then, the friction term in the flow equation was separated and solved explicitly, and then the water level was calculated by simultaneously solving the continuity equation and the remaining terms of the momentum equations.
 
Results: In order to verify the validity of the two-dimensional model, The longitudinal profile of gradually varying flow in a rectangular channel, hydraulic jump, sudden channel opening, and asymmetric dam failure were considered for evaluating the model. The output results of the numerical model show that the developed model is capable of simulating the flow for different conditions with appropriate accuracy.
 
Conclusions: The time splitting algorithms with staggered grid approach are used to handle the coupling of water depth and velocity parameter. In this case, water level variable was defined at the center of the cell and the velocity and discharge were defined at the boundaries of each cell. The proposed model has been tested in several experimental cases of flow. The simulated water depth and velocity in time and space are in good agreement with the measured data

Keywords

Main Subjects


Abbott, M.B., & Basco, D. (1989). Computational Fluid Mechanics. An Introduction for Engineers. Longman Scientific & Technical publications, New York, United States. https://archive.org/details/computationalflu0000abbo
Anastasio, K., & Chan, C.T. (1997). Solution of the 2D shallow water equations using the finite volume method on unstructured triangular meshes. International Journal for Numerical Methods in Fluids, 24(11), 1225–1245. https://doi.org/10.1002/(SICI)1097-0363(19970615)24:11
Busto, S., & Dumbser, M. (2022). A staggered semi-implicit hybrid finite volume / finite element scheme for the shallow water equations at all Froude numbers. Applied Numerical Mathematics, 175, 108-132. https://doi.org/10.1016/j.apnum.2022.02.005
Cebeci, T., & Bradshaw, P. (1977). Momentum transfer in boundary layers. Washington : Hemisphere Publications, Corp, Washington, D.C, United States.                                        https://archive.org/details/momentumtransfer00cebe
Chaudhry, M.H. (2008). Open-Channel Flow. Springer New York, New York, USA.                        https://link.springer.com/book/10.1007/978-0-387-68648-6
Elder, J.W. (1959). The dispersion of marked fluid particles in turbulent shear flow. Journal of Fluid Mechanics, 5(4), 544–560. https://doi.org/10.1017/S0022112059000374
Fennema, R., & Chaudhry, M. (1990). Explicit methods for 2D transient free-surface flows. Journal of Hydraulic Engineering, 116(8), 1013–1034. https://doi.org/10.1061/(ASCE)0733-9429(1990)116:8(1013)
Gharangik, A.M. (1988). Numerical simulation of hydraulic jump. MSc thesis, Washington state university, Pullman, Washington.
Ghobadian, R. (2019). Application of lax and leap-frog schemes in two-dimensional modeling of dam breaking using explicit finite-difference method. Journal of Water and Soil Science, (29)2, 85-99. https://water-soil.tabrizu.ac.ir/article_9273.html (In Persian)
Groosi, F., Cusicahua, A., Shademani, M., & Shakibaeinia, A. (2022). Experimental and numerical investigations of dam break flow over dry and wet beds. International Journal of Mechanical Sciences, 215, 106946. https://doi.org/10.1016/j.ijmecsci.2021.106946
Hou, J., Liang, Q., Zhang, H. & Hinkelmann, R. (2015). An efficient unstructured MUSCL scheme for solving the 2D shallow water equations. Environmental Modeling & Software, 66, 131-152. https://doi.org/10.1016/j.envsoft.2014.12.007
Mehrmousavi, Z. (2019). Development of unsteady flow 2D numerical model and its implementation on dam-break simulation. Ph.D Thesis, Faculty of Agricultural Sciences and Engineering of Razi University, Kermanshah, Iran. (In Persian)
Mehrmousavi, Z., Ghobadian, R., & Javan, M. (2019) Investigation of the effect of reservior shape coefficient on dam-break waves using Leap-Frog and Lax methods in curvilinear coordinates. Journal of Water and Soil Conservation, 26(1), 131-150. https://jwsc.gau.ac.ir/article_4554.html (In Persian)
Molls, T., & Chaudhry, M.H. (1995). Depth-averaged open-channel flow model. Journal of Hydraulic Engineering, 121 (6), 453-465.https://doi.org/10.1061/(ASCE)0733-9429(1995)121:6(453)
Namin, M.M. (2003). A fully three dimensional non hydrostatic free surface flow model for hydro environmental predictions. Ph.D. Thesis, Cardiff School of Engineering of Cardiff University, Cardiff ,United Kingdom.
Namin, M.M., Lin B. & Falconer, R. A. (2004) Modeling estuarine and coastal flows using an unstructured triangular finite volume algorithm. Advances in Water Resources, 27(12), 1179-1197. https://doi.org/10.1016/j.advwatres.2004.08.012
Norouzi, S., Kolahdoozan, M., & Zarrati, A.R. (2022). Development of a 2D depth-averaged model for calculating scouring and deposition in alluvial streams, Scientia Iranica, 29(3), 941-950. https://doi.org/10.24200/sci.2022.50705.1830
Nujic, M. (1995). Efficient implementation of non-oscillatory scheme for the computation of free-surface flows. Journal of Hydraulic Research. 33(1), 101–11. https://doi.org/10.1080/00221689509498687
Papanicolaou, A.N., Elhakeem, M. & Wardman, B. (2011). Calibration and verification of a 2D-hydrodynamic model for simulating flow around bendway weir structures. Journal of Hydraulic Engineering. 137(1), 75–89. https://doi.org/10.1061/(ASCE)HY.1943-7900.0000280.
Patankar, S. (1980). Numerical heat transfer and fluid flow, CRC Press, Boca Raton, USA.
https://doi.org/10.1201/9781482234213
Rajaratnam, N., & Nwachukwu B.A. )1983). Flow near Groin-like structures. Journal of Hydraulics Division. 109, 463-480. https://doi.org/10.1061/(ASCE)0733-9429(1983)109:3(463)
Rodi, W. (1984). Turbulence models and their application in hydraulics: A state of the art review.Balkema Publishers, Leiden, The Netherlands.                                         https://www.researchgate.net/publication/327303543_Turbulence_models_and_their_application_in_hydraulics_A_state-of-the-art_review_third_edition
Sarveram, H., & Shamsaie, A. (2013). A implicit, semi-lagrangian method for 2D numerical modeling of dam break flow. Iranian Water Research Journal, 7(13), 21-30. https://iwrj.sku.ac.ir/article_10953.html (In Persian)
Teymouri yeganeh, M., Heidari, M.M., & Ghobadian, R. (2023) Two-dimensional simulation of flow pattern in open channels using time-splitting method. Journal of Irrigation and Water Engineering, 13, 2131-247. https://www.waterjournal.ir/article_177279.html (In Persian)
Versteeg, H. K., & Malalasekera, W. (1995). An introduction to computational fluid dynamics, the finite volume method, Pearson Education Limited Publications, England, United Kingdom. https://www.researchgate.net/profile/GhassanSmaisim/post/FEM_mesh_generator/attachment/59d655b479197b80779acc78/AS%3A526908706508800%401502636233004/download/110+Versteeg+2007+an+introduction+to+computational+fluid+dynamics+the+finite+volume+method+2nd+edition.pdf
Wu, W. (2004). Depth-averaged two-dimensional numerical modeling of unsteady flow and nonuniform sediment transport in open channels. Journal of Hydraulic Engineering, 130 (10),1013–1024. https://doi.org/10.1061/(ASCE)0733-9429(2004)130:10(1013)
Wu, Y., & Yang, Z. (2024) A depth-averaged SPH-FV landslide dynamic model for evaluating hazard zones. Computers and Geotechnics, 169, 106210.                                https://doi.org/10.1016/j.compgeo.2024.106210
Xie, B.L. (1994). Flow data measured in a channel with sudden expansion. Wuhan University, China. https://www.ncche.olemiss.edu/sites/default/files/files/docs/cche2d/validation.pdf
Zarrati, A.R., Tamai, N., & Jin, Y.C. (2005) Mathematical modeling of meandering channels with a generalized depth averaged model. Journal of Hydraulic Engineering. 131, 467–475. https://doi.org/10.1061/(ASCE)0733-9429(2005)131:6(467)
Zhang, S., Ghidaoui, M., & Gray, W. (2003). A kinetic flux vector splitting scheme for shallow water flows. Advances in Water Resources, 26(6), 635–47. https://doi.org/10.1016/S0309-1708(03)00029-0