• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.5
• /
• pp.1-9
• /
• 2017

The cryogenic nitrogen spray of a swirl injector has been numerically investigated using three dimensional LES turbulence model to analyze the dynamic characteristics under supercritical condition. To predict the precise nitrogen properties under supercritical condition, SRK equation of state, Chung's method for viscosity and thermal conductivity and Takahashi's correlation based on Fuller's theory for diffusion coefficient are implemented. The complex flow structures due to interaction between flow field and acoustic field are observed inside and outside the injector under supercritical condition. FFT, POD, and DMD techniques are employed to understand the coherent structures. By implementing the FFT, the dominant frequencies are identified inside and outside the injector. The coherent flow structures related to the dominant frequencies are visualized using the POD and DMD techniques. In addition, the DMD provides the damping coefficient which is related with the instability prediction.

• 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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.04a
• /
• pp.295-298
• /
• 2007

The effectiveness of passive control techniques for reducing the pressure oscillation generated in a supersonic cavity flow was investigated numerically and experimentally, respectively. The control device includes a sub-cavity installed in the upstream edge of a rectangular cavity. Time-dependent supersonic cavity flow characteristics with turbulent features were examined by using the three-dimensional, mass-averaged Navier-Stokes computation based on a finite volume scheme and large eddy simulation. The results show that the pressure oscillation near the trailing edge dominates overall time-dependent cavity pressure variations. Such an oscillation can be attenuated more significantly in the presence of the sub-cavity compared with the cavity without sub-cavity, and a larger sub-cavity leads to better control performance.

• Nuclear Engineering and Technology
• /
• v.48 no.4
• /
• pp.847-858
• /
• 2016

An overview is given on the work of the Laboratory of Nuclear Energy Systems at ETH, Zurich (ETHZ) and of the Laboratory of Thermal Hydraulics at Paul Scherrer Institute (PSI), Switzerland on tight-lattice bundles. Two-phase flow in subchannels of a tight triangular lattice was studied experimentally and by computational fluid dynamics simulations. Two adiabatic facilities were used: (1) a vertical channel modeling a pair of neighboring sub-channels; and (2) an arrangement of four subchannels with one subchannel in the center. The first geometry was equipped with two electrical film sensors placed on opposing rod surfaces forming the subchannel gap. They recorded 2D liquid film thickness distributions on a domain of $16{\times}64$ measuring points each, with a time resolution of 10 kHz. In the bubbly and slug flow regime, information on the bubble size, shape, and velocity and the residual liquid film thickness underneath the bubbles were obtained. The second channel was investigated using cold neutron tomography, which allowed the measurement of average liquid film profiles showing the effect of spacer grids with vanes. The results were reproduced by large eddy simulation + volume of fluid. In the outlook, a novel nonadiabatic subchannel experiment is introduced that can be driven to steady-state dryout. A refrigerant is heated by a heavy water circuit, which allows the application of cold neutron tomography.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.5
• /
• pp.101-112
• /
• 2013

Large eddy simulation was performed to investigate the effect of linearly increasing wall disturbance on the modification of turbulent characteristics of temperature field in the vicinity of the wall. It was noted that temperature variance increased monotonically whereas temperature dissipation decreased significantly, resulting in a noticeable reduction in both time and length-scales. A sudden drop in turbulent Prandtl number down to around 0.25 in the near-wall region indicated that the similarity between velocity and temperature fields decreases near the wall as a result of linear wall disturbance.

• 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.