• Korean Journal of Agricultural and Forest Meteorology
• /
• v.10 no.3
• /
• pp.82-93
• /
• 2008

In order to ensure a standardized data analysis of the eddy covariance measurements, Hong and Kim's quality control program has been updated and used to process eddy covariance data measured at two levels on the main flux tower at Gwangneung site from January to May in 2005. The updated program was allowed to remove outliers automatically for $CO_2$ and latent heat fluxes. The flag system consists of four quality groups(G, D, B and M). During the study period, the missing data were about 25% of the total records. About 60% of the good quality data were obtained after the quality control. The number of record in G group was larger at 40m than at 20m. It is due that the level of 20m was within the roughness sublayer where the presence of the canopy influences directly on the character of the turbulence. About 60% of the bad data were due to low wind speed. Energy balance closure at this site was about 40% during the study period. Large imbalance is attributed partly to the combined effects of the neglected heat storage terms, inaccuracy of ground heat flux and advection due to local wind system near the surface. The analysis of wind direction indicates that the frequent occurrence of positive momentum flux was closely associated with mountain valley wind system at this site. The negative $CO_2$ flux at night was examined in terms of averaging time. The results show that when averaging time is larger than 10min, the magnitude of calculated $CO_2$ fluxes increases rapidly, suggesting that the 30min $CO_2$ flux is influenced severely by the mesoscale motion or nonstationarity. A proper choice of averaging time needs to be considered to get accurate turbulent fluxes during nighttime.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.3
• /
• pp.180-201
• /
• 2020

It has been widely known that the effect of diffracted waves at the tip of composite breakwater with finite length causes the change of standing wave height along the length of breakwater, the spatial change of wave pressure on caisson, and the occurrence of meandering damage on the different sliding distance in sequence. It is hard to deal with the spatial change of wave force on trunk of breakwater through the two-dimensional experiment and/or numerical analysis. In this study, two and three-dimensional numerical techniques with olaFlow model are used to approach the spatial change of wave force including the impulsive breaking wave pressure applied to trunk of breakwater, the effect of rear region, and the occurrence of diffracted waves at the tip of caisson located on the high crested rubble mound. In addition, it is thoroughly studied the mean wave height, mean horizontal velocity, and mean turbulent kinetic energy through the numerical analysis. In conclusion, it is confirmed that the larger wave pressure occurs at the front wall of caisson around the still water level than the original design conditions when it generates the shock-crushing wave pressure checked by not two-dimensional analysis, but three-dimensional analysis through the change of wave pressure applied to the caisson along the length of breakwater.

• Journal of Korea Water Resources Association
• /
• v.48 no.11
• /
• pp.891-903
• /
• 2015

This paper presents numerical simulations of flow around a circular patch of vegetation using a depth-averaged two-dimensional numerical model which is capable of simulating flow structure in vegetated open channel. In order to account for vegetation effect, drag force terms are included in governing equations. Numerical simulations are conducted with various solid volume fractions (SVF). Flow passes through a circular patch and low velocity region, which is called wake region, is formed downstream of the patch. When SVF is larger than 0.08, a recirculation is observed. The location of the recirculation is moved further downstream as SVF decreases. Von-$K{\acute{a}}rm{\acute{a}}n$ vortex street is developed beyond the wake region due to interaction between two shear layers induced by a circular patch of vegetation. The vortex is developed as SVF is larger than 0.08, and the location of the vortex is consistent with the maximum of turbulence kinetic energy. The location of the peak of turbulence kinetic energy is moved further downstream as SVF decreases.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.8
• /
• pp.537-544
• /
• 2017

The purpose of this study is to analyze the flow characteristics when a staggered hollow fiber membrane module is modeled as a porous medium. The pressure-velocity equation was used for modeling the porous medium, using pressure drop data. In terms of flow characteristics, we compared the case of the "porous medium" when the membrane module was modeled as a porous medium with the case of the "membrane module" when considering the original shape of the membrane module. The difference in pressure drop between the "porous medium" and "membrane module" was less than 0.6%. However, the maximum flow velocity and mean turbulent kinetic energy of the "porous medium" were 2.5 and 95 times larger than those of the "membrane module," respectively. Our results indicate that modeling the hollow fiber module as a porous medium is useful for predicting pressure drop, but not sufficient for predicting the maximum flow velocity and mean turbulent kinetic energy.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.6
• /
• pp.531-552
• /
• 2020

In the previous study, both the wave characteristics at the tip of composite breakwater and on caisson were investigated by applying olaFlow numerical model of three-dimensional regular waves. In this paper, the same numerical model and layout/shape of composite breakwater as applied the previous study under the action of one directional irregular waves were used to analyze two and three-dimensional spatial change of wave force including the impulsive breaking wave pressure applied to trunk of breakwater, the effect of rear region, and the occurrence of diffracted waves at the tip of caisson located on the high crested rubble mound. In addition, the frequency spectrum, mean significant wave height, mean horizontal velocity, and mean turbulent kinetic energy through the numerical analysis were studied. In conclusion, the larger wave pressure occurs at the front wall of caisson around the still water level than the original design conditions when it generates the shock-crushing wave pressure in three-dimensional analysis condition. Which was not occurred by two-dimensional analysis. Furthermore, it was confirmed that the wave pressure distribution at the caisson changes along the length of breakwater when the same significant incident wave was applied to the caisson. Although there is difference in magnitude, but its variation shows the similar tendency with the case of previous study.

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

• The Journal of the Acoustical Society of Korea
• /
• v.34 no.1
• /
• pp.36-45
• /
• 2015

Recently, a lot of studies have been made about the methods used to generate turbulent velocity fields stochastically in order to effectively predict broadband flow noise. Among them, the FRPM (Fast Random Particle Mesh) method which generates turbulence with specific statistical properties using turbulence kinetic energy and dissipation obtained from the steady solution of the RANS (Reynolds Averaged Navier-Stokes) equations has been successfully applied. However, the FRPM method cannot be applied to the flow noise problems involving intrinsic unsteady characteristics such as centrifugal fan. In this paper, to effectively predict the broadband noise generated by centrifugal fan, U-FRPM (unsteady FRPM) method is developed by extending the FRPM method to be combined with the unsteady numerical solutions of the unsteady RANS equations to generate the turbulence considered as broadband noise sources. Firstly, an unsteady flow field is obtained from the unsteady RANS equations through CFD (Computational Fluid Dynamics). Then, noise sources are generated using the U-FRPM method combined with acoustic analogy. Finally, the linear propagation model which is realized through BEM (Boundary Element Method) is combined with the generated sources to predict broadband noise at the listeners' position. The proposed technique is validated to compare its prediction result with the measured data.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.7
• /
• pp.930-938
• /
• 2018

Underwater noise radiated from submarines is directly related to the probability of being detected by the sonar of an enemy vessel. Therefore, minimizing the noise of a submarine is essential for improving survival outcomes. For modern submarines, as the speed and size of a submarine increase and noise reduction technology is developed, interest in flow noise around the hull has been increasing. In this study, a noise analysis technique was developed to predict flow noise generated around a submarine shape considering the free surface effect. When a submarine is operated near a free surface, turbulence-induced noise due to the turbulence of the flow and bubble noise from breaking waves arise. First, to analyze the flow around a submarine, VOF-based incompressible two-phase flow analysis was performed to derive flow field data and the shape of the free surface around the submarine. Turbulence-induced noise was analyzed by applying permeable FW-H, which is an acoustic analogy technique. Bubble noise was derived through a noise model for breaking waves based on the turbulent kinetic energy distribution results obtained from the CFD results. The analysis method developed was verified by comparison with experimental results for a submarine model measured in a Large Cavitation Tunnel (LCT).

• Journal of Power System Engineering
• /
• v.5 no.2
• /
• pp.22-29
• /
• 2001

This paper represents axial mean velocity, turbulent kinetic energy and swirl number based on momentum flux measured in the X-Y plane and Y-Z plane respectively of a cone type gas swirl burner by using X-probe from the hot-wire anemometer system. This experiment is carried out at flow rates 350 and $450{\ell}/min$ respectively, which are equivalent to the combustion air flow rate necessary for heat release 15,000 kcal/hr in gas furnace, in the test section of a subsonic wind tunnel. Axial mean velocities and turbulent kinetic energies show that their maximum values exist centering around narrow slits situated radially on the edge of and in the forefront of a burner until $X/R{\fallingdotseq}1.5$, but they have a peculiar shape like a starfish diffusing and developing into inward and outward of a burner by means of the mixing between flows ejected from narrow slits, an inclination baffle plate and swirl vanes respectively according to downstream regions. Moreover, they show a relatively large value in the inner region of 0.5$S_m$ obtained by integration of velocity profiles shows a characteristic that has an inflection point composing of the maximum and minimum value until X/R<3, but shows close agreement with the geometric swirl number after a distance of X/R=3.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.7
• /
• pp.1675-1683
• /
• 1995

A numerical study of the fluid flow and heat transfer characteristics of the two-dimensional impingement jet with a confinement plate has been carried out. The fluid flow was calculated by solving the full Navier-Stokes equation. In doing that, the well known SIMPLER algorithm was used and the trouble making convection term was discretized according to QUICKER scheme. The energy equation was simply solved by using the SOR method. For the Reynolds number of 10000, two channel heights, say 1.5 and 3.0 times the jet exit width, and two thermal boundary conditions constant wall temperature and constant wall heat flux were considered. Discrete heat sources were flush mounted along the impingement plate at a distance of 0, 2, 3, 4, 5, 6, 10, 12, times the jet exit width from the stagnation point. The length of each heat source is 4 times the jet exit width long. The Nusselt number averaged over each heat source was compared with experiment. Comparison shows that both calculations and experiment have the secondary peak of Nusselt number at downstream of stagnation point, even though there is a little quantitative difference in between. The difference is believed due to abscure thermal boundary condition in experiment and also accuracy of turbulence model used. The secondary peak is shown to be caused by rigorous turbulent flow motion generated as the wall jet flow is retarded and developes into the channel flow without flow reversal.