• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.3
• /
• pp.1015-1027
• /
• 1996

In this paper, a numerical procedure for the calculation of the incompressible Navier-Stokes equations in a generalized nonorthogonal body fitted coordinate system is proposed and is validated through three test problems. Present numerical procedure derives the pressure equation by using the pressure substitution method on the regular grid system, and discretized momentum equations are based on the covariant velocity components. Cavity flow, backward facing step flow, and two dimensional channel flow with a sinusoidal wavy wall are chosen as three test problems. Numerical solutions obtained by present procedure shows a good agreement with previous numerical and/or experimental results. Convergence rate is also satisfactory.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.5
• /
• pp.1613-1623
• /
• 1996

In this paper, two discretization methods, hybrid and QUICK, are tested for the Navier-Stokes equations written in general nonorthogonal body fitted coordinates. Comparison is made by calculating two laminar flows at low Reynolds numbers of 10 - 100. One is a two-dimensional channel of gradually expanding cross section and the other is an axisymmetric flow through a circular tube having a circular constriction. Results show that the QUICK scheme results in a numerical solution more accurate than that of hybrid. The QUICK scheme also shows faster convergence for both test cases. As the number of grid points increases, all numerical solutions converge with more oscillation. The number of grid points in the y-direction(cross stream direction) is also shown to play a significant role in the approximation of convection term within separated flow zone.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.1 no.2
• /
• pp.22-31
• /
• 1997

We developed a general purpose program for the analysis of flows in the combustor with recirculating flow regime and simulated the flows. The program uses non-staggered grids based on finite volume method and the primitive variables are cartesian velocities. The combustion model is irreversible one step reaction with infinite chemistry The Favre averaged governing equations are considered and the clipped gaussian distribution is considered as a probability density function of the conserved scalar. We calculated turbulent diffusion flame with recirculating flow regime. Simulation shows two recirculating regions like experimental results. Velocity, turbulent kinetic energy, temperature and concentration distribution in simulation agree well with experimental data.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.2
• /
• pp.668-676
• /
• 1991

Turbulent flows around tube banks and in the diffuser were studied using a non-orthogonal boundary fitted coordinate system and the modified K-.epsilon. turbulence model. In these cases, many problems emerge which stem from the geometrical complexity of the flow domain and the physical complexity of turbulent flow itself. To treat the complex geometry, governing equations were reformulated in a non-orthogonal coordinate system with Cartesian velocity components and discretised by the finite volume method with a non-staggered variable arrangement. The modified K-.epsilon. model of Hanjalic and Launer was applied to solve above two cases under the condition of strong and mild pressure gradient. The results using the modified K-.epsilon. model results in both test cases.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.73-83
• /
• 1997

이 논문은 최근에 개발된 네가지 종류의 고차 유계 대류항 처리법인 SOUCUP, HLPA, SMARTER 그리고 COPLA 해법들을 비교 분석한다. 모든 해법들을 실제적인 공학적인 문제에의 적용을 위하여 비균일, 비직교 좌표에 공식화 하였다. 해법들의 상대적인 검증을 위하여 여러 가지 종류의 시험 문제들에 적용하여 검증하였다. 수치실험의 결과는 시험한 4종류의 해법들이 유계성을 만족하고 고차 해법의 정확도를 유지하는 것으로 나타났다. HLPA, SMARTER와 COPLA 해법들은 거의 같은 정도의 해의 정확성을 보였으며, SCOCUP 해법은 조금 부정확한 것으로 나타났다.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.21-28
• /
• 1997

This paper presents a two-dimensional horizontal implicit model to general circulation in estuaries and coastal seas. The model is developed in non-orthogonal curvilinear coordinates system, using the Interpolating Matrix Method (IMM), in combination with a technique of multi-domain. In the propose model, the Saint-Venant equations are solved by a splitting-up technique, in the successive steps; convection, diffusion and wave propagation. The ability of the proposed model to deal with full scale nature is illustrated by the interpretation of a dye-tracing experiment in the Gironde estuary.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1149-1164
• /
• 1997

In this study, the applicability of the RNG k-.epsilon. model to the analysis of the complex flows is studied. The governing equations based on a non-orthogonal coordinate formulation with Cartesian velocity components are used and discretized by the finite volume method with non-staggered variable arrangements. The predicted results using the RNG k-.epsilon. model of three complex flows, i.e., the flow over a backward-facing step and a blunt flat plate, the flow around a 2D model car are compared to these from the standard k-.epsilon. model and experimental data. That of the unsteady axisymmetric turbulent flow within a cylinder of reciprocating model engine including port/valve assembly and the spray characteristics within a chamber of direct injection model engine are compared to these from the standard k-.epsilon. model and experimental data. The results of reattachment length, separated eddy size, average surface pressure distribution using the RNG k-.epsilon. model show more reasonable trends comparing with the experimental data than those using the modified k-.epsilon. model. Although the predicted rms velocity using the modified k-.epsilon. model is lower considerably than the experimental data in incylinder flow with poppet valve, predicted axial and radial velocity distributions at the valve exit and in-cylinder region show good agreements with the experimental data. The spray tip penetration predicted using the RNG k-.epsilon. model is more close to the experimental data than that using the modified k-.epsilon. model. The application of the RNG k-.epsilon. model seems to have some potential for the simulations of the unsteady turbulent flow within a port/valve-cylinder assembly and the spray characteristics over the modified k-.epsilon. model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.1
• /
• pp.139-149
• /
• 1994

A numerical simulation of unsteady axisymmetric turbulent flow was performed for a reciprocating engine including port/valve assembly. The governing equations based on a nonorthogonal coordinate formulation with Cartesian velocity components were used and discretised by the finite volume method with non-staggered variable arrangements. The modified $\kappa-\xi$. turbulence model which included the effect of compressibility was used. The results of twodimensional transient calculation for the axisymmetric configuration were compared with the experimental data. Although slightly low rms velocity was predicted compared to the experimental data, predicted velocity distributions at the valve exit and in-cylinder region showed good agreements with the experimental data. The flow at the valve exit was separated at the same valve lift position with the experimental data. Two vortices incylinder region were generated during the initial intake process. The clockwise main vortex became strong and moved upward to the top wall. The counter-clockwise second vortex became weak and stick to the upper left corner of the cylinder. After middle intake process, new vortex adjacent to upper cylinder wall appeared by the piston motion and therefore, the in-cylinder flow was formed into three vortices. The cylinder pressure just before bottom dead center of piston was higher than inlet pressure and then the reverse flow occured at the valve exit. The in-cylinder flow characteristics were strongly dependent on piston motion, but insensitive to valve motion.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1971-1978
• /
• 1992

A turbulent flow in a conical diffuser with total divergence angle of 8.deg. was numerically studied. The low Reynolds number k-.epsilon. model(Launder-Sharma model) was adopted to simulate the turbulence. The continuity and time averaged Navier-Stokes equations in a nonorthogonal coordinate system were solved by a finite volume method based on the fully elliptic formulation. The low Reynolds number k-.epsilon. model reasonably simulates the pressure recovery and the mean velocity components. However, there are also considerable discrepancies between predicted and measured shear stress distribution on the wall and turbulent kinetic energy distributions. It is necessary to investigate the flow structure at the entry of the diffuser, numerically as well as experimentally.