• Wind and Structures
• /
• 제34권1호
• /
• pp.115-125
• /
• 2022

The high-Reynolds number flow around a rectangular cylinder, having streamwise to crossflow length ratio equal to 5 is analyzed in the present paper. The flow is characterized by shear-layer separation from the upstream edges. Vortical structures of different size form from the roll-up of these shear layers, move downstream and interact with the classical vortex shedding further downstream in the wake. The corresponding mean flow is characterized by a recirculation region along the lateral surface of the cylinder, ending by mean flow reattachment close to the trailing edge. The mean flow features on the cylinder side have been shown to be highly sensitive to set-up parameters both in numerical simulations and in experiments. The results of 21 Large Eddy Simulations (LES) are analyzed herein to highlight the impact of the lateral mean recirculation characteristics on the near-wake flow features and on some bulk quantities. The considered simulations have been carried out at Reynolds number Re=DU_∞/ν=40 000, being D the crossflow dimension, U_∞ the freestream velocity and ν the kinematic viscosity of air; the flow is set to have zero angle of attack. Some simulations are carried out with sharp edges (Mariotti et al. 2017), others with different values of the rounding of the upstream edges (Rocchio et al. 2020) and an additional LES is carried out to match the value of the roundness of the upstream edges in the experiments in Pasqualetto et al. (2022). The dimensions of the mean recirculation zone vary considerably in these simulations, allowing us to single out meaningful trends. The streamwise length of the lateral mean recirculation and the streamwise distance from the upstream edge of its center are the parameters controlling the considered quantities. The wake width increases linearly with these parameters, while the vortex-shedding non-dimensional frequency shows a linear decrease. The drag coefficient also linearly decreases with increasing the recirculation length and this is due to a reduction of the suctions on the base. However, the overall variation of C_D is small. Finally, a significant, and once again linear, increase of the fluctuations of the lift coefficient is found for increasing the mean recirculation streamwise length.

• Advances in aircraft and spacecraft science
• /
• 제4권5호
• /
• pp.503-514
• /
• 2017

The attitude control of an aircraft is usually fulfilled by means of thrusters at high altitudes. Therefore, the possibility of using also aerodynamic surfaces would produce the advantage of reducing the amount of fuel for the thrusters to be loaded on board. For this purpose, Zuppardi already considered some aerodynamic problems linked to the use of a wing flap in a previous paper. A NACA 0010 airfoil with a trailing edge flap of 35% of the chord, in the range of angle of attack 0-40 deg and flap deflections up to 30 deg was investigated. Computer tests were carried out in hypersonic, rarefied flow by a direct simulation Monte Carlo code at the altitudes of 65 and 85 km of Earth Atmosphere. The present work continues this subject, considering the same airfoil and free stream conditions but two flap extensions of 45% and 25% of the chord and two flap deflections of 15 and 30 deg. The main purpose is to compare the influence of the flap dimension with that of the flap deflection. The present analysis is carried out in terms of: 1) percentage variation of the global aerodynamic coefficients with respect to the no-flap configuration, 2) increment of pressure and heat flux on the airfoil lower surface due to the Shock Wave-Shock Wave Interaction (SWSWI) with respect to the same quantities with no SWSWI or in no-flap configuration, 3) flap hinge moment. Issues 2) and 3) are important for the design of the mechanical and thermal protection system and of the flap actuator, respectively. Under the above mentioned test and geometrical conditions, the flap deflection is aerodynamically more effective than the flap extension, because it involves higher variation of the aerodynamic coefficients. However, tests verify that a smaller deflection angle involves the advantage of a smaller increment of pressure and heat flux on the airfoil lower surface, due to SWSWI, as well as a smaller hinge moment.

• 한국항공우주학회지
• /
• 제48권10호
• /
• pp.745-752
• /
• 2020

SpaceX Falcon 9에 탑재된 초음속 그리드핀 유동에 관한 수치해석을 통하여 해당 그리드핀의 공력특성을 평가하였다. 3차원 정상 수치해석 계산의 효율성 향상을 위하여 단위 그리드핀 개념을 도입하였다. 본 연구에서는 발사체의 받음각, 외부 프레임, 그리고 커넥팅 로드 등의 영향성을 보정하는 전후 처리과정을 개선하였고, 이에 따라 최종 공력계산의 정확도가 과거 연구에 비하여 향상되었음을 확인하였다. 이후 천음속과 초음속 비행조건에서 다양한 받음각을 갖는 SpaceX Falcon 9 그리드핀의 공력특성을 평가하였다.

• Ocean Systems Engineering
• /
• 제6권3호
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• 수산해양기술연구
• /
• 제27권4호
• /
• pp.286-292
• /
• 1991

현재 사용되고 있는 전개판의 성능을 보다 향상시키기 위한 방안으로 비행기의 익이론에 기초를 두고, 전개판의 경계층을 흡입하거나 제거하는 방법 중에서 역류역을 제거하는 와발생기의 부착을 전개판에 적용하였다. 성능분석을 위해서 단순만곡형, V-만곡형 그리고 슬롯만곡형 전개판의 3종에 와발생기를 부착하지 않은 경우와 부착한 경우에 대하여 성능 및 유체역학적 특성을 성능실험과 가시화실험을 통하여 비교해 보았다. 실험결과를 요약하면 다음과 같다. \circled1 단순만곡형전개판에 와발생기를 부착한 경우 최대전개력계수는 10~15% 증가, 항력계수는 2%미만증가, 양항비는 5~20% 증가하였다. \circled2 V-만곡형전개판에 와발생기를 부착한 경우 양항비가 10~20% 증가하였다. \circled3 슬롯만곡형에 와발생기를 부착한 경우 영각 20$^{\circ}$ 이내에서 최대전개력계수는 약 20%이상 증가, 항력계수는 5~20% 증가, 양항비는 약간 높게 나타났으나, 영각 25$^{\circ}$이상에서는 이들 계수의 값들이 기준형에 비하여 낮게 나타났다. \circled4 와발생기를 부착한 전개판에서 박리점이 기준형에 비하여 다소 후연에 나타났다. \circled5 와발생기에 따른 전.후면의 유속차는 단순만곡형의 경우 모든 영각에서 약 10% 증가, V-만곡형의 경우 전반적으로 증가하였으나 영각 25$^{\circ}$에서는 비슷하게, 슬롯만곡형의 경우 영각 25$^{\circ}$까지는 증가하는 것으로 나타났다. \circled6 와발생기에 따른 박리역의 크기는 단순만곡형에서 약 5~15% 감소, V-만곡형과 슬롯만곡형은 영각 15$^{\circ}$~20$^{\circ}$에서는 약 10% 감소, 영각 25$^{\circ}$ 이상일 때는 비슷하게 나타났다.

• Structural Engineering and Mechanics
• /
• 제75권4호
• /
• pp.487-496
• /
• 2020

Most of the previous works on numerical analysis of galloping of transmission lines are generally based on the quasisteady theory. However, some wind tunnel tests of the rectangular section or hangers of suspension bridges have shown that the galloping phenomenon has a strong unsteady characteristic and the test results are quite different from the quasi-steady calculation results. Therefore, it is necessary to check the applicability of the quasi-static theory in galloping analysis of the ice-covered transmission line. Although some limited unsteady simulation researches have been conducted on the variation of parameters such as aerodynamic damping, aerodynamic coefficients with wind speed or wind attack angle, there is a need to investigate the numerical simulation of unsteady galloping of two-dimensional iced transmission line with comparison to wind tunnel test results. In this paper, it is proposed to conduct a two dimensional (2-D) unsteady numerical analysis of ice-covered transmission line galloping. First, wind tunnel tests of a typical crescent-shapes iced conductor are conducted firstly to check the subsequent quasisteady and unsteady numerical analysis results. Then, a numerical simulation model consistent with the aeroelastic model in the wind tunnel test is established. The weak coupling methodology is used to consider the fluid-structure interaction in investigating a two-dimension numerical simulation of unsteady galloping of the iced conductor. First, the flow field is simulated to obtain the pressure and velocity distribution of the flow field. The fluid action on the iced conduct at the coupling interface is treated as an external load to the conductor. Then, the movement of the conduct is analyzed separately. The software ANSYS FLUENT is employed and redeveloped to numerically analyze the model responses based on fluid-structure interaction theory. The numerical simulation results of unsteady galloping of the iced conduct are compared with the measured responses of wind tunnel tests and the numerical results by the conventional quasi-steady theory, respectively.

• International Journal of Aeronautical and Space Sciences
• /
• 제18권1호
• /
• pp.8-16
• /
• 2017

Most of small air vehicles with moving wing fly at low Reynolds number condition and the reduced frequency of the moving wing ranges from 0.0 to 1.0. The physical phenomena over the wing dramatically vary with the reduced frequency. This study examines experimentally the effect of the reduced frequency at low Reynolds number. The NACA0012 airfoil performs sinusoidal pitching motion with respect to the quarter chord with the four reduced frequencies of 0.1, 0.2, 0.4 and 0.76 at the Reynolds number $2.3{\times}10^4$. Smoke-wire flow visualization, unsteady surface pressure measurement, and unsteady force calculation are conducted. At the reduced frequency of 0.1 and 0.2, various boundary layer events such as reverse flow, discrete vortices, separation and reattachment change the amplitude and the rotation direction of the unsteady force hysteresis. However, the boundary layer events abruptly disappear at the reduced frequency of 0.4 and 0.76. Especially at the reduced frequency of 0.76, the local variation of the unsteady force with respect to the angle of attack completely vanishes. These results lead us to the conclusion that the unsteady aerodynamic characteristics of the reduced frequency of 0.2 and 0.4 are clearly distinguishable and the unsteady aerodynamic characteristics below the reduced frequency of 0.2 are governed by the boundary layer events.

• 항공우주시스템공학회지
• /
• 제12권4호
• /
• pp.57-63
• /
• 2018

적외선 온도 측정 기법은 비접촉식 방법으로 모형의 표면 온도를 가시화할 수 있는 기법이다. 그러나 획득할 수 있는 결과는 2차원 온도 결과로 정량적인 결과를 획득하기에는 한계가 있다. 본 연구는 3차원 매핑 기법을 적외선 온도 측정 기법에 적용하는 것이 목표이다. 풍동 실험은 국방과학연구소에서 보유하고 있는 중형 아음속 풍동에서 수행했으며, 대상 모형은 오자이브 실린더이다. 시험 조건은 유속 20 m/s에서 80 m/s, 받음각은 $0^{\circ}{\sim}90^{\circ}$이다. 3차원 매핑 기법은 마커를 이용하여 실제 모형의 위치정보와 적외선 이미지 상의 위치 정보를 대응시키는 방법을 사용하였다. 그 결과, 모형의 박리점이 이론적인 값과 매우 일치하는 것을 확인하였다.

• 항공우주기술
• /
• 제11권2호
• /
• pp.65-72
• /
• 2012

프로펠러와 고양력장치를 장착한 터보프롭 항공기에 대한 실속 특성 분석을 위해 수치 해석을 수행하였다. 항공기의 실속 특성은 프로펠러와 고양력 장치의 장착 조합에 따른 형상별 전산해석 결과를 통해 정성적으로 분석하였다. 실속 특성 해석은 Spalart-Allmaras 난류 모델을 기반으로 한 3차원 Navier-Stokes 방정식 해법을 이용하였으며 프로펠러의 회전은 슬라이딩 격자기법을 이용하여 모사하였다. 분석 결과 순항 형상의 경우 동체/날개 페어링에서 주요 유동박리가 발생하며 프로펠러 후류로 인해 점차 감소함을 알 수 있었다. 고양력장치를 장착한 경우 나셀 바깥쪽에서 주요 유동박리 현상이 발생하였고 프로펠러가 회전하는 경우에도 상대속도 감소와 유효 받음각 증가로 나셀 바깥쪽 날개 부분은 조기에 실속에 잠김을 알 수 있었다. 프로펠러는 날개의 inboard에서 하강하는 회전 방향이 프로펠러 후류로 인한 실속 지연 측면에서 유리함을 알 수 있었다.

• Wind and Structures
• /
• 제17권6호
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.