Numerical study of thin-film flows and open-channel flows

A numerical procedure is developed capable of simulating gravity-driven film flows in two-dimensions. The moving boundary problem is handled through the ALE formulation. In the case of turbulent fluid flows, the two-equations $k-\epsilon$ closure model is used to model the turbulence. A Chorin-type projection scheme is utilized to decouple the velocity and pressure fields, and the spatial discretization is done using the Finite Element Method. Thin liquid films draining down vertical or inclined planes are susceptible to long wavelength disturbances. An extensive numerical study of the surface wave instability in isothermal thin film flows is done by solving the full-scale nonlinear system. Temporal stability analysis of a spatially periodic disturbance reveals interesting wave dynamics. The transition from nearly sinusoidal supercritical waves to broad-banded solitary waves is found to go through a quasi-periodic regime. In this quasi-periodic state, the fundamental mode and several of its harmonics are in an oscillatory state, with continuous exchange of energy. An extensive parametric search has been done to obtain the phase boundary delineating the quasi-periodic regime. Complex wave interactions such as wave-splitting and wave-merging are discussed. Spatial stability analysis akin to the usual experimental studies is done and comparisons are made with the experiments. For the development of successful theories capable of predicting the formation of bedforms, it is essential to understand the turbulent fluid flow on top of the bed. To this end, mean flow and turbulence characteristics for flow over artificial stream-wise periodic bedforms are obtained. Due to the local accelerations associated with stream-line bending, very large velocities and stresses are found to exist at the tip of the dune The separation wake turbulence is found to completely dominate the wall-generated turbulence, and the maximum turbulence intensity levels occur at a distance, approximately equal to the dune height, away from the bed. The accuracy of the rigid-lid approximation is determined by computing the flow field with and without the rigid-lid approximation. The rigid-lid approximation is found to over-predict the shear at the dune crest. Lastly, the mean flow and turbulence characteristics in hydraulic jump are obtained. In the case of flow with inlet supercritical Froude number 2.0, a small recirculation zone is found to exist at the foot of the jump. The mixing layer turbulence associated with the surface roller and the recirculation zone are found to dominate the wall-generated turbulence

Two-dimensional modeling of flow and sedimentation

A two-dimensional (vertical) flow model is developed to simulate open channel flow. The governing equations are solved without any depth averaging. Free surface movement is handled with the help of Arbitrary Lagrangian Eulerian Method. $k - \epsilon$ closure is employed to determine the turbulent eddy viscosity. Spatial discretization is done using linear 3-noded triangular elements. Galerkin method of weighted residuals is used to obtain weak form formulation. Navier-Stokes equations are solved by marching in time using Fractional Step Scheme. Continuity equation is replaced with a pressure Poisson equation. Supercritical flow and hydraulic jump are simulated to validate the flow code. Suspended sediment transport is modeled using an advection-diffusion transport equation for sediment concentration. Bed load transport is determined with the help of empirical correlations. The boundary condition at the bed is specified as a reference concentration. Sample problem involving bed erosion and deposition is solved