Open access peer-reviewed compact

Program Description of the Diffusion Hydrodynamic Model

By Theodore V. Hromadka II and Chung-Cheng Yen

Reviewed: June 17th 2020Published: September 9th 2020

DOI: 10.5772/intechopen.93209

Downloaded: 259


The numerical algorithm, with a focus on the interface element that is used in the diffusion hydrodynamic model, is presented in this chapter. The source program was written in FORTRAN language, and it can be downloaded from this book companion website. The channel, flood plain, and the interface flow details are discussed.


  • flood plain
  • overflow
  • control volume
  • water elevation
  • channel interface

1. Introduction

A computer program for the two-dimensional diffusion hydrodynamic model which is based on the diffusion form of the St. Venant equations [1, 2, 3, 4, 5] where gravity, friction, and pressure forces are assumed to dominate the flow equation will be discussed in this section.

The DHM consists of a 1-D channel and 2-D flood plain models, and an interface sub-model. The one-dimensional channel element utilizes the following assumptions:

  • Infinite vertical extensions on channel walls (Figure 1)

  • Wetted perimeter is calculated as shown in Figure 1(a)

  • Volumes due to channel skew are ignored (Figure 1(b))

  • All overflow water is assigned to one grid element (Figure 2)

Figure 1.

One-dimensional channel flow element characteristics: (a) element geometrics, (b) element associations to grid elements, and (c) channel element connections.

Figure 2.

Grid element nodal molecule.

The interface model calculates the excess amount of water either from the channel element or from the flood plain element. This excess water is redistributed to the flood plain element or the channel element according to the water surface elevation.

This FORTRAN program has the capabilities to simulate both one-and two-dimensional surface flow problems, such as the one-dimensional open channel flow and two-dimensional dam-break problems illustrated in the preceding pages. Engineering applications of the program will be presented in the next chapter.


2. Interface model

2.1 Introduction

The interface model modifies the water surface elevations of flood plain grids and channel elements at specified time intervals (update intervals). There are three cases of interface situations: (1) channel overflow, (2) grid overflow, and (3) flooding of channel and grid elements.

2.2 Channel overflow

When the channel is overflowing; the excess water is temporarily stored in the vertically extended space (Figure 3(b)). Actually, it is the volume per unit length. This excess water is the product of the depth of water, the width of the channel, and the length of the channel and is subsequently uniformly distributed over the grid elements. In other words, the new grid water surface elevation is equal to the old water surface elevation plus a depth of hw/L, and the channel water surface elevation now matches the parent grid water surface elevation.

Figure 3.

Diffusion hydrodynamic interface model: (a) model interface geometries, (b) channel overflow interface model, (c) grid overflow interface model, (d) grid-channel flooding interface model, and (e) channel-grid flooding interface model.

2.3 Grid overflow

When the water surface elevation of the grid element is greater than a specified surface detention (Figure 3(a)), the excess water drains into the channel element, and the new water surface elevation is changed according to the following two conditions (Figure 3(c)), (a) if v > v’, where v denotes the excess volume of water per unit length and v’ denotes the available volume per unit length, the new water surface of the grid element is ANEW = AOLD − (v − v’)/L, and the new water surface elevation of the channel element is also equal to ANEW, and (b) if v ≤ v’, the new water surface elevation of the grid element is ANEW = AOLD − h and the new water surface elevation of the channel element is BNEW = BOLD + v/w.

2.4 Flooding of channel and grid

When flooding occurs, the water surface elevations of the grid and channel elements are both greater than the specified surface detention elevation. Two cases have to be considered as follows:

  1. If A > B (Figure 3(d)), the new water surface elevation of the grid element is ANEW=BOLD+hLwL, and the new water surface elevation of the channel element is equal to ANEW.

  2. If A < B (Figure 3(e)), the new water surface elevation of the grid element is ANEW=AOLD+h*wL, and the new water surface elevation of the channel element is equal to ANEW.


3. Conclusions

The characteristic features of the diffusion hydrodynamic model are detailed with a focus on its ability to model the interface. The interface component in the model can modify the water elevation in the grids along the flood plain and channel to account for different overflow scenarios, which are also illustrated.

© 2020 The Author(s). Licensee IntechOpen. This compact is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License, which permits use, distribution and reproduction for non-commercial purposes, provided the original is properly cited.

How to cite and reference

Link to this compact Copy to clipboard

Cite this compact Copy to clipboard

Theodore V. Hromadka II and Chung-Cheng Yen (September 9th 2020). Program Description of the Diffusion Hydrodynamic Model, A Diffusion Hydrodynamic Model, Theodore V. Hromadka II, Chung-Cheng Yen and Prasada Rao, IntechOpen, DOI: 10.5772/intechopen.93209. Available from:

compact statistics

259total compact downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next compact

Applications of Diffusion Hydrodynamic Model

By Theodore V. Hromadka II and Chung-Cheng Yen

Related Book

First chapter

Strategies for Testing the Impact of Natural Flood Risk Management Measures

By Barry Hankin, Peter Metcalfe, David Johnson, Nick A. Chappell, Trevor Page, Iain Craigen, Rob Lamb and Keith Beven

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us