• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.11
• /
• pp.1169-1176
• /
• 2011

Liquid film flows are classified into waveless laminar, wavy laminar, and turbulent flows depending on the Reynolds number or the flow stability. Since the wavy motions of the film flows are so intricate and nonlinear, studies on them have largely been experimental. Most numerical approaches have been limited to the waveless flow regime. The various free surface-tracking schemes adopted for this problem were used to more accurately estimate the average film thickness, rather than to capture the unsteady wavy motion. In this study, the wavy motions in laminar wavy liquid film flows with Reynolds numbers of 200-1000 were simulated with various numerical schemes based on the volume of fluid (VOF) method for interface tracking. The results from each numerical scheme were compared with the experimental results in terms of the average film thickness, the wave velocity, and the wave amplitude.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.1
• /
• pp.92-100
• /
• 2018

Vortical systems are considered a main feature to sustain turbulence in a boundary layer through interaction. Such turbulent structures result in frictional drag and erosion or vibration in engineering applications. Research for controlling turbulent flow has been actively carried out, but in order to show the effect of vortices in a turbulent boundary layer, it is necessary to clarify the mechanism by which turbulent energy is transferred. For this purpose, it is convenient to demonstrate and capture phenomena in a laminar boundary layer. Therefore, in this study, the interactions of disturbed flow around a hemisphere on a flat plate in laminar flow were analyzed. In other words, a street of hairpin vortices was generated following a wake region formed after flow separation occurred over a hemisphere. Necklace vortices surrounding the hemisphere also appeared due to a strong adverse pressure gradient that brought high momentum fluid into the wake region thereby leading to an increase in the frequency of hairpin vortices. To mitigate the effect of these necklace vortices, local suction control was applied through a hole in front of the hemisphere. Flow visualization was recorded to qualitatively determine flow modifications, and hot-film measurements quantitatively supported conclusions on how much the power of the hairpin vortices was reduced by local wall suction.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.414-426
• /
• 1997

Numerical computations were conducted to characterize the three-dimensional laminar flow through an injector orifice having an inclined angle of 30 .deg.. For this study, the incompressible Navier-Stokes equations in generalized curvilinear coordinates, using a pseudocompressibility approach for continuity equation, were solved. The computations were performed using the finite difference implicit, approximately factored scheme of Beam and Warming and multi-block grids of complete continuity at block interfaces. The multi-block computations were validated for the steady state using direct comparison of multi-block solutions with equivalent single-block ones, including 2-D 180.deg. TAD and 3-D 90.deg. pipe bend. The comparisons between the numerical solutions and the flow field measurements for a tube with sudden contraction were presented in this work for solution validation. Computational results showed the nature of complex flow fields within the inclined injector orifice, including strong pressure-driven secondary flows in the cross stream induced by the effect of streamline curvature. In addition, asymmetric secondary flows were induced in the Reynolds number range above assumed laminar flow regime considered. However, turbulence calculations and grid dependency studies are needed for more accurate computations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4003-4012
• /
• 1996

Hydro dynamically developing and simultaneously (hydro dynamically and thermally) developing laminar flows of a Bingham plastic in a circular pipe have been investigated numerically. Solutions have been obtained by using a four-step fractional method combined with an equal order bilinear finite element method. For the hydro dynamically developing flow, shorter entrance length is required to reach fully developed velocity field for larger yield stress and non-monotonic pressure drop along the pipe centerline is observed when the yield stress exceeds a certain critical value. For the simultaneously developing flow, the heat transfer characteristics show the same trends as those predicted for the thermally developing flow (Graetz problem).

• Journal of the Korean Society of Mechanical Technology
• /
• v.13 no.4
• /
• pp.1-7
• /
• 2011

The paper concerns numerical study of fully developed laminar flow of a Newtonian water and non-Newtonian fluids, 0.2% aqueous of sodium carboxymethyl cellulose(CMC) solution in eccentric annuli with combined bulk axial flow and inner cylinder rotation. Pressure losses and skin friction coefficients have been measured when the inner cylinder rotates at the speed of 0~200 rpm. A numerical analysis considered mainly the effects of annular eccentricity and inner cylinder rotation. The present analysis has demonstrated the importance of the drill pipe rotation and eccentricity. In eccentricity of 0.7 of a Newtonian water, the flow field is recirculation dominated and unexpected behavior is observed. it generates a strong rotation directed layer, that two opposing effects act to create two local peaks of the axial velocity. The influences of rotation, radius ratio and working fluid on the annular flow field are investigated.

• Journal of Advanced Marine Engineering and Technology
• /
• v.15 no.2
• /
• pp.64-72
• /
• 1991

The unsteady, two-dimensional laminar flow and heat transfer of an incompressible, constant-property fluid flowing around a circular cylinder are numerically analyzed. The Navier-Strokes equation and the energy equation are solved by the finite difference method. The range of the Reynolds number is 10 to 100 and the Prandtl number considered is 0.7. The contours of the flow pattern, equi-vorticity line and isotherm pattern around a circular cylinder are shown. Also, numerical solutions of the surface vorticity, pressure coefficient, drag coefficient, local Nusselt number and mean Nusselt number are obtained. The numerical results for the final time fo calculation are compared with the other available experimental and theoretical results for the steady state and are found to be in good agreement with them.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.3-10
• /
• 1995

The compressible laminar and turbulent viscous flow on a slender body in supersonic speed as well as subsonic speed has been numerically simulated at high angle of attack. The steady and time-accurate compressible thin-layer Navier-Stokes code based on an implicit upwind-biased LU-SGS algorithm has been developed and specifically applied at angles of attack of 20, 30, 40 deg, respectively. The modified eddy-viscosity turbulence model suggested by Degani and Schiff was used to simulate the case of turbulent flow. Any geometric asymmetry and numerical perturbation have not been intentionally or artificially imposed in the process of computation. The purely numerical results for laminar and turbulent cases, however, show clear asymmetric formation of vortices which were observed experimentally. Contrary to the subsonic results, the supersonic case shows the symmetric formation of vortices as indicated by the earlier experiments.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.1-8
• /
• 2010

Laminar separation bubble and transitional flow over the NACA0012 are investigated at a moderate range of Reynolds numbers. A Reynolds-Averaged Navier-Stokes code is coupled with an empirical transition model that can predict transition onset points and the length of transition region. Without solving the boundary layer equations, approximated e-N method is directly applied to the RANS code and iteratively solved together. The computational results are compared with the experimental data for the NACA0012 airfoil. Results of transition onset point and the length are compared well with experimental data and Xfoil prediction. The present RANS results show at high angles of attack better agreement with experimental data than Xfoil results using the boundary layer equations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.9
• /
• pp.1124-1132
• /
• 2004

The inverse problem of determining unknown inlet temperature in thermally developing, hydrodynamically developed two phase laminar flow in a parallel plate duct is considered. The inlet temperature profile is determined by measuring temperature in the flow field. No prior information is needed for the functional form of the inlet temperature profile. The inverse convection problem is solved by minimizing the objective function with regularization method. The conjugate gradient method as iterative method and the Tikhonov regularization method are employed. The effects of the functional form of inlet temperature, the number of measurement points and the measurement errors are investigated. The accuracy and efficiency of these two methods are compared and discussed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1552-1561
• /
• 1997

In this numerical study, it is investigated for the ice-formation phenomena for water flow in a concentric tube. The freezing layers of ice in both the inner and outer wall of a concentric tube are simultaneously considered. In the solution strategy, the complete set of governing equations in both the solid and liquid regions are resolved. Numerical results are obtained by varying the inner/outer wall temperatures and Reynolds number. The results show that the inner/outer wall temperatures have the great effect on the thickness of the solidification layer thereof. The shapes of ice layer in both the inner and outer wall can be expressed as a function of inverse Graetz number. As the wall temperature in inner or outer tube decreases, the heat transfer coefficients in both inner and outer ice layer surfaces increase absolutely.