• Title/Summary/Keyword: Discrete-Element-Method

Search Result 508, Processing Time 0.023 seconds

Finite element modeling of high Deborah number planar contraction flows with rational function interpolation of the Leonov model

  • Youngdon Kwon;Kim, See-Jo;Kim, Seki
    • Korea-Australia Rheology Journal
    • /
    • v.15 no.3
    • /
    • pp.131-150
    • /
    • 2003
  • A new numerical algorithm of finite element methods is presented to solve high Deborah number flow problems with geometric singularities. The steady inertialess planar 4 : 1 contraction flow is chosen for its test. As a viscoelastic constitutive equation, we have applied the globally stable (dissipative and Hadamard stable) Leonov model that can also properly accommodate important nonlinear viscoelastic phenomena. The streamline upwinding method with discrete elastic-viscous stress splitting is incorporated. New interpolation functions classified as rational interpolation, an alternative formalism to enhance numerical convergence at high Deborah number, are implemented not for the whole set of finite elements but for a few elements attached to the entrance comer, where stress singularity seems to exist. The rational interpolation scheme contains one arbitrary parameter b that controls the singular behavior of the rational functions, and its value is specified to yield the best stabilization effect. The new interpolation method raises the limit of Deborah number by 2∼5 times. Therefore on average, we can obtain convergent solution up to the Deborah number of 200 for which the comer vortex size reaches 1.6 times of the half width of the upstream reservoir. Examining spatial violation of the positive definiteness of the elastic strain tensor, we conjecture that the stabilization effect results from the peculiar behavior of rational functions identified as steep gradient on one domain boundary and linear slope on the other. Whereas the rational interpolation of both elastic strain and velocity distorts solutions significantly, it is shown that the variation of solutions incurred by rational interpolation only of the elastic strain is almost negligible. It is also verified that the rational interpolation deteriorates speed of convergence with respect to mesh refinement.

The Multidirectional Random Wave Diffraction in a Partial-Reflecting Harbor due to a Submarine Pit (Pit에 의한 부분반사율을 갖는 항내에서의 다방향 불규칙 파랑회절에 관한 연구)

  • Kim, Sung-Duk;Lee, Hong-Sik
    • Journal of Korean Society of Coastal and Ocean Engineers
    • /
    • v.20 no.3
    • /
    • pp.291-299
    • /
    • 2008
  • The present study is to estimate the effects of diffracted wave fields in a harbor and around harborentrance due to random waves, when a navigation channel is dredged in the vicinity of the a harbor entrance. The cross sections of harbor boundary are considered to be partial or full reflection in this study. The numerical simulation has been performed by the boundary element method, which is to discrete segments of pit- and harbor- boundary with the algorism of auto generated elements. The incident wave conditions are specified using discretized forms of the Mitsuyasu's frequency spectrum and directional function. The results of the present numerical simulation agreed well with those of the published experimental data. It is shown that the ratios of wave height reduction are about 20% for the case of fully reflecting boundary, and 10% for the case of partially reflecting boundary by the effect of placing a pit, respectively.

Estimation of groundwater inflow into an underground oil storage facility in granite

  • Wang, Zhechao;Kwon, Sangki;Qiao, Liping;Bi, Liping;Yu, Liyuan
    • Geomechanics and Engineering
    • /
    • v.12 no.6
    • /
    • pp.1003-1020
    • /
    • 2017
  • Estimation of groundwater inflow into underground opening is of critical importance for the design and construction of underground structures. Groundwater inflow into a pilot underground storage facility in China was estimated using analytical equations, numerical modeling and field measurement. The applicability of analytical and numerical methods was examined by comparing the estimated and measured results. Field geological investigation indicated that in local scale the high groundwater inflows are associated with the appearance of open joints, fractured zone or dykes induced by shear and/or tensile tectonic stresses. It was found that 8 groundwater inflow spots with high inflow rates account for about 82% of the total rate for the 9 caverns. On the prediction of the magnitude of groundwater inflow rate, it was found that could both (Finite Element Method) FEM and (Discrete Element Method) DEM perform better than analytical equations, due to the fact that in analytical equations simplified assumptions were adopted. However, on the prediction of the spatial distribution estimation of groundwater inflow, both analytical and numerical methods failed to predict at the present state. Nevertheless, numerical simulations would prevail over analytical methods to predict the distribution if more details in the simulations were taken into consideration.

Variable Radix-Two Multibit Coding and Its VLSI Implementation of DCT/IDCT (가변길이 다중비트 코딩을 이용한 DCT/IDCT의 설계)

  • 김대원;최준림
    • Journal of the Institute of Electronics Engineers of Korea SD
    • /
    • v.39 no.12
    • /
    • pp.1062-1070
    • /
    • 2002
  • In this paper, variable radix-two multibit coding algorithm is presented and applied in the implementation of discrete cosine transform(DCT) and inverse discrete cosine transform(IDCT). Variable radix-two multibit coding means the 2k SD (signed digit) representation of overlapped multibit scanning with variable shift method. SD represented by 2k generates partial products, which can be easily implemented with shifters and adders. This algorithm is most powerful for the hardware implementation of DCT/IDCT with constant coefficient matrix multiplication. This paper introduces the suggested algorithm, it's proof and the implementation of DCT/IDCT The implemented IDCT chip with 8 PEs(Processing Elements) and one transpose memory runs at a tate of 400 Mpixels/sec at 54MHz frequency for high speed parallel signal processing, and it's verified in HDTV and MPEG decoder.

Image Retrieval using Variable Block Size DCT (가변 블록 DCT를 이용한 영상 검색 기법)

  • 김동우;서은주;윤태승;안재형
    • Journal of Korea Multimedia Society
    • /
    • v.4 no.5
    • /
    • pp.423-429
    • /
    • 2001
  • In this paper, we propose the improved method for retrieving images with DC element of DCT that is used in image compression such as JPEG/MPEG. The existing method retrieves images with DC of fixed block size DCT. In this method, the increase in the block size results in faster retrieving speed, but it lessens the accuracy. The decrease in the block size improves the accuracy, however, it degrades the retrieving speed. In order to solve this problem, the proposed method utilizes the variable block size DCT. This method first determines the existence of object regions within each block, and then creates an image region table. Based on this table, it determines the size of each block, following a simple rule; decrease the block size in the object regions, and increase the block size in the background regions. The proposed method using variable block size DCT improves about 15% in terms of the accuracy. Additionally, when there rarely exist images of same pattern, it is able to retrieve faster only by comparing the image region patterns.

  • PDF

A Study of the Tower Crane Hoisting Time Estimation Simulation Model with Climate Element for the High-Rise Building Construction (기후요소를 고려한 초고층 건설공사의 타워크레인 양중시간 예측 시뮬레이션 모델)

  • Yang, Kanghyeok;Lee, Hyun-Soo;Park, Moonseo;Jung, Minhyeok;Hwang, Sungjoo
    • Korean Journal of Construction Engineering and Management
    • /
    • v.14 no.2
    • /
    • pp.96-107
    • /
    • 2013
  • Tower crane hoisting plan is one of the key element for the success of entire High-Rise Building construction. Hoisting time is the basic factor to appropriate hoisting plan which need to the hoisting load estimate and tower crane selection. With this reason, accurate hoisting time is needed to the proper hoisting plan. The current hoisting time estimation for High-Rise Building focus on the hoisting cycle time estimation with historical data. However, this method underestimated the external influences like environmental factor. Thus, this paper aims to develop the hoisting time estimation model with discrete event simulation which include the wind influences with certain height. According to the simulation result, the hoisting time which applied wind influence is increasing with height growth. Because of the high speed wind, the upper area of building has more operation delay time than the mechanical operation time. Seoul, the research area, has the most fastest wind speed on April and the least on October. Due to these differences of wind speed, the hoisting time is estimated with significant differences between April and October. This hosting time estimation model would be used for estimating the influence of wind. Moreover, this could apply to make the realistic hoisting plan.

Defining the hydraulic excavation damaged zone considering hydraulic aperture change (수리적 간극변화를 고려한 수리적 굴착손상영역의 정의에 관한 연구)

  • Park, Jong-Sung;Ryu, Chang-Ha;Lee, Chung-In;Ryu, Dong-Woo
    • Journal of Korean Tunnelling and Underground Space Association
    • /
    • v.9 no.2
    • /
    • pp.133-141
    • /
    • 2007
  • The excavation damaged zone (EDZ) is an area around an excavation where in situ rock mass properties, stress condition, displacement, groundwater flow conditions have been altered due to the processes induced by the excavation. Various studies have been carried out on EDZ, but most studies have focused on the mechanical bahavior of EDZ by in situ experiment. Even though the EDZ could potentially form a high permeable pathway of groundwater flow, only a few studies were performed on the analysis of groundwater flow in EDZ. In this study, the 'hydraulic EDZ' was defined as the rock zone adjacent to the excavation where the hydraulic aperture has been changed due to the excavation by using H-M coupling analysis. Fundamental principles of distinct element method (DEM) were used in the analysis. In the same groundwater level, the behavior of hydraulic aperture near the cavern was analyzed for different stress ratios, initial apertures, fracture angles and fracture spacings by using a two-dimensional DEM program. We evaluate the excavation induced hydraulic aperture change. Using the results of the study, hydraulic EDZ was defined as an elliptical shape model perpendicular to the joint.

  • PDF

Generation Method of the Rectangular Grid Information for Finite Difference Model (유한차분모형을 위한 직사각형 격자정보 생성기법)

  • 정신택;조범준;김정대
    • Journal of Korean Society of Coastal and Ocean Engineers
    • /
    • v.15 no.3
    • /
    • pp.190-195
    • /
    • 2003
  • For many coastal problems, such as wave transformation, tidal circulation, sediment transports and diffusion phenomena, we resort to numerical techniques. The representative numerical techniques are the method of finite differences and finite elements. The approximate algebraic equations, referred to as finite difference equations(FDEs), are subsequently solved at discrete grid points within the domain of interests. Therefore, a set of grid points within the domain, as well as the boundaries of the domain, must be specified. The generation of grids for FDEs, with uniform spacing, is very simple compared to that of finite elements. However, within a very complex domain, there are few grid generation tools we can use conveniently. Unfortunately, most of the commercial grid generation programs are developed only for finite element method. In this paper, grid generation method using digitizer, with uniform rectangular spacing, are introduced in detail. Didger and Surfer programs by Golden Software are necessary to produce comparatively accurate and simple depth data.

Effect of the Earth Pressure Coefficient on the Support System in Jointed Rock Mass

  • Son, Moorak;Adedokun, Solomon;Hwang, Youngcheol
    • Journal of the Korean GEO-environmental Society
    • /
    • v.16 no.2
    • /
    • pp.33-43
    • /
    • 2015
  • This paper investigated the magnitude and distribution of earth pressure on the support system in jointed rock mass by considering different earth pressure coefficients, rock types and joint inclination angles. The study mainly focused on the effect of the earth pressure coefficients on the earth pressure. Based on a physical model test (Son & Park, 2014), extended studies were conducted considering rock-structure interactions based on the discrete element method, which can consider the joints characteristics of rock mass. The results showed that the earth pressure was highly influenced by the earth pressure coefficients as well as the rock type and joint inclination angles. The effects of the earth pressure coefficients increased when the rock suffered more weathering and has no joint slide. The test results were also compared with Peck's earth pressure for soil ground, and clearly showed that the earth pressure in jointed rock mass can be greatly different from that in soil ground. This study indicated the earth pressure coefficients considering the rock types and joint inclination angles are important parameters influencing the magnitude and distribution of earth pressure, which should be considered when designing the support systems in jointed rock mass.

A Study on Numerical Modeling of Dynamic CPT using Particle Flow Code (입자결합모델을 이용한 동적콘관입시험(DCPT)의 수치해석 모델링에 관한 연구)

  • You, Kwang Ho;Lee, Chang Su;Choi, Jun Seong
    • International Journal of Highway Engineering
    • /
    • v.16 no.2
    • /
    • pp.43-52
    • /
    • 2014
  • PURPOSES : To solve problems in current compaction control DCPT(Dynamic Cone Penetrometer Test), highly correlated with various testing methods, simple, and economic is being applied. However, it、s hard to utilize DCPT results due to the few numerical analyses for DCPT have been performed and the lack of data accumulation. Therefore, this study tried to verify the validation of numerical modeling for DCPT by comparing and analyzing the results of numerical analyses with field tests. METHODS: The ground elastic modulus and PR(Penetration Rate) value were estimated by using PFC(Particle Flow Code) 3D program based on the discrete element method. Those values were compared and analyzed with the result of field tests. Also, back analysis was conducted to describe ground elastic modulus of field tests. RESULTS : Relative errors of PR value between the numerical analyses and field tests were calculated to be comparatively low. Also, the relationship between elastic modulus and PR value turned out to be similar. CONCLUSIONS : Numerical modeling of DCPT is considered to be suitable for describing field tests by carrying out numerical analysis using PFC 3D program.