• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• pp.179-186
• /
• 1999

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67), and to evaluate the performances of k-$\epsilon$ and Baldwin-Lomax turbulence models. A finite volume method is used for spatial discretization. And, the equations are solved implicitly in time with the use of approximate factorization. Upwind difference scheme is used for inviscid terms, but viscous terms are centrally differenced. The flux-difference-splitting of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. And, the results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, we conclude that k-$\epsilon$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost same.

• Journal of the Korean Society of Safety
• /
• v.15 no.3
• /
• pp.96-101
• /
• 2000

The effects of internal fluid motion have to be considered in the analysis of liquid storage containers. Therefore this thesis developed a three-dimensional boundary element-finite element method for the analysis of rectangular liquid storage containers. The irrotational motion of inviscid and incompressible ideal fluid is modeled by using boundary elements and the motion of structure by finite elements. Coupling is performed by using compatibility and equilibrium conditions along the interface. Dynamic response characteristics of rectangular liquid storage containers such as sloshing motion, hydrodynamic pressure, displacement by fluid-structure interaction are investigated.

• 한국가시화정보학회:학술대회논문집
• /
• 2006.12a
• /
• pp.149-152
• /
• 2006

This paper describes the flow visualization around the model of a commercial passenger airplane, Boeing 747-400, which cruises in the transonic speed. The geometry was realized through the reverse engineering based on the photogrammetry. The three-dimensional inviscid steady compressible governing equations are solved in the unstructured grid system under the cruise condition and in a finite volume method. The convective term is processed by the Crank-Nicholson scheme and first order upwind scheme is applied. The lift and drag forces in the wing with engines increase by 1.49% End 3.9%, respectively compared with the wing without engines.

• Advances in aircraft and spacecraft science
• /
• v.10 no.5
• /
• pp.457-471
• /
• 2023

The paper evaluates the aerodynamic coefficients on a blunt-nose re-entry capsule with a conical cross-section followed by a cone-flare body. A computer code is developed to solve three-dimensional compressible inviscid equationsfor flow over a Space Recovery Experiment (SRE) configuration at different flare-cone half-angle at Mach 6 and angle of attack up to 5°, at 1° interval. The surface pressure variation is numerically integrated to obtain the aerodynamic forces and pitching moment. The numerical analysis reveals the influence of flare-cone geometry on the flow characteristics and aerodynamic coefficients. The numerical results agree with wind tunnel results. Increase of cone-flare angle from 25° to 35° results in increase of normal force slope, axial forebody drag, base drag and location of centre of pressure by 62.5%, 56.2% and 33.13%, respectively, from the basic configuration ofthe SRE of 25°.

• Journal of Mechanical Science and Technology
• /
• v.15 no.5
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.3
• /
• pp.296-306
• /
• 2012

A series of numerical simulations are carried out to analyze a supersonic inlet buzz, which is an unsteady pressure oscillation phenomenon around a supersonic inlet. A simple but efficient geometry, experimentally adopted by Nagashima, is chosen for the analysis of unsteady flow physics. Among the two sets of simulations considered in this study, the effects of various throttling conditions are firstly examined. It is seen that the major physical characteristic of the inlet buzz can be obtained by inviscid computations only and the computed flow patterns inside and around the inlet are qualitatively consistent with the experimental observations. The dominant frequency of the inlet buzz increases as throttle area decreases, and the computed frequency is approximately 60Hz or 15% lower than the experimental data, but interestingly, this gap is constant for all the test cases and shock structures are similar. Secondly, inviscid calculations are performed to examine the effect regarding angle of attack. It is found that patterns of pressure oscillation histories and distortion due to asymmetric (or three-dimensional) shock structures are substantially affected by angle of attack. The dominant frequency of the inlet buzz, however, does not change noticeably even in regards to a wide range of angle of attacks.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.3
• /
• pp.160-166
• /
• 2013

Recently submerged objects moving at high speed such as a supercavitating torpedo have been studied for their practical advantage of the dramatic drag reduction. In this study, we are focusing our attention on supercavitating flows around axisymmetric cavitators; a numerical method based on inviscid flow is developed and predicted supercavities around several shapes of 2D and 3D cavitators are presented. The results are validated by comparison of existing theoretical and empirical results. In addition, characteristic features of supercavity shapes and drag forces acting on a real scale torpedo are evaluated according to practically appropriate operating conditions.

• Journal of computational fluids engineering
• /
• v.14 no.1
• /
• pp.70-77
• /
• 2009

In this study, time-dependent numerical analysis was carried out to investigate the plasma jet impingement on a flat plate, and a compressible form of two-dimensional inviscid gas dynamics equations were solved using the flux corrected transport algorithm. The mathematical modeling of Joule heating in the polycarbonate capillary bore and the mass ablation from the bore wall was incorporated in the numerical analysis and the series of computation was performed for three cases depending on the distance of the opposing plate from the capillary exit. The computational results reveal that the presence of the opposing plate does not affect the flow conditions inside the capillary when compared to the case of open-air plasma discharge. In the exterior region, the flow structure shows the typical supersonic underexpanded jet which consists of the strong Mach disk in front of the opposing plate and the barrel shock at the side of the jet. It is found that the shock evolution becomes more quasi-steady when the plate distance decreases. Also, the effects of the distance between the capillary bore exit and the opposing plate on the flow conditions along the opposing plate are investigated and the pressure variation on the plate shows more complicated interaction between the plasma discharge and the opposing plate when the location of plate becomes closer to the capillary exit.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.8
• /
• pp.8-17
• /
• 2005

It is well known that preconditioning methods are efficient for convergence acceleration in the compressible low Mach number flows. In this study, the original Euler equations and three differently nondimensionalized preconditioning methods are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as numerical flux discretization and time integration, respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning allows Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• International Journal of Aeronautical and Space Sciences
• /
• v.11 no.1
• /
• pp.1-9
• /
• 2010

Numerical simulation of unsteady flows around helicopters was conducted to investigate the aerodynamic interaction of main rotor and other components such as fuselage and tail rotor. For this purpose, a three-dimensional inviscid flow solver has been developed based on unstructured meshes. An overset mesh technique was used to describe the relative motion between the main rotor, and other components. As the application of the present method, calculations were made for the rotor-fuselage aerodynamic interaction of the ROBIN (ROtor Body INteraction) configuration and for a complete UH-60 helicopter configuration consisted of main rotor, fuselage, and tail rotor. Comparison of the computational results was made with measured time-averaged and instantaneous fuselage surface pressure distributions for the ROBIN configuration and thrust distribution and available experimental data for the UH-60 configuration. It is demonstrated that the present method is efficient and robust for the simulation of complete rotorcraft configurations.