• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2011.04a
• /
• pp.752-755
• /
• 2011

이종재료의 열전달문제 수치해석시 추가적으로 만족시켜야 하는 계면경계조건들의 존재와 계면경계로 인한 불연속면의 처리는 근사함수의 구성 뿐만 아니라 수치기법의 개발 자체를 어렵게 만든다. 본 논문에서는 계면경계의 불연속성을 모델링하는 특수한 함수를 포함하고 계면경계조건을 항상 만족시킬 수 있는 근사함수를 구성하고, 계면경계문제의 강형식을 직접 이산화하며 고속으로 해를 계산할 수 있는 이동최소제곱 차분법을 제시한다. 계면경계조건이 매입된 이동최소제곱 차분법으로 이종재료의 열전달문제를 해석한 결과, 높은 정확성과 효율성을 갖는 것을 확인할 수 있었다.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2010.04a
• /
• pp.48-51
• /
• 2010

본 논문에서는 자유경계문제 해석을 위해 정확도가 향상된 implicit 이동최소제곱 차분법을 제시한다. 계면경계에 대한 implicit 정의로 인해 비선형 시스템이 구성되고, 매 해석단계마다 절점해와 계면경계의 위치를 반복계산을 통해 찾는다. 계면경계 결정시 속도항을 한 단계 뒤로 지연시켜 explicit하게 근사적으로 계산하던 기존 방법에 비해 계면경계의 위치를 더 정확하게 계산할 수 있고, 결과적으로 해의 정확도가 향상되었다. 계면경계 위치값이 비교적 빠른 속도로 수렴하기 때문에 많은 반복계산이 필요치 않다. 수치예제를 통해 기존의 방법으로 계산한 결과와 비교하여 새롭게 개발한 implicit 방법의 향상된 정확도를 보였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.411-420
• /
• 2009

This study presents an extended finite difference method based on moving least squares(MLS) method for solving potential problems with interfacial boundary. The approximation constructed from the MLS Taylor polynomial is modified by inserting of wedge functions for the interface modeling. Governing equations are node-wisely discretized without involving element or grid; immersion of interfacial condition into the approximation circumvents numerical difficulties owing to geometrical modeling of interface. Interface modeling introduces no additional unknowns in the system of equations but makes the system overdetermined. So, the numbers of unknowns and equations are equalized by the symmetrization of the stiffness matrix. Increase in computational effort is the trade-off for ease of interface modeling. Numerical results clearly show that the developed numerical scheme sharply describes the wedge behavior as well as jumps and efficiently and accurately solves potential problems with interface.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.215-222
• /
• 2019

The heat conduction problem with discontinuous material coefficients generally consists of the conservative equation, boundary condition, and interface condition, which should be additionally satisfied in the solution procedure. This feature often makes the development of new numerical schemes difficult as it induces a layered singularity in the solution fields; thus, a special approximation is required to capture the singular behavior. In addition to the approximation, the construction of a total system of equations is challenging. In this study, a wedge function is devised for enriching the approximation, and the interface condition itself is embedded in the moving least squares(MLS) derivative approximation to consistently satisfy the interface condition. The heat conduction problem is then discretized in a strong form using the developed derivative approximation, which is named as the interface immersed MLS difference method. This method is able to efficiently provide a numerical solution for such interface problems avoiding both numerical quadrature as well as extra difference equations related to the interface condition enforcement. Numerical experiments proved that the developed numerical method was highly accurate and computationally efficient at solving the heat conduction problem with interfacial jump as well as the problem with a geometrically induced interfacial singularity.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.779-787
• /
• 2007

This paper presents a highly efficient moving least squares finite difference method (MLS FDM) for a heat transfer problem of bi-material with interfacial boundary. The MLS FDM directly discretizes governing differential equations based on a node set without a grid structure. In the method, difference equations are constructed by the Taylor polynomial expanded by moving least squares method. The wedge function is designed on the concept of hyperplane function and is embedded in the derivative approximation formula on the moving least squares sense. Thus interfacial singular behavior like normal derivative jump is naturally modeled and the merit of MLS FDM in fast derivative computation is assured. Numerical experiments for heat transfer problem of bi-material with different heat conductivities show that the developed method achieves high efficiency as well as good accuracy in interface problems.

• Journal of the Korea Organic Resources Recycling Association
• /
• v.17 no.1
• /
• pp.73-79
• /
• 2009

The hydraulic conductivity of porous media is the most important property in soil characteristics. The hydraulic conductivity is determined by outdoor and indoor methods. Indoor methods normally use soil columns for flow test. Assumption of the column test is that fluid one-dimensionally flows through the column. However, fluids may flow toward the wall of the column, resulting in "boundary flow". This study investigated the effect of boundary flow on the hydraulic conductivity by using a permeameter excluding boundary flow. The results showed that the hydraulic conductivity excluding boundary flow was much smaller than the hydraulic conductivity employing the conventional determination method. This study also investigated the effects of particle size and surfactant on the hydraulic conductivity. As the particle size increased, the hydraulic conductivity was increased. The hydraulic conductivity was reduced by increasing surfactant concentration. The result showed that the viscosity of fluid significantly affected the determination of hydraulic conductivity.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.315-322
• /
• 2009

This paper presents a novel numerical method based on the extended moving least squares finite difference method(MLS FDM) for solving 1-D Stefan problem. The MLS FDM is employed for easy numerical modelling of the moving boundary and Taylor polynomial is extended using wedge function for accurate capturing of interfacial singularity. Difference equations for the governing equations are constructed by implicit method which makes the numerical method stable. Numerical experiments prove that the extended MLS FDM show high accuracy and efficiency in solving semi-infinite melting, cylindrical solidification problems with moving interfacial boundary.

• Journal of the Microelectronics and Packaging Society
• /
• v.8 no.3
• /
• pp.25-30
• /
• 2001

The stress intensity factors for edge cracks located at the bonding interface between the semiconductor chip and the adhesive layer subjected to a uniform transverse tensile strain are investigated. Such cracks might be generated due to a stress singularity in the vicinity of the free surface. The boundary element method (BEM) is employed to investigate the behavior of interface stresses. The amplitude of complex stress intensity factor depends on the crack length, but it has a constant value at large crack lengths. The rapid propagation of interface crack is expected if the transverse tensile strain reaches a critical value.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.3
• /
• pp.887-897
• /
• 1991

본 연구에서는 이와 같은 배경에서, Fig. 1(f)와 같이 가장 일반적인 이방성 재료가 접합된 경우의 응력확대계수를 명확히 정의하고, 수치해석법으로 구할 수 있는 외삽식을 제안한다. 또한, 탄성문제의 수치해석 방법으로 적은 요소의 분할로써 고 정밀도의 수치해석 결과를 얻을 수 있는 경계요소법(boundary element method:BEM), 특히 저자들이 개발한 복합재료에 대한 2차원 경계요소법 프로그램을 이용하여 이방성 이종재료 접합계면 균열의 응력확대계수를 해석하고, 복합재료내의 섬유방향에 대한 접합계면 균열의 정성적 거동을 고찰하고자 한다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.3
• /
• pp.309-315
• /
• 2001

The Stress intensity factors for edge cracks located at the bonding interface between the elastic semiconductor chip and the viscoelastic adhesive layer have been investigated. Such cracks might be generated due to stress singularity in the vicinity of the free surface. The domain boundary element method(BEM) has been employed to investigate the behavior of interface stresses. The overall stress intensity factor for the case of a small interfacial edge crack has been computed. The magnitude of stress intensity factors decrease with time due to viscoelastic relaxation.