• Wind and Structures
• /
• v.20 no.5
• /
• pp.661-682
• /
• 2015

This study investigates the use of time-frequency coherence analysis for detecting and evaluating coherent "structures" of surface pressures and wind turbulence components, simultaneously on the time-frequency plane. The continuous wavelet transform-based coherence is employed in this time-frequency examination since it enables multi-resolution analysis of non-stationary signals. The wavelet coherence quantity is used to identify highly coherent "events" and the "coherent structure" of both wind turbulence components and surface pressures on rectangular prisms, which are measured experimentally. The study also examines, by proposing a "modified" complex Morlet wavelet function, the influence of the time-frequency resolution and wavelet parameters (i.e., central frequency and bandwidth) on the wavelet coherence of the surface pressures. It is found that the time-frequency resolution may significantly affect the accuracy of the time-frequency coherence; the selection of the central frequency in the modified complex Morlet wavelet is the key parameter for the time-frequency resolution analysis. Furthermore, the concepts of time-averaged wavelet coherence and wavelet coherence ridge are used to better investigate the time-frequency coherence, the coherently dominant events and the time-varying coherence distribution. Experimental data derived from physical measurements of turbulent flow and surface pressures on rectangular prisms with slenderness ratios B/D=1:1 and B/D=5:1, are analyzed.

• 한국연소학회:학술대회논문집
• /
• 1998.10a
• /
• pp.23-30
• /
• 1998

The stabilization characteristics of diffusion flame formed in the wake of a cylindrical flame holder were investigated. Distribution of turbulence intensity, concentration distribution of combustion gas, and ion currents were measured. The turbulence intensity in the wake of cylindrical- game holder is increased with increase of diameter or blockage ratio of grid. If the auxiliary fuel is injected into recirculation zone, the concentration of $C_3H_8$ is high, but the concentration of $CO_2$ is low at the boundary of recirculation zone. The region with highest average value of ion currents in the middle of flame is moved to the upstream side by the turbulent components of main stream. The flame mass with partially active reaction is moved fast for uniform flow and turbulence generator G3, but the flame mass with relatively slow reaction is moved slowly for turbulence generator G1.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.118-131
• /
• 1991

The effects of free surface on vortex shedding phenomena around a bluff body were studied by both numerical simulation and flow visualization experiments. A vortex method, which approximates the vorticity field as the sum of discrete vortices; was used for the numerical simulation. Flow visualization experiments were performed in the KRISO cavitation tunnel. Hydrogen bubble was used as illumination material. Free surface elevation was also measured during experiments. The hydrodynamic drag and lift were predicted by numerical simulation. The predicted period of vortex shedding was compared with the results of experiments.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.4
• /
• pp.391-398
• /
• 2005

In this investigation, turbulent premixed combustion and flame front propagation in a gas turbine combustion chamber is studied. Direct numerical simulation of turbulent reacting flows demands extremely high computational resources, especially in more complicated geometry. The alternative choice may be left for Large Eddy Simulation (LES) by which only large scales are solved directly. In combustion problems, capturing the large scales' behavior without solving the details of small scales is a difficult task. Using a transport equation for description of the flame front propagation and therefore avoiding the calculation of inner flame structure is the basic idea of this study. For this purpose. the so-called G-equation has been used by which any iso-level of the G variable provides the flame location. A comparison with the experiment indicates that the present method can predict a turbulent velocity field and also capture a instantaneous 3-dimensional flame structure.

• Wind and Structures
• /
• v.1 no.1
• /
• pp.111-125
• /
• 1998

This paper deals with the development of variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined element and efficiently used for the construction of a refined mesh without generating distorted elements. A modified Guassian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of the independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems in which the locations to be refined are changed in accordance with the dynamic distribution of velocity gradient, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• Wind and Structures
• /
• v.17 no.4
• /
• pp.361-377
• /
• 2013

Wall functions have been widely used in computational fluid dynamics (CFD) simulations and can save significant computational costs compared to other near-wall flow treatment strategies. However, most of the existing wall functions were based on the asymptotic characteristics of near-wall flow quantities, which are inapplicable in complex and non-equilibrium flows. A modified wall function is thus derived in this study based on flow over a plate at zero-pressure gradient, instead of on the basis of asymptotic formulations. Turbulent kinetic energy generation ($G_P$), dissipation rate (${\varepsilon}$) and shear stress (${\tau}_{\omega}$) are composed together as the near-wall expressions. Performances of the modified wall function combined with the nonlinear realizable k-${\varepsilon}$ turbulence model are investigated in homogeneous equilibrium atmosphere boundary layer (ABL) and flow around a 6 m cube. The computational results and associated comparisons to available full-scale measurements show a clear improvement over the standard wall function, especially in reproducing the boundary layer flow. It is demonstrated through the two case studies that the modified wall function is indeed adaptive and can yield accurate prediction results, in spite of its simplicity.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.5
• /
• pp.671-680
• /
• 2007

The objective of this work is to reduce drag on a bluff body within a viscous flow by applying suction or injection of fluid along the surface of the body. In addition to minimizing drag, the optimal solution tends to reduce boundary layer separation and flow recirculation. When discretized by finite elements, the optimal control problem can be posed as a large-scale nonlinearly-constrained optimization problem. The constraints correspond to the discretized form of the Navier-Stokes equations. Unfortunately, solving such large-scale problems directly is essentially intractable. We developed several Sequential Quadratic Programming methods that are tailored to the structure of the control problem. Example problems of laminar flow around an infinite cylinder in two dimensions are solved to demonstrate the methodology. We use these optimal control techniques to study the influence of number of suction/injection holes and location of holes on the resulting optimized flow. We compare the proposed SQP methods against one another, as well as against available methods from the literature, from the point of view of efficiency and robustness. The most efficient of the proposed methods is two orders of magnitude faster than existing methods.

• Wind and Structures
• /
• v.26 no.5
• /
• pp.279-292
• /
• 2018

Inspired by the energy harvesting eel, a flexible flag behind a D-shape cylinder in a uniform viscous flow was simulated by using the immersed boundary method (IBM) along with low-speed wind tunnel experimentation. The flag in the wake of the cylinder was strongly influenced by the vortices shed from the upstream cylinder under the vortex-vortex and vortex-body interactions. Geometric and flow parameters were optimized for the flexible flag subjected to passive flapping. The influence of length and bending coefficient of the flexible flag, the diameters (D) of the cylinder and the streamwise spacing between the cylinder and the flag, on the energy generation was examined. Constructive and destructive vortex interaction modes, unidirectional and bidirectional bending and the different flapping frequency were found which explained the variations in the energy of the downstream flag. Voltage output and flapping behavior of the flag were also observed experimentally to find a more direct relationship between the bending of the flag and its power generation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.4
• /
• pp.661-674
• /
• 2002

The objective of this study is to develop efficient numerical method to enable solution of optimal control problems of Navier-Stokes flows and to apply these technique to the problem of viscous drag minimization on a bluff body by controlling boundary velocities on the surface of the body. In addition to the industrial importance of the drag reduction problem, it serves as a model for other more complex flow optimization settings, and allows us to study, modify, and improve the behavior of the optimal control methods proposed here. The control is affected by the suction or injection of fluid on portions of the boundary, and the objective function represents the rate at which energy is dissipated in the fluid. This study shows how reduced Hessian successive quadratic programming method, which avoid converging the flow equations at each iteration, can be tailored to these problems.

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

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