• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.42-49
• /
• 2008

A wavelet method is presented in order to improve the computational efficiency of two dimensional unsteady flow problems while maintaining the order of accuracy of conventional CFD schemes. First, by using the interpolating wavelet transformation including decomposition and thresholding, an adaptive dataset to a solution is constructed. Then, inviscid and viscous fluxes are calculated only at the points within an adaptive dataset, which enhances the computational efficiency. Second, thresholding step is modified to maintain the spatial and temporal accuracy of conventional CFD schemes automatically by selecting the threshold value between user-defined value and the magnitude of spatial or temporal truncation error. The wavelet method suggested in this study is successfully applied to various unsteady flow problems and it is shown that the computational efficiency is enhanced with maintaining the computational accuracy of CFD schemes.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.590-598
• /
• 1988

Viscous flows through a screen normal to an uniform flow are numerically simulated. A .kappa.-.epsilon. model is adopted for evaluation of the Reynolds stresses. The existence of a screen is regarded as extra sources in the momentum equations. The amount of extra sources is related to the resistance coefficient and the refraction coefficient of the screen. Flows are numerically simulated for various resistance coefficients and heights of the screen and Reynolds numbers. The present method has been verified to reasonably simulate viscous wakes and shear layers of the screen, for which the inviscid theory is quite limitted. As the fluids approach the screen, the velocity is reduced and the pressure is raised to satisfy the Bernoulli equation. After passing the screen, the velocity shows its minimum value at the down-stream, but static pressure is slowly recovered. A detached separation-bubble from the screen appears as the resistance coefficient is increased to a certain level. Such results are qualitatively in agreement with limitted experimental data available. The turbulent kinetic energy shows its maximum value at further down stream and decrease thereafter.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.42-49
• /
• 2008

A wavelet method is presented in order to improve the computational efficiency of two dimensional unsteady flow problems while maintaining the order of accuracy of conventional CFD schemes. First, by using the interpolating wavelet transformation including decomposition and thresholding, an adaptive dataset to a solution is constructed. Then, inviscid and viscous fluxes are calculated only at the points within an adaptive dataset, which enhances the computational efficiency. Second, thresholding step is modified to maintain the spatial and temporal accuracy of conventional CFD schemes automatically by selecting the threshold value between user-defined value and the magnitude of spatial or temporal truncation error. The wavelet method suggested in this study is successfully applied to various unsteady flow problems and it is shown that the computational efficiency is enhanced with maintaining the computational accuracy of CFD schemes.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.77-84
• /
• 1999

Numerical design optimization techniques are implemented for the improvement of the ram accelerator performance. The design object is to find the minimum ram tube length required to accelerate projectile from initial velocity $V_0$ to target velocity $V_e$. The premixture is composed of $H_2,\;O_2,\;N_2$ and the mole numbers of these species are selected as design variables. The objective function and the constraints are linearized during the optimization process and gradient-based Simplex method and SLP(Sequential Linear Programming) have been employed. With the assumption of two dimensional inviscid flow for internal flow field, the analyses of the nonequilibrium chemical reactions for 8 steps 7 species lave been performed. To determined the tube length, ram tube internal flow field is assumed to be in a quasi-steady state and the flow velocity is divided into several subregions with equal interval. Hence the thrust coefficients and accelerations for corresponding subregions are obtained and integrated for the whole velocity region. With the proposed design optimization techniques, the total ram tube length had been reduced $19\%$ within 7 design iterations. This optimization procedure can be directly applied to the multi-stage, multi-premixture ram accelerator design optimization problems.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.167-173
• /
• 1999

An Euler-based CFD tool has been developed for the performance evaluation of a thrust reverser mounted on a high bypass ratio turbofan engine. The computational domain surrounded by the ground and non-reflection boundary includes the whole nacelle configuration with a deployed thrust reverser. The numerical algorithm is based on the modified Godunovs scheme to allow the second order accuracy in both space and time. The grid system is generated by using eleven multi-blocks, of which the total cell number is 148,400. The thrust reverser is modeled as if it locates at the nacelle simply in all circumferential direction. The existence of a fan and an OGV(Outlet Guide Vane) is simulated by adopting the actuator disk concept, in which predetermined radial distributions of stagnation pressure ratio and adiabatic efficiency coefficient are used for the rotor type disk, and stagnation pressure losses and flow outlet angles for the stator type disk. All boundary conditions including the fan and OGV simulation are treated by Riemann solver. The developed solver is applied to a turbofan engine with a bypass ratio of about 5.7 and the diameter of the fan cowl of 83 inch. The computational results show that the Euler-based inviscid method is very useful and economical to evaluate the performance of thrust reversers.

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

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• 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 the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.94-103
• /
• 1993

This study deals with an investigation on the tip vortex generated by an elliptic wing with section shapes of NACA 0020. The flow structure on the wing surface is investigated by using tufts test as well as observing the cavitation pattern. The surface pressure on a foil surface is measured to complement the visualized flow field. Results show that a strong spanwise pressure gradient is a definite contributor on the formation of tip vorex, and the fluids from both sides contribute to the evolutionary process of tip vortex. On the other hand, a series of experiments are conducted to investigate the detailed structure of tip-vortex at various angles of attack. The tip-vortex formation and development are observed by producing a cavitation, and then by a laser sheet technique in conduction with a dye injection method. The shape of tip-vortex and the distance between a vortex core and the trailing vortex sheet are found to vary with the angle of attack. Overall features of tip flow are evaluated to complement the vortex model based on inviscid theory.