• The Korean Journal of Rheology
• /
• v.8 no.1
• /
• pp.1-10
• /
• 1996

본 논문에서는 협착이 발생된 원관내 뉴턴유체와 비뉴턴유체의 정상 및 맥동유동특 서을 유한요소법으로 해석하였다. 본연구는 맥동유동특성에서 협착관 형상의 변화, 협착이 주기적으로 발생된 협착관에서 협착부사이의 거리와 협착부의 수가유동특성에 미치는 영향 을 검토하였다. 레이놀즈수가 일정할 때 협착이 발생된 원형관내 뉴턴유체와 비뉴턴유체의 맥동유동특성은 정상유동의 경우와 크게 다르게 나타난다. 맥동유동에서는 정상유동보다 관 중앙부에서 속도분포가 훨씬 평탄하고 맥동유동의 속도분포는 감속시에 비하여 가속시에 관 중앙부의 속도분포가 더 평탄하게 나타난다. 정상유동과 맥동유동으 감속시에서는 협착부 하류의 벽면에서 재순환영역이 발생된다. 협착부의 수가 증가하면 각 협착부 주위의 속도장 은 유사하게 나타나지만 전체 압력손실은 크게 증가한다. 협착부사이의 거리가 변화될 경우 맥동유동속도의 국소최대치와 국소최소치의 차이가 가속시에는 거의 없지만 감속시에는 협 착부사이의 거리에 따라 다르게 나타난다.

• Proceedings of the KSR Conference
• /
• 2011.10a
• /
• pp.119-127
• /
• 2011

This paper shows development of 1-3 dimensional hybrid mesh method to analysis flow induced by ultra-high speed vehicle inside a long distance tunnel. For three-dimensional analysis of the tunnel system many meshes are required. However it is not efficient to calculate the whole tunnel system in three-dimension. Therefore in this paper, three-dimension meshes was used to describe stations, shafts and around vehicle, and one-dimension meshes was used to describe the tunnel except these three sections. And unsteady flow analysis of the ultra-high speed vehicle was performed with UDFs in commercial software, Ansys vr. 12.0.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.194-199
• /
• 2007

IB (immersed boundary) method is one of the prominent tool in computational fluid dynamics for the analysis of flows over complex geometries. The IB technique simplyfies the solution procedure by eliminating the requirement of complex body fitted grids and it is also superior in terms of memory requirement. In this study we have developed numerical code (FOTRAN) for the analysis of 2D flow over a cylinder using IB technique. The code is validated by comparing the wake lengths and separation angles given by Guo et. al. We employed fractional-step procedure for solving the Navier-Stokes equations governing the flow and discrete forcing IB technique for imposing boundary conditions. Also we have developed a 3D code for the backward-facing-step flow and flow over a sphere. The reattachment length in backward-facing-step flow was compared with the one given by Nie and Armaly, which has proven the validity of our code.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.41-44
• /
• 2007

A large eddy simulation method with unstructured mesh is presented. Two explicit filtering procedures are adopted for reducing the aliasing error of the nonlinear convective term and measuring the level of subgrid scale velocity fluctuation, respectively. The developed subgrid scale model is basically an eddy viscosity model which depends on both local velocity fluctuation level and local grid scale. As a validation problem, the flows around a sphere of several Reynolds numbers are simulated and some characteristic quantities are compared to experimental data and numerical results in the literature.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.5
• /
• pp.387-395
• /
• 2015

A technique for calculating the additive drag of the inlet in supersonic flow was studied using commercial CFD software, STAR-CCM+, which provides a efficient way of 3 dimensional flow analysis with polyhedron-shaped grid system. Three configurations were chosen and applied to the calculation with various flow conditions of two different free stream Mach No. and some mass flow ratios. Comparisons with results from wind tunnel test gave good agreements. Though computation were carried out with the inviscid and compressible flow around the supersonic inlet for the supercritical condition, ignoring the viscous effects is concluded to give little effects on the accuracy of the additive drag calculation and to make the calculation more efficient owing to less effort and time consumed for grid system build-up and for iteration because of less grid number and simpler boundary condition.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.401-402
• /
• 2006

The turbulent flow around a sphere was investigated using two experimental techniques: smoke-wire flow visualization in wind tunnel at Re=5300, 11000 and PIV measurements in a circulating water channel. The smoke-wire visualization shows flow separation points near an azimuthal angle of $90^{\circ}$, recirculating flow, transition from laminar to turbulent shear layer, evolving vortex roll-up and fully turbulent eddies in the sphere wake. The mean velocity field measured using a PIV technique in x-y center plane demonstrates the detailed near-wake structure such as nearly symmetric recirculation region, two toroidal vortices, laminar separation, transition and turbulent eddies. The PIV measurements of turbulent wake in y-z planes show that a recirculating vortex pair dominates the near-wake region.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.39-44
• /
• 2008

Numerical simulations of laminar flow over a sphere are conducted to investigate the effect of the Reynolds number on the characteristics of vortex shedding. The Reynolds numbers considered are between 300 and 475, covering unsteady planar-symmetric and asymmetric flows. Results show that the unsteady planar-symmetric flow can be categorized into two different regimes: single-frequency regime and multiple-frequency regime. The former has a single frequency component due to periodic shedding of the vortices with the same strength in every shedding cycle, while the latter has multiple frequency components due to cycle-to-cycle variation in the strength of shed vortices with the shedding angle fixed. The multiple-frequency planar-symmetric flow, which is newly found in the present study, occurs at Re=330${\sim}$360 between the single-frequency planar-symmetric flow and the asymmetric flow. On the other hand, the asymmetric flow occurs at Re${\geq}$365, where the vortices shed from the sphere show variation both in strength and shedding angle unlike the planar-symmetric flow. Also, it is shown that the breaking of planar symmetry is closely related to the imbalance of vortical strength between a pair of streamwise vortices.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.810-815
• /
• 2003

Flow-induced noise propagated from flow over a sphere is numerically investigated for laminar flow at Re = 300 and 425, and for turbulent flow at Re = 3700 and $10^4$, where the Reynolds number is based on the freestream velocity and the sphere diameter. The numerical method used for obtaining the flow over a sphere is based on an immersed boundary method in a cylindrical coordinate system. The Curle’s solutions of the Lighthill’s acoustic analogy with and without the far-field and compact-source approximation are used in order to investigate the noise field from flow over a sphere. Since the drag and lift forces change irregularly in time at Re = 425, 3700 and $10^{4}$, the noise propagates in a complicated manner. At Re = 300, 425 and $10^{4}$, the noise from dipole sources is much larger than that from quadrupole sources. On the other hand, at Re = 3700, the quadrupole source becomes dominant. The temporal variation of the flow-induced noise around a sphere is obtained at some observation points, which shows that the peak frequency corresponds to the Strouhal number associated with the wake instability.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1607-1612
• /
• 2004

Numerical simulations are performed to investigate the characteristics of flow past a sphere in uniform shear. The Reynolds numbers considered are Re=300, 425 and 480 based on the inlet center velocity and sphere diameter. The non-dimensional shear rate K of the inlet uniform shear is varied from 0 to 0.15. At Re=300, the head of the hairpin vortex loop always locates on the high-velocity side in uniform shear, and the flow maintains the planar symmetry. At Re=425 and 480, the irregularity in the location and strength of the hairpin vortex appearing in uniform inlet flow is much reduced in uniform shear, but the flows still keep the asymmetry for most inlet shear rates. However, in the cases of K=0.075 and 0.1 at Re=425, the flows become planar symmetric and their characteristics of the evolution of the hairpin vortex loops are different from those of asymmetric flows. A hysteresis phenomenon switching from the planar symmetry to the asymmetry (or vice versa) depending on the initial condition is also observed at Re=425.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1998.04a
• /
• pp.21-21
• /
• 1998

추력방향 제어시스템 설계에 있어서 가장 기본적으로 요구되는 Requirement는 Rocket Motor의 추력, 요구되는 최대 Side Force, Rocket Motor의 외경, System의 총 구동시간, 구동후의 분리여부 등이다. 이러한 Requirement를 만족하기 위해서는 Nozzle 출구의 분출가스 물성치로부터 초음속 유동해석을 통하여 Vane 주위의 속도, 온도, 압력 분포를 구하고, Vane의 받음각 변화에 대한 Aerodynamic Force와 Moment를 계산하고, Side Force를 만족하는 최대 받음각의 결정, Torque를 만족하는 감속기와 Motor의 선정 및 Housing 기본 형상을 설계하였다. 금번 개발에서는 지상 시험용으로서 안전 계수를 Flight Model보다 약간 높게 설계하였으며, 작동 완료 후 System이 Nozzle로부터 떨어져나가는 분리시스템은 포함하지 않았다.