• Wind and Structures
• /
• v.34 no.1
• /
• pp.59-72
• /
• 2022

In this work the numerical results of the flow around a 5:1 rectangular cylinder at Reynolds numbers 3 000 and 40 000, zero angle of attack and smooth incoming flow condition are presented. Implicit Large Eddy Simulations (ILES) have been performed with a high-order accurate spatial scheme and an implicit high-order accurate time integration method. The spatial approximation is based on a discontinuous Galerkin (dG) method, while the time integration exploits a linearly-implicit Rosenbrock-type Runge-Kutta scheme. The aim of this work is to show the feasibility of high-fidelity flow simulations with a moderate number of DOFs and large time step sizes. Moreover, the effect of different parameters, i.e., dimension of the computational domain, mesh type, grid resolution, boundary conditions, time step size and polynomial approximation, on the results accuracy is investigated. Our best dG result at Re=3 000 perfectly agrees with a reference DNS obtained using Nek5000 and about 40 times more degrees of freedom. The Re=40 000 computations, which are strongly under-resolved, show a reasonable correspondence with the experimental data of Mannini et al. (2017) and the LES of Zhang and Xu (2020).

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.107-113
• /
• 2008

The present study numerically simulates the flow and heat transfer characteristics of rib-induced secondary flow in a square cooling channel with staggered V-shaped ribs, extruded on both walls. The rib pitch-to-height ratio (p/h) varies from 2.8 to 10 with the rib-height-to-hydraulic diameter ration (h/$D_h$)of 0.07 and the Reynolds number of 50,000. Shear stress transport (SST) turbulence model is used as a turbulence model. Computational results show that complex secondary flow patterns are generated in the channel due to the snaking flow in the streamwise direction for all tested cases. In the range of p/h=5 to 10 the staggered V-shaped rib gives about 3 times higher heat transfer augmentation than the reference smooth pipe with high heat transfer on both front side and the area around the leading edge of the ribs, while the former cases give about 18 times higher streamwise pressure drop than the latter ones. However, for the thermal performances, based on the equal pumping power condition, the case of p/h=2.8 gives the best result among three cases, mainly due to relatively low streamwise pressure drop, although it gives relatively low heat transfer augmentation.

• Journal of computational fluids engineering
• /
• v.13 no.2
• /
• pp.27-41
• /
• 2008

An unstructured hybrid mesh flow solver has been developed for the simulation of three-dimensional steady and unsteady incompressible flow fields. The incompressible Navier-Stokes equations with an artificial compressibility method were discretized by using a node-based finite-volume method. For the unsteady time-accurate computation, a dual-time stepping method was adopted to satisfy a divergence-free flow field at each physical time step. An implicit time integration method with local time stepping was implemented to accelerate the convergence in the pseudo-time sub-iteration procedure. The one-equation Spalart-Allmaras turbulence model has been adopted to solve high-Reynolds number flow fields. The flow solver was parallelized to minimize the CPU time and to overcome the computational overhead. This method has been applied to calculate steady and unsteady flow fields around submarine configurations and a 3-D infinite cylinder. Validations were made by comparing the predicted results with those of experiments or other numerical results. It was demonstrated that the present method is efficient and robust for the prediction of steady and unsteady incompressible flow fields.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.3
• /
• pp.12-18
• /
• 2012

The Lattice Boltzmann Method has developed for solving the Boltzmann equation in Cartesian domains containing immersed boundaries of arbitrary geometrical complexity moving with prescribed kinematics. When a immersed boundaries are sweeping the fixed fluid node, refilling the node information in a vicinity of fluid nodes is one of the important issues in Lattice Boltzmann Method. In this study, we propose a simple refill algorithm for the particle distribution function based on a proper velocity, density and strain rate to enhance accuracy and stability of the method. The refill scheme based on a asymptotic analysis of LBGK model has improved accuracy than interpolation schemes. The proposed scheme in this study is validated by the simulations of an impulsively started rotating circular cylinder to investigate adaptability for fluid-structure interaction (FSI) problem. This refill scheme has improved stability and accuracy especially at high Reynolds number region.

• Journal of the Korean Society of Visualization
• /
• v.15 no.2
• /
• pp.33-40
• /
• 2017

Abstract In the present study, the characteristics of kinematics and dynamics in the fluttering motion of free-falling plates are investigated at Reynolds number of $10^5$. We record quasi-two-dimensional trajectories of free-falling plates with and without superhydrophobic coating using high-speed camera, and compute the drag and lift forces by trajectory analysis. Translational and angular velocities are modeled as harmonic functions with specific phase differences. In particular, periodic mass elevations near turning points are explained using the suggested models. At each turning point, a sudden drop in lift and a rapid increase in drag occur simultaneously due to fast increase in angle of attack. However, the lift is increased over the buoyancy-corrected weight of plate during gliding flight, resulting in periodic mass elevations near turning points. Superhydrophobicity is shown to increase lift but to reduce drag on a fluttering plate, resulting in the decrease of mean descent speed.

• Wind and Structures
• /
• v.21 no.1
• /
• pp.25-40
• /
• 2015

Oscillating motions of flexible cylinders are associated to some extent with the aerodynamic response of plants. Tandem motions of reeds with flexible stems in a colony are experimentally investigated using an array of flexible cylinders made of polydimethylsiloxane (PDMS). Consecutive images of flexible cylinders subjected to oncoming wind are recorded with a high-speed camera. To quantify oscillating motions, the average bending angle and displacement of flexible cylinders are evaluated using point-tracking method and spectral analysis. The tandem motions of flexible cylinders are closely related to the flow characteristics around the cylinders. Thus, the dynamic motions of a tandem arrangement of flexible cylinders are investigated with varying numbers of cylinders arranged in-line, numbers of cylinders in a group (behaving like a single body), and Reynolds numbers (Re). When the number of cylinders in a group increases, the damping effect caused by the support of downstream cylinders is pronounced. These results would be provide useful information on the tandem-arranged design of complex structures and energy harvesting devices.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.6
• /
• pp.480-488
• /
• 2014

Vortex shedding which is the dominant feature of body wakes and of direct relevance to practical engineering problems, has been intensively studied for flows past a circular cylinder. In contrast, vortex shedding from a hydrofoil trailing edge has been studied to much less extent despite numerous practical applications. The physics of the problem is still poorly understood. The present study deals with $K{\acute{a}}rm{\acute{a}}n$ vortex shedding from a truncated trailing-edge hydrofoil in relatively high Reynolds number flows. The objectives of this paper are twofold. First, we aim to simulate unsteady turbulent flows past a two dimensional hydrofoil through a hybrid particle-mesh method and penalization method. The vortex-in-cell (VIC) method offers a highly efficient particle-mesh algorithm that combines Lagrangian and Eulerian schemes, and the penalization method enables to enforce body boundary conditions by adding a penalty term to the momentum equation. The second purpose is to investigate shedding frequencies of vortices behind a NACA 0009 hydrofoil operating at a zero angle of attack.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2315-2328
• /
• 1993

The turbulence effect of particle deposition on a horizontal free-standing wafer in a vertical flow has been studied numerically by using the low-Reynolds-number k-.epsilon. turbulence model. For both the upper and lower surfaces of the wafer, predictions are made of the averaged particle deposition velocity and its radial distribution. Thus, it is now possible to obtain local information about the particle deposition on a free-standing wafer. The present result indicates that the particle deposition velocity on the lower surface of wafer is comparable to that on the upper one in the diffusion controlled deposition region in which the particle sizes are smaller than $0.1{\mu}m$. And it is found in this region that, compared to the laminar flow case, the averaged deposition velocity under the turbulent flow is about two times higher, and also that the local deposition velocity at the center of wafer is high equivalent to that the wafer edge.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.14 no.12
• /
• pp.996-1003
• /
• 2002

Heat transfer characteristics were investigated for the falling film flowing down the outside of an electrically heated vertical tube. Water was used for the falling film, and its Reynolds number was varied in the range of 70～500. Because water is heated and evaporated as it flows down, both sensible and latent heat transfer should be considered. The effect of the surrounding air movement was investigated by changing the direction of the air injection; without air injection, parallel-flow, and counter-flow. For all cases, sensible teat transfer rate was almost linearly increased with the increasing film flow rate. It was found that the film heat transfer coefficient was hardly influenced by the parallel air flow. However, the counter-flow of air reduced the heat transfer coefficient, which might be caused by the uneven distribution or flooding of the film. At high heat flux, a sudden change of the film heat transfer coefficient was detected as the film flow rate reached the transition value. It is supposed that this phenomenon was caused by the change in the film flow pattern.