• 한국방사성폐기물학회:학술대회논문집
• /
• 한국방사성폐기물학회 2009년도 학술논문요약집
• /
• pp.187-187
• /
• 2009

Opening of fractures induced by shear dilation or normal deformation can be a significant source of fracture permeability change in fractured rock, which is important for the performance assessment of geological repositories for spent nuclear fuel. As the repository generates heat and later cools the fluid-carrying ability of the rocks becomes a dynamic variable during the lifespan of the repository. Heating causes expansion of the rock close to the repository and, at the same time, contraction close to the surface. During the cooling phase of the repository, the opposite takes place. Heating and cooling together with the, virgin stress can induce shear dilation of fractures and deformation zones and change the flow field around the repository. The objectives of this work are to examine the contribution of thermal stress to the shear slip of fracture in mid- and far-field around a KBS-3 type of repository and to investigate the effect of evolution of stress on the rock mass permeability. In the first part of this study, zones of fracture shear slip were examined by conducting a three-dimensional, thermo-mechanical analysis of a spent fuel repository model in the size of 2 km $\times$ 2 km $\times$ 800 m. Stress evolutions of importance for fracture shear slip are: (1) comparatively high horizontal compressive thermal stress at the repository level, (2) generation of vertical tensile thermal stress right above the repository, (3) horizontal tensile stress near the surface, which can induce tensile failure, and generation of shear stresses at the comers of the repository. In the second part of the study, fracture data from Forsmark, Sweden is used to establish fracture network models (DFN). Stress paths obtained from the thermo-mechanical analysis were used as boundary conditions in DFN-DEM (Discrete Element Method) analysis of six DFN models at the repository level. Increases of permeability up to a factor of four were observed during thermal loading history and shear dilation of fractures was not recovered after cooling of the repository. An understanding of the stress path and potential areas of slip induced shear dilation and related permeability changes during the lifetime of a repository for spent nuclear fuel is of utmost importance for analysing long-term safety. The result of this study will assist in identifying critical areas around a repository where fracture shear slip is likely to develop. The presentation also includes a brief introduction to the ongoing site investigation on two candidate sites for geological repository in Sweden.

• Structural Engineering and Mechanics
• /
• 제16권3호
• /
• pp.259-274
• /
• 2003

In this paper, a novel numerical solution technique, the differential cubature method is employed to study the buckling problems of thick plates with arbitrary quadrilateral planforms and non-uniform boundary constraints based on the first order shear deformation theory. By using this method, the governing differential equations at each discrete point are transformed into sets of linear homogeneous algebraic equations. Boundary conditions are implemented through discrete grid points by constraining displacements, bending moments and rotations of the plate. Detailed formulation and implementation of this method are presented. The buckling parameters are calculated through solving a standard eigenvalue problem by subspace iterative method. Convergence and comparison studies are carried out to verify the reliability and accuracy of the numerical solutions. The applicability, efficiency, and simplicity of the present method are demonstrated through solving several sample plate buckling problems with various mixed boundary constraints. It is shown that the differential cubature method yields comparable numerical solutions with 2.77-times less degrees of freedom than the differential quadrature element method and 2-times less degrees of freedom than the energy method. Due to the lack of published solutions for buckling of thick rectangular plates with mixed edge conditions, the present solutions may serve as benchmark values for further studies in the future.

• 한국도로학회논문집
• /
• 제18권3호
• /
• pp.21-32
• /
• 2016

PURPOSES : In this paper, the analytical solutions suggested to simulate the behavior of rheological fluids were rigorously re-derived and investigated for fixed conditions to evaluate the applicability for the solutions on a mini-cone slump test of cement paste. The selected solutions with proper boundary conditions can be used as reference solutions to evaluate the performance of numerical simulation approaches, such as the discrete element method. METHODS : The slump, height, and spread radius for the given boundary and yield stress conditions that are determined by five different analytical solutions are compared. RESULTS : The analytical solution based on fluid mechanics for pure shear flow shows similar results to that for intermediate flow at low yield stresses. The fluid mechanics-based analytical solution resulted in a very similar trend to the geometry-based analytical solution. However, it showed a higher slump at high yield stress and lower slump at low yield stress ranges than the geometry-based analytical model. The analytical solution based on the mini-cone geometry was not significantly affected by the yield criteria, such as von Mises and Tresca. CONCLUSIONS : Even though differences among the analytical solutions in terms of slump and spread radius existed, the difference can be considered insignificant when the solutions were used as reference to evaluate the appropriateness of numerical approaches, such as the discrete element method.

• Geomechanics and Engineering
• /
• 제7권6호
• /
• pp.665-678
• /
• 2014

In discrete element modeling, 2D software has been widely used in order to gain further insights into the fundamental mechanisms with less computational time. The porosities used in 2D DEM studies should be determined with appropriate approaches based on 3D laboratory porosities. This paper summarizes the main approaches for converting porosities from 3D to 2D for DEM studies and theoretical evaluations show that none of the current approaches can be widely used in dealing with soil mechanical problems. Therefore, a parabolic equation and a criterion have been suggested for the determination of 2D porosities in this paper. Moreover, a case study has been used to validate that the 2D porosity obtained from the above suggestion to be rational with both the realistic contact force distribution in the specimen and the good agreement of the DEM simulation results of direct shear tests with the corresponding experimental data. Therefore, the parabolic equation and the criterion are suggested for the determination of 2D porosities in a wide range of polydisperse particle systems, especially in dealing with soil mechanical problems.

• 한국해양공학회지
• /
• 제18권5호
• /
• pp.15-21
• /
• 2004

An edge tone is the discrete tone or narrow-band sound produced by an oscillating free shear layer, impinging on a rigid surface. In this paper, we present a 2-D edge tone to predict the frequency characteristics of the discrete oscillations of a jet-edge feedback cycle, using the finite difference lattice Boltzmann method (FDLBM). We use a modified version of the lattice BGK compressible fluid model, adding an additional term and allowing for longer time increments, compared to a conventional FDLBM, and also use a boundary fitted coordinates system. The jet is chosen long enough in order to guarantee the parabolic velocity profile of the jet at the outlet, and the edge consists of a wedge with an angle of ${\alpha}$ = 23. At a stand-off distance, the edge is inserted along the centerline of the jet, and a sinuous instability wave, with real frequency, is assumed to be created in the vicinity of the nozzle and propagates towards the downstream. We have succeeded in capturing very small pressure fluctuations, resulting from periodical oscillations of a jet around the edge. The pressure fluctuations propagate with the speed of sound. Its interaction with the wedge produces an non-rotational feedback field, which, near the nozzle exit, is a periodic transverse flow, producing the singularities at the nozzle lips.

• 대한기계학회논문집B
• /
• 제21권4호
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• Macromolecular Research
• /
• 제13권3호
• /
• pp.243-252
• /
• 2005

The effect of the vinyl acetate (VA) content on the thermal, viscoelastic, rheological, morphological and mechanical behaviors in various blends of ethylene-vinyl acetate (EVA)/ethylene-$\alpha$-olefin copolymers was investigated using 28, 22 and $15 mol\%$ of VA in EVA. In the DSC melting and crystallization thermograms of all of the EVA systems blended with ethylene-$\alpha$-olefin copolymers, discrete peaks were observed which were related to the constituents. In the dynamic mechanical thermal analysis, the storage modulus increased with increasing content of ethylene-$\alpha$-olefin copolymers. In addition, the transition regions relating to the tan bpeaks varied with the VA content. The crossover point between G' and G" varied depending on the VA contents, and shear-thinning was more prominent in the EVA/EtBC system. In the SEM investigation, a discrete phase morphology was observed in both the EVA/EtBC and EVA/EtOC blends, but the contrast improved with decreasing VA content. However, the tensile strength and modulus improved, but the elongation at break reduced with decreasing VA content, implying that the ethylene-$\alpha$-olefin copolymers play the role of reinforcing materials. Thus, the EVA and ethylene-$\alpha$-olefin components in the copolymers are immiscible in the molten and solid states, but are nevertheless mechanically compatible.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1915-1920
• /
• 2004

An edge tone is the discrete tone or narrow-band sound produced by an oscillating free shear layer impinging on a rigid surface. In this paper we present a two-dimensional edge tone to predict the frequency characteristics of the discrete oscillations of a jet-edge feedback cycle by the finite difference lattice Boltzmann method. We use a new lattice BGK compressible fluid model that has an additional term and allow larger time increment comparing a conventional FDLB model, and also use a boundary fitted coordinates. The jet is chosen long enough in order to guarantee the parabolic velocity profile of the jet at the outlet, and the edge consists of a wedge with an angle of ${\alpha}=23^{\circ}$ . At a stand-off distance ${\omega}$ , the edge is inserted along the centreline of the jet, and a sinuous instability wave with real frequency f is assumed to be created in the vicinity of the nozzle and to propagate towards the downstream. We have succeeded in capturing very small pressure fluctuations result from periodically oscillation of jet around the edge. That pressure fluctuations propagate with the sound speed. Its interaction with the wedge produces an irrotational feedback field which, near the nozzle exit, is a periodic transverse flow producing the singularities at the nozzle lips.

• Journal of Ship and Ocean Technology
• /
• 제12권1호
• /
• pp.1-15
• /
• 2008

An edgetone is the discrete tone or narrow-band sound produced by an oscillating free shear layer, impinging on a rigid surface. In this paper, 2-dimensional edgetone to predict the frequency characteristics of the discrete oscillations of a jet-edge feedback cycle is presented using lattice Boltmznan model with 21 bits, which is introduced a flexible specific heat ratio y to simulate diatomic gases like air. The blown jet is given a parabolic inflow profile for the velocity, and the edges consist of wedges with angle 20 degree (for symmetric wedge) and 23 degree (for inclined wedge), respectively. At a stand-off distance w, the edge is inserted along the centerline of the jet, and a sinuous instability wave with real frequency is assumed to be created in the vicinity of the nozzle exit and to propagate towards the downward. Present results presented have shown in capturing small pressure fluctuating resulting from periodic oscillation of the jet around the edge. The pressure fluctuations propagate with the speed of sound. Their interaction with the wedge produces an irrotational feedback field which, near the nozzle exit, is a periodic transverse flow producing the singularities at the nozzle lips. It is found that, as the numerical example, satisfactory simulation results on the edgetone can be obtained for the complex flow-edge interaction mechanism, demonstrating the capability of the lattice Boltzmann model with flexible specific heat ratio to predict flow-induced noises in the ventilating systems of ship.

• 전산구조공학
• /
• 제7권4호
• /
• pp.85-101
• /
• 1994

적층판의 동적거동에 대한 유한요소해석모델개발을 목적으로 전단변형을 적합하게 고려한 적층판이론에 대한 변분원리를 유도하였다. 유도방법은 Sandhu 등에 의해 개발된 다변수 경계치문제의 변분원리이론을 따랐으며, 지배방정식의 미분연산자 매트릭스를 self-adjoint로 만들기 위하여 convolution을 이중선형사상으로 사용하였다. 유도된 적층판의 범함수에는 경계조건, 초기조건뿐만 아니라 유한요소해석모델에서 생길 수 있는 요소간 불연속조건도 포함시킬 수 있다. 상태변수의 적합함수공간을 확장하거나 특정조건을 적용하므로서 다양한 형태의 범함수를 유도할 수 있으며, 이를 통해 다양한 유한요소해석모델의 개발이 가능함을 논하였다.