• 한국군사과학기술학회지
• /
• 제3권2호
• /
• pp.23-34
• /
• 2000

Three dimensional, compressible, mass weighted averaging of Favre, Navier-Stokes system with k-$\varepsilon$ turbulence, is numerically discretized to compute three dimensional multiple jet interaction flow fields for a hybrid projectile containing three rocket motors in the ogive section. Numerical flow field computations have been made for angled nose jets and rockets at supersonic speed using multiblock structured grid. The jet conditions include very high jet to free stream pressure ratio and high temperature. It is shown that the strength of nozzle stagnation pressure affects the flow field near the side nozzle and the high stagnation pressure increases total amount of drag by a few percent. However, minor drag loss due to the pressure drag might be fully overcomed by an additional axial thrust. The results of present study can be applied for the design of future hybrid projectile.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2002년도 춘계 학술대회논문집
• /
• pp.29-34
• /
• 2002

TIn this study, the counter-jet flows which designed for improvement of aerodynamic performance of the blunt body vehicle have been analyzed. The variations of the drag force and jet penetration depth due to changes in the stagnation properties of counter jet new such as total pressure, mach number, and total temperature. The counter jet flow, which is injected toward incoming supersonic freestream at stagnation region of blunt cone-cylinder vehicle, have been studied by using upwind flux difference splitting navier-stokes method. The changes in the stagnation pressure and Mach number resulted in large effects on the wall pressure and drag force, on the other hand tile total temperature changes did not.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1999년도 유체기계 연구개발 발표회 논문집
• /
• pp.134-140
• /
• 1999

A three-dimensional CSCM upwind flux difference splitting Navier-stokes code with two-equation turbulence models was developed to predict the transonic flows in centrifugal compressor diffuser. The k-$\epsilon$ model of Abe et al. performed well in predicting the pressure distribution in the shock wave/turbulent boundary-layer interaction. Three turbulence models predicted the similar distribution of static pressure through the diffuser and showed a good agreement with the experimental results. The secondary flows in the corner were predicted well by these turbulence models. The pressure increase before the throat of the diffuser vane is important for the overall pressure recovery. As the mass flow rate increased the blockage decreased at the throat. The pressure coefficient distribution through the diffuser depended on the throat blockage not on the rotational speed of the impeller.

• International Journal of Aeronautical and Space Sciences
• /
• 제3권1호
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• 대한기계학회논문집B
• /
• 제25권10호
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. 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, it is concluded that Abid'k-$\omega$ 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 the same.

• 한국전산유체공학회지
• /
• 제19권3호
• /
• pp.91-97
• /
• 2014

A high resolution numerical method aimed at solving cavitating flow was proposed and applied to gas-liquid two-phase shock tube problem with arbitrary void fraction. The present method with compressibility effects employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. The Jacobian matrix from the inviscid flux of constitute equation is diagonalized analytically and the speed of sound for the two-phase media is derived by eigenvalues. So that the present method is appropriate for the extension of high order upwind schemes based on the characteristic theory. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results of high speed flow phenomena such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and solutions at isothermal condition are provided and discussed.

• 한국수자원학회논문집
• /
• 제44권6호
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.