• Journal of Mechanical Science and Technology
• /
• v.18 no.12
• /
• pp.2190-2203
• /
• 2004

A combined procedure for two-dimensional Delaunay mesh generation algorithm and an adaptive remeshing technique with higher-order compressible flow solver is presented. A pseudo-code procedure is described for the adaptive remeshing technique. The flux-difference splitting scheme with a modified multidimensional dissipation for high-speed compressible flow analysis on unstructured meshes is proposed. The scheme eliminates nonphysical flow solutions such as the spurious bump of the carbuncle phenomenon observed from the bow shock of the flow over a blunt body and the oscillation in the odd-even grid perturbation in a straight duct for the Quirk's odd-even decoupling test. The proposed scheme is further extended to achieve higher-order spatial and temporal solution accuracy. The performance of the combined procedure is evaluated on unstructured triangular meshes by solving several steady-state and transient high-speed compressible flow problems.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.60-68
• /
• 2005

For study on the unsteady wall interference effect, flows around a forced oscillating airfoil in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The Spalart-Allmaras one-equation model is employed for the turbulence effect. The computed results of the oscillating airfoil having a thin wake showed that the lift curve slope is increased and the magnitude of hysteresis loop is reduced by the interference effects. Since the vortex around the airfoil is generated and convected downstream faster than the free-air condition, the phase of lift, drag and pitching moment coefficients was shifted. The pressure on the test section wall shows harmonic terms having the oscillating frequency contained in the wail effect.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.3 no.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.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.137-143
• /
• 1999

The unsteady transonic viscous flow has been analyzed over an oscillatory airfoil. The CSCM(Conservative Supra Characteristic Method) upwind flux difference splitting method and the iterative time marching scheme having first order accuracy in time and second to third order accuracy in space was applied on dynamic meshes. A steady flow field of Mach number 0.7 has been calculated for the verification of unsteady algorithm. The time-accurate unsteady calculations have been done for NACA 0012 airfoil oscillating around quarter chord about freestream Mach number 0.6 on dynamic meshes. The results have been compared with the AGARD Case 3 experimental data. The periodic characteristics have been compared with the experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• 유체기계공업학회:학술대회논문집
• /
• 1999.12a
• /
• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.249-257
• /
• 2007

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results 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 exact solutions are provided and discussed.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• Journal of computational fluids engineering
• /
• v.20 no.1
• /
• pp.1-9
• /
• 2015

A developed code is applied to simulate relative motion of floating bodies in a side-by-side arrangement, including effects of a fender and a hawser. The developed code is based on the flux-difference splitting scheme for immiscible incompressible fluids and the hybrid Cartesian/immersed boundary method. To validate the developed code for free surface flows around deforming boundaries, the water wave generation is simulated, which is caused by bed movement. The computed wave profile and time histories of wave elevation are compared with other experimental and computational results. The effects of a fender and a hawser are modeled by asymmetric force acting on the floating bodies according to a relative displacement with the bounds, in which the fender and the hawser exert no force on the bodies. It has been observed that the floating body can be accelerated by a gap flow due to a phase difference caused by the free surface. Grid independency is established for the computed time history of the body velocity, based on three different size grids.