• Journal of the Korean Geotechnical Society
• /
• v.34 no.3
• /
• pp.13-27
• /
• 2018

To study how stable the rubble mound breakwaters are, one can look to the research of wave induced seepage flow through the pores of the rubble mound. Seepage flow is generally generated by the difference between the water level around the breakwater during a typhoon. The existing stability analysis method of the rubble mound is the static analysis which simply considers the force equilibrium taking into account the horizontal force acting on the concrete block induced by a wave (calculated by Goda equation) and the vertical force induced by the weight inclusive of the concrete block, quarry run, filter, and armor layer above the slipping plane. However, this static method does not consider the wave-induced seepage flow in the rubble mound. Such seepage may decrease the stability of the rubble mound. The stability of a rubble mound breakwater under the action of seepage was studied based on the results of CFD software (OpenFOAM) and Limit Equilibrium Method (GeoStudio). The numerical analysis result showed that the seepage flow decreased the stability of the rubble mound breakwaters. The results of the numerical analyses also revealed the stability of the rubble mound was varied with time. Especially, the most critical state happened at the condition of overtopping the concrete block, acting strong uplift pressure raising along side and underneath the concrete block, and generating high pore pressure inside rubble mound due to seepage flow. Therefore, it may be necessary to conduct a dynamic analysis considering the effect of wave-induce seepage flow together with the static analysis.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.6
• /
• pp.207-216
• /
• 2003

To evaluate the smear zone caused by the drain installation, 10 piezometers were installed in the typical soft ground in the western coastal area of Korea. The dynamic pore water pressure developed during the drain installation was monitored using piezometers installed at a distance of 10cm, 20cm, 30cm, 40cm and 50cm from the location of the drain. The decay of pore pressure with time after pushing piezometers to depths of 5 meters and 7 meters during the drain installation was monitored to assess flow and consolidation characteristics of the soil after disturbance of the soil due to the drain installation. The drain installation results in shear strain and displacement of the soil and it decreases the permeability of the soil. Hence, the comparison between dissipation of the pore water pressure process in 10 pieszometers before as well as after installation of the drain indicated the diameter of disturbance zone and smear zone, which is related to the cross-sectional dimension of the mandrel. In addition, Cone-pressuremeter(CPM) tests were performed to obtain rigidity index of the soil for an interpretation of the dissipation processes. It has been evaluated by in-situ tests that the smear zone is from 3.0 to 3.6 times of the cross-sectional dimension of the mandrel. The hydraulic conductivity expressed in terms of the coefficient of consolidation after the drain installation was calculated from 3 to 8 times decrease evaluated by Teh & Houlsby equation and CPM test results.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.1
• /
• pp.82-89
• /
• 2016

In this study, we proposed a pier of pontoon-type, "Very Large Floating Structure" (VLFS), with the length of 500m, breadth of 200 m and height of 2 m in Yeosu domestic port. Since this structure ought to endure wave loads for long periods at sea, it is essential to analyze the wave response characteristics. Direct-method is used to analyze the fluid-structure problem and the coupled motion of equation is used to obtain response results. The structural part is calculated by using finite element method (FEM) and the fluid part is analyzed by using boundary element method (BEM). Dynamic responses caused by the elastic deformation and rigid motion of structure are analyzed by numerical calculation. To investigate response characteristics of the pier in regular waves, several factors such as the wavelength, water depth, wave direction and flexural rigidity of structure are considered. As a result, wave response of pier changed at the point of $L/{\lambda}$ 1.5 and represented the torsional phenomenon according to the various incident waves. And the responses showed increasing tendency as the water depths increase at the incident point in case of $L/{\lambda}=8.0$ and peak point of vertical displacement amplitude moved from side to side as the flexural rigidity of structure changes.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.4
• /
• pp.225-243
• /
• 2014

As the economic matters are involved, applying the WEC, which is used for controlling waves as well as utilizing the wave energy on existing breakwater, is preferred rather than installing exclusive WEC. This study examines the OWC mounted on a pressurized air chamber floating breakwater regarding the functionality of both breakwater and wave-power generation. In order to verify the performance as a WEC, the velocity of air flow from pressurized air chamber to WEC has to be evaluated properly. Therefore, numerical simulation was implemented based on BEM from linear velocity potential theory as well as Boyle's law with the state equation to analyze pressurized air flow. The validity of the obtained values can be determined by comparing the previous results from numerical analysis and empirically obtained values of different shapes. In the actual numerical analysis, properties of wave deformation around OWC system mounted on fixed type and floating type breakwaters, motions of the structure with air flow velocities are investigated. Since, the wind power generating system can be hybridized on the structure, it is expected to be applied on complex power generation system which generates both wind and wave power energy.

• Journal of Korea Soil Environment Society
• /
• v.1 no.1
• /
• pp.89-102
• /
• 1996

An integrated model is presented to describe underground flow and mass transport, using a multicomponent multiphase approach. The comprehensive governing equation is derived considering mass and force balances of chemical species over four phases(water, oil, air, and soil) in a schematic elementary volume. Compact and systemati notations of relevant variables and equations are introduced to facilitate the inclusion of complex migration and transformation processes, and variable spatial dimensions. The resulting nonlinear system is solved by a multidimensional finite element code. The developed code with dynamic array allocation, is sufficiently flexible to work across a wide spectrum of computers, including an IBM ES 9000/900 vector facility, SP2 cluster machine, Unix workstations and PCs, for one-, two and three-dimensional problems. To reduce the computation time and storage requirements, the system equations are decoupled and solved using a banded global matrix solver, with the vector and parallel processing on the IBM 9000. To avoide the numerical oscillations of the nonlinear problems in the case of convective dominant transport, the techniques of upstream weighting, mass lumping, and elementary-wise parameter evaluation are applied. The instability and convergence criteria of the nonlinear problems are studied for the one-dimensional analogue of FEM and FDM. Modeling capacity is presented in the simulation of three dimensional composite multiphase TCE migration. Comprehesive simulation feature of the code is presented in a companion paper of this issue for the specific groundwater or flow and contamination problems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.3
• /
• pp.61-71
• /
• 1988

The stage-discharge relation curve(rating curve) is the basic formula in hydrologic analysis. It plays an important role in converting to the discharge from available flood water level data including the daily mean stage. However, the river induces a cross section change at the gauging station because of the composed material of the river bed and three processes of the stream flow; i.e., erosion, transportation, and sedimentation. Rating curve has to be revised according to the temporal variation of the river bed due to the those factors. In this study, the basic rating curve is developed with respect to the current river bed to convert the existing rating curves and also to seize the hydraulic and geometric characteristics for the temporal variation of the river bed, relationships among the basic rating curve and the existing rating curves, water level, cross sectional area, and flow velocity are analyzed. Indogyo station, which is not only the key station of the Han river but also greatly changed the river bed after completion of the Han river development plan during the year 1983 to 1986, was chosen for the study. In this study, the river bed is assumed in a dynamic equilibrium condition. The basic rating curve is developed using hydrologic data of the physical year of 1987. For a given discharge, relationships for conversion of previous data, stage and velocity, the current one are formulated. To verify the usefulness of the relationships, stage-cross sectional area and stage velocity formula are also derived. Both hydrologic method using continuity equation and statistical method by the rating curve are compared and checked, then the validation of the both are positively shown.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.343-351
• /
• 2016

Cryogenic spray characteristics of a nitrogen swirl injector operating in supercritical environment have been numerically investigated. By comparing the equation of states(EOS) used for supercritical condition, SRK EOS was applied to predict the nitrogen thermodynamic property under supercritical environment. A Chung's method was implemented for the calculation of viscosity and conductivity and Takahashi's correlation based on Fuller's Theorem was implemented for the calculation of diffusion coefficient. By injecting the nitrogen with 5 bar differential pressure into 50 bar chamber filled with nitrogen, numerical simulation has been conducted. The dynamic Smagorinsky sub-grid scale (SGS) model has been compared with the algebraic Smagorinsky SGS model using FFT frequency analysis. The instability at the liquid film and gas core inside injector and the propagation of pressure oscillation into the injector has been investigated. The spreading angle of swirl injector obtained by numerical calculation has been validated with experimental result.

• Journal of the Korean GEO-environmental Society
• /
• v.16 no.12
• /
• pp.5-12
• /
• 2015

Since blasting noise is impact noise, it will give a sudden shock to the human. In the case, such as the blast vibration, it has given aging buildings and livestock great damage to move the vibration along ground in nearby regions. In this study, the influence of the blasting generated during excavation was analyzed for effects on regional. A couple of field and laboratory surveys about geological were carried out to figure out the geological ratio in the study-performed area. Blast vibration noise was compared to the domestic and abroad case studies and each of the institutions permissible standards established the most appropriate criteria in site condition. The vibration velocity of blasting vibration exploits the values which were measured from test blasting on the ground in order to examine blasting effect. Considering the blasting point as the shortest distance from safety facilities (farms, private houses, etc.), the examination of the influence range, which was derived from the vibration velocity of blasting vibration, was performed to figure out how the point affected the safety facilities. Three-dimensional numerical analysis was performed a time history analysis in order to analyze the behavior of the structure for a dynamic blast load, which was determined in three directions of the blast vibration value. The results of three-dimensional numerical analysis and the blasting effect of blasting vibration estimation equation blasting vibration of impact circle with accompanying test blasting were compared. And the analysis confirmed similar results figures.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.3
• /
• pp.35-47
• /
• 1996

A surface panel method treating a boundary-value problem of the Dirichlet type is presented to design a three dimensional body with free surface corresponding to a prescribed pressure distribution. An integral equation is derived from Green's theorem, giving a relation between total potential of known strength and the unknown local flux. Upon discretization, a system of linear simultaneous equations is formed including free surface boundary condition and is solved for an assumed geometry. The pseudo local flux, present due to the incorrect positioning of the assumed geometry, plays a role f the geometry corrector, with which the new geometry is computed for the next iteration. Sample designs for submerged spheroids and Wigley hull and carried out to demonstrate the stable convergence, the effectiveness and the robustness of the method. For the calculation of the wave resistance, normal dipoles and Rankine sources are distributed on the body surface and Rankine sources on the free surface. The free surface boundary condition is linearized with respect to the oncoming flow. Four-points upwind finite difference scheme is used to compute the free surface boundary condition. A hyperboloidal panel is adopted to represent the hull surface, which can compensate the defects of the low-order panel method. The design of a 5500TEU container carrier is performed with respect to reduction of the wave resistance. To reduce the wave resistance, calculated pressure on the hull surface is modified to have the lower fluctuation, and is applied as a Dirichlet type dynamic boundary condition on the hull surface. The designed hull form is verified to have the lower wave resistance than the initial one not only by computation but by experiment.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.137-145
• /
• 2014

We present a design sensitivity analysis(DSA) method for multiscale problems based on bridging scale decomposition. In this paper, we utilize a bridging scale method for the coupled system analysis. Since the analysis of full MD systems requires huge amount of computational costs, a coupled system of MD-level and continuum-level simulation is usually preferred. The information exchange between the MD and continuum levels is taken place at the MD-continuum boundary. In the bridging scale method, a generalized Langevin equation(GLE) is introduced for the reduced MD system and the GLE force using a time history kernel is applied at the boundary atoms in the MD system. Therefore, we can separately analyze the MD and continuum level simulations, which can accelerate the computing process. Once the simulation of coupled problems is successful, the need for the DSA is naturally arising for the optimization of macro-scale design, where the macro scale performance of the system is maximized considering the micro scale effects. The finite difference sensitivity is impractical for the gradient based optimization of large scale problems due to the restriction of computing costs but the analytical sensitivity for the coupled system is always accurate. In this study, we derive the analytical design sensitivity to verify the accuracy and applicability to the design optimization of the coupled system.