• Wind and Structures
• /
• v.30 no.4
• /
• pp.339-366
• /
• 2020

This study aims at modeling boundary layers (BLs) encountered in sparse and built environments (i.e. open, suburban and urban) at the subsonic Wind Tunnel (WT) at Ryerson University (RU). This WT has an insignificant turbulence intensity and requires a flow-conditioning system consisting of turbulence generating elements (i.e., spires, roughness blocks, barriers) to achieve proper turbulent characteristics. This system was developed and validated in the current study in three phases. In phase I, several Computational Fluid Dynamic (CFD) simulations of the tunnel with generating elements were conducted to understand the effect of each element on the flow. This led to a preliminary design of the system, in which horizontal barriers (slats) are added to the spires to introduce turbulence at higher levels of the tunnel. This design was revisited in phase II, to specify slat dimensions leading to target BLs encountered by tall buildings. It was found that rougher BLs require deeper slats and, therefore, two-layer slats (one fixed and one movable) were implemented to provide the required range of slat depth to model most BLs. This system only involves slat movement to change the BL, which is very useful for automatic wind tunnel testing of tall buildings. The system was validated in phase III by conducting experimental wind tunnel testingof the system and comparing the resulting flow field with the target BL fields considering two length scales typically used for wind tunnel testing. A very good match was obtained for all wind field characteristics which confirms accuracy of the system.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.10
• /
• pp.925-931
• /
• 2010

The paper describes a Large Eddy Simulation (LES) study of turbulent flow around a cavity. A series of three-dimensional cavities placed in a turbulent boundary layer are simulated at a Reynolds number of $1.0{\times}10^5$ by considering U and h, which represent the velocity at the top and the depth of the cavity, respectively. In order to obtain the appropriate solution for the filtered Navier-Stokes equation for incompressible flow, the computational mesh forms dense close to the wall of the cavity but relatively coarse away from the wall; this helps reduce computation cost and ensure rapid convergence. The Boussinesq hypothesis is employed in the subgrid-scale turbulence model. In order to determine the subgrid-scale turbulent viscosity, the Smagorinsky-Lilly SGS model is applied and the CFL number for time marching is set as 1.0. The results show the flow variations inside cavities of different sizes and shapes.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.05a
• /
• pp.210-213
• /
• 2010

Recent experimental data shows that the noticeable feature of irregular roughened spots on the fuel surface occurs during the combustion test. The generation of these unexpected patterns is likely to be resulted from the disturbed boundary layer due caused by wall blowing which is intended to simulate the process of fuel vaporization. LES without chemical reaction was conducted to investigate the flow characteristics at the near-fuel surface and the behavior of turbulent structures which is evolved by the wall blowing at the Reynolds number of 23,000. Cylindrical geometry was considered to get the most reality of the calculation results because real hybrid rocket motor is circular grain configuration. It was shown that the wall blowing pushed turbulent structures upwards making them tilted and this skewed displacement, in effect, left the foot prints of the structures on the surface. This change of kinematics may explain the formation of irregular isolated spots on the fuel surface observed in the experiment.

• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.354-359
• /
• 2013

일반적으로 건축물의 설계시 풍동 실험을 통한 풍환경의 평가를 수행하고 있으며, 이는 환경 영향 평가법에서 정한 건축 사업 시행 시 수반되어야 할 자연환경, 생활환경 그리고 사회경제환경의 영향 평가의 일환으로 실시되고 있다. 그러나, 풍동 실험의 경우 여러 가지 현실적 제약조건으로 설계와 실험의 피드백 (Feedback)이 원활하지 못하며, 특히 대상 건축물이 공장과 같이 대기 오염원이 되는 경우 실험은 더욱 어려운 형편이다. 이에 대한 보완책으로 전산 유체 역학을 이용한 건축물의 풍압 해석에 의한 풍하중 추정이나 인접 지형-지물의 영향을 고려한 건축물 주위의 풍환경 평가가 있다. 전산 모사에 의해 풍동 실험의 미비점을 보완하고, 보다 상세한 정보를 확보함으로써 건축물의 구조적 안전성의 증대와 환경 피해 감소를 기할 수 있다. 그러나 복잡한 지형-지물이나 건축물 주위의 풍환경에 대한 전산 모사는 주로 두 가지의 기술적 어려움을 수반하게 되다. 그 중 하나는 고정 경계면을 이루는 형상의 복잡성으로 인해 기존에 많이 이용하고 있는 Body-fitted 격자계를 이용하는 경우, 격자 생성 과정이 매우 복잡하고 어려울 뿐 만 아니라 생성된 격자가 주로 비정렬 (unstructured) 특성을 갖게 되어 수치해석 과정의 효율을 저하시키는 요인이 되며, 격자의 형상도 수치해석의 수렴성을 저하시키는 예가 많다. 다른 어려움으로 풍환경은 전형적인 난류 유동장으로서 난류의 전산 해석은 아직도 해결하지 못한 부분이 많다는 점이다. 이에 본 논문에서는 복잡한 지형-지물이나 건축물의 풍하중과 풍환경의 전산 모사 기술 확보를 위하여 수행중인 연구의 일환으로 물체 형상의 기하학적 복잡성의 극복을 위한 가상경계법 (Immersed Boundary Method)과 난류 유동장의 물리적 엄밀성을 높이기 위한 다와동 모사 (Large Eddy Simulation)을 이용한 물체 형상과 무관한 유동장 해석 기술 개발에 대하여 다루고자 한다. 먼저 최근에 유동 해석에 이용되는 방법인 가상경계법(IBM)은 물체를 포함한 전체 전산 영역을 직교 좌표계에 의해 이산화하고, 유동장내 존재하는 물체의 표면에서의 점착 조건을 만족시키기 위하여 지배 방정식에 적절한 외력을 추가로 고려하는 방법이다. 본 연구에서는 가상경계법을 이용하여 경계층에 위치한 건물 형상의 각진 물체 주위 사이에 형성되는 공동 내부의 비정상 유속 및 압력에 대한 전산 해석을 수행하고, 풍상측 전면에 형성되는 경계층에 의한 영향을 분석하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.3
• /
• pp.193-204
• /
• 2011

To develop more robust oil boom which is vulnerable to various failure mode under severe weather condition, highly accurate wave model is developed using Spatially filtered Navier-Stokes Eq., LDS (Lagrangian Dynamic Smagorinsky model) for residual stresses, SPH (Smoothed Particle Hydrodynamics). To clarify the hydraulic characteristics of floating type oil boom, we numerically simulate the behavior of oil spill around oil boom under very energetic progressive waves. At the first stage, we firmly anchored the oil boom, and then, allowed the excursion of the oil boom. It turns out that oil boom with skirt of enough length (longer than 30% of depth) effectively confines the oil spill even against very energetic waves. We can also observe obliquely descending vertical eddies between y = 1~2 m as horizontal vortices shedding at the interface of oil spill and water are diffused toward the bottom, which is believed to be the birth, growing and break-down of Kelvin-Helmholz wave.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.558-558
• /
• 2012

다수의 자연 하천은 유사 이동과 하안의 침식으로 인하여 사행 하천이 발생한다. 사행 하천에서의 이차류는 원심력, 편수위로 발생하는 중력에 의한 압력차, 그리고 난류로 인하여 발생하는 응력으로 인하여 형성되며 나선형 구조 형태로 표면 유속은 외부로 향하고 하천 바닥의 유속은 안으로 향하게 된다. 이러한 2차류의 형성은 주 흐름의 특성을 변형시킨다. 자연하천에서 2차류는 주 흐름의 15-25%의 크기를 가지고 있으나, 하상의 변화, 유사의 이동 등과 연관되므로 2차류의 영향을 정확히 해석하는 것은 수리학적으로 매우 중요하다. 본 연구에서는 사행수로에서 발생하는 2차류 거동을 수치모의를 통하여 수행하였다. 우선 2차류의 분석을 위하여 실험을 통한 결과물을 비교하였다. 자연 하천의 특성을 반영할 수 있도록 서일원(2006)이 수행한 S-자 형태의 실험 수로의 실험 결과를 분석하였다. 수치 모의를 위하여 3차원 전산유체역학 프로그램을 사용하여 사행수로의 2차원 유속 구조를 모의할 수 있도록 하였다. FLOW-3D 프로그램을 이용하여 실험 결과와 모의 결과를 비교할 수 있도록 하였으며 비교 후 보정을 실시하였다. 모의는 주로 LES (Large Eddy Simulation) 모형을 통해 이루어졌으며, 이를 통하여 실험에서 획득한 결과와 비슷한 유속구조 분포를 확인할 수 있었다. 보정 및 검증 후 수치 모의를 통한 유속 데이터를 이용하여 민감도 분석을 실행하였다. 이후로는 수로의 만곡부, 조도, 수심 등 인공수로의 조건을 변경하여 수치 모의를 수행하였다. 보정된 결과를 이용하여 추가적인 모의를 통한 유속 분포 구조의 비교가 이루어졌다. 이를 통하여 각 조건이 이차류의 크기에 미치는 영향을 확인할 수 있었으며, 모의를 통한 유속분포 결과는 대체적으로 실험을 통한 이차류의 연직분포 구조와 일치하였다.

• Journal of Aerospace System Engineering
• /
• v.13 no.1
• /
• pp.1-10
• /
• 2019

The flow characteristics in a stirred tank are very useful in a wide variety of industrial applications. Generally, the flow pattern, power consumption and mixing time in stirred vessels depend not only on the design of the impeller, but also on the tanks' geometry and internal structure. In this study, the analysis of an unstable and unsteady complicated flow characteristics generated by the interaction between the baffle shape and impeller were performed using the ANSYS FLUENT LES Turbulence Model. The study compared the predictions of CFD with the interaction between two types of rotating impellers (axial and radial flows) and the shapes of three baffles. The results of the comparison verified that the design model showed a relatively efficient trend in the mixing flow fields and characteristics around the impeller and baffles during agitation.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.6
• /
• pp.56-63
• /
• 2018

In low swirl combustors the flame is lifted above the nozzle to achieve balance between the flame speed and velocity field at the exit of the nozzle. Characterization of the flame liftoff height is important because it affects the stability of the combustor and degradation of the nozzle material. In experiments, a counter-intuitive trend of flame liftoff heights with respect to inlet velocities was observed. To elucidate the complicated flow field in a low swirl combustor having swirl vanes and a turbulence generator, a series of numerical simulations of non-reacting flows was conducted by varying the inlet velocity. The flow structures at the exit of the nozzle with respect to the inlet velocities are investigated to support the observation in the experiments.

• Wind and Structures
• /
• v.28 no.5
• /
• pp.285-298
• /
• 2019

For super long-span bridges, the aerodynamic forces induced by the flow passing the box girder should be considered carefully. And the Reynolds number sensitively of aerodynamic characteristics is one of considerable issue. In the study, a numerical study on the Reynolds number sensitivity of aerodynamic characteristic (flow pattern, pressure distribution and aerodynamic forces) of a twin-box girder were carried out using large eddy simulation (LES) with the dynamic Smagorinsky-Lilly subgrid model. The results show that the aerodynamic characteristics have strong correlation with the Reynolds number. At the leading edge, the flow experiences attachment, departure, and reattachment stages accompanying by the laminar transition into turbulence, causing pressure plateaus to form on the surface, and the pressure plateaus gradually shrinks. Around the gap, attributing that the flow experiences stages of laminar cavity flow, the wake with alternate shedding vortices, and turbulent cavity flow in sequence with an increase in the Reynolds number, the pressures around the gap vary greatly with the Reynold number. At the trailing edge, the pressure gradually recovers as the flow transits to turbulence (the flow undergoes wake instability, shear layer transition-reattachment station), In addition, at relative high Reynolds numbers, the drag force almost does not change, however, the lift force coefficient gradually decreases with an increase in Reynolds number.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.43-58
• /
• 2022

The BARC flow is studied via Direct Numerical Simulation at a relatively low turbulent Reynolds number, with focus on the geometrical representation of the leading-edge (LE) corners. The study contributes to further our understanding of the discrepancies between existing numerical and experimental BARC data. In a first part, rounded LE corners with small curvature radii are considered. Results show that a small amount of rounding does not lead to abrupt changes of the mean fields, but that the effects increase with the curvature radius. The shear layer separates from the rounded LE at a lower angle, which reduces the size of the main recirculating region over the cylinder side. In contrast, the longitudinal size of the recirculating region behind the trailing edge (TE) increases, as the TE shear layer is accelerated. The effect of the curvature radii on the turbulent kinetic energy and on its production, dissipation and transport are addressed. The present results should be contrasted with the recent work of Rocchio et al. (2020), who found via implicit Large-Eddy Simulations at larger Reynolds numbers that even a small curvature radius leads to significant changes of the mean flow. In a second part, the LE corners are fully sharp and the exact analytical solution of the Stokes problem in the neighbourhood of the corners is used to locally restore the solution accuracy degraded by the singularity. Changes in the mean flow reveal that the analytical correction leads to streamlines that better follow the corners. The flow separates from the LE with a lower angle, resulting in a slightly smaller recirculating region. The corner-correction approach is valuable in general, and is expected to help developing high-quality numerical simulations at the high Reynolds numbers typical of the experiments with reasonable meshing requirements.