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

We calculate the coordinates of an axisymmetric nozzle with a central body. This nozzle ensures a transonic flow with a plane sound surface, which is orthogonal to the symmetry axis and has a wall kink at the sonic point, The Chaplygin transformation in the subsonic part of the flow leads the Dirichlet problem for a system of nonlinear equations. The definition domain of the solution in the velocity-hodograph plane is taken as a rectangle. This enables one to obtain the nozzle with a monotonic distribution of velocity along its subsonic part. In the nonlinear differential equation, the linear Chaplygin operator for plane flows is separated, which allows the iterative calculation of the solution. The supersonic part of the nozzle is calculated under the assumption that the flow at the nozzle exit is uniform and parallel to the symmetry axis; i.e., the supersonic jet outflows to the submerged space with the same pressure. The calculation is performed by the characteristic method. The exact solution of Tricomi equation for near-sonic flows with the straight sonic line is used to 'move away' the sound plane. The velocity distribution alone the supersonic part of the nozzle is also monotonic, which ensures the absence of the boundary-layer separation and, therefore, the adequacy of the ideal-gas model. calculations show that the flow in the supersonic part of the nozzle is continuous (compression shocks are absent)

• International Journal of Aeronautical and Space Sciences
• /
• 제18권3호
• /
• pp.403-413
• /
• 2017

Dual-point aerodynamic design optimization is conducted for DLR-F6 wing-body-nacelle-pylon configuration adopting an efficient surface mesh movement method for complex junction geometries. A three-dimensional unstructured Euler solver and its discrete adjoint code are utilized for flow and sensitivity analysis, respectively. Considered design conditions are a low-lift condition and a cruise condition in a transonic regime. Design objective is to minimize drag and reduce shock strength at both flow conditions. Shape deformation is made by variation of the section shapes of inboard wing and pylon, nacelle vertical location and nacelle pitch angle. Hicks-Henne shape functions are employed for deformation of the section shapes of wing and pylon. By the design optimization, drag coefficients were remarkably reduced at both design conditions retaining specified lift coefficient and satisfying other constraints. Two-point design results show mixed features of the one-point design results at low-lift condition and cruise conditions.

• 한국군사과학기술학회지
• /
• 제25권1호
• /
• pp.45-54
• /
• 2022

There has been necessity to supplement the fin database to improve the accuracy of low-fidelity aerodynamic solver for missile configuration. In this study, fin database is expanded by in-house solver, utilized in the triservice data the previously established into regions beyond means of CFD. Fin alone data of CFD analysis results in the original region is matched well with triservice data originated from the wind tunnel tests. Extensive fin aerodynamic data from CFD analysis is added to the existing database of the low-fidelity solver. For confirmation, aerodynamic characteristics of body-tail and body-canard-tail missile configurations is computed using upgraded low-fidelity solver at transonic region. The result using improved solver shows good agreements with wind tunnel test and CFD analysis results, which implies that it becomes more accurate.

• 동력기계공학회지
• /
• 제3권4호
• /
• pp.16-21
• /
• 1999

Numerical simulation on two-dimensional turbine cascade flow has been performed using compressible Navier-Stokes equations. The flow equations are written in a cartesian coordinate system, then mapped into a generalized body-fitted ones. All direction of viscous terms are incoporated and turbulent effects are modeled using the extended ${\kappa}-{\epsilon}$ model. Equations are discretized using control volume SIMPLE algorithm on the nonstaggered grid sysetem. Applications are made at a VKI turbine cascade flow in atransonic wind-tunnel and compared to experimental data. Present numerical results are shown to be in good agreement with the experimental results and simulate the compressible viscous flow characteristics inside the turbine blade passage.

• 대한기계학회논문집
• /
• 제16권2호
• /
• pp.313-324
• /
• 1992

Two-dimensional Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system, then mapped into a generalized body-fitted coordinate system. All direction of viscous terms are incorporated and turbulent effects are modeled using the Baldwin-Lomax algebraic model. Equation are discretized using finite difference method on the C-type grids and solved using implicit LU-ADI decomposition scheme. Calculations are made at a VKI turbine cascade flow in a transonic wind-tunnel and compared to experimental data. Present numerical scheme is shown to be in good agreement with the previous experimental results and simulates the two-dimensional viscous flow phenomena.

• Advances in aircraft and spacecraft science
• /
• 제6권4호
• /
• pp.283-296
• /
• 2019

The influence of different planform parameters on the aerodynamic performance of large-scale subsonic and transonic Blended Wing Body (BWB) aircraft have gained comprehensive research in the recent years, however, it is not the case for small-size low subsonic speed Unmanned Aerial Vehicles (UAVs). The present work numerically investigates aerodynamics governing four different trailing edge geometries characterizing BWB configurations in standard flight conditions at angles of attack from $-4^{\circ}$ to $22^{\circ}$ to provide generic information that can be essential for making well-informed decisions during BWB UAV conceptual design phase. Simulation results are discussed and comparatively analyzed with useful implications for formulation of proper mission profile specific to every BWB configuration.

• 한국전산유체공학회지
• /
• 제20권1호
• /
• pp.47-56
• /
• 2015

In this paper, we study the effect of various turbulence models by comparing the aerodynamic characteristics and the flow patterns computed for aircraft models. An in-house CFD solver, MSAPv, that solves the three dimensional RANS equations with the turbulence model equations is used. The turbulence models used in this study are the Spalart-Allmaras model, Menter's $k-{\omega}$ SST model, Coakley's $q-{\omega}$ model, and Huang and Coakley's $k-{\varepsilon}$ model. DLR-F6 WB and WBNP configurations are selected for the study. We concentrate on the separated flow pattern variations with the turbulence models at the wing-body junction and the wing-pylon junction as well as drag polar curves.

• International Journal of Aeronautical and Space Sciences
• /
• 제8권1호
• /
• pp.95-104
• /
• 2007

A new design approach of complex geometries such as wing/body configuration is arranged by using overset mesh techniques under large scale computing environment. For an in-depth study of the flow physics and highly accurate design, several special overlapped structured blocks such as collar grid, tip-cap grid, and etc. which are commonly used in refined drag prediction are adopted to consider the applicability of the present design tools to practical problems. Various pre- and post-processing techniques for overset flow analysis and sensitivity analysis are devised or implemented to resolve overset mesh techniques into the design optimization problem based on Gradient Based Optimization Method (GBOM). In the pre-processing, the convergence characteristics of the flow solver and sensitivity analysis are improved by overlap optimization method. Moreover, a new post-processing method, Spline-Boundary Intersecting Grid (S-BIG) scheme, is proposed by considering the ratio of cell area for more refined prediction of aerodynamic coefficients and efficient evaluation of their sensitivities under parallel computing environment. With respect to the sensitivity analysis, discrete adjoint formulations for overset boundary conditions are derived by a full hand-differentiation. A smooth geometric modification on the overlapped surface boundaries and evaluation of grid sensitivities can be performed by mapping from planform coordinate to the surface meshes with Hicks-Henne function. Careful design works for the drag minimization problems of a transonic wing and a wing/body configuration are performed by using the newly-developed and -applied overset mesh techniques. The results from design applications demonstrate the capability of the present design approach successfully.

• 대한기계학회논문집B
• /
• 제21권2호
• /
• pp.316-328
• /
• 1997

There were appreciable progresses on the study of shock wave / boundary layer interaction control in the transonic flow without nonequilibrium condensation. But in general, the actual flows associated with those of the airfoil of high speed flight body, the cascade of steam turbine and so on accompany the nonequilibrium condensation, and under a certain circumstance condensation shock wave occurs. Condensation shock wave / boundary layer interaction control is quite different from that of case without condensation, because the droplets generated by the result of nonequilibrium condensation may clog the holes of the porous wall for passive control and the flow interaction mechanism between the droplets and the porous system is concerned in the flow with nonequilibrium condensation. In these connections, it is necessary to study the condensation shock wave / boundary layer interaction control by passive cavity in the flow accompanying nonequilibrium condensation with condensation shock wave. In the present study, experiments were made on a roof mounted half circular arc in an indraft type supersonic wind tunnel to evaluate the effects of the porosity, the porous wall area and the depth of cavity on the pressure distribution around condensation shock wave. It was found that the porosity of 12% which was larger than the case of without nonequilibrium condensation produced the largest reduction of pressure fluctuations in the vicinity of condensation shock wave. The results also showed that wider porous area, deeper cavity for the same porosity of 12% are more favourable "passive" effect than the cases of its opposite. opposite.

• 한국음향학회지
• /
• 제37권4호
• /
• pp.163-173
• /
• 2018

본 논문에서는 비행 중 비행체 표면에 작용하는 음향하중 예측을 수행하였다. 비행 중 음향하중은 비행체 표면의 압력 변동에 의해 발생한다. 기존의 비행 중 음향하중 예측방법은 반경험적 방법으로 이론과 실험 결과를 기반으로 도출한 경험식을 활용한다. 하지만 경험식의 입력 값으로 사용되는 비행체 주변 유동특성 및 경계층 파라미터를 매번 실험을 통해 얻는 것에는 한계가 있다. 따라서 본 논문에서는 전산유체해석(Computational Fluid Dynamics, CFD) 결과를 반경험적 방법과 혼합하는 하이브리드 방법을 이용하여 비행 중 비행체에 작용하는 음향하중을 예측하였다. Cone-cylinder-flare 형상 비행체에 대해 아음속, 천음속, 초음속, 최대동압도달(Maximum dynamic pressure, Max-q) 시점의 비행 환경에 대한 음향하중 예측을 수행하였다. 하이브리드 방법 적용 시 전산유체해석결과를 기반으로 한 경계층 끝단 영역 판단 방법에 대해 비교하였고 여러 연구자에 의해 제시된 경험식에 따른 음향하중 예측결과를 비교하였다.