• Journal of Wetlands Research
• /
• v.13 no.3
• /
• pp.395-409
• /
• 2011

The objective of this study is to develop an accurate and robust two-dimensional finite element method for turbulence simulation in open channels. The model is based on Streamline Upwind/Petrov-Galerkin finite element method and Boussinesq's eddy viscosity theory. The method developed in the study is depth-averaged mixing length model which assumes anisotropic and local equilibrium state of turbulence. The model calibration and validation were performed by comparing with analytical solutions and observed data. Several numerical simulations were carried out, which examined the performance of the turbulence model for the purpose of sensitivity analysis. The uniform channels that appear horizontal flow and vertical flow were carried out. The model was also applied to the Han river was in for the applicability test. The results were compared with the observed data. The suggested model displayed reasonable flow distribution compare to the observed data in natural river flow. As a result of this study, the two-dimensional finite element model provides a reliable results for flow distribution based on the turbulence simulation in open channels.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.321-324
• /
• 2011

A computational analysis of nozzle flow characteristics and plume structure using Reynolds-averaged Navier-Stokes equations with $k-{\omega}$ SST turbulence model was conducted to examine performance of the supersonic nozzle employed in a small liquid-rocket engine for ground firing test. Computed results and experimental outcome of 2-D converging-diverging nozzle flow were compared for verifying the computational capability as well as the turbulence model validity. Numerical computations of 2-D axisymmetric nozzle flow was carried out with the selected model. As a result, flow separation with backflow appeared around the nozzle exit. This investigation was reported as a background data for the optimal nozzle design of small liquid-propellant rocket engine for ground test.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.1
• /
• pp.9-16
• /
• 2012

This paper presents the numerical results of a simulation of the free surface flow around a blunt bow ship model and focuses on the validation of the proposed method with a brief investigation of the relation between the resistance and free surface behavior. A finite volume method based on the Reynolds Averaged Navier-Stokes (RANS) approach is used to solve the governing flow equations, where the free surface, including wave breaking,is captured by using a two-phase Level-Set (LS) method. For turbulence closure, a two equation k-${\varepsilon}$ model with the standard wall function technique is used. Finally, the numerical results are compared with the available experimental data, showing good agreement.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.8-14
• /
• 2017

In this paper, separated flow characteristics is studied using the RANS(Reynold-averaged Navier-Stokes) modeling. The analysis is performed for the NASA's hump configuration which is the combination of a flat plate and a hump. This configuration was used in NASA's flow control workshop and it was one of validation cases for RANS and LES simulations. The separation occurs at the 65% of model length where a slot is positioned for the flow control. No flow control case and steady suction case are studied using RANS modeling. The Spalart-Allmaras model and the SST(Shear Stress Transport) model are applied and their accuracy are compared. To correlate CFD analysis with experimental data, the optimal boundary condition was investigated and the effect of a cavity around the slot is studied for the no flow case.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2001

Wave-induced currents drive nearshore transport processes, and hence an accurate understanding of wave-current interaction is required for proper management of coastal zone. This paper presents details of an adaptive quadtree grid based numerical model of the coupled wave climate and depth-averaged current field. The model accounts for wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, mixing processes, bottom friction effects, and movement of land-water interface at the shoreline. The wave period- and depth-averaged governing equations arc discrctized explicitly by means of an Adarns¬Bashforth second-order finite difference technique on adaptive hierarchical staggered quadtree grids. Results from the numerical model are in reasonable agreement with the laboratory data of longshore current generated by oblique waves on a plane beach (Visser 1980, 1991).

• Atmosphere
• /
• v.22 no.4
• /
• pp.429-436
• /
• 2012

The performances of the Seoul National University Urban Canopy Model (SNUUCM) under different meteorological conditions (clear, cloudy, and rainy conditions) in summertime are compared using observation dataset obtained at an urban site. The daily-averaged net radiation, sensible heat flux, and storage heat flux are largest in clear days and smallest in rainy days, but the daily-averaged latent heat flux is similar among clear, cloudy, and rainy days. That is, the ratio of latent heat flux to net radiation increases in order of clear, cloudy, and rainy conditions. In general, the performance of the SNUUCM is better in clear days than in cloudy or rainy days. However, the performance in simulating sensible heat flux in clear days is as poor as that in rainy days. For all the meteorological conditions, the performance in simulating latent heat flux is worst among the performances in simulating net radiation, sensible heat flux, and latent heat flux. The normalized mean error for latent heat flux is largest in rainy days in which the relative importance of latent heat flux in the surface energy balance becomes greatest among the three conditions. This study suggests that improvements to the parameterization of processes that are related to latent heat flux are particularly needed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.2
• /
• pp.245-252
• /
• 2002

Numerical analysis was conducted to characterize particle deposition on a horizontal rotating disk with thermophorectic effect under laminar flow field. The particle transport mechanisms considered were convection, Brownian diffusion, gravitational settling and thermophoresis. The averaged particle deposition velocities and their radial distributions for the upper surface of the disk were calculated from the particle concentration equation in a Eulerian frame of reference for rotating speeds of 0∼1000rpm and temperature differences of 0∼5K. It was observed from the numerical results that the rotation effect of disk increased the averaged deposition velocities, and enhanced the uniformity of local deposition velocities on the upper surface compared with those of the disk at rest. It was also shown that the heating of the disk with ΔT=5K decreased deposition velocity over a fairly broad range of particle sizes. Finally, an approximate deposition velocity model for the rotating disk was suggested. The comparison of the present numerical results with the results of the approximate model and the available experimental results showed relatively good agreement between them.

• Journal of the Korea Concrete Institute
• /
• v.21 no.3
• /
• pp.377-386
• /
• 2009

A reference diffusion coefficient model from ACI life-365 is drawn from test results by NT build 443. This test method gives a time-averaged diffusion coefficient during immersion period, thus the ACI model uses the time-averaged diffusion coefficient as a reference value. ACI model needs to be revised, considering the difference between the time-average value and reference value at specified time. In this study, firstly the analytic solutions of diffusion equation are derived considering the initiation time and period of exposure to chloride, and secondly the time-averaged diffusion coefficient from NT build 443 is converted into the diffusion coefficient at reference time. From this study, the reference diffusion coefficient of ACI model should be modified to be about 10% larger values than those of present ACI model. For convenient design of service life, previous relationship between the chloride diffusion coefficient from NT build 443 and that from NT build 492 is also modified. To compare the chloride diffusion coefficients of ACI and JCI models, the reference chloride diffusion coefficient with respect to the JCI model is drawn in the similar form of ACI model's, and service life prediction by ACI life-365 method is confirmed to give a conservative result.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2002.10a
• /
• pp.123-128
• /
• 2002

In this paper, an averaged constitutive model of concrete and reinforcing bars for RC structure and path-dependent Ohsaki's model for soil are applied, and an elasto-plastic interface model having thickness is preposed for seismic analysis of underground RC subway station structure. A finite element analysis technique is developed by applying aforementioned constitutive equations and verified through seismic analysis of underground RC subway station. Then, failure mechanisms of the RC subway station structure under seismic action are numerically derived. Then, failure modes and damage levels of the station are also analytically evaluated for the cases of several designs of the underground RC station.