• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.6
• /
• pp.374-385
• /
• 2013

Sediment particle motions have been numerically simulated over a sinusoidal ripple. Turbulent boundary layer flows are generated by Large Eddy Simulation, and the sediment particle motions are simulated using Lagrangian particle tracking method. Two unsteady flow conditions are used in the experiment by employing two different wave amplitudes while keeping other conditions such as wave period same. As expected, the amount of suspended sediment particles is clearly dependent on the wave amplitude as it is increasing with increasing flow intensity. However, it is also observed that the pattern of suspension may be different as well due to the only different condition caused by wave amplitude. Specially, the time of maximum sediment suspension within the wave period is not coincident between the two cases because sediment suspension is strongly affected by the existence of turbulent eddies that are formed at different times over the ripple between the two cases as well. The role of these turbulent eddies on sediment suspension is important as it is also confirmed in previous researches. However, it is also found the time of these eddies' formation may also dependent on the wave amplitude over rippled beds. Therefore, it has been proved that various flow as well as geometric conditions under waves has to be considered in order to have better understanding on the sediment suspension process over ripples. In addition, it is found that high turbulent energy and strong upward flow velocities occur during the time of eddy formation, which also supports high suspension rate at these time steps. The results indicate that the relationship between the structure of flows and bedforms has to be carefully examined in studying sediment suspension at coastal regions.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.670-691
• /
• 2018

This study presents the research results of the BMT(Benchmark Model Test) simulations of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to predict fault reactivation and the coupled hydro-mechanical behavior of fault. BMT scenario simulations of Task B were conducted to improve each numerical model of participating group by demonstrating the feasibility of reproducing the fault behavior induced by water injection. The BMT simulations consist of seven different conditions depending on injection pressure, fault properties and the hydro-mechanical coupling relations. TOUGH-FLAC simulator was used to reproduce the coupled hydro-mechanical process of fault slip. A coupling module to update the changes in hydrological properties and geometric features of the numerical mesh in the present study. We made modifications to the numerical model developed in Task B Step 1 to consider the changes in compressibility, Permeability and geometric features with hydraulic aperture of fault due to mechanical deformation. The effects of the storativity and transmissivity of the fault on the hydro-mechanical behavior such as the pressure distribution, injection rate, displacement and stress of the fault were examined, and the results of the previous step 1 simulation were updated using the modified numerical model. The simulation results indicate that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing interaction and collaboration with other research teams of DECOVALEX-2019 Task B and validated using the field experiment data in a further study.

• Journal of Korean Society of Disaster and Security
• /
• v.14 no.4
• /
• pp.47-60
• /
• 2021

Recently, as the occurrence frequency of sudden floods due to climate variability increased, the damage of aging chuteway slabs of spillway are on the rise. Accordingly, a wide array of field survey, hydraulic experiment and numerical simulation have been conducted to find the cause of damage on chuteway slabs. However, these studies generally reviewed the flow characteristics and distribution of pressure on chuteway slabs. Therefore the derivation of damage on chuteway slabs was relatively insufficient in the literature. In this study, the cavitation erosion and hydraulic jacking were assumed to be the causes of damage on chuteway slabs, and the phenomena were reproduced using 3D numerical models, FLOW-3D and COMSOL Multiphysics. In addition, the cavitation index was calculated and the von Mises stress by uplift pressure distribution was compared with tensile and bending strength of concrete to evaluate the possibility of cavitation erosion and hydraulic jacking. As a result of numerical simulation on cavitation erosion and hydraulic jacking under various flow conditions with complete opening gate, the cavitation index in the downstream of spillway was less than 0.3, and the von Mises stress on concrete was 4.6 to 5.0 MPa. When von Mises stress was compared with tensile and bending strength of concrete, the fatigue failure caused by continuous pressure fluctuation occurred on chuteway slabs. Therefore, the cavitation erosion and hydraulic jacking caused by high speed flow were one of the main causes of damage to the chuteway slabs in spillway. However, this study has limitations in that the various shape conditions of damage(cavity and crack) and flow conditions were not considered and Fluid-Structure Interaction (FSI) was not simulated. If these limitations are supplemented and reviewed, it is expected to derive more efficient utilization of the maintenance plan on spillway in the future.

• Journal of the Korea Society for Simulation
• /
• v.18 no.3
• /
• pp.73-82
• /
• 2009

Due to the many advantages including low price, low power consumption, and miniaturization, the CMOS camera has been utilized in many applications, including mobile phones, the automotive industry, medical sciences and sensoring, robotic controls, and research in the security field. In particular, the 360 degree omni-directional camera when utilized in multi-camera applications has displayed issues of software nature, interface communication management, delays, and a complicated image display control. Other issues include energy management problems, and miniaturization of a multi-camera in the hardware field. Traditional CMOS camera systems are comprised of an embedded system that consists of a high-performance MCU enabling a camera to send and receive images and a multi-layer system similar to an individual control system that consists of the camera's high performance Micro Controller Unit. We proposed the SL-AVS (Small Size/Low power Around-View System) to be able to control a camera while collecting image data using a high speed synchronization technique on the foundation of a single layer low performance MCU. It is an initial model of the omni-directional camera that takes images from a 360 view drawing from several CMOS camera utilizing a 110 degree view. We then connected a single MCU with four low-power CMOS cameras and implemented controls that include synchronization, controlling, and transmit/receive functions of individual camera compared with the traditional system. The synchronization of the respective cameras were controlled and then memorized by handling each interrupt through the MCU. We were able to improve the efficiency of data transmission that minimizes re-synchronization amongst a target, the CMOS camera, and the MCU. Further, depending on the choice of users, respective or groups of images divided into 4 domains were then provided with a target. We finally analyzed and compared the performance of the developed camera system including the synchronization and time of data transfer and image data loss, etc.

• Progress in Medical Physics
• /
• v.20 no.4
• /
• pp.290-297
• /
• 2009

In case of radiation treatment using small field high-energy photon beams, an accurate dosimetry is a challenging task because of dosimetrically unfavorable phenomena such as dramatic changes of the dose at the field boundaries, dis-equilibrium of the electrons, and non-uniformity between the detector and the phantom materials. In this study, the absorbed dose in the phantom was measured by using an ion chamber and a diode detector widely used in clinics. $GAFCHROMIC^{(R)}$ EBT films composed of water equivalent materials was also evaluated as a small field detector and compared with ionchamber and diode detectors. The output factors at 10 cm depth of a solid phantom located 100 cm from the 6 MV linear accelerator (Varian, 6 EX) source were measured for 6 field sizes ($5{\times}5\;cm^2$, $2{\times}2\;cm^2$, $1.5{\times}1.5\;cm^2$, $1{\times}1\;cm^2$, $0.7{\times}0.7\;cm^2$ and $0.5{\times}0.5\;cm^2$). As a result, from $5{\times}5\;cm^2$ to $1.5{\times}1.5\;cm^2$ field sizes, absorbed doses from three detectors were accurately identified within 1%. Wheres, the ion chamber underestimated dose compared to other detectors in the field sizes less than $1{\times}1\;cm^2$. In order to correct the observed underestimation, a convolution method was employed to eliminate the volume averaging effect of an ion chamber. Finally, in $1{\times}1\;cm^2$ field the absorbed dose with a diode detector was about 3% higher than that with the EBT film while the dose with the ion chamber after volume correction was 1% lower. For $0.5{\times}0.5\;cm^2$ field, the dose with the diode detector was 1% larger than that with the EBT film while dose with volume corrected ionization chamber was 7% lower. In conclusion, the possibility of $GAFCHROMIC^{(R)}$ EBT film as an small field dosimeter was tested and further investigation will be proceed using Monte Calro simulation.

• Proceedings of the Korea Society of Poultry Science Conference
• /
• 2004.11a
• /
• pp.80-90
• /
• 2004

Three experiments were performed to test the assumption that imbalanced dietary amino acid mixtures must lead to increased heat production (HP). The first experiment was based on diets formulated to have a wide range of crude protein (CP) concentrations but a fixed concentration of lysine, formulated to be the first-limiting amino acid. In the second (converse) experiment, lysine concentration was varied over a wide range while CP content was kept constant. To prevent the masking of dietary effects by thermoregulatory demands, the third experiment was performed at 30 $^{\circ}C$ with the diets similar to the diets used in the second experiment. The detailed relationships among amino acid balance, nitrogen (N) metabolism and energy (E) metabolism were investigated in a computer-controlled chamber calorimetry system. The results of experiments were compared with the predictions of a computerised simulation model of E metabolism. In experiment 1. with constant lysine and varying CP, there was a 75 % increase in N intake as CP concentration increased. This led to a 150 % increase in N excretion. with no significant change in HP. Simulated HP agreed with the empirically determined results in not showing a trend with dietary CP. In experiment 2, with varying lysine but constant CP, there was a 3-fold difference in daily weight gain between the lowest and highest lysine diets. HP per bird increased significantly with dietary lysine concentration. There was still an effect when HP was adjusted for body weight differences, but it failed to maintain statistical significance. Simulated HP results agreed in showing little effect of varying lysine concentration and growth rate on HP. Based on the results of these two experiments, the third experiment was designed to test the response of birds to dietary lysine in high ambient temperature. In experiment 3 which performed at high ambient temperature (30 $^{\circ}C$), HP per bird increased significantly with dietary lysine content, whether or not adjusted for body-weight. The trend was greater than in the previous experiment (20 $^{\circ}C$).

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.10 no.2
• /
• pp.67-85
• /
• 2007

A three-dimensional hydrodynamic model with the fine grid is applied to simulate the barotropic tides, tidal currents, residual currents and salinity dispersions in the Yellow Sea and the East China Sea. Data inputs include seasonal hydrography, mean wind and river input, and oceanic tides. Computed tidal distributions of four major tides($M_2,\;S_2,\;K_1$ and $O_1$) are presented and results are in good agreement with the observations in the domain. The model reproduces well the tidal charts. The tidal residual current is relatively strong around west coast of Korea including the Cheju Island and southern coast of China. The current by $M_2$ has a maximum speed of 10 cm/s in the vicinity of Cheju Island with a anti-clockwise circulation in the Yellow Sea. General tendency of the current, however, is to flow eastward in the South Sea. Surface residual current simulated with $M_2$ and with $M_2+S_2+K_1+O_1$ tidal forcing shows slightly different patterns in the East China Sea. The model shows that the southerly wind reduces the southward current created by freshwater discharge. In summer during high runoff(mean discharge about $50,000\;m^3/s$ of Yangtze), low salinity plume-like structure(with S < 30.0 psu) extending some 160 km toward the northeast and Changjiang Diluted Water(CDW), below salinity 26 psu, was found within about 95 km. The offshore dispersion of the Changjiang outflow water is enhanced by the prevailing southerly wind. It is estimated that the inertia of the river discharge cannot exclusively reach the around sea of Cheju Island. It is noted that spatial and temporal distribution of salinity and the other materials are controlled by mixture of Changjiang discharge, prevailing wind, advection by flowing warm current and tidal current.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.11 no.4
• /
• pp.181-190
• /
• 2008

To response climate change and Kyoto protocol and to reduce greenhouse gas emissions, marine geological storage of $CO_2$ is regarded as one of the most promising option. Marine geological storage of $CO_2$ is to capture $CO_2$ from major point sources(eg. power plant), to transport to the storage sites and to store $CO_2$ into the marine geological structure such as deep sea saline aquifer. To design a reliable $CO_2$ marine geological storage system, it is necessary to perform numerical process simulation using thermodynamic equation of state. The purpose of this paper is to compare and analyse the relevant equations of state including ideal, BWRS, PR, PRBM and SRK equation of state. To evaluate the predictive accuracy of the equation of the state, we compared numerical calculation results with reference experimental data. Ideal and SRK equation of state did not predict the density behavior above $29.85^{\circ}C$, 60 bar. Especially, they showed maximum 100% error in supercritical state. BWRS equation of state did not predict the density behavior between $60{\sim}80\;bar$ and near critical temperature. On the other hand, PR and PRBM equation of state showed good predictive capability in supercritical state. Since the thermodynamic conditions of $CO_2$ reservoir sites correspond to supercritical state(above $31.1^{\circ}C$ and 73.9 bar), we conclude that it is recommended to use PR and PRBM equation of state in designing of $CO_2$ marine geological storage process.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.11 no.4
• /
• pp.191-198
• /
• 2008

To design a reliable $CO_2$ marine geological storage system, it is necessary to perform numerical process simulation using thermodynamic equation of state. $CO_2$ capture process from the major point sources such as power plants, transport process from the capture sites to storage sites and storage process to inject $CO_2$ into the deep marine geological structure can be simulate with numerical modeling. The purpose of this paper is to compare and analyse the relevant equations of state including ideal, BWRS, PR, PRBM and SRK equation of state. We also studied the effect of thermodynamic equation of state in designing the compression and transport process. As a results of comparison of numerical calculations, all relevant equation of state excluding ideal equation of state showed similar compression behavior in pure $CO_2$. On the other hand, calculation results of BWRS, PR and PRBM showed totally different behavior in compression and transport process of captured $CO_2$ mixture from the oxy-fuel combustion coal-fired plants. It is recommended to use PR or PRBM in designing of compression and transport process of $CO_2$ mixture containing NO, Ar and $O_2$.