NUMERICAL COMPUTATION OF ONE- AND TWO-LAYER SHALLOW FLOW MODEL

  • Komang Dharmawan Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Udayana, Indonesia https://orcid.org/0000-0002-7021-1386
  • Putu Veri Swastika Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Udayana, Indonesia
  • G K Gandhiadi Department of Mathematics, Faculty of Mathematics and Natural Sciences, Universitas Udayana, Indonesia
Keywords: Saint Venant, Dam-Break, Sub-Maximal Exchange, Momentum-Conserving Staggered-Grid

Abstract

In this research, we study a proficient computational model designed to simulate shallow flows involving one- and two-layer shallow flow. This numerical model is built upon the Saint Venant equations, which are widely used in hydraulics to depict the behavior of shallow water flow. The numerical scheme used here is constructed based on the conventional leapfrog technique implemented on a staggered grid framework, referred to as MCS. The primary objective of this research is to re-examine and implement the MCS in accurately modelling the free surface and interface waves produced by different flows passing through irregular geometries. Unlike the conventional MCS, we modify the momentum conservation principle to be more general, accommodating a non-negative wet cross-sectional area due to irregular geometry. We successfully conduct numerous numerical simulations by examining various scenarios involving one-layer and two-layer flow through irregularly shaped channels or structures. Our results show that the correct surface wave profile generated by a one-dimensional dam break through the triangular obstacle in the open channel can be simulated very well. Comparison with the existing experimental data seems promising although some disparities are being found due to dispersive phenomena with RMSE less than 5%. Furthermore, our scheme is successfully extended to simulate the steady sub-maximal exchange in two-layer flows using specific boundary conditions. The alignment between the submaximal numerical results with exchange flow theory is noticeable in the interface profile, characteristics of flow conditions and the flux values achieved when the steady situation occurs. These satisfying results indicate that our proposed numerical model can be used for practical needs involving various flow situations both one and two-layer cases

Downloads

Download data is not yet available.

References

D. Bourgault, P. S. Galbraith, and C. Chavanne, “Generation of internal solitary waves by frontally forced intrusions in geophysical flows,” Nat. Commun., vol. 7, 2016, doi: 10.1038/ncomms13606.

M. J. Castro, J. A. García-Rodríguez, J. M. González-Vida, J. Marías, C. Parés, and M. E. Vázquez-Cendón, “Numerical simulation of two-layer shallow water flows through channels with irregular geometry,” J. Comput. Phys., vol. 195, no. 1, pp. 202–235, 2004, doi: 10.1016/j.jcp.2003.08.035.

M. J. Castro, J. A. García-Rodríguez, J. M. González-Vida, J. Macías, and C. Parés, “Improved FVM for two-layer shallow-water models: Application to the Strait of Gibraltar,” Adv. Eng. Softw., vol. 38, no. 6, pp. 386–398, 2007, doi: 10.1016/j.advengsoft.2006.09.012.

S. R. Pudjaprasetya, V. M. Risriani, and Iryanto, “Numerical simulation of propagation and run-up of long waves in U-shaped bays,” Fluids, vol. 6, no. 4, pp. 1–13, 2021, doi: 10.3390/fluids6040146.

K. R. Helfrich, “Time Dependent Two-Layer Hydraulic Exchange Flows,” J. Phys. Oceanogr., vol. 25, no. 3, pp. 359–372, 1995, doi: 10.1175/1520-0485(1995)025<0359:TDTLHE>2.0.CO;2.

P. Brandt, W. Alpers, and J. O. Backhaus, “Study of the generation and propagation of internal waves in the Strait of Gibraltar using a numerical model and synthetic aperture radar images of the European ERS 1 satellite,” J. Geophys. Res. C Oceans, vol. 101, no. C6, pp. 14237–14252, 1996, doi: 10.1029/96JC00540.

P. Brandt, A. Rubino, W. Alpers, and J. O. Backhaus, “Internal waves in the Strait of Messina studied by a numerical model and synthetic aperture radar images from the ERS 1/2 satellites,” J. Phys. Oceanogr., vol. 27, no. 5, pp. 648–663, 1997, doi: Numerical study of cascading dam-break characteristics using SWEs.

S. Dai, Y. He, J. Yang, Y. Ma, S. Jin, and C. Liang, “Numerical study of cascading dam-break characteristics using SWEs and RANS,” Water Supply, vol. 20, no. 1, pp. 348–360, Feb. 2020, doi: 10.2166/ws.2019.168.

Y. Zhang and P. Lin, “An improved SWE model for simulation of dam-break flows,” Proc. Inst. Civ. Eng. Water Manag., vol. 169, no. 6, pp. 260–274, 2016, doi: Dam-break flows during initial stage using SWE and RANS approaches.

H. Ozmen-Cagatay and S. Kocaman, “Dam-break flows during initial stage using SWE and RANS approaches,” J. Hydraul. Res., vol. 48, no. 5, pp. 603–611, Oct. 2010, doi: 10.1080/00221686.2010.507342.

P. V. Swastika, M. Fakhruddin, S. Al Hazmy, S. Fatimah, and A. De Souza, “A novel technique for implementing the finite element method in a shallow water equation,” MethodsX, vol. 11, p. 102425, Dec. 2023, doi: 10.1016/j.mex.2023.102425.

R. Iacono, “Analytic solutions to the shallow water equations,” Phys. Rev. E, vol. 72, no. 1, p. 017302, Jul. 2005, doi: 10.1103/PhysRevE.72.017302.

M. J. Castro et al., “The numerical treatment of wet/dry fronts in shallow flows: Application to one-layer and two-layer systems,” Math. Comput. Model., vol. 42, no. 3–4, pp. 419–439, 2005, doi: 10.1016/j.mcm.2004.01.016.

B. Arry Sanjoyo, M. Hariadi, and M. H. Purnomo, “Stable Algorithm Based On Lax-Friedrichs Scheme for Visual Simulation of Shallow Water,” Emit. Int. J. Eng. Technol., vol. 8, no. 1, pp. 19–34, Jun. 2020, doi: 10.24003/emitter.v8i1.479.

A. Arakawa, “Computational Design for Long-Term Numerical Integration of the Equations of Fluid Motion: Two-Dimensional Incompressible Flow. Part I,” Journal of Computational Physics, vol. 1, no. 1, pp. 119–143, 1997.

G. S. Stelling and S. P. A. Duinmeijer, “A staggered conservative scheme for every Froude number in rapidly varied shallow water flows,” Int. J. Numer. Methods Fluids, vol. 43, no. 12, pp. 1329–1354, 2003, doi: 10.1002/fld.537.

P. V. Swastika, S. R. Pudjaprasetya, L. H. Wiryanto, and R. N. Hadiarti, “A momentum-conserving scheme for flow simulation in 1D channel with obstacle and contraction,” Fluids, vol. 6, no. 1, p. 26, 2021, doi: 10.3390/fluids6010026.

P. V. Swastika and S. R. Pudjaprasetya, “The Momentum Conserving Scheme for Two-Layer Shallow Flows,” Fluids, vol. 6, p. 346, 2021, doi: 10.3390/fluids6100346.

N. Erwina, D. Adytia, S. R. Pudjaprasetya, and T. Nuryaman, “Staggered Conservative Scheme for 2-Dimensional Shallow Water Flows,” Fluids, vol. 5, no. 3, pp. 1–18, 2020, doi: 10.3390/fluids5030149.

S. R. Pudjaprasetya and R. Sulvianuri, “The momentum conservative scheme for simulating nonlinear wave evolution and run-up in U-shaped bays,” Jpn. J. Ind. Appl. Math., vol. 40, no. 1, pp. 737–754, Jan. 2023, doi: 10.1007/s13160-022-00549-4.

P. V. Swastika, S. R. Pudjaprasetya, and N. Subasita, “Numerical simulation of two-layer shallow water flows; Exchange Flow in Lombok Strait,” East Asian J. Appl. Math., vol. Accepted and to be appear 2025 (in Private Communication).

S. Soares-Frazão, “Experiments of dam-break wave over a triangular bottom sill,” J. Hydraul. Res., vol. 45, no. sup1, pp. 19–26, Dec. 2007, doi: 10.1080/00221686.2007.9521829.

D. M. Farmer and L. Armi, “Maximal two-layer exchange over a sill and through the combination of a sill and contraction with barotropic flow,” J. Fluid Mech., vol. 164, no. 10, pp. 53–76, 1986, doi: 10.1017/S002211208600246X.

I. Magdalena, A. A. A. Hariz, M. Farid, and M. S. B. Kusuma, “Numerical studies using staggered finite volume for dam break flow with an obstacle through different geometries,” Results Appl. Math., vol. 12, p. 100193, 2021, doi: 10.1016/j.rinam.2021.100193.

S. R. Pudjaprasetya and I. Magdalena, “Momentum conservative schemes for shallow water flows,” East Asian J. Appl. Math., vol. 4, no. 2, pp. 152–165, 2014, doi: 10.4208/eajam.290913.170314a.

G. S. Stelling, A. K. Wiersma, and J. B. T. M. Willemse, “Practical Aspects of Accurate Tidal Computations,” J. Hydraul. Eng., vol. 112, no. 9, pp. 802–816, 1986, doi: 10.1061/(asce)0733-9429(1986)112:9(802).

A. Balzano, “Evaluation of methods for numerical simulation of wetting and drying in shallow water flow models,” Coast. Eng., vol. 34, no. 1–2, pp. 83–107, 1998, doi: 10.1016/S0378-3839(98)00015-5.

K. Dharmawan, P. V. Swastika, G. K. Gandhiadi, and S. R. Pudjaprasetya, “A Non-hydrostatic Model for Simulating Dam-Break Flow Through Various Obstacles,” MENDEL, vol. 30, no. 1, pp. 33–42, Jun. 2024, doi: 10.13164/mendel.2024.1.033.

Published
2024-07-31
How to Cite
[1]
K. Dharmawan, P. Swastika, and G. Gandhiadi, “NUMERICAL COMPUTATION OF ONE- AND TWO-LAYER SHALLOW FLOW MODEL”, BAREKENG: J. Math. & App., vol. 18, no. 3, pp. 1509-1518, Jul. 2024.