• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.4
• /
• pp.313-327
• /
• 2008

Thermo-mechanical coupled behavior of an underground structure during a fire accident have not been fully understood yet. Moreover, when such a thermo-mechanical coupled behavior is not considered in numerical analyses based on conventional heat transfer theory, fire-induced damage zone in an underground structure can be considerably underestimated. This study aims to develop a FEM-based numerical technique to simulate the thermo-mechanical coupled behavior of an underground structure in a fire accident. Especially, an element elimination model is newly proposed to simulate fire-induced structural loss together with a convective boundary condition. In the proposed model, an element where the maximum temperature calculated from heat transfer analysis is over a prescribed critical temperature is eliminated. Then, the proposed numerical technique is verified by comparing numerical results with experimental results from real fire model tests. From a series of parametric studies, the key parameters such as critical temperature, element size and temperature-dependent convection coefficients are optimized for the RABT and the RWS fire scenarios.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.6
• /
• pp.461-468
• /
• 2012

In this study, a Cartesian-grid method based on finite volume approach is applied to simulate the ship motions in large amplitude waves. Fractional step method is applied for pressure-velocity coupling and TVD limiter is used to interpolate the cell face value for the discretization of convective term. Water, air, and solid phases are identified by using the concept of volume-fraction function for each phase. In order to capture the interface between air and water, the tangent of hyperbola for interface capturing (THINC) scheme is used with weighed line interface calculation (WLIC) method which considers multidimensional information. The volume fraction of solid body embedded in the Cartesian grid system is calculated using a level-set based algorithm, and the body boundary condition is imposed by a volume weighted formula. Numerical simulations for the two-dimensional barge type model and Wigley hull in linear waves have been carried out to validate the newly developed code. To demonstrate the applicability for highly nonlinear wave-body interactions such as green water on the deck, numerical analysis on the large-amplitude motion of S175 containership is conducted and all computational results are compared with experimental data.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.2
• /
• pp.67-75
• /
• 2008

It is interesting to note that distinct element method has been used extensively to model the response of micro and discontinuous behavior in geomechanics. Impressive advances related to response of distinct particles have been conducted and there were difficulties in considering fluid effect simultaneously. Current distinct element methods are progressively developed to solve particle-fluid coupling focused on fluid flow through soil, rock or porous medium. In this research, numerical simulations of fluid injection into particulate materials were conducted to observe cavity initiation and propagation using distinct element method. After generation of initial particles and wall elements, confining stress was applied by servo-control method. The fluid scheme solves the continuity and Navior-Stokes equations numerically, then derives pressure and velocity vectors for fixed grid by considering the existence of particles within the fluid cell. Fluid was injected as 7-step into the assembly in the x-direction from the inlet located at the center of the left boundary under confining stress condition, $0.1MP{\alpha}\;and\;0.5MP{\alpha}$, respectively. For each simulation, movement of particles, flow rate, fluid velocity, pressure history, wall stress including cavity initiation and propagation by interaction of flulid-paricles were analyzed.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.3 no.1
• /
• pp.16-24
• /
• 1998

In this study, Ocean General Circulation Model (OGCM) has been developed as a counterpart of Atmospheric General Circulation (AGCM) for the study of coupled ocean-atmosphere climate system. The oceanic responses to given atmospheric boundary conditions have been investigated using the OGCM. In an integration carried out over 100 simulated years with climatological monthly mean data (EXP 1), most parts of the model reached a quasi-equilibrium climate reproducing many of the observed large-scale oceanic features remarkably well. Some observed narrow currents, however, such as North Equatorial Counter Current, were inevitably distorted due to the model's relatively coarse resolution. The seasonal changes in sea ice cover over the southern oceans around Antarctica were also simulated. In an experiment (EXP 2) under boundary condition of 10-year monthly data (1982-1991) from NCEP/NCAR Reanalysis Project model properly reproduced major oceanic changes during the period, including El Ni$\tilde{n}$os of 1982-1983 and 1986-87. During the ENSO periods, the experiment showed eastward expansion of warm surface waters and a negative vertical velocity anomalies along' the equator in response to expansion of westerly current velocity anomalies as westerly wind anomalies propagated eastward. Simulated anomalous distribution and the time behavior in response to El Ni$\tilde{n}$o events is consistent with that of the observations. These experiments showed that the model has an ability to reproduce major mean and anomalous oceanic features and can be effectively used for the study of ocean-atmosphere coupling system.

• Journal of the Korean Electrochemical Society
• /
• v.14 no.4
• /
• pp.231-238
• /
• 2011

In this work, the numerical analysis for the zincate behavior at a zinc electrode with an electrolyte flow was carried out for a ZAFC. The Nernst-Planck equation with a boundary condition of Butler-Volmer type was adopted to describe electrochemical effects of mass transfer, migration, kinetics of electrode. The Navier-Stokes equation, coupling to the Nernst-Planck equation, is also applied to describe the internal electrolyte flow fields. The validity of the numerical model is proved through the comparative analysis between numerical and experimental results. The concentration of zincate and the current density were also investigated at a zinc anode according to various electrolyte velocities. We have found the concentration of zincate decreased and the current density increased with an increase in the electrolyte velocity.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.549-560
• /
• 2021

Forces exerted on a shield TBM (tunnel boring machine) such as cutter head torque, thrust force, chamber pressure, and upward force are key factors determining TBM performance. However, the forces acting on the TBM when tunnelling the mixed ground have different tendencies compared to that of the uniform ground, which could impair TBM performance. In this study, the effect of mixed ground tunnelling was numerically investigated with torque, thrust force, chamber pressure, and upward force. A coupled discrete element method (DEM) and finite difference method (FDM) model for TBM driving model was used. This numerical study simulates TBM tunnelling in mixed ground composed of upper weathered granite soil and lower weathered rock. The effect on the force acting on the TBM according to the location and slope of the boundary of the mixed ground was numerically examined.

• The Journal of the Acoustical Society of Korea
• /
• v.35 no.4
• /
• pp.261-271
• /
• 2016

In the prediction of response of a pile in soil, numerical approaches such as a finite element method are generally applied due to complicate nonlinear behaviors of soils. However, the numerical methods based on the finite elements require heavy efforts in pile and soil modelling and also take long computing time. So their usage is limited especially in the early design stage in which principal dimensions and properties are not specified and tend to vary. On the contrary, theoretical approaches adopting linear approximations for soils are relatively simple and easy to model and take short computing time. Therefore, if they are validated to be reliable, they would be applicable in predicting responses of a pile in soil, particularly in early design stage. In case of wind turbines regarded in this study, it is required to assess their natural frequencies in early stages, and in this simulation the supporting pile inserted in soil could be replaced with a simplified elastic boundary condition at the bottom end of the wind turbine tower. To do this, analysis for a pile in soil is performed in this study to extract the spring constants at the top end of the pile. The pile in soil can be modelled as a beam on elastic spring by assuming that the soils deform within an elastic range. In this study, it is attempted to predict pile deformations and influence factors for lateral loads by means of the beam-on-spring model. As two example supporting structures for wind turbines, mono pile and suction pile models with different diameters are examined by evaluating their influence factors and validated by comparing them with those reported in literature. In addition, the deflection profiles along the depth and spring constants at the top end of the piles are compared to assess their supporting features.