A mapping technique for numerical computations of bed evolutions

Free surface flow is one of the most difficult problems in engineering to be solved, since velocity and pressure fields depend on the free surface. On the other hand, the position of the free surface is unknown previously. Furthermore, the boundary condition on the free surface is expressed by a complicated equation. In an alluvial stream, where the boundaries of the domain are not fixed, addition of free surface at the bed will increase this difficulty. A domain mapping technique is developed in this paper to study the bed evolutions. The flow is considered 2D, choosing two coordinates in streamwise and upward directions. With a proper transformation, the hydrodynamics and sediment transport governing equations in irregular domain will be mapped into a simple rectangular one. The new domain can be discretize by finite elements. The transformed governing equations are solved to obtain desired variables in the mapped domain. With a proper transformation, there is no need of inverse mapping to obtain the free water surface profile and bedform evolution and migration in the actual domain. The model has been applied to streams with movable bed and the results show a good agreement with the experimental experiences.