• 한국태양에너지학회 논문집
• /
• 제29권3호
• /
• pp.73-82
• /
• 2009

포화상태 핵비등과 저 Re수의 흐름비등에서 얻어진 실험결과를 바탕으로 하여, 기포가 성장하는 동안의 등가 기포 직경과 열전단율의 거동에 대한 기포 형상 가정의 효과를 제시하기 위한 해석적인 연구를 수행하였다. 이러한 목적을 달성하기 위하여, 등가 기포 반경이 기포가 성장하는 동안 촬영된 기포의 이미지로부터 얻어질 수 있는 형상 가정을 이용하여 계산되었다. 그리고 열전달율을 포화상태 핵비등 동안 미세크기의 히터와 휘스톤브리지 회로를 이용하여 측정하였다. 그리고, 기포 형상 가정의 효과를 실험결과와 비교하였고, 이를 통해 단일 기포의 성장 거동을 분석하기 위한 기포 형상 가정이 매우 중요함을 보였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.133-140
• /
• 1998

A Navier-Stokes code based on a unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number $k-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for a domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. Unsteadiness inside boundary layer is entrained when a unsteady vortex impinge on the blade surface. It shoves that local peak value inside the boundary layer and also local minimum near the edge of boundary layer as it developes along the blade surface. The unsteadiness inside the boundary layer is almost isolated from the free stream unsteadiness and being convected at local boundary layer speed, less than the free stream value.

• International Journal of Aeronautical and Space Sciences
• /
• 제14권3호
• /
• pp.201-209
• /
• 2013

Much numerical and experimental research has been done for the flow around an oscillating airfoil. The main research topics are vortex shedding, dynamic stall phenomenon, MAV's lift and thrust generation. Until now, researches mainly have been concentrated on analyzing the wake flow for the variation of frequency and amplitude at a low angle of attack. In this study, wake structures and acoustic wave propagation characteristics were studied for a plunging airfoil at high angle of attack. The governing equations are the Navier-Stokes equation with LES turbulence model. OHOC (Optimized High-Order Compact) scheme and 4th order Runge-Kutta method were used. The Mach number is 0.3, the Reynolds number is, and the angle of attack is from $20^{\circ}$ to $50^{\circ}$. The plunging frequency and the amplitude are from 0.05 to 0.15, and from 0.1 to 0.2, respectively. Due to the high resolution numerical method, wake vortex shedding and pressure wave propagation process, as well as the propagation characteristics of acoustic waves can be simulated. The results of frequency analysis show that the flow has the mixed characteristics of the forced plunging frequency and the vortex shedding frequency at high angle of attack.

• 대한의용생체공학회:의공학회지
• /
• 제33권4호
• /
• pp.177-183
• /
• 2012

The inspiratory flow rate of a human is changed with the amount of the workload. The flow characteristic is affected by the inspiratory flow rate. In the flow field of airway, the both of turbulence intensity and secondary flow affect the deposition pattern of particles which is important for the drug-aerosol targeting. Thus the analysis of the flow characteristic in a human airway is important. The purpose of this study is to investigate the effects of the inspiratory flow rate on the flow characteristics in a human airway. The tubular airway is consistent with the oral cavity, pharynx, larynx and trachea. The relatively inspiratory flow rate is used at each case of human states regarding the workload. By the effect of geometric airway changes, transition to turbulent airflow after the larynx can occur with relaminarization further downstream. The low Reynolds number k-${\omega}$ turbulence model is used for analysis with flow regime. As the inspiratory flow rate is larger, the turbulence kinetic energy and secondary flow intensity increase in airway. On the other hand, the area of recirculation zone is smaller.

• 한국가시화정보학회:학술대회논문집
• /
• 한국가시화정보학회 2006년도 추계학술대회 논문집
• /
• pp.7-10
• /
• 2006

An experimental study has been conducted to quantitatively characterize the motion of neutrally buoyant particles in 2-dimensional Poiseuille flow through the micron-sized circular capillaries in the range of Re (Reynolds number) $\approx0.1\sim100$. $A{\mu}-PTV$ (Particle Tracking Velocimetry) system is adopted, which consists of a double-headed Nd:YAG laser, an epi-fluorescence microscope and a cooled CCD camera. Since high shear rate can be induced due to the scale effect even at low Re, it is shown that in micro scale neutrally buoyant particles in Poiseuille flow drift away from the wall and away from the center of the capillary. Consequently, particles accumulate at the equilibrium position of $0.52\sim0.64R$ with R being the radius of the capillary, which is analogous to that of tube flow in macro scale. There is a plateau in equilibrium position at small Re, while equilibrium position starts increasing at $Re\approx30$. The outermost edge of particle cluster is closer to the center of the capillary than that in previous studies due to low Re effect. The present study quantitatively presents characteristics of particle motion in circular capillaries. Furthermore, it is expected to give optimum factors for designing microfluidic systems that are to be used fur plasma separation from the blood.

• 대한기계학회논문집B
• /
• 제20권12호
• /
• pp.4013-4026
• /
• 1996

A numerical simulation has been carried out for the jet impinging on a flat plate and a semi-circular concave surface. In this computation finite volume method was employed to solve the full Navier-Stokes equation based on a non-orthogonal coordinate with non staggered variable arrangement. The standard k-.epsilon. turbulent model and low Reynolds number k-.epsilon. model(Launder-Sharmar model) with Yap's correction were adapted. The accuracy of the numerical calculations were compared with various experimental data reported in the literature and showed good predictions of centerline velocity decay, wall pressure distribution and skin friction. For the jet impingement on a semi-circular concave surface, potential core length was calculated for two different nozzle(round edged nozzle and rectangular edged nozzle) to consider effects of the nozzle shape. The result showed that round edged nozzle had longer potential core length than rectangular edged nozzle for the same condition. Heat transfer rate along the concave surface with constant heat flux was calculated for various nozzle exit to surface distance(H/B) in the condition of same jet velocity. The maximum local Nusselt number at the stagnation point occurred at H/B = 8 where the centerline turbulent intensity had maximum value. The predicted Nusselt number showed good agreement with the experimental data at the stagnation point. However heat transfer predictions along the downstream were underestimated. This results suggest that the improved turbulence modeling is required.

• 한국추진공학회지
• /
• 제17권5호
• /
• pp.54-59
• /
• 2013

현재 Micro Shock Tube는 다양한 공학응용분야에 적용되고 있으며, 특히 우주항공 및 연소기술 그리고 약물전달 등의 분야에서 광범위한 잠재력을 가진 장치 중 하나이다. 그러나 Micro Shock Tube에서의 유동 특성은 작은 직경으로 인해 형성되는 매우 낮은 Reynolds Number와 높은 Knudsen Number의 영향으로 일반적으로 잘 알려진 Macro Shock Tube의 유동 특성과 상이하게 나타난다. 본 연구에서는 이러한 Micro Shock Tube의 유동 특성을 상세히 연구하기 위해 직경이 다른 두 가지 Micro Shock Tube의 실험을 수행하였다. 충격파 전파를 측정하기 위해 고압관의 파막압력 그리고 저압관의 세 지점에 센서를 설치하여 압력을 측정하고 분석하였다. 본 연구로부터, 동일한 파막압력에서 Micro Shock Tube 직경의 증가에 따라 충격파 전파속도가 증가하였고, 반사파의 영향도 더 크게 받았다.

• 한국전산유체공학회지
• /
• 제19권4호
• /
• pp.1-7
• /
• 2014

A numerical study of transient conjugate heat transfer on micro heater in a micro-channel substrate under a sinusoidal heat load was conducted. It was found that the time constant is not affected by the maximum heating magnitude of the sinusoidal heat load. However, the time constant increases with low duration of the sinusoidal heating period and low Reynolds number. Moreover, there is a threshold where a heater temperature do not reach to time constant at low thermal diffusivity, low flow rate, and low pulse duration of the sinusoidal heating. The time constant should be considered for transient convective heat transfer under transient sinusoidal heat load in a micro heat sink.

• 한국전산유체공학회지
• /
• 제18권2호
• /
• pp.49-55
• /
• 2013

The goal of the research is to evaluate the open source code of OpenFOAM for the use of nuclear plant flow simulation objectively. Of the various incompressible flow solvers, simpleFoam, pimpelFoam are then tested under three validated cases (backward facing step, flow over circular cylinder and turbulent round jet flow). For the evaluation of steady state incompressible laminar flow simulation, low reynolds number of backward facing step flow was solved by simpleFoam. The resultant of the reattached lengths turned out to be similar with the other experimental and simulation results. For transient flow simulation, flow over circular cylinder and turbulent round jet flow were solved by pimpleFoam. The simulation accuracy was evaluated by comparing the resultant flow patterns with the description of the characteristics of the flow over the circular cylinder. The quantitative accuracy was evaluated for no more than 85% by comparing it to the decaying constants of the turbulent round jet velocity.

• 대기
• /
• 제27권4호
• /
• pp.445-453
• /
• 2017

This study investigates turbulent flow around a square cylinder mounted on a flat surface at high Reynolds number using a large-eddy simulation (LES) model, particularly focusing on vortex streets behind the square cylinder. Total 9 simulation cases with different inflow wind directions, inflow wind speeds, and cylinder widths in the x- and y-directions are considered to examine the effects of inflow wind direction, speed, and cylinder widths on turbulent flow and vortex streets. In the control case, the inflow wind parallel to the x-direction has a maximum speed of $5m\;s^{-1}$ and the width and height of the cylinder are 50 m and 200 m, respectively. In all cases, down-drafts in front of the cylinder and updrafts, wakes, and vortex streets behind the cylinder appear. Low-speed flow below the cylinder height and high-speed flow above it are mixed behind the cylinder, resulting in strong negative vertical turbulent momentum flux at the boundary. Accordingly, the magnitude of the vertical turbulent momentum flux is the largest near the cylinder top. In the case of an inflow wind direction of $45^{\circ}$, the height of the boundary is lower than in other cases. As the inflow wind speed increases, the magnitude of the peak in the vertical profile of mean turbulent momentum flux increases due to the increase in speed difference between the low-speed and high-speed flows. As the cylinder width in the y-direction increases, the height of the boundary increases due to the enhanced updrafts near the top of the cylinder. In addition, the magnitude of the peak of the mean turbulent momentum flux increases because the low-speed flow region expands. Spectral analysis shows that the non-dimensional vortex generation frequency in the control case is 0.2 and that the cylinder width in the y-direction and the inflow wind direction affect the non-dimensional vortex generation frequency. The non-dimensional vortex generation frequency increases as the projected width of the cylinder normal to the inflow direction increases.