• Geomechanics and Engineering
• /
• v.23 no.3
• /
• pp.189-206
• /
• 2020

The accurate modeling of boundary conditions is important in simulations of the discrete element method (DEM) and can affect the numerical results significantly. In conventional triaxial compression (CTC) tests, the specimens are wrapped by flexible membranes allowing to deform freely. To accurately model the boundary conditions of CTC, new flexible boundary algorithms for 2D and 3D DEM simulations are proposed. The new algorithms are computationally efficient and easy to implement. Moreover, both horizontal and vertical component of confining pressure are considered in the 2D and 3D algorithms, which can ensure that the directions of confining pressure are always perpendicular to the specimen surfaces. Furthermore, the boundaries are continuous and closed in the new algorithms, which can prevent the escape of particles from the specimens. The effectiveness of the proposed algorithms is validated by biaxial and triaxial simulations of granular materials. The results show that the algorithms allow the boundaries to deform non-uniformly on the premise of maintaining high control accuracy of confining pressure. Meanwhile, the influences of different lateral boundary conditions on the numerical results are discussed. It is indicated that the flexible boundary is more appropriate for the models with large strain or significant localization than rigid boundary.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.34 no.6
• /
• pp.258-265
• /
• 2022

In Korea, the hydraulic model tests for measuring the wave overtopping have been almost conducted with no bottom slope or single slope condition in Korea. In this study, the bottom seabed for the coastal road area was fabricated at the wave flume and the wave overtopping was measured. The overtopping rate was also measured with the numerical modelling by OLAFoam. The measuring data were compared with EurOtop manual. It could be known the the influence of the foreslope in front of the vertical wall was significant and the these effects should be concerned when designing the coastal structures. And also it could be known that OLAFoam could be used to predict the wave overtopping rate for the complex bottom topography.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.1
• /
• pp.55-59
• /
• 2005

The present data concerns the wave height distribution around the Ieodo underwater rocks and it was obtained from a 3D hydraulic model experiment which was performed in 1999 by Konkuk University and Korea Ocean Research and Development Institute. The experiment was separately undertaken for 4 different wave directions (NNW, SE, S, NNW) under which wave heights were measured at every 1m interval within the preset grid area, 16m×18m. It was observed that the wave breaking occurred on the top of the Ieodo model for all wave directions. This data may be effectively used for improving or verifying the performance of numerical wave propagation models in the area with the local breaking wave zones.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.4
• /
• pp.3-11
• /
• 2004

A finite volume model is developed to simulate the surface irrigation at a paddy field. The model's capabilities are validated through comparison with the simulafed results and the observed data obtained by various experimental tests, and the simulated results are in good agreement with the observed pending depth. The result of surface irrigation simulation shows that the longer the paddy field's the length of long-sided becomes, the longer the advance and storage time is taken. To analyze surface irrigation performance with variable inflow rate, three patterns of flow variation-constant rate, initially high then low, and initially low then high-were studied. The results show that at the pattern with initially high followed by low during the latter half of the irrigation the advance time is shortest, but the pending depth of irrigation completion and irrigation effiency are the little difference between irrigation patterns.

• Geomechanics and Engineering
• /
• v.9 no.6
• /
• pp.815-827
• /
• 2015

The mechanical behavior of soil and soil-rock mixture is investigated via the discrete element method. A non-overlapping combination method of spheres is used to model convex polyhedron rock blocks of soil-rock mixture in the DEM simulations. The meso-mechanical parameters of soil and soil-rock interface in DEM simulations are obtained from the in-situ tests. Based on the Voronoi cell, a method representing volumtric strain of the sample at the particle scale is proposed. The numerical results indicate that the particle rotation, occlusion, dilatation and self-organizing force chains are a remarkable phenomena of the localization band for the soil and soil-rock mixture samples. The localization band in a soil-rock mixture is wider than that in the soil sample. The current research shows that the 3D discrete element method can effectively simulate the mechanical behavior of soil and soil-rock mixture.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.1
• /
• pp.114-125
• /
• 2017

This paper addresses a study of inner-tank sloshing effect on motion responses of a Floating Liquefied Natural Gas (FLNG) system, through experimental analysis and numerical modeling. To investigate hydrodynamic characteristics of FLNG under the conditions of with and without LNG-tank sloshing, a series of numerical simulations were carried out using potential flow solver SESAM. To validate the numerical simulations, model tests on the FLNG system was conducted in both liquid and solid ballast conditions with 75% tank filling level in height. Good correlations were observed between the measured and predicted results, proving the feasibility of the numerical modeling technique. On the verified numerical model, Response Amplitude Operators (RAOs) of the FLNG with 25% and 50% tank filling levels were calculated in six degrees of freedom. The influence of tank sloshing with varying tank filling levels on the RAOs has been presented and analyzed. The results showed that LNG-tank sloshing has a noticeable impact on the roll motion response of the FLNG and a moderate tank filling level is less helpful in reducing the roll motion response.

• Structural Engineering and Mechanics
• /
• v.44 no.6
• /
• pp.753-762
• /
• 2012

The ratcheting characteristics of cylindrical shell under cyclic axial loading are investigated. The specimens are subjected to stress-controlled cycling with non-zero mean stress, which causes the accumulation of plastic strain or ratcheting behavior in continuous cycles. Also, cylindrical shell shows softening behavior under symmetric axial strain-controlled loading and due to the localized buckling, which occurs in the compressive stress-strain curve of the shell; it has more residual plastic strain in comparison to the tensile stress-strain hysteresis curve. The numerical analysis was carried out by ABAQUS software using hardening models. The nonlinear isotropic/kinematic hardening model accurately simulates the ratcheting behavior of shell. Although hardening models are incapable of simulating the softening behavior of the shell, this model analyzes the softening behavior well. Moreover, the model calculates the residual plastic strain close to the experimental data. Experimental tests were performed using an INSTRON 8802 servo-hydraulic machine. Simulations show good agreement between numerical and experimental results. The results reveal that the rate of plastic strain accumulation increases for the first few cycles and then reduces in the subsequent cycles. This reduction is more rapid for numerical results in comparison to experiments.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.1
• /
• pp.37-47
• /
• 2008

Design and construction of tunnels require understanding the influence of groundwater. Particularly, it is essential to know how the drainage conditions at the tunnel boundary affect flow behavior of ground adjacent to the tunnels. In this study flow behavior of a leaking tunnel was investigated using physical model tests for tunnel depths and various hydraulic boundary conditions. Particular concerns were given to flow lines toward tunnels. Test results showed that the boundary conditions hardly influence on flow patterns and time required to reach steady state conditions. It is revealed that with an increase in water depth, flow lines concentrated to the drain holes. The physical tests were numerically simulated. Numerical results showed that the flow behavior was represented appropriately by considering filter-drain hole drainage rather than boundary drainage all over the lining.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.6
• /
• pp.111-120
• /
• 1999

Leachate flow behavior due to intermediate cover soil of low hydraulic conductivity and the applicability of pumping method for reducing the leachate level in the landfill are analyzed with the numerical flow model, MODFLOW. Using the hydraulic conductivity and storativity data obtained from the field pumping and slug tests(Jang and Cho, 1999), the hydraulic condition within the landfill is validated. The optimum rate of pumping, the radius of influence, and the efficiency of horizontal drain are analyzed for reducing the leachate level in the landfill. From the results of the analyses, the barrier effect that the buried cover soil of low hydraulic conductivity prevents the vertical movement of leachate flow through the cover soil, which is found from the in-situ geotechnical studies(Jang and Cho, 1999), is identified again. Also, the installation of horizontal drains to the pumping well can increase the pumping rate from 120 ton/day per a well to 300 ton/day. The length of horizontal drain did not influence significantly on the drawdown-time curve of leachate in the landfill.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.1
• /
• pp.38-49
• /
• 2015

Past study of tsunami heavily relied on the numerical modelling using 2D Boussinesq Eq. and Solitary wave. Lately, based on the fact that numerically simulated run up heights based on solitary wave are somewhat smaller than the measured one, Leading Depression N (LDN) Wave has been elaborated, which can account the advancement of a shore line before tsunami strikes a shore. Thereafter it is reported that more accurate simulation can be possible once LDN is deployed. On the other hand, there were numerous reports indicating that stable LDN wave can't be sustained in the hydraulic model test. These conflicts between the hydraulic model tests and numerical results have their roots on the assumption made in the derivation of Boussinesq type wave model such as that wave nonlinearity is equally balanced with wave dispersiveness. Hence, in the numerical simulation based on the Boussinesq type wave model, wave dispersiveness is inevitably underestimated, especially in deep water. Based on this rationale, we developed the modified methodology for the generation of stable LDN wave in the 3D numerical wave flume, and proceeded to numerically analyze the depression effect of Hybrid Breaker on the run up height due to tsunami using the Navier Stoke Equation. The verification of newly proposed wave model in this study was carried out using the run up height from the hydraulic model test. It was shown that Hybrid Breaker consisting of three water chamber and slope at its front can reduce 13% of run up height for H = 5m, and 10% of run up height for H = 6m.