• Journal of the Korean Society for Marine Environment & Energy
• /
• v.15 no.3
• /
• pp.186-197
• /
• 2012

In order to understand hydrodynamic and water quality changes on the Saemangeum reservoir in accordance to inner development plan, intensive numerical simulations using EFDC have been done. Due to inner dike construction and proposed dredging plans, stratification might occur and yield flow field change. It should be noticed that very conditional gate operation schedule adjusting target water elevation of -1.5 meter causes severe stratification and hence plays an important role in poor water qualities. By using random walk particle tracking residence simulations, it is found that hydrodynamic characteristics depends greatly on riverine inflow conditions. It is also inferred that the northern part of the Mangyeong reservoir behaves as a dead zone and acts as major reasoning of water quality deterioration owing to benthic flux from long-term residing settled sediment.

• The Journal of Engineering Geology
• /
• v.10 no.1
• /
• pp.13-25
• /
• 2000

A three-dimensional discrete fracture network model based on probabilistic characteristics of fracture geometry and transmissivity was designed to calculate the conductivity tensor and to estimate theanisotropy of conductivity. The conductivities, $K_p$, obtained from the numerical simulation of single-holepacker test corresponded well to those from the field tests. From this, it can be concluded that thefracture network model designed in this study can represent hydraulic characteristics of in-situ fractured rock mass. Block-scale conductivities, $K_b$, estimated from the modelling of steady-state flow through the REV-scale block were ranged between the arithmetic mean and harmonic mean of theconductivity estimates from packer tests. The conductivity along north-south direction was 1.4 timesgreater than that along the east-west direction. It was concluded that the anisotropy of conductivitywas insignificant. It was also found that there was a little correlation between $K_b$ and $K_p$. This would be to that the conductivities from the packer test simulation was strongly dependent on thetransmissivity and the number of fractures within the packer test intervals.

• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.70-80
• /
• 2015

Construction waste is known to include a large part of coarse and fine aggregates, which can be recirculated in the industry. Separating those aggregates economically from the waste has been thus considered to be one of the most important issues in this field. In particular, paste mixed in the waste causes significant complain from the inhabitants living near the place where waste-processing equipments are built and operated. In this study, we investigate the operational principle of a newly developed paste separator by using theoretical (in this first part) and CFD (in the second part) analysis. The separator consists of a rotor which turned out to play a significant role in separating those pastes from the aggregates. Under suitable assumptions regarding the air flow velocity as well as the particle velocity, we show that particles can be stagnant at the outlet of the roto channel for a wide range of parameter values, which allow the particles to get enough time to settle down via the gravitation. We also demonstrate such phenomenon by using a simple numerical simulation.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.4
• /
• pp.353-362
• /
• 2009

In the process of turbine modernizations, the investigation of the influences of water passage roughness on radial flow machine performance is crucial and validates the efficiency step up between reduced scale model and prototype. This study presents the specific losses per component of a Francis turbine, which are estimated by CFD simulation. Simulations are performed for different water passage surface roughness heights, which represents the equivalent sand grain roughness height. As a result, the boundary layer logarithmic velocity profile still exists for rough walls, but moves closer to the wall. Consequently, the wall friction depends not only on roughness height but also on its shape and distribution. The specific losses are determined by CFD numerical simulations for each component of the prototype, taking into account its own specific sand grain roughness height. The model efficiency step up between reduced scale model and prototype value is finally computed by the assessment of specific losses on prototype and by evaluating specific losses for a reduced scale model with smooth walls. Furthermore, surveys of rough walls of each component were performed during the geometry recovery on the prototype and comparisons are made with experimental data from the EPFL Laboratory for Hydraulic Machines reduced scale model measurements. This study underlines that if rough walls are considered, the CFD approach estimates well the local friction loss coefficient. It is clear that by considering sand grain roughness heights in CFD simulations, its forms a significant part of the global performance estimation. The availability of the efficiency field measurements provides an unique opportunity to assess the CFD method in view of a systematic approach for turbine modernization step up evaluation. Moreover, this paper states that CFD is a very promising tool for future evaluation of turbine performance transposition from the scale model to the prototype.

• Tunnel and Underground Space
• /
• v.23 no.6
• /
• pp.470-479
• /
• 2013

Hydraulic fracturing is used as a method for promoting the fluid flow in the rock and, in the energy field such as geothermal development and the development of sales gas, many studies has been actively conducted. In many cases, hydraulic fracturing is not performed in isotropic rock and especially in the case of sedimentary rocks, hydraulic fracturing is conducted in the transverse isotropic rock. The direction of the crack growth on hydraulic fracturing does not necessarily coincides with the direction of maximum principal stress in the transverse isotropic rock. Therefore, in this study, bonded particle model with hydro-mechanical coupling analysis was adopted for analyzing the characteristics of hydraulic fracturing in transverse isotropic rock. In addition, experiments of hydraulic fracturing were conducted in laboratory-scale to verify the validity of numerical analysis. In this study, the crack growth and crack patterns showed significant differences depending on the viscosity of injection fluid, the angle of bedding plane and the influence of anisotropy. In the case of transverse isotropic model, the shear crack growth due to hydraulic fracturing appeared prominently.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.6
• /
• pp.383-396
• /
• 2016

In this study, wave model was conducted on the presence or absence of water level changes and currents effects in coastal waters coexisting with waves and currents, then the results were compared. The flow field applied the results of the RIAMOM model and the wave model applied the SWAN model. Among ECMWF, NCEP and JMA, wind data applied JMA data sets which agreed well with the observed data comparatively. Numerical simulation was conducted for 8 months from January to August 2016. For each case, the deviation of wave height was calculated for the high wave of more than 2.5 m for comparison with observed data. As a result, the deviation of wave height was not significant both considering water level changes and currents effects or not at wave observation stations installed in deep waters. However, a significant deviation of wave height of 5~10% was obtained depending on water level changes and currents effects at the comparison point in shallow waters.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.112-122
• /
• 1993

A potential-based panel method is presented for the analysis of a partially or supercavitating two-dimensional hydrofoil at a finite submergence beneath a free surface, treating without approximation the effects of the finite Froude number and the hydrostatic pressure. Free surface sources and normal dipoles are distributed on the foil and cavity surfaces, their strength being determined by satisfying the kinematic and dynamic boundary conditions on the foil-cavity boundary. The cavity surface is determined iteratively as a part of the solution. Numerical results show that the wave profile is altered significantly due to the presence of the cavity. The buoyancy effect due to the hydrostatic pressure, which has usually been neglected in most of the cavitating flow analysis, is found playing an important role, especially for the supercavitating hydrofoil; the gravity field increases the cavity size in shallow submergence, but decreases it when deeply submerged, while the lift reduces at all submergence depth.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.2 s.42
• /
• pp.71-80
• /
• 2005

The solution to a liquid sloshing problem is challenge to the field of engineering. This is not only because the dynamic boundary condition at the free surface is nonlinear, but also because the position of the free surface varies with time in a manner not known a priori. Therefore, this nonlinear phenomenon, which is characterized by the oscillation of the unrestrained free surface of the fluid, is a difficult mathematical problem to solve numerically and analytically. In this study, three-dimensional boundary element method(BEM), which is based on the so-called an arbitrary Lagrangian-Eulerian(ALE) approach for the fluid flow problems with a free surface, was formulated to solve the behavior of the nonlinear free surface motion. An ALE-BEM has the advantage to track the free surface along any prescribed paths by using only one displacement variable, even for a three-dimensional problem. Also, some numerical examples were presented to demonstrate the validity and the applicability of the developed procedure.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.9 no.3
• /
• pp.85-91
• /
• 2005

Radiative heating of a liquid rocket base plane due to plume emission is numerically investigated. Calculation of flow and temperature fields around rocket nozzle precedes and thereby realistic plume shape and temperature distribution inside the plume are obtained. Based on the calculated temperature field, radiative transfer equation is solved by discrete ordinate method. With the sample rocket plume, the averaged radiative heat flux reaching the base plane is calculated about 5 kw/m$^{2}$ at the flight altitude of 10.9 km. This value is small compared with radiative heat flux caused by constant-temperature (1500 K) plume emission, but it is not negligibly small. At higher. altitude (29.8km), view factor between the base plane and the exhaust plume is increased due to the increased expansion angle of the plume. Nevertheless, the radiative heating disappears since the base plane is heated to high temperature (above 1000 K due to convective heat transfer.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.133-138
• /
• 2002

A 3-D axisymmetric computer program is developed to predict the NO behavior in SNCR system for the stoker incinerator with the waste treatment capacity, 200ton/day. To this end a turbulent reacting flow field calculation is made using proper assumption and empiricism. The stoker bed is assumed to be a homogeneous waste-volatilized gaseous state. The initial composition or reactants are assumed based on the data of the ultimate analysis. Turbulent is resolved by k-e model and turbulent reaction is handled by eddy-breakup model harmonized with empirical chemistry data for gaseous combustion, NO and urea reaction. The liquid droplet is traced by Lagrangian method incorporated by aerodynamic drag, Coriolis and crntrifugal forces. Radiation is treated by sensible heat loss model. Calculation results are in good agreement with experimental data at the outlet of post combustion chamber in Daejon 4th industrial complex. The flue gas shows the temperature range of $900\sim1000^{\circ}C$, velocity of 5m/s and NO concentration of 140ppm at the exit while the measured temperature, flue gas velocity and NO concentration are $967^{\circ}C$, $3\sim4m/s$ and $100\sim200ppm$respectively. Using the developed computer program a parametric study has been made with the variation of heat content of waste, castable length and SNCR variables for the determination of proper injector location. In general, the calculated results are consistent and physically acceptable.