• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.5
• /
• pp.101-108
• /
• 2002

A numerical simulation of an incompressible cavity flow is conducted. Two dimensional Navier-Stokes equations are integrated using staggered grid and a finite volume method with C-QUICK scheme for spatial derivatives and fully implicit scheme for the time derivatives. SIMPLE-C algorithm is employed to solve the pressure field. Computational results show that the third eddy is generated in the shear layer mode but not in the steady mode. This signifies that the third eddy plays an important role in cavity flow stability. As a means to control the flow, jet blowing is applied to a position below the cavity upstream edge. Effects of flow control parameters on the stability such as the frequency, the phase, and the velocity magnitude are reported.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.3
• /
• pp.375-385
• /
• 1996

Turbulent flow and heat transfer characteristics around staggered tube banks were studied using the 3-D Navier-Stokes equations and energy equation governing a steady incompressible flow, which were reformulated in a non-orthogonal coordinate system with cartesian velocity components and discretized by the finite volume method with a non-staggered variable arrangement. The predicted turbulent kinetic energy using RNG $k-{\varepsilon}$ model was lower than that of standard $k-{\varepsilon}$ model but showed same result for mean flow field quantities. The prediction of the skin friction coefficient using RNG $k-{\varepsilon}$ model showed better trend with experimental data than standard $k-{\varepsilon}$ model result. The inclined flow showed higher velocity and skin friction coefficient than transverse flow because of extra strain rate ($\frac{{\partial}w}{{\partial}y}$). Also, this was why the inclined flow showed higher local heat transfer coefficient than the transverse flow.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1023-1031
• /
• 1989

An experimental study is conducted to investigate the turbulent wake flow around tubes in vertical single row tube banks. All measurements are performed at Reynolds number(Re$_{max}$) 4.2*10$_{3}$ - 2.5*10$_{4}$ with varying tube spacings from the wide pitch ratio(H/D=3.07) to the very narrow one(H/D=1.23). Flow patterns are visualized using the smoke-wire method. Mean static pressures, velocity components, and various statistical quantities of turbulence are obtained by the computer on-line technique. In the case of wide tube spacings, the near wakes of tube show similar trends to those of a single tube, and their flow indicats an anisotropic turbulence. However, as the pitch ratio decreases, wide and narrow wakes appear alternately behind adjacent tubes due to the deflected flow. Also, in the case of H/D .leq. 1.54, Karman vortex is not formed at the side of relatively wide wake.e.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.2
• /
• pp.1-9
• /
• 1999

Applicability of the pressure gradient method which is formulated based on pressure gradient is verified against turbulent flow analysis. In the pressure gradient method, pressure gradient instead of pressure itself is obtained using continuity constraint. Since correct pressure gradient is found only when mass conservation is satisfied, pressure gradient method can reflect physics of flow field properly The pressure gradient method is formulated with semi-staggered grid system which locates each primitive variables on the same grid point but evaluates pressure gradient in-between. This grid system ensures easy programming and reflection of correct physics in analysis. For verifying applicability of this method, the pressure gradient method is applied to turbulent flow analysis with low Reynolds number $\kappa$-$\varepsilon$ model. Turbulent flows include fully developed channel flow, backward-facing step flow, and conical diffuser flow. Prediction results show that the pressure gradient method can be applied to turbulent flow analysis. However, the pressure gradient method requires somewhat long computation time. Proper way to find optimum under-relaxation factor, $\gamma$, is also need to be developed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.10
• /
• pp.961-969
• /
• 2013

In this study, the numerical analysis of a turbulent flow around both a staggered and an inline tube bundle was conducted using ANSYS CFX V.13, a commercial CFD software. The flow was assumed to be steady, incompressible, and isothermal. According to the CFD Best Practice Guideline, the sensitivity study for grid size, accuracy of the discretization scheme for convection term, and turbulence model was conducted, and its result was compared with the experimental data to estimate the applicability of the CFD Best Practice Guideline. It was concluded that the CFD Best Practice Guideline did not always guarantee an improvement in the prediction performance of the commercial CFD software in the field of tube bundle flow.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.1
• /
• pp.111-117
• /
• 2010

A three-dimensional numerical model is employed to investigate wave deformation due to a submerged structure. The three-dimensional numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES(large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS(sub-grid scale) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF(volume-of-fluid) method is used to track the distorted and broken free surface. A simple linear wave is generated on a constant depth and compared with analytical solutions. The model is then applied to study wave deformation due to a submerged structure and the predicted results are compared with available laboratory measurements.

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

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.