• Geomechanics and Engineering
• /
• 제15권4호
• /
• pp.1017-1027
• /
• 2018

To explore the influence of coal thickness on the mechanical behavior and the failure characteristics of rock-coal-rock (RCR) mass, the experimental investigation of uniaxial compressive tests was conducted first and then a systematic numerical simulation by particle flow code (PFC2D) was performed to deeply analyze the failure mechanical behavior of RCR specimens with different coal thicknesses in conventional compression tests. The overall elastic modulus and peak stress of RCR specimens lie between the rock and the coal. Inter-particle properties were calibrated to match the physical sample strength and the stiffness response. Numerical simulation results show that the deformation and strength behaviors of RCR specimens depend not only on the coal thickness, but also on the confining pressure. Under low confining pressures, the overall failure mechanism of RCR specimen is the serious damage of coal section when the coal thickness is smaller than 30 mm, but it is shear failure of coal section when the coal thickness is larger than 30 mm. Whereas under high confining pressures, obvious shear bands exist in both the coal section and the rock section when the coal thickness is larger than 30 mm, but when the coal thickness is smaller than 30mm, the failure mechanism is serious damage of coal section and shear failure of rock section.

• 해양환경안전학회지
• /
• 제24권5호
• /
• pp.611-618
• /
• 2018

본 연구에서는 이수의 유동 특성을 분석하기 위한 기초 연구로서 상용 코드인 ANSYS CFX 14.5를 이용하여 고체-액체 2상 유동에 대한 수치해석적 연구를 수행하였다. 고체-액체 2상 유동 현상을 모사하기 위해서 균질류 모델과 분리류 모델을 사용하였다. 분리류 모델에서는 Gidaspow의 항력모델을 적용하였으며, 고체 입자에 운동 이론 모델을 적용하였다. 기존의 실험 결과를 기반으로 본 연구에서 사용한 수치해석 모델의 유효성을 검토하였으며, 수치해석은 직경 54.9 mm, 길이 3 m의 수평관에서 체적 분율 0.1~0.5, 속도 1~5 m/s 범위에서 수행되었다. 그리고 압력강하와 고체 입자의 체적 분율 분포를 확인하였으며, 압력강하는 균질류 모델과 분리류 모델이 각각 MAE 17.04 %, 8.98 % 이내에서 실험결과를 잘 예측하였다. 관의 하부에서 높은 체적 분율이 나타나며, 상부로 갈수록 체적 분율은 감소하였다. 그리고 속도가 증가할수록 높이 변화에 따른 체적 분율 분포의 변화는 감소하였으며, 수치해석 결과는 이러한 유동 특성을 잘 예측하였다.

• 천문학회지
• /
• 제35권4호
• /
• pp.159-174
• /
• 2002

Cosmological hydrodynamic simulations of large scale structure in the universe have shown that accretion shocks and merger shocks form due to flow motions associated with the gravitational collapse of nonlinear structures. Estimated speed and curvature radius of these shocks could be as large as a few 1000 km/s and several Mpc, respectively. According to the diffusive shock acceleration theory, populations of cosmic-ray particles can be injected and accelerated to very high energy by astrophysical shocks in tenuous plasmas. In order to explore the cosmic ray acceleration at the cosmic shocks, we have performed nonlinear numerical simulations of cosmic ray (CR) modified shocks with the newly developed CRASH (Cosmic Ray Amr SHock) numerical code. We adopted the Bohm diffusion model for CRs, based on the hypothesis that strong Alfven waves are self-generated by streaming CRs. The shock formation simulation includes a plasma-physics-based 'injection' model that transfers a small proportion of the thermal proton flux through the shock into low energy CRs for acceleration there. We found that, for strong accretion shocks, CRs can absorb most of shock kinetic energy and the accretion shock speed is reduced up to $20\%$, compared to pure gas dynamic shocks. For merger shocks with small Mach numbers, however, the energy transfer to CRs is only about $10-20\%$ with an associated CR particle fraction of $10^{-3}$. Nonlinear feedback due to the CR pressure is insignificant in the latter shocks. Although detailed results depend on models for the particle diffusion and injection, these calculations show that cosmic shocks in large scale structure could provide acceleration sites of extragalactic cosmic rays of the highest energy.

• Structural Engineering and Mechanics
• /
• 제69권1호
• /
• pp.11-20
• /
• 2019

The biaxial failure mechanism of transversally bedding concrete layers was numerically simulated using a sophisticated two-dimensional discrete element method (DEM) implemented in the particle flow code (PFC2D). This numerical modelling code was first calibrated by uniaxial compression and Brazilian testing results to ensure the conformity of the simulated numerical model's response. Secondly, 21 rectangular models with dimension of $54mm{\times}108mm$ were built. Each model contains two transversely bedding layers. The first bedding layer has low mechanical properties, less than mechanical properties of intact material, and second bedding layer has high mechanical properties, more than mechanical properties of intact material. The angle of first bedding layer, with weak mechanical properties, related to loading direction was $0^{\circ}$, $15^{\circ}$, $30^{\circ}$, $45^{\circ}$, $60^{\circ}$, $75^{\circ}$ and $90^{\circ}$ while the angle of second layer, with high mechanical properties, related to loading direction was $90^{\circ}$, $105^{\circ}$, $120^{\circ}$, $135^{\circ}$, $150^{\circ}$, $160^{\circ}$ and $180^{\circ}$. Is to be note that the angle between bedding layer was $90^{\circ}$ in all bedding configurations. Also, three different pairs of the thickness were chosen in models, i.e., 5 mm/10 mm, 10 mm/10 mm and 20 mm/10 mm. The result shows that in all configurations, shear cracks develop between the weaker bedding layers. Shear cracks angel related to normal load change from $0^{\circ}$ to $90^{\circ}$ with increment of $15^{\circ}$. Numbers of shear cracks are constant by increasing the bedding thickness. It's to be noted that in some configuration, tensile cracks develop through the intact area of material model. There is not any failure in direction of bedding plane interface with higher strength.

• 한국터널지하공간학회 논문집
• /
• 제12권1호
• /
• pp.31-41
• /
• 2010

본 연구에서는 실험 및 수치해석적인 접근방법을 통하여 대심도 과지압 구간에서 발생하는 암반의 스폴링(spalling) 및 팽창모드에 대한 모델링 기법을 연구하였다. 이에 대한 실험적 접근 방법으로서 축 방향 압축을 받는 직사각형 인공암석보(beam)에 4점 휨 시험을 결합한 축방향 압축 휨 시험을 수행하여 대심도 지하공간의 응력모드와 유사한 조건 하에서의 암석의 균열 팽창 및 스폴링 과정을 고찰하였다. 또한, 수치해석적 접근방법으로서 기존의 연속체 해석으로는 모사하기 힘든 암석의 균열과정 및 팽창특성을 개별 입자해석 프로그램인 PFC2D를 이용하여 모델링 하였다. 본 연구 결과 휨 실험에서 구한 팽창시점은 스폴링에 요구되는 응력수준을 평가하는데 중요한 지표가 됨을 알 수 있었으며, 또한 수치해석 결과도 유사한 결과를 모사할 수 있음을 보여주었다.

• 한국신재생에너지학회:학술대회논문집
• /
• 한국신재생에너지학회 2010년도 추계학술대회 초록집
• /
• 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의 열전도도는 접촉된 두 입자의 열전도도와 접촉면의 평균 열전도도를 고려하여 적용하였다.

• 한국해양환경ㆍ에너지학회지
• /
• 제16권4호
• /
• pp.255-262
• /
• 2013

본 연구에서는 연속체 이론을 배경으로 하며 일반적으로 많이 사용되는 Navire-Stokes방정식이 아닌 입자의 확률분포를 배경으로 하는 Boltzmann 방정식을 이용하여 자유수면을 포함하는 유동을 해석하는 전산시뮬레이션 코드를 개발하였다. 댐 붕괴시뮬레이션에 적용하여 코드의 검증을 수행하였으며, 기존의 실험 및 계산결과와 비교함과 동시에 격자볼츠만 시뮬레이션의 특성을 분석하였다. 공학적 응용을 위해서 임의 형상의 물체가 존재시에 자유수면 시뮬레이션도 수행하였다.

• Structural Engineering and Mechanics
• /
• 제4권3호
• /
• pp.227-242
• /
• 1996

An important class of problems in the field of geotechnical engineering may be analyzed with the aid of a simple integro-differential equation. Behavior of "rigid" piles(say concrete piles), "deformable" piles(say gravel piles), pile groups, pile-raft foundations, heavily reinforced earth, flow within circular silos and down drag on cylindrical structures (for example the crusher unit of a mineral processing complex) are the type of situations that can be handled by this type of equation. The equation under consideration has the form; $$\frac{{\partial}w(r,\;z)}{{\partial}z}+f(z){\int}^z_0g({\xi})(\frac{{\partial}^2w(r,\;{\xi})}{{\partial}r^2}+\frac{1}{r}\frac{{\partial}w(r,\;{\xi})}{{\partial}r})d{\xi}+h(r,\;z)=0$$ where w(r, z) is the vertical displacement of a soil particle expressed as a function of the polar cylindrical space coordinates (r, z) and the symbols f, g and h represent soil properties and the loading conditions. The merit of the analysis is its simplicity (both in concept and in application) and the ease with which it can be expressed in a computer code. In the present paper the analysis is applied to investigate the behavior of a single rigid pile to bedrock. The emphasis, however, is placed on developing the equation, the numerical techique used in its evaluation and validation of the technique, hereafter called the ID technique, against a formal program, CRISP, which uses the FEM.

• 한국전산유체공학회지
• /
• 제15권3호
• /
• pp.39-45
• /
• 2010

The turbulent flows in a tunnel mock-up($10L{\times}0.5W{\times}0.25H$ m3 : scale reduction 1/20) with rectangular cross section were investigated. The instantaneous velocity fields of Re = 49,029, 89,571 were measured by the 2-D PIV system which is consisted of double pulsed Nd:Yag laser and the tracer particles in the straight-duct mock-up where the flows were fully developed. The mean velocity profiles were taken from the ensemble averages of 1,000 instantaneous velocity fields. Simultaneously, numerical simulations(RANS) were performed to compare with experimental data using STREAM code. Non-linear eddy viscosity model (NLEVM : Abe-Jang-Leschziner Eddy Viscosity Model) was employed to resolve the turbulent flows in the duct. The calculated mean velocity profiles were well compared with PIV results. In the log-law profiles, the experimental data were in good agreement with numerical simulations all the way to the wake region except the viscous sub-layer (near wall region).

• Computers and Concrete
• /
• 제20권1호
• /
• pp.39-47
• /
• 2017

In this paper, mode I fracture toughness of rock was determined by direct and indirect methods using Particle Flow Code simulation. Direct methods are compaction tension (CT) test and hollow centre cracked quadratic sample (HCCQS). Indirect methods are notched Brazilian disk (NBD) specimen, the semi-circular bend (SCB) specimen, hollow centre cracked disc (HCCD), the single edge-notched round bar in bending (SENRBB) specimen and edge notched disk (END). It was determined that which one of indirect fracture toughness values is close to direct one. For this purpose, initially calibration of PFC was undertaken with respect to data obtained from Brazilian laboratory tests to ensure the conformity of the simulated numerical models response. Furthermore, the simulated models in five introduced indirect tests were cross checked with the results from direct tests. By using numerical testing, the failure process was visually observed. Discrete element simulations demonstrated that the macro fractures in models are caused by microscopic tensile breakages on large numbers of bonded discs. Mode I fracture toughness of rock in direct test was less than other tests results. Fracture toughness resulted from semi-circular bend specimen test was close to direct test results. Therefore semi-circular bend specimen can be a proper test for determination of Mode I fracture toughness of rock in absence of direct test.