• Journal of computational fluids engineering
• /
• v.19 no.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.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.37-45
• /
• 2002

A two-dimensional unsteady, inviscid flow solver has been developed for the simulation of airfoil-vortex interactions on unstructured dynamically adapted meshes. The Euler solver is based on a second-order accurate implicit time integration using a point Gauss-Seidel relaxation scheme and a dual time-step subiteration. A vertex-centered, finite-volume discretization is used in conjunction with the Roe's flux-difference splitting. An unsteady solution-adaptive dynamic mesh scheme is used by adding and deleting mesh points to take account of both spatial and temporal variations of the flow field. The effect of vortex interaction on the aeroelastic displacement of an airfoil attached to the idealized two degree-of-freedom spring system is investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.168-173
• /
• 1998

A numerical analysis of axisymetric backward facing step main nozzle flow in air jet loom has been accomplished. To obtain basic design data for an optimum main nozzle for an air-jet loom and to predict the transonic/supersonic flow, a characteristic based upwind flux difference splitting compressible Navier-Stokes method has been used. The wall static pressure of the main nozzle and the flow velocity changes in the nozzle tube were analyzed by changing air tank pressures and acceleration tube lengths. The flow inside the nozzle experiences double choking one at the needle tip and the other at the acceleration tube exit at tank pressures over $4kg_f/cm^2$. The tank pressure $P_t$ leading to the critical condition depends on the acceleration tube length; i.e, $P_t$ is higher for longer acceleration tubes. The $P_t$ value required to bring the acceleration tube exit to the critical condition is nearly constant regardless of acceleration tube length. The round needle tip shape might lead to less total pressure loss when compared with step shape.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.3
• /
• pp.1-13
• /
• 2001

RANS solver has been developed to solve the flows past a ship with a drift angle. The solver employs a finite volume method for the spatial discretization and Euler implicit method for the time integration. Turbulent flows are simulated by Spalart-Allmaras one-equation model. Developed solver is applied to analyze the hydrodynamic forces and flows of two tankers with a same forebody but different afterbodies. The computed flows and hydrodynamic forces are compared with the measured flows and captive model test data. The computed results show good agreements with experimental data and show clearly the effects of stern hull form on the hydrodynamic forces and the flows.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• Journal of Korea Water Resources Association
• /
• v.44 no.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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.6
• /
• pp.495-500
• /
• 2017

In the boundary layer of supersonic or hypersonic vehicles, there is the conversion from kinetic energy to thermal energy, called aerodynamic heating. Aerodynamic heating has to be considered to design supersonic vehicles, because it induces severe heat flux to surface. Transient heat transfer analysis with CFD is used to predict thermal response of vehicles, however transient heat transfer analysis needs excessive computing powers. Loosely coupled method is widely used for evaluating thermal response, however it needs to be revised for overestimated heat flux. In this research, quasi-transient method, which is combined loosely coupled method and conjugate heat transfer analysis, is proposed for evaluating thermal response with efficiency and reliability. Defining reference time of splitting flight scenario for transient simulation is important on accuracy of quasi-transient method, however there is no algorithm to determine. Therefore the research suggests the algorithm with various flow conditions to define reference time. Supersonic flow field of blunt body with constant acceleration is calculated to evaluate quasi-transient method. Temperature difference between transient and quasi-transient method is about 11.4%, and calculation time reduces 28 times for using quasi-transient method.