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

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

• 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.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.24 no.5
• /
• pp.577-588
• /
• 2011

This paper provides a novel extended Moving Least Squares(MLS) difference method for the potential problem with weak and strong discontinuities. The conventional MLS difference method is enhanced with jump functions such as step function, wedge function and scissors function to model discontinuities in the solution and the derivative fields. When discretizing the governing equations, additional unknowns are not yielded because the jump functions are decided from the known interface condition. The Poisson type PDE's are discretized by the difference equations constructed on nodes. The system of equations built up by assembling the difference equations are directly solved, which is very efficient. Numerical examples show the excellence of the proposed numerical method. The method is expected to be applied to various discontinuity related problems such as crack problem, moving boundary problem and interaction problems.

• Journal of Energy Engineering
• /
• v.24 no.1
• /
• pp.114-122
• /
• 2015

In the high convective gas flow condition, irregular shaped water waves from which droplet entrainment occurs are generated under horizontally stratified two-phase flow condition. KAERI proposed a new mechanistic droplet entrainment model based on the momentum balance equation consisting of the shear stress, surface tension, and gravity forces. However, this model requires correlation or experimental data of several physical parameters related to the wave characteristics. In the present study, we tried to measure the physical parameters such as wave slope, wave hypotenuse length, wave velocity, wave frequency, and wavelength experimentally. For this, an experiment was conducted in the horizontal rectangular channel of which width, height, and length are, respectively, 40 mm, 50 mm, and 4.2 m. In the present test, the working fluids are chosen as air and water. The PIV technique was applied not only to obtain images for phase interface waves but also to measure the velocity field of the water flow. Additionally, we developed the parallel wire conductance probe for the confirmation of wave height from PIV image. Finally, we measured the physical parameters to be used in the validation of new droplet entrainment model.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1991.07a
• /
• pp.93-97
• /
• 1991

수표면에 방류되는 온배수등의 흐름과 같은 표층방류밀도분류는 자유난류의 전단류 효과와, 방류수와 주위수의 밀도차에 기인하는 부력효과를 동시에 받는 흐름장을 형성한다. 또한, 이 흐름은 수표면 및 밀도계면에 의해 2 개의 자유경계에 둘러싸인 특이한 경계조건때문에 개수로 흐름으로 대표되는 자유전단류와 구별된다.(중략)

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.6
• /
• pp.1082-1091
• /
• 1989

본 논문에서는 응고계면을 기준으로 하여 액상역, 고상역 및 로울등을 각각 계산 가능한 영역으로 좌표변환하는 경계공정법을 사용하여 로울두께 방향의 온도분 포와 고상역과 액상역의 속도 분포를 고려한 2차원 응고해석을 하여 모델재료를 이용 한 실험결과와 비교 검토하여 이론해석의 타당성을 검정한다. 그리고 열전도율이 연 강보다 적어 박판제조가 어려운 재료인 스테인리스강을 용탕으로부터 직접 생산하기 위한 압연조건을 정량적으로 제시하며, 압연로울의 냉각특성을 밝힌다.

• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.185-193
• /
• 1997

The order of the stress singularity that occurs at the termination of an interface between materials exhibiting bilinear stress-strain response under plane strain conditions has been calculated, The governing equation of elasticity together with traction-free boundary condition and interface continuity condition defines a two-point boundary value problem. The stress components near the free edge are assumed to be proportional to r/sup s-1/, with solutions existing only for certain values of s. Finding these values entails the solution of an eigenvalue problem. Because it has been impossible to integrate the differential equations analytically, the integration has been performed numerically with a shooting method coupled with a Newton improvement scheme.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.1
• /
• pp.39-48
• /
• 2013

This paper develops a 2-D moving least squares(MLS) difference method for Stefan problem by extending the 1-D version of the conventional method. Unlike to 1-D interfacial modeling, the complex topology change in 2-D domain due to arbitrarily moving boundary is successfully modelled. The MLS derivative approximation that drives the kinetics of moving boundary is derived while the strong merit of MLS Difference Method that utilizes only nodal computation is effectively conserved. The governing equations are differentiated by an implicit scheme for achieving numerical stability and the moving boundary is updated by an explicit scheme for maximizing numerical efficiency. Numerical experiments prove that the MLS Difference Method shows very good accuracy and efficiency in solving complex 2-D Stefan problems.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2015.10a
• /
• pp.100-103
• /
• 2015

본 연구에서는 선박의 내항 및 조종성능평가를 위해 사용되는 2차원부가질량계수와 파랑감쇠계쉬를 desingularized indirect boundary integral method를 이용하여 주파수영역에서 계산하였다. 본 방법을 Rankine source를 유체경계면 전체에 분포하는 방법으로 수학적으로 간단하고 경계조건의 변경 및 적용이 용이하다는 장점이 있다. 하지만 물체표면, 자유수면 그리고 방사경계면의 계산영역 및 요소의 배치에 따라 계산정확도에 차이가 발생한다. 따라서 본 연구에서는 수치시험을 통해, 자유수면과 방사계면의 적절한 계산영역과 요소의 개수를 결정하였다. 계산의 정도와 효율성을 확보하기 위하여, 자유수면과 방사경계면의 계산영역과 요소의 분포를 파의주파수에 따라 달리 적용하였다. 본 연구에서 제안된 수치방법을 활용하여 계산된 동유체력과 실험에서 계측된 동유체력을 비교하여 본 방법의 정도를 확인하였다.