• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.335-340
• /
• 2010

In this paper development of an automatic grid generation program for flow field calculation around 3D wing is described and its application is also introduced. The program is developed by using JAVA programming language and a graphic library, JOGL, and it can be usee either as an application program on a local computer or as a applet in the network environment. Currently, The program provides NACA series 4-digit airfoils as the wing cross-section shape and it offers a non-complicated GUI program which can easily generate structured grids for wings based on user's parameter input. Grid generated by the program can be selected as one of two types; O-type and C-type. In this research advancing layer method(ALM) augmented by elliptic smoothing method is used for the FLUENT. It is shown that by using current program high-quality structured grids around 3D wings can be easily generated, and typical grid generation results and flow solutions are demonstrated. Study on effects of geometric parameters on flow field is also tried by changing major wing parameters such as incidence angle type of wing-tip and sweepback angle.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.166-170
• /
• 2008

The changing process from hovering mode to transition one is of importance to determine a stability of tilt-rotor aircraft, which is utilized in UAV(Unmaned Aerial Vehicle). The analysis on fluid flows and aerodynamic characteristics according to variation of tilting angle of rotor is essential step in development of tilt-rotor. In the present study, the computation domain is divided into the rotating and stationary regions in order to consider the rotating blades. For the convenient realization of various tilting angle as well as application of boundary condition, the whole computation region is constructed into sphere domain. The near farfield boundary condition is adopted. The airfoil used in computation is NACA 0012. The computation results for the hovering mode are validated by comparing with previously conducted experimental results. From the results, the flow fields around rotor blade and the aerodynamic characteristics in transition mode are observed. The computational result will provide the basis for development and performance evaluation of tilt-type aircraft.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.11 s.254
• /
• pp.1084-1092
• /
• 2006

Shape optimization of airfoil in WIG craft has been performed by considering the ground effect. The WIG craft should satisfy various aerodynamic characteristics such as lift, lift to drag ratio, and static height stability. However, they show a strong trade-off phenomenon so that it is difficult to satisfy aerodynamic properties simultaneously. Optimization is carried out through the multi-objective genetic algorithm. A multi-objective optimization means that each objective is considered separately instead of weighting. Due to the trade-off, pareto sets and non-dominated solutions can be obtained instead of the unique solution. NACA0015 airfoil is considered as a baseline model, shapes of airfoil are parameterized and rebuilt with four-Bezier curves. There are eighteen design variables and three objective functions. The range of design variables and their resolutions are two primary keys for the successful optimization. By two preliminary optimizations, the variation can be reduced effectively. After thirty evolutions, the non-dominated pareto individuals of twenty seven are obtained. Pareto sets are all the set of possible and excellent solution across the design space. At any selections of the pareto set, these are no better solutions in all design space.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.115-119
• /
• 2004

The effect of a continuous blowing or suction on an oscillating 2-D NACA0012 airfoil was investigated numerically for the dynamic stall control. The influence of control parameter variation was also studied in the view point of aerodynamic characteristics. The result showed that the blowing control kept a higher lift drag ratio before stall angle but the dynamic stall angle was not exceed to without control result. As the slot position was closer to leading edge, the positive control effect becomes greater. The stronger jet and the smaller jet angel made more favorable roles on the control performance. In the cases of the suction, the overall control features were similar to those of the blowing, but dynamic stall angle was increased, i.e. suction was more effective to control dynamic stall. It was also founded that the suction control was showed better control effect as the slot position moves to trail edge within thirty percentage of chord length. In the simulation for the jet strength and the jet angle control, the same tendencies were observed to those of blowing cases.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.5
• /
• pp.685-694
• /
• 2002

Effects of wake-passing orientation and frequency on the wake-induced boundary layer transition on a NACA0012 airfoil are investigated. The wakes are generated by rotating cylinders clockwise (CW) and counterclockwise (CCW) around the airfoil. Time- and phase-averaged streamwise mean velocities and turbulent fluctuations are measured with a single hot-wire probe. Wall skin frictions are estimated by the Computational Preston Tube Method (CPM). The pressure distribution on the airfoil is different according to the wake-passing orientation and frequency. Turbulent patches are generated in the laminar boundary layer due to the passing wake and the boundary layer becomes temporarily transitional. The transition process is significantly affected by the pressure gradient and the turbulent patches. For the receding wake, the turbulent patches propagate more rapidly than those for the approaching wake because adverse pressure gradient becomes larger. As the frequency increases, onset location of transition moles upstream and the boundary layer near the trailing edge becomes more transitional.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.6
• /
• pp.776-785
• /
• 2001

Hot-wire measurements are performed in boundary layers developing on a NACA0012 airfoil over which wakes pass periodically. The Reynolds number based on chord length of the airfoil is 2$\times$10(sup)5 and the wakes are generated by circular cylinders rotating clockwise and counterclockwise around the airfoil. This paper and its companion Part II describe the phenomena of wake-induced transition of the boundary layers on the airfoil using measured data; phase-and time-averaged streamwise mean velocities, turbulent fluctuations, integral parameters and wall skin frictions. This paper describes the background and facility together with results of time-averaged quantities. Due to the passing wake with mean velocity defects and high turbulence intensities, the laminar boundary layer is periodically disturbed at the upstream station and becomes steady-state transitional boundary layer at the downstream station. The velocity defect in the passing wake changes the local pressure at the leading of the airfoil, significantly affects the time-mean pressure distribution on the airfoil and eventually, has influence on the transition process of the boundary layer.

• Advances in aircraft and spacecraft science
• /
• v.2 no.4
• /
• pp.427-443
• /
• 2015

The present paper is the follow-on of a former work in which the influence of the gas-surface interaction models was evaluated on the aerodynamic coefficients of an aero-space-plane and on a section of its wing. The models by Maxwell and by Cercignani-Lampis-Lord were compared by means of Direct Simulation Monte Carlo (DSMC) codes. In that paper the diffusive, fully accommodated, semi-specular and specular accommodation coefficients were considered. The results pointed out that the influence of the interaction models, considering the above mentioned accommodation coefficients, is pretty strong while the Cercignani-Lampis-Lord and the Maxwell models are practically equivalent. In the present paper, the comparison of the same models is carried out considering the dependence of the accommodation coefficients on the angle of incidence (or partial accommodation coefficients). More specifically, the normal and the tangential momentum partial accommodation coefficients, obtained experimentally by Knetchel and Pitts, have been implemented. Computer tests on a NACA-0012 airfoil have been carried out by the DSMC code DS2V-64 bits. The airfoil, of 2 m chord, has been tested both in clean and flapped configurations. The simulated conditions were those at an altitude of 100 km where the airfoil is in transitional regime. The results confirmed that the two interaction models are practically equivalent and verified that the use of the Knetchel and Pitts coefficients involves results very close to those computed considering a diffusive, fully accommodated interaction both in clean and flapped configurations.