• Wind and Structures
• /
• v.18 no.2
• /
• pp.155-180
• /
• 2014

The safety of road vehicles on the ground in crosswind has been investigated for many years. One of the most important fundamentals in the safety analysis is aerodynamic characteristics of a vehicle in crosswind. The most common way to study the aerodynamic characteristics of a vehicle in crosswind is wind tunnel tests to measure the aerodynamic coefficients and/or pressure coefficients of the vehicle. Due to the complexity of wind tunnel test equipment and procedure, the features of flow field around the vehicle are seldom explored in a wind tunnel, particularly for the vehicle moving on the ground. As a complementary to wind tunnel tests, the numerical method using computational fluid dynamics (CFD) can be employed as an effective tool to explore the aerodynamic characteristics of as well as flow features around the vehicle. This study explores crosswind effects on a high-sided lorry on the ground with and without movement through CFD simulations together with wind tunnel tests. Firstly, the aerodynamic forces on a stationary lorry model are measured in a wind tunnel, and the results are compared with the previous measurement results. The CFD with unsteady RANS method is then employed to simulate wind flow around and wind pressures on the stationary lorry. The numerical aerodynamic forces are compared with the wind tunnel test results. Furthermore, the same CFD method is extended to investigate the moving vehicle on the ground in crosswind. The results show that the CFD results match with wind tunnel test results and the current way using aerodynamic coefficients from a stationary vehicle in crosswind is acceptable. The CFD simulation can provide more insights on flow field and pressure distribution which are difficult to be obtained by wind tunnel tests.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.10
• /
• pp.939-946
• /
• 2008

Euler codes or Navier-Stokes codes for compressible flows suffer severe degradation in convergence as Mach number approaches zero. The convergence problem arose from the wide disparity in characteristic speeds can be solved using preconditioning methods without large modifications. In this paper, a preconditioned RANS(Reynolds Averaged Navier-Stokes) solver is developed for analysis of low Mach number flows. In order to validate the method, computational examples are chosen and the results are compared with the experimental data and the existing computed results showing a good accuracy and convergence characteristics for steady inviscid, laminar and turbulent flows at low Mach number.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.54-54
• /
• 2011

기존의 착저식 방파제를 보완하기 위하여 부유식 방파제가 개발되었고, 많은 연구가 선행되어 왔다. 부유식 방파제의 최대 장점은 경제성과 친환경성이다. 그러나 부유식 방파제는 소파성능이 떨어진다는 단점이 있으며, 이를 개선하기 위해 잠재와 혼용, 배열형에 관한 연구등이 선행되어왔다. 그러나 이것은 경제성이라는 강점을 고려하지 못하였다. 그래서 본 연구에서는 부유식 방파제의 중요한 장점중 하나인 경제성을 고려해, 단면현상 변화만을 이용하여 부유식 방파제의 소파성능 개선하고자 하였다. RANS(Reynolds averaged Navier-Stokes) 방정식에 기초하여 VOF법과 $k-{\varepsilon}$ 난류모델을 결합한 수치모델인 CADMAS-SURF를 이용하였으며, 구조물 단면형상 변화를 이용해 와의 상호 간섭을 유도하였고, 이에 따른 투과율 변화를 관찰 하였다. 결과를 살펴보면 요철1 단면에서는 구조물 전면 하단부분과 구조물 후면 하단부분에서 와의 간섭이 일어났으며 가장 아래 요철 부분에서 유속의 전달현상이 보인다. 투과계수는 일반적인 부유식 방파제와 마찬가지로 L/B가 1~4사이 값인, 비교적 단주기에서는 0.3~0.4의 투과율을 보였으나 L/B가 5를 넘어가면서 0.45~0.55의 투과율을 보였고, 요철2 단면에서는 전면과 후면에서 발달한 와가 전, 후면 돌출부에 의해 바닥까지 전파되지 못하는 양상을 보였으며, 돌출부 사이 중앙부분에서 가장 활발한 와의 간섭을 관찰 할 수 있었다. 돌출부 아래에서 역시 강력한 와의 간섭을 보이고 있다. 투과율 역시 가장 낮은 값을 보였으며 비교적 단주기 구간인 B/L 1~4 에서는 0.2~0.35 사이의 값을 가졌으며 5~10사이구간에서는 0.35~0.34의 값을 보이고 있다. 이 같은 결과는 와의 간섭이 가장 활발하게 나타난 결과로 보인다. 그리고 요철 3단면에서는 전면 돌출부 끝단에서의 활발한 와의 간섭을 관찰 할 수 있었다. 투과율은 세 단면 중 가장 높은 값의 투과율을 보이지만 B/L 3~4 구간에서 요철1 경우보다 낮은 값의 투과율을 보이고 있다. 결과에서 보듯이 도출부의 적절한 조합과 배치를 통해 언급한 연구목표(와의 생성과 간섭, 방파효율 개선)를 달성하였고 추후에 돌출부의 크기와 배치, 흘수의 영향, 수심의 영향 등을 고려한 연구가 진행된다면 더욱 우수한 단면형상을 개발 할 것이라 예상된다.

• Wind and Structures
• /
• v.23 no.6
• /
• pp.505-521
• /
• 2016

The objective of this study is to investigate numerically the effect of building roof shaps on wind flow and pollutant dispersion in a street canyon with one row of trees of pore volume, $P_{vol}=96%$. A three-dimensional computational fluid dynamics (CFD) model is used to evaluate air flow and pollutant dispersion within an urban street canyon using Reynolds-averaged Navier-Stokes (RANS) equations and the Explicit Algebraic Reynolds Stress Models (EARSM) based on k-${\varepsilon}$ turbulence model to close the equation system. The numerical model is performed with ANSYS-CFX code. Vehicle emissions were simulated as double line sources along the street. The numerical model was validated by the wind tunnel experiment results. Having established this, the wind flow and pollutant dispersion in urban street canyons (with six roof shapes buildings) are simulated. The numerical simulation results agree reasonably with the wind tunnel data. The results obtained in this work, indicate that the flow in 3D domain is more complicated; this complexity is increased with the presence of trees and variability of the roof shapes. The results also indicated that the largest pollutant concentration level for two walls (leeward and windward wall) is observed with the upwind wedge-shaped roof. But the smallest pollutant concentration level is observed with the dome roof-shaped.

• The KSFM Journal of Fluid Machinery
• /
• v.17 no.6
• /
• pp.52-58
• /
• 2014

In the present work, design optimization of film-cooling hole array on the pressure side of high pressure turbine nozzle was conducted. There are four rows of fan-shaped film cooling holes on the nozzle pressure side surface and each row has a straight array of holes in the spanwise direction for baseline model. For design optimization, hole distributions in streamwise and spanwise directions for three rows of holes except first row are parameterized as a 2nd-order shape function. Three-dimensional compressible RANS equations are used for flow and thermal analysis around the nozzle surface and optimization technique using Design of Experiment, Kriging surrogate model and Genetic Algorithm is used. The results shows that averaged adiabatic wall temperature at the whole nozzle surface decreases about 2.7% and averaged film cooling effectiveness at the pressure side of nozzle increased about 8.2%.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.2
• /
• pp.75-83
• /
• 2016

Since a turbomachine has complex flow characteristics, which are caused by adverse pressure gradient and high speed motion, an elaborate turbulence model is needed to accurately predict the flow. Some turbulence models such as an algebraic or a two-equation eddy viscosity model have been used for in-house RANS-code, but it is difficult to obtain good result for several complex flows. In this study, Durbin's V2-F turbulence model, which has been known for better prediction for severe flow separation, is applied to T-Flow. It was validated for simple cases such as channel and compressor cascade, and its applicability to turbomachinery was shown by analyzing internal flow of a single rotor. As a result, the V2-F turbulence model shows better blade surface pressure distribution than the one-and-two equation turbulence model.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• 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.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.20 no.6
• /
• pp.46-53
• /
• 2016

It is required to predict a limit cycle amplitude controlled by system's nonlinear behavior as well as an eigen-frequency and initial growth rate of instabilities under the linear motions, in order to fully understand combustion instabilities in a lean premixed gas turbine combustor. Special focus of the current work is placed on the limit cycle amplitude prediction using flame describing function(FDF) where the ratio of a heat release fluctuation to a given flow perturbation is expressed as a function of frequency and amplitude. In this study, the CFD modeling work based on RANS is carried out to obtain FDF, which makes that the nonlinear thermo-acoustic model is successfully developed for predicting the limit cycle amplitude of the combustion instability.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.5
• /
• pp.355-363
• /
• 2013

Panel methods are essential tools for analyzing a fluid-flow problem around complex three dimensional bodies, but they lack ability to solve viscous effects. On the other hand, CFD methods are considered as powerful tools for analyzing fluid-flow characteristics including viscosity. However, they also have short falls, requiring more computing time and showing different results depending on the selection of turbulence models and grid systems. In this paper a hybrid scheme combining a panel method and a CFD method is suggested. The scheme adopts a panel method for far-field solution where viscous effects are negligible and a CFD method for the solution of RANS equations in near-field where viscous effects are relatively strong. The intermediate region between the far-field and near-field is introduced where the calculated field point velocities by the panel method are given as boundary velocities for the CFD method. To verify the scheme, calculated results, by a panel method, a CFD method and the hybrid scheme, for a two dimensional foil section are compared. The suggested hybrid scheme gives reasonable results, while computation time and memory can be dramatically reduced. By using the hybrid scheme efforts can be concentrated for the local flow near the leading and trailing edges, by providing more dense grid system, where detailed flow characteristics are required.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.2
• /
• pp.147-157
• /
• 2011

A RANS(Reynolds averaged Navier-Stokes) based numerical method is developed for the evaluation of ship resistance and self-propulsion performances. In the usability aspect of CFD for the hull form design, the field grid around practical hull forms is generated by solving a grid Poisson equation based on the hull surface grid generated from station offsets and centerline profile. A body force technique is introduced to model the effects of the propeller in which the propeller loads are obtained from potential flow analysis using an unsteady lifting surface method. The free surface is captured by using a two-phase level-set method and the realizable $k-{\varepsilon}$ model is used for turbulence closure. The hull attitude in vertical plane, i.e., trim and sinkage, is calculated by using a quasi-steady method and then considered in the computation by translating and rotating the grid system according to the values. For the validation of the proposed method, the numerical results of resistance tests for KCS, KLNG, and KVLCC1 and of self-propulsion test for KCS are compared with experimental data.