• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.2
• /
• pp.224-232
• /
• 2000

Four turbulent $k-{\varepsilon}$ models (i.e., standard model, modified models with streamline curvature modification and/or preferential dissipation modification) are applied in order to analyze the turbulent flow of wall-attaching offset jet. For numerical convergence, this paper develops a method of slowly increasing the convective effect induced by skew-velocity in skew-upwind scheme (hereafter called Partial Skewupwind Scheme). Even though the method was simple, it was efficient in view of convergent speed, computer memory storage, programming, etc. The numerical results of all models show good prediction in first order calculations (i.e., reattachment length, mean velocity, pressure), while they show some deviations in ·second order (i.e., kinetic energy and its dissipation rate). Like the previous results obtained by upwind scheme, the streamline curvature modification results in better prediction, while the preferential dissipation modification does not.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.1
• /
• pp.10-17
• /
• 2008

In this paper an accurate and stable gridless method that can be applied to multi-dimensional convection problems is developed on a flow directional local grid. A two dimensional pure convection problem is calculated and more accurate and stable solution is obtained compared with other schemes in grid method. The tested numerical schemes include 1st-order upwind scheme, 2nd-order Leith scheme, 3rd-order MUSCL, and QUICK scheme. It is seen that more accurate results are expected when the schemes combined with a MMT control limiter.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.771-779
• /
• 2004

Unsteady behavior of the early wake in the viscous flow field past an impulsively started semicircular cylinder is studied numerically. In this paper, we propose the hybrid diffusion scheme to simulate dynamic characteristics of wake such as a fishtail-like flapping and an alternate vortex-shedding more accurately. This diffusion scheme based on particle strength exchange is mixed with the stochastic nature of random walk method. Also, the viscous splitting algorithm which calculates convective and diffusion terms successively is applied in order to handle random walk method effectively. Consequently, the early behavior of wake due to the breakdown of symmetrical vortici balance is more practically simulated with the vortex particle method.

• Journal of Computational Design and Engineering
• /
• v.3 no.4
• /
• pp.330-336
• /
• 2016

The present work is concerned with heat and mass transfer of an electrically conducting second grade MHD fluid past a semi-infinite stretching sheet with convective surface heat flux. The analysis accounts for thermophoresis and thermal radiation. A similarity transformations is used to reduce the governing equations into a dimensionless form. The local similarity equations are derived and solved using Nachtsheim-Swigert shooting iteration technique together with Runge-Kutta sixth order integration scheme. Results for various flow characteristics are presented through graphs and tables delineating the effect of various parameters characterizing the flow. Our analysis explores that the rate of heat transfer enhances with increasing the values of the surface convection parameter. Also the fluid velocity and temperature in the boundary layer region rise significantly for increasing the values of thermal radiation parameter.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.4
• /
• pp.429-438
• /
• 2017

Cranke-Nicolson scheme in native OpenFOAM source libraries was not able to provide 2nd order temporal accuracy of velocity and pressure since the volume flux of convective nonlinear terms was 1st accurate in time. In the present study the simplest way of getting the volume flux with 2nd order accuracy was proposed by using old fluxes. A possible numerical instability originated from an explicit estimation of volume fluxes could be handled by introducing a weighting factor which was determined by observing the ratio of the finally corrected volume flux to the intermediate volume flux at the previous step. The new calculation of volume fluxes was able to provide temporally accurate velocity and pressure with 2nd order. The improvement of temporal accuracy was validated by performing numerical simulations of 2D Taylor-Green vortex of which an exact solution was known and 2D vortex shedding from a circular cylinder.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.61-62
• /
• 2003

In this paper, the new higher order wall boundary conditions are proposed for solving the incompressible flows. The square driven cavity flows are simulated by using the variable order method of lines with the present wall boundary conditions. The variable order method of lines is constructed by the spatial discretization, i.e., the variable order proper convective scheme for convective terms and the modified differential quadrature method for diffusive terms, and time integration. The 2nd, 4th, and 6th order solutions are presented and these results show this higher order boundary conditions are very promising for the incompressible flow simulations.

• Journal of Mechanical Science and Technology
• /
• v.15 no.8
• /
• pp.1188-1195
• /
• 2001

A numerical analysis has been carried out to investigate the influences of thermosolutal convection on the heat and mass transfer and solute segregation in crystals grown by the vertical Bridgman technique. The governing equations are solved by a finite-volume method using the power law scheme and the SIMPLE algorithm in which body-fitted coordinate system has been used. A primary convective cell driven by thermal gradients forms in the bulk of the domain, while a secondary convective cell driven by solutal gradients forms near interface. As the solutal Rayleigh number increases, secondary cell becomes to be stronger and has a great influence on the radial concentration along the interface.

• Journal of computational fluids engineering
• /
• v.3 no.2
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1158-1174
• /
• 1988

A numerical technique is developed for the solution of fully developed turbulent recirculating flow in the passage of variable area using the non-orthogonal coordinate transformation. In the numerical analysis, primitive pressure-velocity finite difference equations were solved by SIMPLER algorithm with 2-equation turbulence model and algebraic stress model (ASM). QUICK scheme on the differencing of convective terms which is free from the inaccuracies of numerical diffusion has been applied to the variable grids and the results compared with those from HYBRID scheme. In order to test the effect of streamline curvatures on turbulent diffusion Lee and Choi streamline curvature correction model which has been obtained by modifying the Leschziner and Rodi's model is testes. The ASM was also employed and the results are compared to those from another turbulence model. The results show that difference of convective differencing schemes and turbulence models give significant differences in the prediction of velocity fields in the expansion region and outlet region of the combustor, however show little differences in the parallel flow region.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1996.11a
• /
• pp.107-116
• /
• 1996

Numerical simulations have been conducted to characterize unsteady flow structures in an axisymmetric supersonic inlet diffuser with sinusoidal pressure oscillations at the diffuser exit. The formulation is based on the unsteady Navier-Stokes equations and turbulence closure is achieved using a two-layer model with a too-Reynolds-number scheme for the near-wall treatment. The governing equations are formulated in an integral form, and are discretized by the four-stage Runge-Kutta scheme for temporal terms and the Harten-Yee upwind TVD scheme for convective terms. Results indicated that the inlet shock characteristics are significantly modified by acoustic oscillations originating from the combustor. The characteristics of shock/boundarv-layer interactions (such as the size of separation bubble, terminal shock shape, and vorticity intensity) are also greatly iufluenced by the shock oscillation due to acoustic waves.