Abstract:A two-dimensional (2D) finite-difference shallow water model based on a second-order hybrid type of total variation diminishing (TVD) approximate solver with a MUSCL limiter function was developed to model flooding and inundation problems where the evolution of the drying and wetting interface is numerically challenging. Both a minimum positive depth (MPD) scheme and a non-MPD scheme were employed to handle the advancement of drying and wetting fronts. We used several model problems to verify the model, including a dam break in a slope channel, a dam break flooding over a triangular obstacle, an idealized circular dam-break, and a tide flow over a mound. Computed results agreed well with the experiment data and other numerical results available. The model was then applied to simulate the dam breaking and flooding of Hsindien Creek, Taiwan, with the detailed river basin topography. Computed flooding scenarios show reasonable flow characteristics. Though the average speed of flooding is 6-7 m s 1 , which corresponds to the subcritical flow condition (Fr < 1), the local maximum speed of flooding is 14Ð12 m s 1 , which corresponds to the supercritical flow condition (Fr ³ 1Ð31). It is necessary to conduct some kind of comparison of the numerical results with measurements/experiments in further studies. Nevertheless, the model exhibits its capability to capture the essential features of dam-break flows with drying and wetting fronts. It also exhibits the potential to provide the basis for computationally efficient flood routing and warning information.

A direct-forcing immersed boundary method (DFIB) with both virtual force and heat source is developed here to solve Navier-Stokes and the associated energy transport equations to study some thermal flow problems caused by a moving rigid solid object within. The key point of this novel numerical method is that the solid object, stationary or moving, is first treated as fluid governed by Navier-Stokes equations for velocity and pressure, and by energy transport equation for temperature in every time step. An additional virtual force term is then introduced on the right hand side of momentum equations in the solid object region to make it act exactly as if it were a solid rigid body immersed in the fluid. Likewise, an additional virtual heat source term is applied to the right hand side of energy equation at the solid object region to maintain the solid object at the prescribed temperature all the time. The current method was validated by some benchmark forced and natural convection problems such as a uniform flow past a heated circular cylinder, and a heated circular cylinder inside a square enclosure. We further demonstrated this method by studying a mixed convection problem involving a heated circular cylinder moving inside a square enclosure. Our current method avoids the otherwise requested dynamic grid generation in traditional method and shows great efficiency in the computation of thermal and flow fields caused by fluid-structure interaction.

No abstract

No abstract

Vortex-induced vibration (VIV) is an important physical phenomenon as one design a riser or a cylindrical structure in ocean. As the riser or the cylindrical structure is adjacent to a seabed, the boundary effect on VIV is not fully understood yet. The direct-forcing immersed boundary (DFIB) method is used to investigate a two-degree-of-freedom VIV of a flexible supported circular cylinder adjacent to a plane boundary in this study. Furthermore, the effect of the VIV of cylinder on skin friction of the plane boundary is investigated. The effects of varying reduced velocity and gap ratio on VIV are discussed. Only a single vortex street is found when the cylinder is close to plane boundary. Hydrodynamic coefficients of the freely vibrating cylinder are analyzed in time and spectral domains. Furthermore, nearly round oval-shaped motion is observed as the so-called lock-in phenomenon occurs. The skin friction of the plane boundary is predicted by the DFIB model. Results show that the vibrating cylinder in the boundary layer flow can reduce the friction effectively. This proposed DFIB model can be useful for the investigation of VIV of the structures under the plane boundary effect even for a small gap between the cylinder and the boundary.