• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1541-1546
• /
• 1993

To extend the widely used Gibson and Launder's second order closure model to the low-Reynolds-number region near a wall, modifications have been made for velocity pressure-gradient interaction and dissipation terms in the stress equations, and also for the dissipation rate equation. From the computation of fully developed plane channel flow, it is found that the results with present model agree well with the data of direct numerical simulation in the predictions of stress components. And, the computed mean velocity profile coincides with the universal velocity law.

• Journal of Mechanical Science and Technology
• /
• v.16 no.9
• /
• pp.1121-1136
• /
• 2002

The flow characteristics of film coolant issuing into turbulent boundary layer developing on a convex surface have been investigated by means of flow visualization and three-dimensional velocity measurement. The Schlieren optical system with a spark light source was adopted to visualize the jet trajectory injected at 35° and 90° inclination angles. A five-hole directional pressure probe was used to measure three-dimensional mean velocity components at the injection angle of 35°. Flow visualization shows that at the 90° injection, the jet flow is greatly changed near the jet exit due to strong interaction with the crossflow. On the other hand, the balance between radial pressure gradient and centrifugal force plays an important role to govern the jet flow at the 35° injection. The velocity measurement shows that at a velocity ratio of 0.5, the curvature stabilizes downstream flow, which results in weakening of the bound vortex structure. However, the injectant flow is separated from the convex wall gradually, and the bound vortex maintains its structure far downstream at a velocity ratio of 1.98 with two pairs of counter rotating vortices.

• 한국연소학회:학술대회논문집
• /
• 2002.11a
• /
• pp.49-56
• /
• 2002

A two-dimensional direct numerical simulation is performed to investigate the formation characteristics of a single vortex interacting with $CH_4/N_2$-Air counterflow nonpremixed flame. The numerical method was based on a predictor-corrector scheme for a low Mach number flow. The detailed transport properties and a 16-step augmented reduced mechanism are adopted in this calculation. The budgets of the vorticity transport equation arc examined to reveal the mechanisms leading to the formation, evolution and dissipation of a single vortex interacting with counterflow nonpremixed flame. It is found that the stretching term, which depends on the azimuthal component of vorticity, and radial velocity, mainly generates vortieitv in non-reacting and reacting flows. The viscous and baroclinic torque term destroy the vorticity in non-reacting flow. In addition, the baroclinic torque term due to density and pressure gradient generates vorticity, while viscous and the volumetric expansion terms due to density gradient destroy vorticity in reacting flow.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.4
• /
• pp.243-252
• /
• 2021

The present study aims to investigate the interaction of a wire-type turbulence stimulator and the laminar boundary layer on a flat plate by flow field measurement. For the towing tank tests, a one-dimensional Laser Doppler Velocimetry (LDV) attached on a two-axis traverse was used to measure the streamwise velocity component of the boundary layer flow in zero pressure gradient, disturbed by a turbulence stimulator. The wire diameter was 0.5 and 1.0 mm according to the recommended procedures and guidelines suggested by the International Towing Tank Conference. Turbulence development by the stimulator was identified by the skin friction coefficient, mean and Root Mean Square (RMS) of the streamwise velocity. The laminar boundary layer with the absence of the wire-type stimulator was similar to the Blasius solution and previous experimental results. By the stimulator, the mean and RMS of the streamwise velocity were increased near the wall, showing typical features of the fully developed turbulent boundary layer. The critical Reynolds number was reduced from 2.7×10⁵ to 1.0×10⁵ by the disturbances caused by the wire. As the wire diameter and the roughness Reynolds number (Rek) increased, the disturbances by the stimulator increased RMS of the streamwise velocity than turbulent boundary layer.

• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

• Wind and Structures
• /
• v.29 no.1
• /
• pp.65-75
• /
• 2019

A heated square cylinder (with height $A^*$) is kept parallel to the cold wall at a fixed gap height $0.5A^*$ from the wall. Another adiabatic rectangular cylinder (of same height $A^*$ and width $0.5A^*$) is placed upstream in an inline tandem arrangement. The spacing between the two cylinders is fixed at $3.0A^*$. The inlet flow is taken as Couette-Poiseuille flow based non-linear velocity profile. The conventional fluid (also known as base fluid) is chosen as water (W) whereas the nanoparticle material is selected as $Al_2O_3$. Numerical simulations are performed by using SIMPLE algorithm based Finite Volume approach with staggered grid arrangement. The dependencies of hydrodynamic and heat transfer characteristics of the cylinder on non-dimensional parameters governing the nanofluids and the fluid flow are explored here. A critical discussion is made on the mechanism of improvement/reduction (due to the presence of the upstream cylinder) of heat transfer and drag coefficient, in comparison to those of an isolated cylinder. It is observed that the heat transfer increases with the increase in the non-linearity in the incident velocity profile at the inlet. For the present range studied, particle concentration has a negligible effect on heat transfer.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.419-428
• /
• 2014

This paper presents a numerical simulation method for the flow around advancing ships in regular waves by using a rectilinear grid system. Because the grid lines do not consist with body surface in the rectilinear grid system, the body geometries are defined by the interaction points of those grid lines and the body surface. For the satisfaction of body boundary conditions, no-slip and divergence free conditions are imposed on the body surface and body boundary cells, respectively. Meanwhile, free surface is defined with the modified marker density method. The pressure on the free surface is determined to make the pressure gradient terms of the governing equations continuous, and the velocity around the free surface is calculated with the pressure on the free surface. To validate the present numerical method, a vortex induced vibration (VIV) phenomenon and flows around an advancing Wigley III ship model in various regular waves are simulated, and the results are compared with existing and corresponding research data. Also, to check the applicability to practical ship model, flows around KRISO Container Ship (KCS) model advancing in calm water are numerically simulated. On the simulations, the trim and the sinkage are set free to compare the running attitude with some other experimental data. Moreover, flows around the KCS model in regular waves are also simulated.