• 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)은 물체를 포함한 전체 전산 영역을 직교 좌표계에 의해 이산화하고, 유동장내 존재하는 물체의 표면에서의 점착 조건을 만족시키기 위하여 지배 방정식에 적절한 외력을 추가로 고려하는 방법이다. 본 연구에서는 가상경계법을 이용하여 경계층에 위치한 건물 형상의 각진 물체 주위 사이에 형성되는 공동 내부의 비정상 유속 및 압력에 대한 전산 해석을 수행하고, 풍상측 전면에 형성되는 경계층에 의한 영향을 분석하였다.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2386-2396
• /
• 1996

A numerical simulation has been performed for isothermal and reacting flows in an exisymmetric, bluff-body research combustor. The present formulation is based on the density-weighted averaged Navier-Stokes equations together with a k-epsilon. turbulence model and a modified eddy-breakup combustion model. The PISO algorithm is employed for solution of thel Navier-Stokes system. Comparison between measurements and predictions are made for a centerline axial velocities, location of stagnation points, strength of recirculation zone, and temperature profile. Even though the numerical simulation gives acceptable agreement with experimental data in many respects, the present model is defictient in predicting the recoveryt rate of a central near-wake region, the non-isotropic turbulence effects, and variation of turbulent Schmidt number. Several possible explanations for these discrepancies have been discussed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.8
• /
• pp.1-7
• /
• 2003

The numerical simulation of flow-filed around a square cylinder near a wall with $\varepsilon$-SST turbulence model is carried out in this study. The newly suggested $\varepsilon$-SST turbulence model that modifies the original SST turbulence model is proved to yield more accurate results than the other 2-equation turbulence models in large separation region around a bluff body. Therefore, $\varepsilon$-SST turbulence model can be effectively applied for predicting the flow-fields with large separation. And it is found that vortex shedding is suppressed below the critical gap height, the Strouhal number is affected by the gap height and the wall boundary layer thickness.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.69-69
• /
• 2011

본 논문에서는 비정상 상태의 비압축성 유동장을 해석하기 위하여 물체맞춤격자방법이 아닌 가상경계법을 사용하였다. 가상경계법은 구조격자를 사용하여 구조물 경계면에서 Momentum Forceing을 사용하여 가상의 경계를 만들어 유동장을 해석하는 방법이다. Navier-Stoke 방정식의 수치 이산화 방법으로 Kim et al(1985)이 사용한 Fractional Step Method(FSM)을 사용하였다. 시간에 대하여 semi-implicit FSM를 사용하였고, 확산항에 대해서는 2차 정확도의 Crank-Nicolson Method를 대류항은 3차 정확도의 Runge-Kutta Method를 사용 하였다. 본 연구에서는 가상경계법을 이용한 유동장 해석이 교량 단면에 대하여 수치해석이 가능한지 검토하였다. 가상경계법은 현재 많은 연구가 유선형의 구조물에 대하여 수행되어 오고 있다. 교량 단면과 같은 각 진 구조물에 대한 검토는 아직 미비한 실정이다. 가상경계법에서 다루고 있는 구조물 경계면에서의 Momentum Forcing 방법이 유선형의 구조물에 맞추어 연구가 진행되었기 때문이다. 먼저 본 연구의 프로그램을 검증하기 위하여 원형 실린더에 대하여 가상경계법을 적용한 결과 Re 수 200에서 Strouhal Number, 양력계수, 항력계수를 이전 연구 결과와 비교하였다. Williamson(1988)과 Zhang(1995)의 연구결과와 유사한 결과를 얻을 수 있다. 그리고 교량의 단면과 같은 각진 구조물(Bluff Body)에 대하여 가상경계법 적용하였다. 본 논문의 연구에서 평가 대상으로 하고 있는 2차원 교량 단면에 대하여 유동장 해석을 하였다. 본 논문에서 정량적인 유체력과 유동장에 대한 비교 및 검토가 이루어지지 못했지만 압력장과 유선의 형태가 이론적인 값을 벗어나지 않고 있는 것으로 확인 되었다. Re 수 2700에서 전산 해석을 수행하였으며, 교량 단면 주위의 압력계수와 박리현상 그리고 후류에서의 Vortex shedding 현상이 모두 적절한 분포가 나타나는 것을 확인할 수 있었다. 따라서 가상경계법을 이용하여 각진 구조물에 대한 주위 유동장해석에 대한 가능성을 확인하였으며, 풍동실험과의 결과비교를 통하여 가상경계법을 이용하여 교량 단면 주위의 유동장 해석 결과를 정량적으로 비교할 것이다.

• Wind and Structures
• /
• v.21 no.2
• /
• pp.223-239
• /
• 2015

Open decks are a widely used deck configuration in long-span cable-stayed bridges; however, incorporating aerodynamic countermeasures are advisable to achieve better aerodynamic performance than a bluff body deck alone. A sectional model of an open deck cable-stayed bridge with a main span of 400 m was selected to conduct a series of wind tunnel tests. The influences of five practical aerodynamic countermeasures on flutter and vortex-induced vibration (VIV) performance were investigated and are presented in this paper. The results show that an aerodynamic shape selection procedure can be used to evaluate the flutter stability of decks with respect to different terrain types and structural parameters. In addition, the VIV performance of $\prod$-shaped girders for driving comfortableness and safety requirements were evaluated. Among these aerodynamic countermeasures, apron boards and wind fairings can improve the aerodynamic performance to some extent, while horizontal guide plates with 5% of the total deck width show a significant influence on the flutter stability and VIV. A wind fairing with an angle of $55^{\circ}C$ showed the best overall control effect but led to more lock-in regions of VIV. The combination of vertical stabilisers and airflow-depressing boards was found to be superior to other countermeasures and effectively boosted aerodynamic performance; specifically, vertical stabilisers significantly contribute to improving flutter stability and suppressing vertical VIV, while airflow-depressing boards are helpful in reducing torsional VIV.

• Korea-Australia Rheology Journal
• /
• v.21 no.1
• /
• pp.27-37
• /
• 2009

In this work, we numerically investigate the effect of viscoelasticity on 2D laminar vortex dynamics in flows past a single rotating cylinder for rotational rates $0{\leq}{\alpha}{\leq}5$ (the rotational rate ex is defined by the ratio of the circumferential rotating velocity to free stream velocity) at Re=100, in which the vortex shedding has been predicted to occur in literature for Newtonian fluids. The objective of the present research is to develop a promising technique to fully suppress the vortex shedding past a bluff body by rotating a cylinder and controlling fluid elasticity. The predicted vortex dynamics with the present method is consistent with the previous works for Newtonian flows past a rotating cylinder. We also verified our method by comparing our data with the literature in the case of viscoelastic flow past a non-rotating cylinder. For $0{\leq}{\alpha}{\leq}1.8$, the frequency of vortex shedding slightly decreases but the fluctuation of drag and lift coefficient significantly decreases with increasing fluid elasticity. We observe that the vortex shedding of viscoelastic flow disappears at lower ${\alpha}$ than the Newtonian case. At ${\alpha}$=5, the relationship between the frequency of vortex shedding and Weissenberg number (Wi) is predicted to be non-monotonic and have a minimum around Wi=0.25. The vortex shedding finally disappears over critical Wi number. The present results suggest that the vortex shedding in the flow around a rotating cylinder can be more effectively suppressed for viscoelastic fluids than Newtonian fluids.

• Wind and Structures
• /
• v.23 no.2
• /
• pp.127-142
• /
• 2016

Aerodynamic effects, such as drag force and flow-induced vibration (FIV), on civil engineering structures can be minimized by optimally modifying the structure shape. This work investigates the turbulent wake of a square prism with its faces modified into a sinusoidal wave along the spanwise direction using three-dimensional large eddy simulation (LES) and particle image velocimetry (PIV) techniques at Reynolds number $Re_{Dm}$ = 16,500-22,000, based on the nominal width ($D_m$) of the prism and free-stream velocity ($U_{\infty}$). Two arrangements are considered: (i) the top and bottom faces of the prism are shaped into the sinusoidal waves (termed as WSP-A), and (ii) the front and rear faces are modified into the sinusoidal waves (WSP-B). The sinusoidal waves have a wavelength of $6D_m$ and an amplitude of $0.15D_m$. It has been found that the wavy faces lead to more three-dimensional free shear layers in the near wake than the flat faces (smooth square prism). As a result, the roll-up of shear layers is postponed. Furthermore, the near-wake vortical structures exhibit dominant periodic variations along the spanwise direction; the minimum (i.e., saddle) and maximum (i.e., node) cross-sections of the modified prisms have narrow and wide wakes, respectively. The wake recirculation bubble of the modified prism is wider and longer, compared with its smooth counterpart, thus resulting in a significant drag reduction and fluctuating lift suppression (up to 8.7% and 78.2%, respectively, for the case of WSP-A). Multiple dominant frequencies of vortex shedding, which are distinct from that of the smooth prism, are detected in the near wake of the wavy prisms. The present study may shed light on the understanding of the underlying physical mechanisms of FIV control, in terms of passive modification of the bluff-body shape.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.296-307
• /
• 2020

Permanent magnet sodium flowmeter (PMSF) have been used to measure the sodium flow in fast breeder reactors. Due to the effects of irradiation, thermal cycling, time lapse, etc., the magnetic flux density of the PMSF will decrease after being used in the reactor for a period of time. Therefore, it must be calibrated regularly. But some flowmeters that immersed in sodium cannot be removed for an off-line calibration, so the on-line calibration is required. However, the best online calibration accuracy of PMSF using cross-correlation analysis method was 2.0-level without considering the repeatability. In order to further improve this work, the operational principle of the transducer in PMSF is analyzed and the design principle of the transducer is proposed. The transducers were tested on the sodium flow loop to collect the experimental data. The signal characteristics are analyzed from the time and frequency domains, respectively. The cross-correlation analysis method based on biased estimation is adopted to obtain the flow rate. The verification experimental results showed that the measurement accuracy is 1.0-level when the flow velocity is above 0.5 m/s, and the measurement accuracy is 3.0-level when the flow velocity is in the range of 0.2 m/s to 0.5 m/s.

• Wind and Structures
• /
• v.30 no.2
• /
• pp.181-198
• /
• 2020

In the traditional buffeting response analysis method, the spanwise incomplete correlation of buffeting forces is always assumed to be same as that of the incident wind turbulence and the action of the signature turbulence is ignored. In this paper, three typical bridge decks usually adopted in the real bridge engineering, a single flat box deck, a central slotted box deck and a two-separated paralleled box deck, were employed as the investigated objects. The wind induced pressure on these bridge decks were measured via a series of wind tunnel pressure tests of the sectional models. The influences of the wind speed in the tests, the angle of attack, the turbulence intensity and the characteristic distance were taken into account and discussed. The spanwise root coherence of buffeting forces was also compared with that of the incidence turbulence. The signature turbulence effect on the spanwise root coherence function was decomposed and explained by a new empirical method with a double-variable model. Finally, the formula of a sum of rational fractions that accounted for the signature turbulence effect was proposed in order to fit the results of the spanwise root coherence function. The results show that, the spanwise root coherence of the drag force agrees with that of incidence turbulence in some range of the reduced frequency but disagree in the mostly reduced frequency. The spanwise root coherence of the lift force and the torsional moment is much larger than that of the incidence turbulence. The influences of the wind speed and the angle of attack are slight, and they can be ignored in the wind tunnel test. The spanwise coherence function often involves several narrow peaks due to the signature turbulence effect in the high reduced frequency zone. The spanwise coherence function is related to the spanwise separation distance and the spanwise integral length scales, and the signature turbulence effect is related to the deck-width-related reduced frequency.