• Journal of Mechanical Science and Technology
• /
• v.20 no.2
• /
• pp.262-270
• /
• 2006

Direct numerical simulations (DNS) of turbulent channel flows up to $Re_{\tau}=1270$ are performed to investigate an elliptic feature and strain rate field on cross sections of coherent fine scale eddies (CFSEs) in wall turbulence. From DNS results, the CFSEs are educed and the strain rate field around the eddy is analyzed statistically. The principal strain rates (i.e. eigenvalues of the strain rate tensor) at the CFSE centers are scaled by the Kolmogorov length $\eta$ and velocity $U_k$. The most expected maximum (stretching) and minimum (compressing) eigenvalues at the CFSE centers are independent of the Reynolds number in each $y^+$ region (i. e. near-wall, logarithmic and wake regions). The elliptic feature of the CFSE is observed in the distribution of phase-averaged azimuthal velocity on a plane perpendicular to the rotating axis of the CFSE $(\omega_c)$. Except near the wall, phase-averaged maximum $(\gamma^{\ast}/\gamma_c^{\ast})$ and minimum $(\alpha^{\ast}/\alpha_c^{\ast})$ an eigenvalues show maxima on the major axis around the CFSE and minima on the minor axis near the CFSE center. This results in high energy dissipation rate around the CFSE.

• Ocean Systems Engineering
• /
• v.12 no.3
• /
• pp.335-358
• /
• 2022

Hydrodynamic field alteration around a cylindrical pier using a curved vane is numerically investigated. The curved vane with various angles ranged from 10 to 220 degree is placed at the upstream of the cylindrical pier. Laminar flow is adopted in order to perform the steady-state analysis. It is found that the flow separation leads to the formation of four bubbles depending on the value of the curved vane angle. Two bubbles are located in the region between the rear of the curved vane and the leading surface of the cylindrical pier, while the remaining two bubbles are located at the wake zone behind the cylindrical pier. Numerical analysis is performed to reveal the hydrodynamic field and influence of curved vane on the formation and evolution of the bubbles. It is found that the center and size of the bubble depend mainly on the value of the curved vane angle. It is observed that the flow velocity vector shows clearly the alteration in the flow velocity direction especially at the leading surface and rear surface of the curved vane owing to the occurrence of flow separation and flow dissipation along the circumference of the vane.

• Journal of Environmental Science International
• /
• v.17 no.6
• /
• pp.699-709
• /
• 2008

To understand the development mechanism of the aerosols in the surface boundary layer, the variation in the aerosol number concentration due to the divergence and convergence of the wind fields was investigated. The aerosol number concentration was measured in the size ranges of $0.3{\sim}10.0{\mu}m$ using a laser particle counter(LPC) from 0000 LST on 03 Feb. to 0600 LST on 07 Feb. 2004 at Mokpo in Korea during snowfall. The Velocity Azimuth Display(VAD) technique was used to retrieve the radar wind fields such as the horizontal wind field, divergence, and deformations including the vertical air velocity from a single Doppler radar. As a result, the distribution of the aerosol number concentration is apparently different for particles larger than $1{\mu}m$ during snowfall, and it has a tendency to increase at the beginning of the snowfall. The increase and decrease in the aerosol concentration due to the convergence and divergence of the wind fields corresponded to the particles with diameters greater than $1{\mu}m$. It is found that the fluctuations in the aerosol number concentration are well correlated with the development and dissipation of snowfall radar echoes due to the convergence and divergence of horizontal wind fields near the surface boundary layer in the inland during the snowfall.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.1
• /
• pp.1-12
• /
• 1990

The study is concerned with the analysis of combined forging of non-axisymmetric shapes with inclined protrusions by UBET technique. Work hardening is considered for the given range of strain rate during the forging process. A complex shape with inclined cavities is analyzed by subdividing the workpiece into finite UBET elements for which simple velocity fields are applicable. An experimental set-up was designed and manufactured for the experiment, and experiments are carried out with lead billets. The devised set-up can be used for closed-die forging of complex shapes with protrusions in which the dies can be separated automatically for easy removal of the forged products. Based on the derived kinematically admissible velocity fields for corresponding UBET elements, general computer programs have been developed. Since the energy dissipation rate for each elemental region is provided by subprograms (Subroutine or Function), the developed program can be applied to the forging problems of various shapes. The present study has shown that the method developed can be effectively applied to forging of non-axisymmetric shapes with complicated protrusions.

• Journal of the Korean Society for Precision Engineering
• /
• v.27 no.11
• /
• pp.46-56
• /
• 2010

The present research works investigated into the low velocity impact characteristics of DP 780 high strength steel sheet with 1.7 mm in thickness subjected to free boundary condition using three-dimensional finite element analysis. Finite element analysis was carried out via ABAQUS explicit code. Hyper-elastic model and the damping factor were introduced to improve an accuracy of the FE analysis. An appropriate FE model was obtained via the comparison of the results of the FE analyses and those of the impact tests. The influence of the impact energy and nose diameter of the impact head on the force-deflection curves, impact time, absorption characteristics of the impact energy, deformation behaviours, and stress-strain distributions was quantitatively examined using the results of FE analysis. The results of the FE analysis showed that the absorption rate of impact energy lies in the range of the 70.7-77.5 %. In addition, it was noted that the absorption rate of impact energy decreases when the impact energy increases and the nose diameter of the impact head decreases. The local deformation of the impacted region was rapidly increased when the impact energy was larger than 76.2 J and the nose diameter was 20 mm. A critical impact energy, which occur the instability of the DP780, was estimated using the relationship between the plastic strain and the impact energy. Finally, characteristics of the plastic energy dissipation and the strain energy density were discussed.

• Journal of Korea Water Resources Association
• /
• v.50 no.9
• /
• pp.637-646
• /
• 2017

This study presents a numerical model for simulating dense interflows. The governing equations are provided and the finite difference method is used with the $k-{\varepsilon}$ turbulence model. The model is used to simulate a dense interflow established in a deep ambient water, resulting velocity and excess density profiles. It is observed that velocity decreases in the longitudinal direction due to water entrainment in the vicinity of the outlet and rarely changes for increased Richardson number. Similarity collapses of velocity and excess density are obtained, but those of turbulent kinetic energy and dissipation rate are not. A shape factor for the dense interflow is obtained from the simulated profiles. The value of this shape factor can be used in the layer-averaged modeling of dense interflows. In addition, a buoyancy-related parameter ($c_{3{\varepsilon}}$) for the $k-{\varepsilon}$ model and the volume expansion coefficient (${\beta}_0$) are obtained from the simulated results.

• Journal of Korea Water Resources Association
• /
• v.44 no.6
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.2
• /
• pp.253-261
• /
• 1987

Thermoacoustic oscillation of an air column induced by heated wires is investigated analytically and experimentally. Acoustic power generation from a single heater wire is estimated based on the result of heat transfer analysis and expressed in terms of the efficiency factor indicating the conversion efficiency from heat to acoustic energy. It is shown that the efficiency factor becomes maximum when the wire radius is the order of the coustic boundary layer thickness and the flow velocity is close to the thermal diffusion velocity. Onset condition of the column oscillation is obtained by equating the acoustic power generation at the heater to the power loss due to thermoviscous dissipation at the tube wall and the convection and radiationloss at the open ends of the tube. In estimating the acoustic power generation, the heater is treated as a stretched single wire by correcting the flow velocity to take into account the interactions between adjacent heater wires. Experiment is performed by using a spiral heater of 1mm diameter in an air column of 37mm diameter. The heat input to drive the oscillation is measured and compared with the theoretical prediction. A good agreement is found between the theory and experiment, which is regarded as a substantial verification of the present analysis.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.10a
• /
• pp.407-414
• /
• 2000

The friction pendulum type seismic isolation system (FPS) has been developed to provide a simple and effective way to achieve earthquake resistance for buildings . The major advantages are: the isolation frequency can be easily achieved by designing a curvature of the surface and does not depend on the supported weight of a structure. The function of carrying vertical load is separated to the function of providing horizontal stiffness. Next the friction provides sufficient energy dissipation to protect the structure from earthquake response and resistance to the weak external disturbances such as wind load and ground vibrations due to traffic. In this paper, the friction coefficients are evaluated from number of experiments on the FPS test specimens. The relations between friction coefficient and the test waveform, velocity, and pressure are reviewed and further works are discussed.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.766-771
• /
• 2001

An experimental study was undertaken to investigate the effect of coriolis force for the turbulent flow at low Reynolds numbers in a rotating straight square duct. The study was carried out using a hot-wire anemometer. The flow Reynolds number based on the hydraulic diameter ranged from 4,000 to 18,000 and Rotation number ranged from 0 to 0.196. At Re=9000, developing turbulent flow was calculated for mean velocity and Reynolds stress. Pressure coefficient and energy dissipation spectrum were also calculated.