• 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.35 no.6
• /
• pp.561-567
• /
• 2011

The inverse radiation analysis of a two-dimensional irregular configuration using unstructured triangular meshes is presented. In this study, an enclosure filled with an absorbing, emitting and scattering medium with diffusely emitting and reflecting opaque boundaries is considered. The finite volume method is applied to solve the radiative transfer equation in order to simulate the measured incident radiation values which are used as input data for the inverse analysis. The conjugate gradient method is adopted for the estimation of wall emissivities by minimizing the objective function at each iteration step. To verify the performance of the unstructured grid system, we compare the results with those using a structured grid system for the two-dimensional lopsided shape. The effect of measurement errors on the estimation accuracy is also investigated.

• Journal of computational fluids engineering
• /
• v.9 no.3
• /
• pp.73-80
• /
• 2004

A LS(Level Set) formulation is developed for computing two-phase flows on non- orthogonal meshes. Compared with the VOF(Volume-of-Fluid) method based on a non-smooth volume-fraction function, the LS method can calculate an interfacial curvature more accurately by using a smooth distance function. Also, it is quite straightforward to implement for 3-D irregular meshes compared with the VOF method requiring much more complicated geometric calculations. The LS formulation is implemented into a general purpose program for 3-D flows and verified through several test problems.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.58-65
• /
• 1994

It is very important to understand characteristics of solution due to the variation of computational grid sizes, especially when turbulence model not incorporating wall-function is used. The present paper performs numerical investigation on the grid dependency of numerical solution for three dimensional turbulent flow field around a ship. In the present study a finite volume method with a modified sub-grid scale turbulence model and a numerically constructed non-orthogonal curvilinear coordinate system capable of conforming complex ship geometries are used. Numerical studies are then performed for a mathematical Wigley hull and the Series 60, $C_B=0.8$ hull forms. The results for various grid sizes are compared with each other and with measured data to show grid dependencies of numerical solutions.

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