• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.2
• /
• pp.15-22
• /
• 2003

Computational study using the 2-dimensional, compressible, Navier-Stokes equations is performed to predict the discharge coefficient of air flow through a film-cooling hole. In order to investigate the effect of internal/external flows on discharge coefficient, the present computational results which are obtained for three flow cases, only external flow, only internal flow, and no flow, are compared with experimental ones. It is found that the computational results predict the discharge coefficient of the film-cooling hole in a reasonable accuracy and the external crossflow reduces the discharge coefficient, while the internal crossflow increases the discharge coefficient in a range of momentum flux ratio $I_{c-jet}$ > 1 due to the total pressure loss and boundary layer effect.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.9
• /
• pp.805-815
• /
• 2011

Nonlinear Parabolized Stability Equations(NSPE) can be effectively used to study more throughly the transition process. NPSE can efficiently analyze the stability of a nonlinear region in transition process with low computational cost compared to Direct Numerical Simulation(DNS). In this study, NPSE in general coordinate system is formulated and a computer code to solve numerically the equations is developed. Benchmark problems for incompressible and compressible boundary layers over a flat plate are analyzed to validate the present code. It is confirmed that the NPSE methodology constructed in this study is an efficient and effective tool for nonlinear stability analysis.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.5
• /
• pp.323-334
• /
• 2020

Since boundary layer transition has a significant impact on the aero-thermodynamic performance of hypersonic flight vehicles, capability of accurate prediction of transition location is essential for design and performance analysis. In this study, γ-Re_θt model is improved to predict transition of hypersonic boundary layers and validated. A coefficient in the production term of the intermittency transport equation that affects the transition onset location is constructed and applied as a function of Mach number, wall temperature, and freestream stagnation temperature based on the similarity numerical solution of compressible boundary layer. To take into account a Mach number dependency of transition onset momentum thickness Reynolds number and transition length, additional correlation equations are determined as function of Mach number and applied to Re_θc and F_length correlations of the baseline model. The suggested model is implemented to a commercial CFD code in consideration of practical use. Analysis of hypersonic flat plate and circular cone boundary layers is carried out by using the model for validation purpose. An improvement of prediction capability with respect to variation of Mach number and unit Reynolds number is identified from the comparison with experimental data.

• Steel and Composite Structures
• /
• v.26 no.6
• /
• pp.771-791
• /
• 2018

This paper deals with the low velocity impact response and dynamic stresses of composite sandwich truncated conical shells (STCS) with compressible or incompressible core. Impacts are assumed to occur normally over the top face-sheet and the interaction between the impactor and the structure is simulated using a new equivalent three-degree-of-freedom (TDOF) spring-mass-damper (SMD) model. The displacement fields of core and face sheets are considered by higher order and first order shear deformation theory (FSDT), respectively. Considering continuity boundary conditions between the layers, the motion equations are derived based on Hamilton's principal incorporating the curvature, in-plane stress of the core and the structural damping effects based on Kelvin-Voigt model. In order to obtain the contact force, the displacement histories and the dynamic stresses, the differential quadrature method (DQM) is used. The effects of different parameters such as number of the layers of the face sheets, boundary conditions, semi vertex angle of the cone, impact velocity of impactor, trapezoidal shape and in-plane stresses of the core are examined on the low velocity impact response of STCS. Comparison of the present results with those reported by other researchers, confirms the accuracy of the present method. Numerical results show that increasing the impact velocity of the impactor yields to increases in the maximum contact force and deflection, while the contact duration is decreased. In addition, the normal stresses induced in top layer are higher than bottom layer since the top layer is subjected to impact load. Furthermore, with considering structural damping, the contact force and dynamic deflection decrees.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1144-1153
• /
• 2002

A hybrid performance prediction method is proposed in the present study. A channel diffuser is divided into four subregions: vaneless space, semi-vaneless space, channel, and channel exit region. One-dimensional compressible core flow and boundary layer calculation of each region with an incidence loss model and empirical correlation of residuary pressure recovery coefficient of a channel predict the performance of diffusers. Three channel diffusers are designed and tested for validating the developed prediction method. The pressure distributions from an impeller exit to the channel diffuser exit are measured and discussed for various operating conditions from choke to nearly surge conditions. The strong non-uniform pressure distribution which is caused by impeller-diffuser interaction is obtained over the vaneless and semi-vaneless spaces. The predicted performance shows good agreement with the measured performance of diffusers at a design condition as well as at off-design conditions.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.1-7
• /
• 2003

The V2F turbulence model, which has shown very good performance in several test cases at low speeds, has been applied to supersonic ramp flow with 20. corner angle at the free stream Mach number of 2.79. The flow is known to manifest strong shock wave/turbulent boundary layer interactions. As a comparative study, low-Reynolds k-ε models are also considered. While the V2F model predicts wall-pressure distribution well, it relatively predicts larger separation bubble and higher skin-friction after the reattachment than the experimental data. Although the ellpticity of f equation is the characteristics of incompressible flows, the converged solutions are acquired in the compressible flow with shock waves. The effect of the realizability constraints used in the model is also examined. In contrast to the result of impinging jet flows, the realizability bounds proposed by Durbin deterioate the overall solutions of the supersonic ramp flow.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.6
• /
• pp.1282-1291
• /
• 1989

본 연구에서는 이러한 연구의 시작단계로 노즐내의 열 경계층 방정식을 정리하고 압축성 와점도모델을 도입하여 유한 차분법을 이용한 전산프로그램을 작성 하였다. 한편 실제로 추진모테에서는 고온 고속의 열기류가 분출되기 때문에 추진 하고 있는 동안 노즐벽을 형성하고 있는 내열재가 심하게 손상되어 표면의 상태가 매우 거칠게 된다. 더구나 경우에 따라서는 내열재가 용발(ablation)하게 된다. 이러한 상태를 감안하여 마찰계수와 열전달 계수를 합리적으로 추정해야만 노즐의 설계와 주변장치를 합리적으로 수행할 수 있다. 따라서 본 연구에는 주로 경계층 내의 압력구배, 압축성의 효과, 물성치의 변화를 고려한 기존 난류모델에 근거하여 프로그램을 작성하고, 이것을 토대로 노즐표면 조도의 영향 및 분출(blowing)의 영 향을 중점적으로 고려하여 그 특성을 연구하였다. 노즐벽에서 분출을 고려한 이유는 표면이 용발할 때 표면의 온도가 거의 일정하게 유지된 상태로 노즐표면이 화학작용을 수반하면서 가스화됨을 초보적으로 고려해보기 위함이다.

• Journal of Industrial Technology
• /
• v.15
• /
• pp.153-168
• /
• 1995

The purpose of this study is to develop a method for predicting the aerodynamic performance of the subsonic airfoils in the 2-dimensional, steady and viscous flow. For this study, the airfoil geometry is specified by adopting the longest chord line system and by considering local surface curvature. In case of the inviscid-incompressible flow, the analysis is accomplished by the linearly varying strength vortex panel method and the Karman-Tsien correction law is applied for the inviscid-compressible flow analysis. The Goradia's integral method and the Truckenbrodt integral method are adopted for the boundary layer analysis of the laminar flow and the turbulent flow respectively. Viscous and inviscid solutions are converged by the Lockheed iterative calculating method using the equivalent airfoil geometry. And the analysis of the seperated flow is performed using the Dvorak and Maskew's method as the basic method. The wake effect is also considered and its geometry expressed by the formula of Summey & Smith when no seperation occurs. A computational efficiency is verified by the comparison of the computational results with experimental data and by the shorter execution time.

• Journal of Ship and Ocean Technology
• /
• v.6 no.2
• /
• pp.37-50
• /
• 2002

In the present numerical simulation of viscous free surface flow around a ship, two-fluids in-compressible Reynolds-averaged Navier-Stokes equations with the standard $\textsc{k}-\varepsilon$turbulence model are discretized on a regular grid by using a finite volume method. A local level-set method is introduced for capturing the free surface movement and the influence of the viscous layer and dynamic boundary condition of the free surface are implicitly considered. Partial differential equations in the level-set method are discretized with second order ENO scheme and explicit Euler scheme in the space and time integration, respectively. The computational results for the Series-60 model with $C_B=0.6$ show a good agreement with the experimental data, but more validation studies for commercial complicated hull forms are necessary.