• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.172-177
• /
• 2004

In the present work, we have interests on the modification of parallel implemented with MPI(Message Passing Interface) programming method, 3-Dimensional, unsteady, incompressible Navier-Stokes equation solver to analyze the low-Reynolds number flow In order to accurate calculation aerodynamic coefficients in low-Reynolds number flow field, we modified the two-equation turbulence model. This paper describes the development and validation of a new two-equation model for the prediction of flow transition. It is based on Mentor's low Reynolds $\kappa-\omega$ model with modifications to include Total Stresses Limitation (TSL) and Separation Transition Trigger (STT)

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.15-25
• /
• 2001

Numerical prediction of the diffusion controlled transition in a turbine gas pass is important because it can change the local heat transfer rate over a turbine blade as much as three times. In this study, the gas flow over turbine blade is simplified to the flat plate boundary layer, and an adaptive grid scheme redistributing grid points within the computation domain is proposed with a great emphasis on the construction of the grid control function. The function is sensitized to the second invariant of the mean strain tensor, its spatial gradient, and the interaction of pressure gradient and flow deformation. The transition process is assumed to be described with a κ-ε turbulence model. An elliptic solver is employed to integrate governing equations. Numerical results show that the proposed adaptive grid scheme is very effective in obtaining grid independent numerical solution with a very low grid number. It is expected that present scheme is helpful in predicting actual flow within a turbine to improve computation efficiency.

• Journal of KSNVE
• /
• v.10 no.4
• /
• pp.602-609
• /
• 2000

When an axisymmetric body moves through air the boundary layer near the stagnation region remains laminar and subsequently it goes through transition to turbulent. The experimental investigation described in this paper concerns the characteristics of wall pressure fluctuations at the initial stage of boundary layer flow including transition. Flush-mounted microphones are used to measure the wall pressure fluctuations at the transition and turbulent boundary layer region of a blunt axisymmetric body in the low noise wind tunnel. It if found from this study that the wall pressure fluctuations in the transition region is higher than that in the turbulent region.

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

In the present paper, some difficulties encountered in predicting airfoil characteristics are described and solutions for those problems are discussed Since drag is determined by the amounts of pressure and, especially, shear stress, accurate estimation of shear stress is very crucial. However shear stress computation is dependent on the grid density and turbulence model, it should be consistent in preparing grid and turbulence model. When the transition from laminar to turbulent happen at the middle of airfoil, CFD solver should divide the region into laminar and turbulent region based on the transition location.

• Proceeding of EDISON Challenge
• /
• 2014.03a
• /
• pp.521-526
• /
• 2014

현대에 이르러 초경량 무인 비행기에 대한 많은 연구가 진행되고 있다. 이러한 비행체는 저레이놀즈수 영역에서 사용되는 특성으로 인해, 경계층 내에서 박리현상과 난류영역으로의 천이 등과 같은 여러 복합적인 현상을 발생시킴으로써 비행체의 공력특성에 큰 영향을 미친다. Bumpy Airfoil은 저레이놀즈수 유동에서의 이와 같은 문제를 해결하기 위해 제안된 익형이다. 따라서 본 논문은 전산열유체해석 프로그램인 EDISON_전산열유체를 이용하여 Bumpy Airfoil 형상에 대한 공력특성을 연구하였고, 발생하는 양항비를 원 익형과 비교하였다. 비압축성 조건 내에서, 공력 성능 향상을 위한 Bumpy Airfoil의 형상 변수로 Bump 개수와 높이를 선정하여 받음각에 따른 유동장을 분석하고 양항비를 수치해석 및 고찰하였다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.19-23
• /
• 2010

Injection characteristics of a kerosene swirl injector of liquid rocket engine operating at supercritical environment have been investigated. Kerosene surrogate models are proposed to model the kerosene properties. Turbulent numerical model is based on large eddy simulation and contains Soave modification of Redlich-Kwong equation of state and Chung's model. Numerical analysis results at supercritical environment are compared with the one at transcritical condition. Differences of density and viscosity are analyzed at both liquid film and core gas in the swirl injector.

• Tribology and Lubricants
• /
• v.9 no.1
• /
• pp.22-31
• /
• 1993

The thermohydrodynamic performance of large tilting pad journal bearings is analyzed, taking into account the three dimensional variation of lubricant viscosity. The eddy viscosity model based on wall formula is applied. The effects of temperature rise and turbulence on the bearing performance are studied in comparision with the isothermal or the laminar analysis. It is shown that these effects have significant influence on temperature distribution, load capacity and power loss of the bearing.

• Journal of Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1179-1183
• /
• 2007

관망 내에서 흐름의 연속 방정식과 운동량 방정식을 상 미분으로 전개하여 해석한 특성선 방법은 주로 가압 관망체계(Pressurized Pipeline System)에서의 부정류 해석(Unsteady Analysis)에 사용 된다. 그러나 이특성선 방법은 천이류 해석을 위한 관망 재구성 과정에서 Courant수 조건의 만족을 위한 관의 재배열에 천문학적인 계산용량과 시간이 필요하다는 단점이 있다. 이는 현장 적용 시 압력파 전파속도의 불확실성과 연계되어 상당한 장해요소가 되고 있다. 이에 대안적인 방법으로서 임펄스응답법이 개발되었다. 이는 경계지점에서 복소수 유량에 대한 복소수 수두의 비율로써 정의된 관망에서의 수리임피던스를 역퓨리에 변환에 적용하여, 주파수 영역의 수치를 시간 영역으로 변환하여 응답함수를 산출한 후, 산출된 응답함수와 구해진 경계지점에서의 유량과의 적분을 통하여 임의의 지점에서의 수두 및 유량을 계산하는 방법이다. 임펄스 응답법은 관 부속물관의 특성을 기술하는 수학적 표현의 난해함으로 인해 지금까지는 단일관에 대한 연구에만 국한되어 왔다. 본 연구에서는 임펄스응답법을 수리구조물이 부착된 관망에 적용하여 다양한 조건에서 천이류 분석을 시행하였다. 즉, 에어챔버 및 서지탱크와 같은 수리구조물을 각각에 대한 수리임피던스를 구하고, 가지관 및 통합 관성항으로 취급하여 수리구조물을 처리하였다. 그리고 이러한 결과를 특성선방법과 비교하여 그 적절성을 검증하였는데, 특성선 방법에 의한 모의 결과와 비교하였을 때, 일치하는 결과를 나타내었다. 임펄스응답법에 의한 모의 결과에서 감쇄효과를 과대평가하는 경향이 관찰되었다. 이는 임펄스 응답법의 가정에 기인한 것으로써 난류 상태의 흐름에서 상당한 불일치를 가져올 수 있으나, 수리 구조물에 의한 수격압이 감쇄되는 과정에서 대부분 흐름이 층류 상태로 전환된다고 가정 할 때는 상당한 적용성이 있다. 본 연구는 수리구조물이 부착된 관망의 해석함에 있어서 임펄스응답법의 적용이 가능함을 보였고, 이는 보다 복잡한 관망에서의 천이류 해석이 가능함을 시사한다.