• 대한조선학회논문집
• /
• 제30권4호
• /
• pp.61-73
• /
• 1993

본 연구에서는 탱크 내에 부분적으로 적재되어 있는 액체화물의 유동문제를 유한차분법을 이용해 해석하였다. 자유표면의 변화를 추적하고 유체장 내의 연속방정식과 Navier-Stokes 방정식을 만족시키기 위해 SOLA-SURF 기법이 적용되었다. 특히, 심한 유동으로 인한 탱크 상부에서의 충격압력을 예측하기 위해 충격완화영역의 개념을 도입하여 갑작스런 경계조건의 변화를 억제함으로써, 보다 현실적인 압력값을 예측할 수 있도록 하였다. 수치계산에서는 실험결과가 알려진 세가지 모델들에 대한 계산을 수행하였는데, 계산결과가 실험결과에 비교적 잘 일치하고 있었다. 그리고, 이 기법의 실선적용에 대한 가능성을 살펴보기 위해 30만톤급 초대형 유조선에 대한 실선계산도 수행하였다. 이러한 적용을 바탕으로 본 연구에서 적용된 수치기법이 sloshing 충격압력의 추정에 실용적으로 적용할 수 있음들 확인하였다.

• 한국해양공학회지
• /
• 제31권1호
• /
• pp.14-21
• /
• 2017

The erosion of solid particles in a pipe elbow was numerically investigated. A numerical procedure to estimate the sand erosion rate, as well as the particle motion, in the pipe elbow flow was introduced. This procedure was performed based on the combined empirical erosion model and computational fluid dynamics (CFD) analysis to consider the interaction between the particle motion and the eroded surface. The underlying turbulent flow on an Eulerian frame is described by the Reynolds averaged Navier-Stokes (RANS) equations with a $k-{\epsilon}$ turbulent model. The one-way coupled Eulerian-Lagrangian motion of the air flow and sand particles is employed to simulate the particle trajectories and particle-wall interactions on the pipe surfaces. The predicted CFD erosion magnitudes are compared with experimental data from pipe elbows. The erosion rate results do not reveal a good accordance between the simulation and experimental results. It seems that the CFD shows a slightly over-predicted erosion ratio.

• 한국전산유체공학회지
• /
• 제20권4호
• /
• pp.93-101
• /
• 2015

This paper provides numerical results of the simulation for the flow around the hull and the propeller of KCS model ship advancing in shallow water conditions. A finite volume method is used to solve the unsteady Reynolds averaged Navier-Stokes(RANS) equations, where the wave-making problem is solved by using a volume-of-fluid(VOF) method. The wave formed near the hull surface in shallow water conditions shows a deep trough dominant pattern that causes the loss of buoyancy followed by hull squat. The flow past the hull increases as the depth of water decreases. However, the axial flow velocity around the stern shows a reduction in magnitude by the effect of shallow water accompanied by the hull-propeller interaction. As a results, the thrust and torque coefficient increase about 8.3% and 6.2%, respectively for a depth of h/T=3.0 corresponding to a depth Froude number of $F_h=0.693$. The resistance coefficient increases about 11.6% at this Froude number condition.

• 한국전산유체공학회지
• /
• 제17권3호
• /
• pp.11-17
• /
• 2012

In this study, aerodynamic analyses have been conducted for 4-Bladed Vertical-Axis Wind Turbine (VAWT) configuration and the results are compared with experimental data. Reynolds-averaged Navier-Stokes equation with LES turbulence model is solved for unsteady flow problems. In addition, the computation results by standard k-${\omega}$ and SST k-${\omega}$ turbulence models are also presented and compared. An experiment model of 4-Bladed VAWT model has been designed and constructed herein. Experimental tests for aerodynamic performance of the present VAWT model are practically conducted using the vehicle mounted testing system. Comparison results between the experiment and the computational fluid dynamics (CFD) analyses are presented in order to show the accuracy of CFD analyses using the different turbulent models.

• 한국전산유체공학회지
• /
• 제14권3호
• /
• pp.105-114
• /
• 2009

A numerical study of the flow characteristics inside a U-type circular tube is carried out in this paper. The numerical simulations carried out by using a Navier-Stokes code which is commercially available. Before detailed numerical simulations, validation of present numerical approach is made by comparing numerical solutions with experimental data. Numerical simulations are performed to study the effect of curvature on the flow characteristics inside a U-type tube. Numerical solutions show that a significant effect on the secondary flow structure in the cross section of the tube, especially in the curved section is shown when the curvature ratio, ratio of curvature to tube diameter, is smaller than about 3.5. As the curvature ratio decreases below 3.5, a counter rotating vortex is found below the primary vortex in the cross section of the tube. Another dramatic change of the flow structure is the formation of streamwise separation zone when the curvature ratio is decreased below 1.25.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 춘계학술대회논문집
• /
• pp.113-116
• /
• 2008

This work presents numerical optimization for design of staggered arrays of dimples printed on opposite surfaces of a cooling channel with a fast and elitist Non-Dominated Sorting of Genetic Algorithm (NSGA-II) of multi-objective optimization. As Pareto optimal front produces a set of optimal solutions, the trends of objective functions with design variables are predicted by hybrid multi-objective evolutionary algorithm. The problem is defined by three non-dimensional geometric design variables composed of dimpled channel height, dimple print diameter, dimple spacing and dimple depth to maximize heat transfer rate compromising with pressure drop. Twenty designs generated by Latin hypercube sampling were evaluated by Reynolds-averaged Navier-Stokes solver and the evaluated objectives were used to construct Pareto optimal front through hybrid multi-objective evolutionary algorithm. The optimum designs were grouped by k-mean clustering technique and some of the clustered points were evaluated by flow analysis. With increase in dimple depth, heat transfer rate increases and at the same time pressure drop also increases, while opposite behavior is obtained for the dimple spacing. The heat transfer performance is related to the vertical motion of the flow and the reattachment length in the dimple.

• Journal of Mechanical Science and Technology
• /
• 제15권3호
• /
• pp.366-373
• /
• 2001

A numerical analysis of shock wave/boundary layer interaction in transonic/supersonic axial flow compressor cascade has been performed by using a characteristics upwind Navier-Stokes method with various turbulence models. Two equation turbulence models were applied to transonic/supersonic flows over a NACA 0012 airfoil. The results are superion to those from an algebraic turbulence model. High order TVD schemes predicted shock wave/boundary layer interactions reasonably well. However, the prediction of SWBLI depends more on turbulence models than high order schemes. In a supersonic axial flow cascade at M=1.59 and exit/inlet static pressure ratio of 2.21, k-$\omega$ and Shear Stress Transport (SST) models were numerically stables. However, the k-$\omega$ model predicted thicker shock waves in the flow passage. Losses due to shock/shock and shock/boundary layer interactions in transonic/supersonic compressor flowfields can be higher losses than viscous losses due to flow separation and viscous dissipation.

• 대한기계학회논문집B
• /
• 제21권4호
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• 대한기계학회논문집B
• /
• 제24권6호
• /
• pp.860-869
• /
• 2000

Numerical and experimental investigations are performed for the molecular transition and slip flows in pumping channels of a disk-type drag pump. The flow occurring in the pumping channel develops from the molecular transition to the slip flow traveling downstream. Two different numerical methods are used in this analysis: the first one is a continuum approach in solving the Navier-Stokes equations with slip boundary conditions, and the second one is a stochastic approach through the use of the direct simulation Monte Carlo method. In the experimental study, the inlet pressures are measured for various outlet pressures in the range of 0.1{\sim}4Torr. From the present study, the numerical results of predicting the performance, obtained by both methods, agree well with the experimental data for the range of Knudsen number $Kn{\leq}0.1$ (i.e., the slip flow regime). But the results from the second method only agree with the experimental data for Kn>0.1(i.e., the molecular transition regime)

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2007년도 춘계 학술대회논문집
• /
• pp.64-67
• /
• 2007

In the present study, the wind pressure transients on platform screen door in side platforms caused by passing trains have been investigated numerically. The transient compressible 3-D full Navier-Stokes solution is obtained with actual operational condition of subway train and the moving mesh technique is adopted considering the train movement. To achieve more accurate analysis, the entrance and exit tunnel of platform are included in a computational domain and detailed shape of train is also modeled Numerical analyses were conducted on five operational condition which are different velocity variation of subway train, existence of stationary train and passing each other trains. The results show that pressure load on platform screen door is maximized when the two trains are passing each other. It is also seen from the computational results that the maximum pressure variation was found to be satisfactory to various foreign standards.