• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.1
• /
• pp.37-50
• /
• 2016

A large eddy simulation (LES) of a turbulent channel flow is performed by using the third order low-storage Runge-Kutta method in time and second order finite difference formulation in space with staggered grid at a Reynolds number, $Re_{\tau}=590$ based on the channel half width, ${\delta}$ and wall shear velocity, $u_{\tau}$. To reduce the calculation cost of LES, algebraic wall model (AWM) is applied to approximate the near-wall region. The computation is performed in a domain of $2{\pi}{\delta}{\times}2{\delta}{\times}{\pi}{\delta}$ with $32{\times}20{\times}32$ grid points. Standard Smagorinsky model is used for subgrid-scale (SGS) modeling. Essential turbulence statistics of the flow field are computed and compared with Direct Numerical Simulation (DNS) data and LES data using no wall model. Agreements as well as discrepancies are discussed. The flow structures in the computed flow field have also been discussed and compared with LES data using no wall model.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.1636-1640
• /
• 2006

SMAC method, one of the computational fluid dynamics techniques, is modified from the original MAC method for the time-dependent variation of flow analysis. The Navier-Stokes equations for incompressible time-dependent viscous flow are applied, and also marker particles that present the visualization of flow analysis are used. In this study SMAC technique is used to analyze the flow behavior in the water-filling of virtual channel. Then by changes of diameter of two ingates, the change of velocity and discharge when two ingates are filled the water to virtual channel are simulated. As a result, water-filling flow pattern in the virtual channel is simulated very well. Therefore, this numerical simulation will also be applied for the design of structures as open flume and porous breakwater.

• Journal of Ocean Engineering and Technology
• /
• v.8 no.1
• /
• pp.105-113
• /
• 1994

Double-diffusive convection induced by simultaneously-imposed lateral temperature and concentration gradients in a rectangular enclosure with aspect retio, 2.0 has been studied experimentally for adiabatic and isothermal horizontal boundary conditions. Visual observations show two distinct flow structures depending on the buoyancy ratio. The unicell flow structure is observed for a lower buoyancy ratio while the layered flow structure appears for a higher buoyancy ratio. There exists an unstable flow regime between two buoyancy ratios.

• Korea-Australia Rheology Journal
• /
• v.12 no.3_4
• /
• pp.143-149
• /
• 2000

In this article, we considered shear-induced microstructure and rheological behavior of micellar solutions of cationic surfactant, cetylpyridinium chloride (CPC) in the presence of a structure-forming additive, sodium salicylate (NaSal). Shear viscosity, shear moduli and flow birefringence were measured as functions of the surfactant and additive concentrations. In the presence of NaSal, the micellar solution exhibited the non-linear rheological behavior due to the formation of supramolecular structures when the molar ratio of NaSal to CPC exceeded a certain threshold value. Flow birefringence probed the change in micelle alignment under shear flow. At low shear rates, the flow birefringence increased as the shear rate increased. On the other hand, fluctuation of flow birefringence appeared from the shear rate near the onset of shear thickening, which was caused by shear-induced coagulation or aggregation. These results were confirmed by the SEM images of in situ gelified micelle structure through sol-gel route.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.100-109
• /
• 2011

Large eddy simulation(LES) of fully developed turbulent pipe flow has been performed to investigate the effect of Reynolds number on the flow field at $Re_{\tau}$=180, 395, 590 based on friction velocity and pipe radius. A dynamic subgrid-scale model for the turbulent subgrid-scale stresses was employed to close the governing equations. The mean flow properties, mean velocity profiles and turbulent intensities obtained from the present LES are in good agreement with the previous numerical and experimental results currently available. The Reynolds number effects were observed in the higher-order statistics(Skewness and Flatness factor). Furthermore, the budgets of the Reynolds stresses and turbulent kinetic energy were computed and analyzed to elucidate the effect of Reynolds number on the turbulent structures.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.86-103
• /
• 2006

Using a combination of force transducer measurement to quantify net lift force, high frame rate camera to quantify and subtract inertial contributions, and Digital Particle Image Velocimetry (DPIV) to calculate aerodynamic contributions in the spanwise plane, the contribution of spanwise flow to the generation of lift force in wings undergoing a pure flapping motion in hover is shown as a function of flapping angle throughout the flapping cycle. These experiments were repeated at various flapping frequencies and for various wing planform sizes for flat plate and span wise cambered wings. Despite the previous identification of the importance of span wise fluid structures in the generation of lift force in flapping wings throughout the existing body of literature, the direct contribution of spanwise flow to lift force generated has not previously been quantified. Therefore, in the same manner as commonly applied to investigate the chordwise lift distribution across an airfoil in flapping wings, spanwise flow due to bulk flow and rotational fluid dynamic mechanisms will be investigated to validate the existence of a direct component of the lift force originating from the flapping motion in the spanwise plane instead.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.359-368
• /
• 2002

This paper describes a CFD investigation into the flow over the cab of a bluff-fronted lorry. Several different simulations were undertaken, using the commercial codes: CFX, Fluent and PowerFLOW. Using the $k-{\varepsilon}$ turbulence model, the flow over the cab was symmetric, however, using more accurate turbulence models such as the RNG $k-{\varepsilon}$ model or the Reynolds Stress Model, the flow was asymmetric. The paper discusses whether this phenomenon is a real effect or whether it is a solver artefact and the study is supported by experimental evidence. The findings are preliminary, but suggest that it has a physical origin and that it may be aspect ratio-dependent.

• Journal of the Korean Society of Combustion
• /
• v.13 no.3
• /
• pp.33-40
• /
• 2008

Two- and Three- dimensional numerical simulations are carried out to understand the combustion characteristics of LNG-fueled gas turbine combustor for power generation using imported and domestic natural gases. Reacting flow characteristics of the swirl stabilized natural gas combustor were understood from the numerical results with the flow conditions selected from the gas turbine operation data. The thermal influences of different natural gases were very small and the fuel composition and flow rate were considered to be tuned well. The flow structures of the recirculation and combustion region was understood from the comparison of the two- and three-dimensional results. The complexity of the three-dimensional swirl flows inside the gas turbine combustor with multiple swirlers was understood those resulting from the interactions of the stage and pilot burners.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.11-16
• /
• 2002

In this paper, we present two successful results from active controls of flows over a circular cylinder and a sphere for drag reduction. The Reynolds number range considered for the flow over a circular cylinder is 40-3900 based on the free-stream velocity and cylinder diameter, whereas for the flow over a sphere it is $10^{5}$ based on the free-stream velocity and sphere diameter. The successful active control methods are a distributed (spatially periodic) forcing and a high-frequency (time periodic) forcing. With these control methods, the mean drag and lift fluctuations decrease and vortical structures are significantly modified. For example, the time-periodic forcing at a high frequency (larger than 20 times the vortex shedding frequency) produces $50{\%}$ drag reduction for the flow over a sphere at $Re=10^{5}$. The distributed forcing applied to the flow over a circular cylinder results in a significant drag reduction at all the Reynolds numbers investigated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.11
• /
• pp.827-833
• /
• 2009

The present study numerically investigates three-dimensional laminar flow past a rotating circular cylinder placed in a uniform stream. For the purpose of a careful analysis of the modification of flow by the effect of the rotation on the flow, numerical simulations are performed at a various range of rotational coefficients ($0{\leq}{\alpha}{\leq}2.5$) at one Reynolds number of 300. As ${\alpha}$ increases, flow becomes stabilized and finally a steady state beyond the critical rotational coefficient. The 3D (three dimensional) wake mode of the stationary cylinder defined at this Reynolds number has been disorganized according to ${\alpha}$, which were observed by the visualization of 3D vortical structures. The variation of the Strouhal number is very weak when the wake pattern is changed according to the rotational coefficient. As ${\alpha}$ increases, the lift increases, whereas the drag decreases.