• Journal of the Korean Society for Geothermal and Hydrothermal Energy
• /
• v.17 no.2
• /
• pp.1-10
• /
• 2021

Lithium-ion batteries with high energy density, long cycle life and other advantages, have been widely used to energy storage systems(ESS). But as ESS fires frequently occur, the safety concern has become the main obstacle that hinders the large-scale applications of lithium-ion batteries. Especially, thermal runaway is the key scientific problem in battery safety research. Therefore, in this study, we performed a numerical analysis on the thermal runaway phenomenon of NCM111, NCM523 and NCM622 batteries using a two-dimensional analysis model. The results show that the two-dimensional simulation results are generally matched with three-dimensional simulation. Also, In the case of NCM111 with a low Ni content in the temperature range used in this study, thermal runaway phenomenon does occurred very slowly, but as the Ni content is increased, the thermal runaway phenomenon occurs rapidly and the thermal stability tends to be decreased. And, in NCM523 and NCM622 batteries, chain reactions occur almost simultaneously, but in the case of NCM111 battery, it is found that after the SEI(Solid Electrolyte Interface) layer decomposition reaction, the cathode-electrolyte reaction is appeared sequentially. After that, the anodic decomposition reaction is increased and leads to the thermal runaway reaction.

• Geomechanics and Engineering
• /
• v.17 no.1
• /
• pp.1-9
• /
• 2019

This paper proposes a new numerical approach to model geocell reinforced soils, where the geocell is described as membrane elements and the complex interaction between geocell and soil is realized by coupling their degrees of freedom. The effectiveness and robustness of this approach are demonstrated using two examples, i.e., a geocell-reinforced foundation and a large scale retaining wall project. The first example validates the approach against established solutions through a comprehensive parametrical study to understand the influence of geocell on the improvement of bearing capacity of foundations. The study results show that reducing the geocell pocket size has a strong effect on improving the bearing capacity. In addition, when the aspect ratio maintains the same value, the bearing capacity improvement with increasing geocell height is insignificant. Comparing with the field monitoring and measurement in the project, the second example investigates the application of the approach to practical engineering projects. This paper provides a practically feasible and efficient modelling approach, where no explicit interface or contact is required. This allows geocell reinforced soils in large scale project can be effectively modelled where the mechanism for complex geocell-soil interaction can be explicitly observed.

• Structural Engineering and Mechanics
• /
• v.84 no.6
• /
• pp.743-765
• /
• 2022

Near-surface excavations may cause the tilting and destruction of the adjacent superstructures in big cities. The stability of a huge excavation and its nearby superstructures was studied in this paper. Some test instruments monitored the deformation and loads at the designed location. Then the numerical models of the excavation were made in FLAC3D (a three-dimensional finite difference code) and Plaxis-3D (a three-dimensional finite element code). The effects of different supporting and reinforcement tools such as nails, piles, and shotcretes on the stability and bearing capacity of the foundation were analyzed through different numerical models. The numerically approximated results were compared with the corresponding in-field monitored results and reasonable compatibility was obtained. It was concluded that the displacement in excavation and the settlement of the nearby superstructure increases gradually as the depth of excavation rises. The effects of support and reinforcements were also observed and modeled in this study. The settlement of the structure gradually decreased as the supports were installed. These analyses showed that the pile significantly increased the bearing capacity and decreased the settlement of the superstructure. As a whole, the monitoring and numerical simulation results were in good consistency with one another in this practically important project.

• Journal of the Korean Geotechnical Society
• /
• v.37 no.4
• /
• pp.5-17
• /
• 2021

In this paper, three-dimensional analysis of the curved levee was performed according to curvature angle, and radius of curvature to investigate the property of seepage. The hydraulic gradients in the curved parts of levees decreased in the outer levee and increased in the inner levee, compared to the two-dimensional analysis. The smaller the curvature angle and the radius of curvature, the larger the change of the hydraulic gradient, compared to the two-dimensional analysis. The effect of curvature radius on the hydraulic gradient was greater than the curvature angle. As a result of evaluating the piping safety factor for the critical hydraulic gradient, the safety factor was increased by 2~5% in the outer levee and decreased by 4~12% in the inner levee, compared to the two-dimensional analysis. Considering this reduction, if the two-dimensional analysis is performed on the curved part of the levee, and if designed the safety factor for piping is 0.1~0.3 greater than allowable FS=2.0, the safety factor of the curved part is slightly reduced, but there is no difficulty in securing stability.

• Tribology and Lubricants
• /
• v.36 no.5
• /
• pp.262-266
• /
• 2020

Herein, we present a numerical investigation of wear analysis of sliding systems with a constant speed subjected to Archard's wear law. For this investigation, we compared two methods: eigen-wear analysis and adaptive meshing technique. The eigen-wear analysis is advantageous to predict the evolution of contact pressure due to wear using the initial contact pressure and contact stiffness. The adaptive meshing technique in finite element analysis is employed to obtain transient wear behavior, which needs significant computational resources. From the eigen-wear analysis, we can determine the appropriate element size required for finite element analysis and the time increment required for wear evolution by a dimensionless variable above a certain value. Since the prediction of wear depends on the maximum contact pressure, the finite element model should have a reasonable representation of the maximum contact pressure. The maximum contact pressure and wear amount according to this dimensionless variable shows that the number of fine meshes in the contact area contributes more to the accuracy of the wear analysis, and the time increment is less sensitive when the number of contact nodes is significantly larger. The results derived from a two-dimensional wear model can be applied to a three-dimensional wear model.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.3
• /
• pp.180-198
• /
• 2019

In recent years, countries like Europe and Japan have been involved in many researches on the Low-Crested Structure (LCS) which is the method to protect beach erosion and it is regarded as an alternative to the submerged breakwaters, and compiled its results and released the design manual. In the past, studies on LCS have focused on two-dimensional wave transmission and calculating required weight of armor units, and these were mainly examined and discussed based on experiments. In this study, three-dimensional numerical analysis is performed on permeable LCS. The open-source CFD code olaFlow based on the Navier-Stokes momentum equations is applied to the numerical analysis, which is a strongly nonlinear analysis method that enables breaking and turbulence analysis. As a result, the distribution characteristics of the LCS such as water level, water flow, and turbulent kinetic energy were examined and discussed, then they were carefully compared and examined in the case of submerged breakwaters. The study results indicate that there is a difference between the flow patterns of longshore current near the shoreline, the spatial distribution of longshore and on-offshore directions of mean turbulent kinetic energy in case of submerged breakwaters and LCS. It is predicted that the difference in these results leads to the difference in sand movement.

• 한국정보디스플레이학회:학술대회논문집
• /
• 2007.08a
• /
• pp.517-519
• /
• 2007

We report our study on the implementation of Q tensor approach into three-dimensional finite element method (FEM) numerical solver. The comparative simulation results demonstrated the possibility of a different director configuration in between Q tensor method and vector method. The comparative study confirmed that Q Tensor implementation is more appropriate for OCB analysis than the vector method.

• Journal of the Korean Geosynthetics Society
• /
• v.18 no.3
• /
• pp.23-32
• /
• 2019

This study investigated the scour characteristics of artificial reef according to unreinforced and reinforced artificial reefs through laboratory test and numerical analysis. In this study, geogrid was used as a reinforcement of seabed. Three different reinforced areas were selected; one time, three times and five times bottom area of artificial reef. Two-dimensional water channel test was carried out to investigate scour patten of sand and silt grounds. Numerical analysis was also performed to see the velocity vector and scouring patten around artificial reef which was same condition with experiment. The numerical results were in good agreements with the experimental test results. It was found that the geogrid-installed soft ground under artificial reef tends to reduce scour, compared to unreinforced soft ground.

• Journal of Navigation and Port Research
• /
• v.31 no.8
• /
• pp.675-681
• /
• 2007

A numerical procedure is described for estimating the effects of the multi-directional irregular waves on the structural responses of the Tension Leg Platform (TLP). The numerical approach is based on a three dimensional source distribution method for hydrodynamic forces, a three dimensional frame analysis method for structural responses, in which the superstructure of TLP is assumed to be flexible instead of rigid. Hydrodynamic and hydrostatic forces on the submerged surface of a TLP have been accurately calculated by excluding the assumption of the slender body theory. The hydrodynamic interactions among TLP members, such as columns and pontoons, and the structural damping are included in structural analysis. The spectral description used in spectral analysis of directional waves for the linear system of a TLP in the frequency domain is sufficient to completely define the structural responses. This is due to both the wave inputs and responses are stationary Gaussian random process of which the statistical properties in the amplitude domain are well known. The numerical results for the linear motion responses and tension variations in regular waves are compared with the experimental and numerical ones, which are obtained in Yoshida et al.(1983). The results of comparison confirmed the validity of the proposed approach.

• Structural Engineering and Mechanics
• /
• v.3 no.4
• /
• pp.325-340
• /
• 1995

The iterative procedure to use the nonconforming elements in elasto-plastic problems is established and applied to the variable node transition solid element developed for the automated three-dimensional local mesh refinement. Through numerical tests, the validity and performance of the element are examined. As the nonlinear iterative procedure presented in this paper is accomplished for the general three-dimensional case, it can also be easily applied to the two-dimensional elements such as membranes, plates and shells.