• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2002년도 추계 학술대회논문집
• /
• pp.17-22
• /
• 2002

A numerical solution procedure has been developed to analyze the flow field resulted from the injection of a transverse jet through the divergent flap of a 2DCD nozzle for thrust vector control. The formulation employs the compressible Navier-Stokes equations in conservation law form and a two equation $\kappa-\omega$ turbulence model. Detailed numerical results are presented in this paper for the 2DCD nozzle over a range of secondary to primary injection mass flow ratios and nozzle pressure ratios.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.13-18
• /
• 2009

Flight vehicles such as wheel wells and bomb bays have many cavities. The flow around a cavity is characterized as an unsteady flow because of the formation and dissipation of vortices brought about by the interaction between the free stream shear layer and the internal flow of the cavity. The resonance phenomena can damage the structures around the cavity and negatively affect the aerodynamic performance and stability of the vehicle. In this study, a numerical analysis was performed for the cavity flows using the unsteady compressible three-dimensional Reynolds-Averaged Navier-Stokes (RANS) equation with Wilcox's turbulence model. The Message Passing Interface (MPI) parallelized code was used for the calculations by PC-cluster. The cavity has aspect ratios (L/D) of 2.5 ~ 7.5 with width ratios (W/D) of 2 ~ 4. The Mach and Reynolds numbers are 0.4 ~ 0.6 and $1.6{\times}106$, respectively. The occurrence of oscillation is observed in the "shear layer and transient mode" with a feedback mechanism. Based on the Sound Pressure Level (SPL) analysis of the pressure variation at the cavity trailing edge, the dominant frequencies are analyzed and compared with the results of Rossiter's formula. The dominant frequencies are very similar to the result of Rossiter's formula and other experimental data in the low aspect ratio cavity (L/D = ~ 4.5). In the large aspect ratio cavity, however, there are other low dominant frequencies due to the leading edge shear layer with the dominant frequencies of the feedback mechanism. The characteristics of the acoustic wave propagation are analyzed using the Correlation of Pressure Distribution (CPD).

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

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

The heat transfer characteristics of a turbulent flow in a ribbed two-dimensional channel have been investigated numerically. The fully elliptic governing equations, coupled with a four-equation turbulence model, $\kappa-\omega-\bar{t^2}-\epsilon_t$, are solved by a finite volume method of SIMPLE type. Calculations have been carried out for three rib cross-sections : square, triangular, and semicircular, with various rib pitches and Reynolds numbers. The procedure appears to be satisfactory as the results for the square rib compare favorably with available experimental data and earlier calculation. The optimal rib pitch that yields the maximum heat transfer has been identified. It is also found that the square rib is most effective in enhancing the heat transfer. The semicircular rib, on the other hand, incurs the least amount of pressure drop but the improvement in heat transfer is substantially lower.

• 한국항공우주학회지
• /
• 제39권4호
• /
• pp.289-296
• /
• 2011

본 논문에서는 항공기의 항력 감소를 위해 페어링 형상을 설계하고 설계 형상에 대한 공력해석을 수행하였다. 페어링은 날개 끝과 날개-동체 접합부에 대해 적용되었다. 날개끝 페어링의 경우 날개 전체에 대한 유동 해석을 수행하였고 날개-동체 페어링에 대해서는 전기체에 대한 해석을 수행하였다. Menter의 ${\kappa}-{\omega}$ SST 난류 모델이 적용된 3차원 RANS 코드를 이용하였으며 유동 해석을 실시하였다. 페어링의 형상과 유무에 따른 해석 결과를 비교하여 각각의 페어링 날개와 항공기의 항력에 미치는 영향을 확인하였다.

• 대한조선학회논문집
• /
• 제51권1호
• /
• pp.67-77
• /
• 2014

In this study, Explicit Algebraic Reynolds Stress Model (EARSM) which is based on the existing ${\kappa}-{\omega}$ model has been applied to the flow field analysis around ship hulls. Existing transport equations for the turbulent kinetic energy and the dissipation rate are used in almost the same form and anisotropy terms of Reynolds stresses are newly considered. The well-known KVLCC2 and KCS hull forms are selected as validation cases, which were also used in 2010 Workshop on CFD in Ship Hydrodynamics. In case of KVLCC2 double model, comparison of mean velocity distribution, turbulent kinetic energy, and Reynolds stresses near the propeller plane has been carried out and wave elevation and wave profiles have been additionally studied for KCS and KVLCC2 with free surface models. Some improved results for mean velocity distribution at the propeller plane have been obtained while there is little change in free surface wave profiles.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.592-595
• /
• 2011

A numerical investigation on a suction vortices, free vortices and subsurface vortices behavior in the model sump system with multi-intakes is performed A test model sump and piping system were designed based on Froude similitude for the prototype of the recommended structure layout by HI-9.B Standard for Pump Intake Design of the Hydraulic Institute. A numerical analysis of three dimensional multiphase flows through the model sump is performed by using the finite volume method of the CFX code with multi-block structured grid systems. A ${\kappa}-{\omega]$ ShearStressTransportturbulencemodelandthe Rayleigh-Plesset cavitation model are used for solving turbulence cavitating flow. From the numerical analysis, several types of vortices are reproduced and their formation, growing shedding and detailed vortex structures are investigated. To reduce abnormal vortices, an anti-vortex device is considered and its effect is investigated and discussed.

• 대한기계학회논문집B
• /
• 제30권10호
• /
• pp.966-972
• /
• 2006

A numerical simulation is performed fur developing turbulent flow in a strongly curved 180 deg pipe and its downstream tangent by a new solution code(PowerCFD) which adopts an unstructured cell-centered method. The governing equations are discretized as the full elliptic from of the equations of motion. Three typical two-equation turbulence models of low-Reynolds-number form are used to approximate the turbulent stress field. Solutions fur both streamwise and circumferential velocity components are compared with the experimental data by Azzola et at.(1986). The ${\kappa}-{\omega}$ model by Wilcox(1988) is found to give better prediction performance than the other two. Predicted secondary velocities and streamwise velocity component contours at sequential longitudinal stations are also presented in order to enable a detailed description of the complete flow. It is also found that, in the bend both mean streamwise and secondary velocities never achieve a fully-developed state and the code is capable of producing very well the complex nature of steady flow in a strongly curved pipe.

• 한국전산유체공학회지
• /
• 제21권2호
• /
• pp.70-80
• /
• 2016

The comparison of two commercial codes(FLUENT and STAR-CCM+) and an open-source code(OpenFOAM) are carried out for the aerodynamic analysis of flight vehicles at low speeds. Tailless blended-wing-body UCAV, main wing and propeller of HALE UAV(EAV-3) are chosen as geometries for the investigation. Using the same mesh, incompressible flow simulations are carried out and the results from three different codes are compared. In the linear region, the maximum difference of lift and drag coefficients of UCAV are found to be less than 2% and 5 counts, respectively and shows good agreement with wind tunnel test data. In a stall region, however, the reliability of RANS simulation is found to become poor and the uncertainty according to code also increases. The effect of turbulence models and meshes generated from different tools are also examined. The transition model yields better results in terms of drag which are much closer to the test data. The pitching moment is confirmed to be sensitive to the existence and the location of transition. For the case of EAV-3 wing, the difference of results with ${\kappa}-{\omega}$ SST model is increased when Reynolds number becomes low. The results for the propeller show good agreement within 1% difference of thrust. The reliability and uncertainty of three codes is found to be reasonable for the purpose of engineering use. However, the physical validity and reliability of results seem to be carefully examined when ${\kappa}-{\omega}$ SST model is used for aerodynamic simulation at low speeds or low Reynolds number conditions.

• 한국항공운항학회:학술대회논문집
• /
• 한국항공운항학회 2004년도 추계학술발표대회 논문집
• /
• pp.240-246
• /
• 2004

For the prediction on the onset of oblique shock wave-induced vortex breakdown, computational studies on the Oblique Shock wave/Vortex Interaction (OSVI) are conducted and compared with both experimental results and analytic model. A Shock-stable numerical scheme, the Roe scheme with Mach number-based function (RoeM), and a two-equation eddy viscosity-transport approach are used for three-dimensional turbulent flow computations. The computational configuration is identical to available experiment, and we attempt to ascertain the effect of parameters such as a vertex strength, streamwise velocity deficit, and shock strength at a freestream Mach number of 2.49. Numerical simulations using the ${\kappa}-{\omega}SST$ turbulence model and suitably modeled vortex profiles are able to accurately reproduce many fine features through a direct comparison with experimental observations. The present computational approach to determine the criterion on the onset of oblique shock wave-induced vortex breakdown is found to be in good agreement with both the experimental result and the analytic prediction.