• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 PARALLEL CFD 2006
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 대한조선학회논문집
• /
• 제48권3호
• /
• pp.222-232
• /
• 2011

This article examines the scale effect of the flow characteristics, resistance and propulsion performance on a 317k VLCC. The turbulent flows around a ship in both towing and self-propulsion conditions are analyzed by solving the Reynolds-averaged Navier-Stokes equation together with the application of Reynolds stress turbulence model. The computations are carried out in both model- and full-scale. A double-body model is applied for the treatment of free surface. An asymmetric body-force propeller is used. The speed performances including resistance and propulsion factors are obtained from two kinds of methods. One is to analyze the computational results in model scale through the revised ITTC' 78 method. The other is directly to analyze the computational results in full scale. Based on the computational predictions, scale effects of the resistance and the self-propulsion factors including form factor, thrust deduction fraction, effective wake fraction and various efficiencies are investigated. Scale effects of the streamline pattern, hull pressure and local flow characteristics including x-constant sections, propeller and center plane, and transom region are also investigated. This study presents a useful tool to hull-form and propeller designers, and towing-tank experimenters to take the scale effect into consideration.

• 대한조선학회논문집
• /
• 제47권2호
• /
• pp.163-177
• /
• 2010

This article examines hydrodynamic characteristics and speed performances of a ship attached with a full spade and a twisted rudder based on a computational method. For this study, a 13,100 TEU container carrier is selected. The turbulent flows around a ship are analyzed by solving the Reynolds-averaged Navier-Stokes equation together with the application of Reynolds stress turbulence model. The computations are carried out at the conditions of rudder, bare hull, hull-rudder and hull-propeller-rudder. An asymmetric body-force propeller is applied. The speed performance is predicted by the model-ship performance analysis method of the revised ITTC'78 method. The hydrodynamic forces are compared in both rudder-open-water and self-propulsion conditions. The flow characteristics, the speed performance including propulsion factors and the rudder-cavitation performance are also compared. The model tests are conducted at a deep-water towing tank to validate the computational predictions. The computational predictions show that the twisted rudder is superior to the full spade rudder in the respect of the speed and the cavitation performances.

• 대한조선학회논문집
• /
• 제39권3호
• /
• pp.28-40
• /
• 2002

지면효과를 받는 2차원 날개 주위의 난류유동을 비압축성 RANS(Reynolds Averaged Navier Stokes) 방정식과 유한차분법(Finite Difference Method)을 이용하여 해석하였다. 높은 레이놀즈수에 효과적인 Baldwin-Lomax 난류모델을 사용하였다. 본 연구의 목적은 지면효과를 받는 2차원 날개단면에서의 각기 다른 두 바닥 경계조건(이동지면, 고정지면)에 따른 유동의 특성을 파악하는 것이다. Clark-Y(t/C 11.7%)날개단면의 계산 결과와 발표된 계산결과 및 실험 값과의 비교를 통해 본 수치해석 프로그램의 정확성을 검증하였다. NACA4412 날개단면에 대해 지면과의 높이변화에 대해서 두 바닥 경계조건에 대해서 유동해석을 수행하였다 계산결과에 의하면 이동지면과 고정지면에 대해서 양력과 모멘트는 별 차이가 없으나 항력은 고정지면의 경우가 이동지면의 경우보다 다소 작았다. 따라서 풍동시험에서 고정지면의 결과는 이동지면에 비해 상대적으로 저항이 낮게 평가될 가능성이 있다고 본다.

• Journal of Mechanical Science and Technology
• /
• 제20권8호
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• 한국항공우주학회지
• /
• 제36권1호
• /
• pp.72-78
• /
• 2008

본 연구에서는 블레이드 구조 변형 효과를 고려하여 스테이터-로터 상호간섭 케스케이드 모델의 성능평가를 위한 유체-구조 연계해석 시스템을 개발하였다. 고정된 스테이터와 회전하는 로터는 상호간섭 영향이 유동해석에 고려되었으며, 레이놀즈-평균화 난류 방정식인 Spalart-Allmaras 모델과 k-ω SST 난류 모델이 압축성 유동박리 효과를 고려한 유동하중을 예측하기 위해 적용되었다. 정적인 유체-구조 연계해석과 수렴율 증진을 효과적으로 수행하기 위하여 큰 인공 감쇠를 가지는 연계 Newmark 시간적분 기법을 적용하였다. 수치실험을 통해 탄성축 위치에 따른 구조변형 효과가 케스케이드 성능에 미치는 영향을 파악하였다. 구조변형 효과가 고려된 경우 일반적인 강체 블레이드 모델에 대한 성능예측 결과와 다소 차이가 유발될 수 있음을 보였으며 공력탄성학적 영향을 고찰하였다.

• 대한기계학회논문집
• /
• 제19권7호
• /
• pp.1675-1683
• /
• 1995

A numerical study of the fluid flow and heat transfer characteristics of the two-dimensional impingement jet with a confinement plate has been carried out. The fluid flow was calculated by solving the full Navier-Stokes equation. In doing that, the well known SIMPLER algorithm was used and the trouble making convection term was discretized according to QUICKER scheme. The energy equation was simply solved by using the SOR method. For the Reynolds number of 10000, two channel heights, say 1.5 and 3.0 times the jet exit width, and two thermal boundary conditions constant wall temperature and constant wall heat flux were considered. Discrete heat sources were flush mounted along the impingement plate at a distance of 0, 2, 3, 4, 5, 6, 10, 12, times the jet exit width from the stagnation point. The length of each heat source is 4 times the jet exit width long. The Nusselt number averaged over each heat source was compared with experiment. Comparison shows that both calculations and experiment have the secondary peak of Nusselt number at downstream of stagnation point, even though there is a little quantitative difference in between. The difference is believed due to abscure thermal boundary condition in experiment and also accuracy of turbulence model used. The secondary peak is shown to be caused by rigorous turbulent flow motion generated as the wall jet flow is retarded and developes into the channel flow without flow reversal.

• 한국수자원학회논문집
• /
• 제38권3호
• /
• pp.189-198
• /
• 2005

본 연구에서는 조력발전소 배수갑문의 형상과 배치에 따른 방류능력을 해석하는데 3차원 수치모의가 효과적으로 이용될 수 있음을 보였다. 3차원 수치모형은 RANS를 지배방정식으로 하는 FLOW-3D 모형을 이용하였다. 본 연구결과 배수갑문의 방류능력은 물받이길이와 도류벽의 접근각도에 큰 영향을 받는 것으로 나타났다. 그리고 이의 개선 여부에 따라 $10\%$ 이상의 방류량 차이가 발생하였다. 또한 방류량은 배수문과 수차구조물을 연결하는 구조물의 형상과 물받이 끝 사면경사의 영향을 받는 것으로 나타났다. 본 연구에서는 배수갑문의 설계시 방류능력 개선을 위해서는 수리학적 검토가 필요하며, 수치모형실험이 수리모형실험과 더불어 유용한 해석도구로 이용될 수 있음을 보였다.

• 한국해안·해양공학회논문집
• /
• 제28권6호
• /
• pp.332-349
• /
• 2016

OLAFOAM은 파동역학의 시뮬레이션을 위하여 $OpenFOAM^{(R)}$을 확장한 강력한 CFD코드이며, $OpenFOAM^{(R)}$은 다양한 분야에서 각각 수치계산의 목적에 대응할 수 있도록 많은 Solver를 제공하고 있다. OLAFOAM의 기본방정식은 VARANS식에 기초하고, 수치기법으로는 유한체적법을 적용하며, 프로그램은 C++로 코딩되어 Linux운영체제에서 실행된다. 본 연구는 OLAFOAM을 이용하여 먼저 1) 단파와 규칙파하 투과성구조물에서 파의 변형, 2) 규칙파하 잠제에 의한 파의 변형 및 3) 흐름하 규칙파의 변형과 연직유속분포에 대해 기존의 각 실험결과와 비교 검토하여 OLAFOAM의 타당성을 검증하였다. 이로부터 지금까지 거의 검토되지 않은 규칙파와 흐름의 공존장에 설치된 투과성잠제에 대해 배후경사면을 불투과성 혹은 투과성으로 고려한 경우 흐름방향 등의 변화에 따른 잠제 주변에서 수위, 파고, 주파수스펙트럼, 쇄파, 평균유속 및 난류운동에너지 등의 변동특성을 면밀히 검토하였다. 결과로부터 흐름방향(순방향과 역방향)에 따른 파고변화는 난류운동에너지와 밀접한 관계를 가지는 것 등을 알 수 있었다.