• Journal of applied mathematics & informatics
• /
• v.11 no.1_2
• /
• pp.357-364
• /
• 2003

A fully discrete $H^1-mixed$ finite element approximation for the single-phase Stefan problem is introduced and the unique existence of the approximation is established. And some numerical experiments are given.

• Journal of applied mathematics & informatics
• /
• v.28 no.3_4
• /
• pp.953-966
• /
• 2010

In this paper, we develop discontinuous Galerkin methods with penalty terms, namaly symmetric interior penalty Galerkin methods to solve nonlinear parabolic equations. By introducing an appropriate projection of u onto finite element spaces, we prove the optimal convergence of the fully discrete discontinuous Galerkin approximations in ${\ell}^2(L^2)$ normed space.

• 한국신재생에너지학회:학술대회논문집
• /
• 2010.11a
• /
• pp.123.1-123.1
• /
• 2010

수평형 지중열교환기의 최적설계를 위해서는 되메움재의 광물특성 및 입자크기, 열전도도(thermal conductivity), 열용량(heat capacity)등과 같은 열적 특성을 파악 하는 것은 중요하다. 수평형 지중 열교환기용 되메움재의 열전도도를 파악하기 위해 비정상 열선법을 적용한 QTM-500을 사용하여 포화도에 따른 천연규사-물-공기 혼합물의 열전도도를 측정하였다. 측정된 열전도도를 개별요소법(Discrete Element Mothod)에 근거한 2차원 수치해석 프로그램인 PFC2D(Particle Flow Code in 2 Dimension)를 이용하여 비교 분석하였다. 수치해석에서는 혼합물의 건조밀도를 일정하게 유지한 상태에서 포화도에 따라 가상의 물 입자 개수를 변화시켰다. 개별요소법을 이용한 열전달 수치해석에서는 입자의 접촉을 통해 발생한 thermal pipe에 의해 열전달이 이루어진다. 이러한 thermal pipe의 열전도도는 접촉된 두 입자의 열전도도와 접촉면의 평균 열전도도를 고려하여 적용하였다.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.16 no.2
• /
• pp.149-154
• /
• 2017

Numerous studies of the dynamic behavior of powders have been performed by Discrete Element Method (DEM). The behavior of powders can be analyzed using the DEM assuming that the powder is composed of spherical particles. Moreover, the assumption of spherical particle reduces the computing time significantly. However, the biggest problem with this assumption is the real shape of the particles. Some types of particles, such as calcium carbonate and colloidal copper, are needle shaped. Thus, analysis based on spherical particles can produce errors because of the incorrect assumption. In this research, we developed a new model to simulate needle-shaped particles using the DEM. In the model, a series of particles are connected and regarded as a rod. There is no relative motion among the particles. Thus, the behavior of the rod is rigid motion. To validate the developed model, we carried out the drop-and-bounce test with different initial angles. The results showed negligible error of less than 2%.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.17 no.3
• /
• pp.22-27
• /
• 2018

One problem of the Discrete Element Method is the assumption of a spherical particle shape, which reduces the computing time but significantly limits the application of the DEM to analysis. This limitation can be overcome by a recently developed rigid-rod model. However, the rigid-rod model has an essential problem related with friction: it always contains friction error because of the bumpy surface. To overcome this issue, we suggest a superposed rigid-rod model in this paper. The superposed rigid-rod model is notably consistent with the theoretical value in terms of the velocity and angular velocity after the collision. The estimated error is negligible(less than 2%). Then, the developed model is applied to hopper discharging. The developed model shows no problem in the discharging flow from the hopper.

• Computers and Concrete
• /
• v.5 no.4
• /
• pp.343-358
• /
• 2008

This article focuses on concrete structures submitted to impact loading and is aimed at predicting local damage in the vicinity of an impact zone as well as the global response of the structure. The Discrete Element Method (DEM) seems particularly well suited in this context for modeling fractures. An identification process of DEM material parameters from macroscopic data (Young's modulus, compressive and tensile strength, fracture energy, etc.) will first be presented for the purpose of enhancing reproducibility and reliability of the simulation results with DE samples of various sizes. The modeling of a large structure by means of DEM may lead to prohibitive computation times. A refined discretization becomes required in the vicinity of the impact, while the structure may be modeled using a coarse FE mesh further from the impact area, where the material behaves elastically. A coupled discrete-finite element approach is thus proposed: the impact zone is modeled by means of DE and elastic FE are used on the rest of the structure. The proposed approach is then applied to a rock impact on a concrete slab in order to validate the coupled method and compare computation times.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.15 no.4
• /
• pp.1-7
• /
• 2016

The pick cutter, which directly contacts and crushes the rock, is the expendable part of a roadheader. The arrangement and angle of attachment of the pick cutter are important factors that determine excavator performance. It is necessary to numerically calculate the contact between the pick cutter and rock. The rock is defined as a set of particles using the discrete element method. The parallel bond model is used to define the bonds between particles. The properties of granite that are measured by the uniaxial compressive test are applied to the numerical rock model. The pick cutter is defined by the polygon elements. The linear cutting simulation is considered to simulate the contact between the pick cutter and rock. The results of the simulation show the rock breaking due to contact with the pick cutter.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.2
• /
• pp.87-95
• /
• 2003

In this paper, the discrete element method was applied to evaluate the stability of composite cover and lining system of landfill. This method is capable of estimating the distribution of tensile force and shear stress mobilized in each liner component and its interfaces, based on a relationship of force and displacement. It was assumed that the cover soil and geomembrane were comprised of slices connected with elastoplastic Winkler springs and tensile spring respectively. Parametric study using this method was performed and compared with other techniques based on limit equilibrium method fur the example analysis.

• Geomechanics and Engineering
• /
• v.13 no.6
• /
• pp.977-995
• /
• 2017

It is well accepted that rock failure mechanism influence the cutting efficiency and determination of optimum cutting parameters. In this paper, an attempt was made to research the factors that affect the failure mechanism based on discrete element method (DEM). The influences of cutting depth, hydrostatic pressure, cutting velocity, back rake angle and joint set on failure mechanism in rock-cutting are researched by PFC2D. The results show that: the ductile failure occurs at shallow cutting depths, the brittle failure occurs as the depth of cut increases beyond a threshold value. The mean cutting forces have a linear related to the cutting depth if the cutting action is dominated by the ductile mode, however, the mean cutting forces are deviate from the linear relationship while the cutting action is dominated by the brittle mode. The failure mechanism changes from brittle mode with larger chips under atmospheric conditions, to ductile mode with crushed chips under hydrostatic conditions. As the cutting velocity increases, a grow number of micro-cracks are initiated around the cutter and the volume of the chipped fragmentation is decreasing correspondingly. The crack initiates and propagates parallel to the free surface with a smaller rake angle, but with the rake angle increases, the direction of crack initiation and propagation is changed to towards the intact rock. The existence of joint set have significant influence on crack initiation and propagation, it makes the crack prone to propagate along the joint.

• 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.