• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.9
• /
• pp.1073-1084
• /
• 1999

The effects of secondary flow on the structure of a turbulent wake generated by a flat plate was investigated experimentally. The secondary flow was induced In a $90^{\circ}$ curved duct in which the flat plate wake generator was installed. The wake generator was installed in such a way that the wake velocity gradient exists in the span wise direction of the curved duct. Measurements were made in the plane containing the mean radius of curvature where pressure gradient and curvature effects were small compared with the secondary flow effect. All six components of the Reynolds stresses were measured in the curved duct. Turbulence intensities in the curved wake are higher than those in the straight wake due to an increase of the turbulent kinetic energy production by the secondary flow. In the inner wake region, shear stress and strain in the plane containing the velocity gradient of the wake show opposite signs with respect to each other, so that eddy viscosity Is negative in this region. This indicates that gradient-diffusion type turbulence models are not appropriate to simulate this type of flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1340-1346
• /
• 1999

The turbulent shear flow around a surface-mounted vertical fence was investigated using the two-frame PTV system. The Reynolds number based on the fence height(H) was 2950. From this study, it is revealed that at least 400 instantaneous velocity field data are required for ensemble average to get reliable turbulence statistics, but only 100 field data are sufficient for the time-averaged mean velocity information. Various turbulence statistics such as turbulent intensities, turbulence kinetic energy and Reynolds shear stress were calculated from 700 instantaneous velocity vector fields. The fence flow has an unsteady recirculation region behind the fence, followed by a slow relaxation to the flat-plate boundary layer flow. The time-averaged reattachment length estimated from the streamline distribution is about 11.2H. There exists a region of negative Reynolds shear stress near the fence top due to the highly convex (stabilizing) streamline-curvature of the upstream flow. The large eddy structure in the separated shear layer seems to have significant influence on the development of the separated shear layer and the reattachment process.

• Nuclear Engineering and Technology
• /
• v.23 no.3
• /
• pp.275-284
• /
• 1991

It is important to predict the main feature of fully developed turbulent secondary flow through infinite triangular arrays of parallel rod bundles. One-equation turbulence model which include anisotropic eddy viscosity model was applied to predict the exact velocity field. For a constant properties, Reynolds equations were solved by the finite element method. Mean axial velocity near the wall is simulated by the law of the wall. The numerical results showed good agreement with avaiable experimental data. The strength of the secondary flow increased with Reynolds number but decreased with rod spacing, P/D (pitch-to-diameter). The secondary flow affects remarkably the distribution of the axial velocity, wall shear stress and turbulent kinetic energy in the closely packed rod array bundles.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.13 no.6
• /
• pp.1034-1042
• /
• 2010

This study investigates effects of three different parameterizations of vertical mixing scheme on upper ocean simulation of the East Sea, focusing on the seasonal variations of the sea surface temperature(SST) and the mixed layer depth(MLD) using an ocean general circulation model(GFDL MOM1.1). The considered vertical mixing schemes are the Laplacian scheme(L scheme) that use a constant eddy coefficient, the Mellor-Yamada scheme(MY scheme), and a new scheme(Noh scheme). The Noh scheme, a second-order turbulence closure, was developed considering recent observational evidences such as the enhancement of turbulent kinetic energy near the sea surface. During summer L scheme underestimates the SST, while MY scheme overestimates the SST, compared to climatological SST. Noh scheme produces the SST in better agreement with climatological one. During winter all schemes overestimate the SST up to $4^{\circ}C$ compared to climatological SST. Vertical profiles of the basin-mean temperature show that L scheme produces higher temperature below the thermocline than those of other schemes. The winter MLD simulated from L scheme is rather large compared to that from other schemes, but the differences in MLD during summer are not significant.

• Nuclear Engineering and Technology
• /
• v.51 no.2
• /
• pp.573-578
• /
• 2019

Anisotropic distribution of the turbulent kinetic energy and the near-field excitations are the main causes of the steady state Flow-Induced Vibration (FIV) which could lead to fretting wear damage in vertically arranged supported slender rods. In this article, a combined Computational Fluid Dynamics (CFD) and Computational Structural Mechanic (CSM) approach named two-way Fluid-Structure Interaction (FSI) is used to investigate the modal characteristics of a typical rod's vibration. Performance of an Unsteady Reynolds-Average Navier-Stokes (URANS) and Large Eddy Simulation (LES) turbulence models on asymmetric fluctuations of the flow field are investigated. Using the LES turbulence model, any large deformation damps into a weak oscillation which remains in the system. However, it is challenging to use LES in two-way FSI problems from fluid domain discretization point of view which is investigated in this article as the innovation. It is concluded that the near-wall meshes whiten the viscous sub-layer is of great importance to estimate the Root Mean Square (RMS) of FIV amplitude correctly as a significant fretting wear parameter otherwise it merely computes the frequency of FIV.

• Wind and Structures
• /
• v.34 no.1
• /
• pp.43-58
• /
• 2022

The BARC flow is studied via Direct Numerical Simulation at a relatively low turbulent Reynolds number, with focus on the geometrical representation of the leading-edge (LE) corners. The study contributes to further our understanding of the discrepancies between existing numerical and experimental BARC data. In a first part, rounded LE corners with small curvature radii are considered. Results show that a small amount of rounding does not lead to abrupt changes of the mean fields, but that the effects increase with the curvature radius. The shear layer separates from the rounded LE at a lower angle, which reduces the size of the main recirculating region over the cylinder side. In contrast, the longitudinal size of the recirculating region behind the trailing edge (TE) increases, as the TE shear layer is accelerated. The effect of the curvature radii on the turbulent kinetic energy and on its production, dissipation and transport are addressed. The present results should be contrasted with the recent work of Rocchio et al. (2020), who found via implicit Large-Eddy Simulations at larger Reynolds numbers that even a small curvature radius leads to significant changes of the mean flow. In a second part, the LE corners are fully sharp and the exact analytical solution of the Stokes problem in the neighbourhood of the corners is used to locally restore the solution accuracy degraded by the singularity. Changes in the mean flow reveal that the analytical correction leads to streamlines that better follow the corners. The flow separates from the LE with a lower angle, resulting in a slightly smaller recirculating region. The corner-correction approach is valuable in general, and is expected to help developing high-quality numerical simulations at the high Reynolds numbers typical of the experiments with reasonable meshing requirements.

• Wind and Structures
• /
• v.11 no.4
• /
• pp.275-289
• /
• 2008

Large eddy simulations have been performed within and over different types of urban building arrays. This paper adopted three dimensionless parameters, building frontal area density (${\lambda}_f$) the variation degree of building height (${\sigma}_h$), and the staggered degree of building range ($r_s$), to study the systematic influence of building spacing, height and layout on wind and turbulent characteristics. The following results have been achieved: (1) As ${\lambda}_f$ decrease from 0.25 to 0.18, the mean flow patterns transfer from "skimming" flow to "wake interference" flow, and as ${\lambda}_f$ decrease from 0.06 to 0.04, the mean flow patterns transfer from "wake interference" flow to "isolated roughness" flow. With increasing ${\lambda}_f$, wind velocity within arrays increases, and the vortexes in front of low buildings would break, even disappear, whereas the vortexes in front of tall buildings would strengthen and expand. Tall buildings have greater disturbance on wind than low buildings do. (2) All the wind velocity profiles and the upstream profile converge at the height of 2.5H approximately. The decay of wind velocity within the building canopy was in positive correlation with ${\lambda}_f$ and $r_s$. If the height of building arrays is variable, Macdonald's wind velocity model should be modified through introducing ${\sigma}_h$, because wind velocity decreases at the upper layers of the canopy and increases at the lower layers of the canopy. (3) The maximum of turbulence kinetic energy (TKE) always locates at 1.2 times as high as the buildings. TKE within the canopy decreases with increasing ${\lambda}_f$ and $r_s$ but the maximum of TKE are very close though ${\sigma}_h$ varies. (4) Wind velocity profile follows the logarithmic law approximately above the building canopy. The Zero-plane displacement $z_d$ heighten with increasing ${\lambda}_f$, whereas the maximum of and Roughness length $z_0$ occurs when ${\lambda}_f$ is about 0.14. $z_d$ and $z_0$ heighten linearly with ${\sigma}_h$ and $r_s$, If ${\sigma}_h$ is large enough, $z_d$ may become higher than the average height of buildings.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.25 no.3
• /
• pp.55-66
• /
• 2020

Near-inertial waves (NIW) which are primarily generated by wind can contribute to vertical mixing in the ocean. The energetic NIW can be generated by typhoon due to its strong wind and preferable wind direction changes especially on the right-hand side of the typhoon. Here we investigate the generation and distribution of NIW using the output of a real-time ocean forecasting system. Five-year model outputs during 2013-2017 are analyzed with a focus on cases of energetic NIW generation by the passage of three typhoons (Halong, Goni, and Chaba) over the East Sea. Calculations of wind energy input (${\bar{W}}_I$), and horizontal kinetic energy in the mixed layer (${\bar{HKE}}_{MLD}$) reveal that the spatial distribution of ${\bar{HKE}}_{MLD}$, which is strengthened at the right-hand side of typhoon tracks, is closely related with ${\bar{W}}_I$. Horizontal kinetic energy in the deep layer (${\bar{HKE}}_{DEEP}$) shows patch-shaped distribution mainly located at the southern side of the East Sea. Spatial distribution of ${\bar{HKE}}_{DEEP}$ shows a close relationship with negative relative vorticity regions caused by warm eddies in the upper layer. Monthly-mean ${\bar{HKE}}_{MLD}$ and ${\bar{HKE}}_{DEEP}$ during a typhoon passing over the East Sea shows about 2.5-5.7 times and 1.2-1.6 times larger values than those during summer with no typhoons, respectively. In addition, their magnitudes are respectively about 0.4-1.0 and 0.8-1.0 times from those during winter, suggesting that the typhoon-induced NIW can provide a significant energy to enhance vertical mixing at both the mixed and deep layers during summer.

• Ocean Science Journal
• /
• v.42 no.2
• /
• pp.89-102
• /
• 2007

The surface circulation of northern South China Sea (hereafter SCS) for the period 1987-2005 was studied using the data of more than 500 satellite-tracked drifters and wind data from QuikSCAT. The mean flow directions in the northern SCS except the Luzon Strait (here after LS) during the periods October_March was southwestward, and $April{\sim}September$ northeastward. A strong northwestward intrusion of the Kuroshio through the LS appears during the $October{\sim}March$ period of northeasterly wind, but the intrusion became weak between April and September. When the strong intrusion occurred, the eddy kinetic energy (EKE) in the LS was $388cm^2/s^2$ which was almost 2 times higher than that during the weak-intrusion season. The volume transport of the Kuroshio in the east of the Philippines shows an inverse relationship to that of the LS. There is a six-month phase shift between the two seasonal phenomena. The volume transport in the east of the Philippines shows its peak sis-month earlier faster than that of the LS. The strong Kuroshio intrusion is found to be also related to the seasonal variation of the wind stress curl generated by the north easterly wind. The negative wind stress curl in the northern part of LS induces an anticyclonic flow, while the positive wind stress curl in the southern part of LS induces a cyclonic flow. The northwestward Kuroshio intrusion in the northern part of LS happened with larger negative wind stress curl, while the westward intrusion along $20.5^{\circ}N$ in the center of the LS occurred with weaker negative wind stress curl.

• Journal of Korea Water Resources Association
• /
• v.33 no.5
• /
• pp.581-592
• /
• 2000

Turbulence structure and suspended sediment transport capacity in vegetated open-channel flows are investigated numerically in the present paper. The $\textsc{k}-\;\varepsilon$ model is employed for the turbulence closure. Mean velocity and turbulence characteristics including turbulence intensity, Reynolds stress, and production and dissipation of turbulence kinetic energy are evaluated and compared with measurement data available in the literature. The numerical results show that mean velocity is diminished due to the drag provided by vegetation, which results in the reduction of turbulence intensity and Reynolds stress. For submerged vegetation, the shear at the top of vegetation dominates turbulence production, and the turbulence production within vegetation is characterized by wakes. For emergent condition, it is observed that the turbulence generation is dominated by wakes within vegetation. In general, simulated profiles compares favorably to measured data. Computed values of eddy viscosity are used to solve the conservation equation for suspended sediment, yielding sediment concentration more uniform over the depth compared with the one in the plain channel. The simulation reveals that the suspended load decreases as the vegetation density increases and the suspended load increases as the particle diameter decreases for the same vegetation density.