We developed a two-dimensional cellular automata based model, which couples the dynamics of overland flow, infiltration processes and surface evolution through sediment transport. The natural hill slopes are represented by their topographic elevation and spatially varying soil properties infiltration rates and surface roughness coefficients. This model allows modeling of Hortonian overland flow and infiltration during complex rainfall events. An advantage of the cellular automata approach over the kinematic wave equations is that wet/dry interfaces that often appear with rainfall overland flows can be accurately captured and are not a source of numerical instabilities. An adaptive explicit time stepping scheme allows for rainfall events to be adequately resolved in time, while large time steps are taken during dry periods to provide for simulation run time efficiency. The time step is constrained by the CFL condition and mass conservation considerations. The spatial discretization is shown to be first order accurate. For validation purposes, hydrographs for non-infiltrating and infiltrating plates are compared to the kinematic wave analytic solutions and data taken from literature. Results show that our cellular automata model quantitatively accurately reproduces hydrograph patterns. However, recent works have showed that even through the hydrograph is satisfyingly reproduced, the flow field within the plot might be inaccurate. For a more stringent validation, we compare steady state velocity, water flux, and water depth fields to rainfall simulation experiments. Comparisons show that our model is able to accurately capture these flow properties. Currently, a sediment transport and deposition module is being implemented and tested.

[1] Jozefik, Z., Schmidt, H., Hinz, Ch. (2015). A cellular automata approach for modeling surface water runoff. European Geosciences Union General Assembly, Vienna, Austria 


The correct description of turbulent buoyant multi-phase flow, e.g. flow of CO2 at very high pressures, is very challenging, mostly due to the amount of different length and time schales. Goal of the project was to capture a broad range of those scales by using a stochastic approach.

[1] Schulz, F. T., Glawe, C., Schmidt, H., Kerstein, A. R. (2013). Toward modeling of CO2 multi-phase flow patterns using a stochastic multi-scale approach Environmental Earth Science, DOI 10.1007/s12665-013-2461-5

[2] Glawe, C., Gonzalez-Juez, E. D., Schmidt, H., Kerstein, A.R. (2013). ODTLES simulations of turbulent flows through heated channels and ducts, 8th International Symposium on Turbulence and Shear Flow Phenomena, Poitiers, France

[3] Schulz, F., Glawe, C., Schmidt, H., Kerstein, A.R. (2013). Liquid jet simulation using one-dimensional turbulence, 8th International Symposium on Turbulence and Shear Flow Phenomena, Poitiers, France

This website uses cookies. Those have two functions: On the one hand they are providing basic functionality for this website. On the other hand they allow us to improve our content for you by saving and analyzing anonymized user data. You can redraw your consent to using these cookies at any time. Find more information regarding cookies on our Data Protection Declaration and regarding us on the Imprint.


These cookies are needed for a smooth operation of our website.


Name Purpose Lifetime Type Provider
_pk_id Used to store a few details about the user such as the unique visitor ID. 13 months HTML Matomo
_pk_ref Used to store the attribution information, the referrer initially used to visit the website. 6 months HTML Matomo
_pk_ses Short lived cookie used to temporarily store data for the visit. 30 minutes HTML Matomo