A HIGH RESOLUTION FINITE VOLUME METHOD FOR SOLVING SHALLOW WATER EQUATIONS(2)

6. applications of dam-break computation

through the computation of 1d dam-break waves in a horizontal and frictionless channel and the comparison with stoker's theoretical solution, it is shown that steep and nonoscillatory numerical solutions could be obtained using the hybrid type of tvd scheme . two typical examples of 2d dam-break problems are solved and discussed by solving the shallow water equations using above finite volume tvd scheme.

6.1 rectangular dam-break

consider a 2d partial dam-break model with a non-symmetrical breach. it is assumed that in the center of a 200m × 200m channel, a partial dam breaking takes place instantaneously. the breach is 75m in length, which has distances of 30m from the left bank and 95m from the right. the initial water height is 10m and 5m respectively. no slope and friction are considered. the results displaying the views of the water surface elevation, contour of the surface elevation and velocity field are shown in figure3 at time t=7.2s after the dam failure. at the instant of breaking of the dam, water is released through the breach, forming a positive wave propagating downstream and a negative wave spreading upstream. these results agree quite well with the results of using finite difference hybrid type of tvd scheme and those in ref. .

fig. 3(a) water surface elevation for a rectangular dam-break

fig. 3(b) contour of surface elevation for a rectangular dam-break

6.2 circular dam-break

another typical example is based on the hypothetical test case studied by alcrudo and garcia-navarro [7], which involves the breaking of a circular dam. it is an important test example for the analysis of the algorithm performance and solving a complex shallow water problem. the physical model is that two regions of still water are separated by a cylindrical wall of radius 11m. the water depth inside the dam is 10m, whilst outside the dam is 1m. at the instant of dam failure the circular wall is assumed to be removed completely and no slope and friction is considered, then the circular dam-break waves will spread and propagate radially and symmetrically. the results with above method at time t=0.69s are shown in figures 4 (a), (b) and (c) which denote the water surface elevation, contour of surface elevation and velocity field respectively. it can be clearly seen that the waves spread uniformly and symmetrically. these results agree quite well with those given by alcrudo and garcia-navarro

, zhao et al. , alastansiou and chan and they can be tested each other. it demonstrates that the present method is reliable and fine.

fig. 3(c) velocity field for a rectangular dam-break

fig. 4(a) water surface elevation for a circular dam-break circular dam-break

fig. 4(b) contour of surface elevation for a circular dam-break

fig. 4(c) velocity field for a circular dam-break

7. summary and conclusions

tvd scheme is playing an important role in gas dynamics because of its high accuracy, good shock-capturing ability and nonoscillatory numerical performance. but it is constructed based on finite difference method. in this paper a new geometry and topology is defined for the extension of nodes to elements. with the conservative type of the shallow water equations, a hybrid type second order tvd scheme is applied and two-step runge –kutta method is adopted in time, then a finite volume tvd scheme for the shallow water equations on arbitrary quadrilateral elements is developed. the numerical results of two types of dam-break problem show that the method is sufficiently robust and can handle discontinuities and complex flow problems efficiently. the results presented in this paper are in excellent agree with those reported recently and even display sharper discontinuities and the maximum values attenuate more slowly. it can be foreseen that this method has much broader application foreground. as for further studies, such as in the cases of a channel having bend, bifurcation and inner islands, will discuss in another paper.

references

1. a. harten, 1983: high resolution schemes for hyperbolic conservation laws, journal of computational physics, 49, 357-393.

2. s. y. hu, w. y. tan, 1990: numerical modeling of bores due to dam-break, journal of hydrodynamics, ser. a., 5(2), 90~98 (in chinese).

3. j. h. tao, w. d. zhang, 1993: the simulation of one and two dimensional dam-breaking waves by tvni scheme, journal of tian jin university, (1), 7~15 (in chinese).

4. j. y. yang, c. a. hsu, and s. h. chang, 1993: computations of free surface flows, part 1: 1-d dam-break flow, journal of hydraulic research, 31(1).

5. j. s. wang, h. g. ni, s. jin and j. c. li, 1998: simulation of 1d dam-break flood wave routing and reflection by using tvd schemes, journal of hydraulic engineering, (5), 7~11 (in chinese).

6. j. s. wang, h. g. ni, and s. jin, 1998: a high accurate numerical simulation of the propagation and diffraction for 2d dam-break bores, journal of hydraulic engineering, (10), 1~6 (in chinese).

7. f. alcrudo, p. garcia-navarro, 1993: a high resolution godunov-type scheme in finite volumes for the 2d shallow water equation, international journal for numerical method in fluids, 16, 489-505 1993.

8. d. h. zhao, h. w. shen, j. s. lai, and g. q. tabios ⅲ , 1996: approximate riemann solvers in fvm for 2d hydraulic shock wave modeling, journal of hydraulic engineering, 692-702.

9. k. alastansiou, c. t. chan, 1997: solution of the 2d shallow water equations using the finite volume method on unstructured triangular meshes, international journal for numerical method in fluids, 24, 1225-1245.

10. j. s. wang, 1998: a study of numerical simulation of dam-break bores by applying tvd schemes, ph.d. thesis, dalian university of technology (in chinese).