Implementation of physical boundary conditions into computational domain in modelling of oscillatory bottom boundary layers

2010-11-30
Tiğrek, Şahnaz
Yılmaz, Bilgi
This paper discusses the importance of realistic implementation of the physical boundary conditions into computational domain for the simulation of the oscillatory turbulent boundary layer flow over smooth and rough flat beds. A mathematical model composed of the Reynolds averaged Navier-Stokes equation, turbulent kinetic energy (k) and dissipation rate of the turbulent kinetic energy (epsilon) has been developed. Control-volume approach is used to discretize the governing equations to facilitate the numerical solution. Non-slip condition is imposed on the bottom surface, and irrotational main flow properties are applied to the upper boundary. The turbulent kinetic energy is zero at the bottom, whereas the dissipation rate is approaching to a constant value, which is proportional to the kinematic viscosity times the second derivative of the turbulent kinetic energy. The output of the model is compared with the available experimental studies conducted in oscillatory tunnels and wave flume. It is observed that the irrotational flow assumption at the upper boundary is not realistic in case of water tunnels. Therefore, new upper boundary conditions are proposed for oscillatory tunnels. The data of wave flume show good agreement with the proposed numerical model. Additionally, several factors such as grid aspect ratio, staggered grid arrangement, time-marching scheme and convergence criteria that are important to obtain a robust, realistic and stable code are discussed. Copyright (C) 2009 John Wiley & Sons, Ltd.
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS

Suggestions

Boundary element solution of unsteady magnetohydrodynamic duct flow with differential quadrature time integration scheme
Bozkaya, Canan; Tezer, Münevver (Wiley, 2006-06-20)
A numerical scheme which is a combination of the dual reciprocity boundary element method (DRBEM) and the differential quadrature method (DQM), is proposed for the solution of unsteady magnetohydro-dynamic (MHD) flow problem in a rectangular duct with insulating walls. The coupled MHD equations in velocity and induced magnetic field are transformed first into the decoupled time-dependent convection-diffusion-type equations. These equations are solved by using DRBEM which treats the time and the space deriva...
Global volume conservation in unsteady free surface flows with energy absorbing far-end boundaries
Demirel, Ender; Aydın, İsmail (Wiley, 2010-10-30)
A wave absorption filter for the far-end boundary of semi-infinite large reservoirs is developed for numerical simulation of unsteady free surface flows. Mathematical model is based on finite volume solution of the Navier-Stokes equations and depth-integrated continuity equation to track the free surface. The Sommerfeld boundary condition is applied at the far-end of the truncated computational domain. A dissipation zone is formed by applying artificial pressure on water surface to dissipate the kinetic ene...
Solution to transient Navier-Stokes equations by the coupling of differential quadrature time integration scheme with dual reciprocity boundary element method
Bozkaya, Canan; Tezer, Münevver (Wiley, 2009-01-20)
The two-dimensional time-dependent Navier-Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equati...
On efficient use of simulated annealing in complex structural optimization problems
Hasançebi, Oğuzhan (Springer Science and Business Media LLC, 2002-01-01)
The paper is concerned with the efficient use of simulated annealing (SA) in structural optimization problems of high complexity. A reformulation of the working mechanism of the Boltzmann parameter is introduced to accelerate and enhance the general productivity of SA in terms of convergence reliability. Two general and complementary parameters, referred as "weighted Boltzmann parameter" and "critical Boltzmann parameter," are proposed, Several alternative methodologies are suggested for these two parameter...
The DRBEM solution of incompressible MHD flow equations
Bozkaya, Nuray; Tezer, Münevver (Wiley, 2011-12-10)
This paper presents a dual reciprocity boundary element method (DRBEM) formulation coupled with an implicit backward difference time integration scheme for the solution of the incompressible magnetohydrodynamic (MHD) flow equations. The governing equations are the coupled system of Navier-Stokes equations and Maxwell's equations of electromagnetics through Ohm's law. We are concerned with a stream function-vorticity-magnetic induction-current density formulation of the full MHD equations in 2D. The stream f...
Citation Formats
Ş. Tiğrek and B. Yılmaz, “Implementation of physical boundary conditions into computational domain in modelling of oscillatory bottom boundary layers,” INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, pp. 973–991, 2010, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/57129.