• Journal of computational fluids engineering
• /
• v.4 no.3
• /
• pp.28-34
• /
• 1999

The objective of this study is to present a numerical treatment of the pressure gradient when control volumes are sharing the interface between a homogeneous fluid and a porous medium. Two possible approaches, e.g. linear interpolation and extrapolation, are considered, and they are applied to the case of a steady and two-dimensional curved channel flow which is partially filled with a porous medium. It was found that the linear extrapolation produces a continuous velocity-field at the interface and thus is recommended. On the contrary, the linear interpolation entails a discontinuous velocity field at the interface, thereby warning its use in connection with the Brinkman-Forchheimer-extended Darcy flow model.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.386-386
• /
• 2015

담수가 해수에서 흘러드는 하구에서는 성층이 관측되며 이것은 난류의 미세구조를 변화시키는 주요 원인으로 작용한다. 이러한 성층화 현상은 하구 내 부유사의 군집인 하구 최대혼탁수(Estuarine Turbidity Maximum, ETM)의 형성에 영향을 주게 된다. 본 연구는 성층의 하구 최대 혼탁수 생성 메커니즘에 관심을 두고 수치모델링을 활용한 미세 난류의 부유사 거동 분석에 초점을 두었다. 성층과 전단응력 사이의 난류 혼합을 대표하는 유동인 켈빈-헬름홀츠 불안정성(Kelvin-Helmholtz Instability)을 도입하고 성층 경계면 근처에서 부유사의 이송을 높은 레이놀즈수(Reynolds number) 유동에서 RANS(Reynolds-averaged Navier-Stokes Simulation)보다 다양한 규모의 에너지 획득이 가능하여 미세 난류 구조 재현에 장점을 갖는 Large-eddy Simulation(LES)를 활용하여 모사하였다. 여기에서, 부유사는 주위 유동의 물리적 특성 변화에 영향을 미치지 않는 Passive scalar로 가정하였으며 $6^{th}$-order Lagrangian 다항식 보간법을 적용하여 입자의 이동 속도를 계산하고 이를 시간에 대해 적분함으로써 이동 궤적을 추적하였다. 수치 모델 결과 Lock-exchange 유동 내에서 켈빈-헬름홀츠 불안정성이 발생함에 따라 경계면 주위에 위치한 부유사가 billow 내에서 트랩핑(trapping)되는 것을 보여주어 KH-billow 혹은 braids 내의 미세 난류에 의한 영향이 확인되었다. 본 연구에서는 LES를 활용하여 성층류 및 성층류 내의 부유사 혼합을 모사하여 난류의 정도에 따른 이동 궤적의 차이에 대해서 분석함으로써 성층의 난류 강도 저하에 따른 부유사의 군집으로의 영향에 대해 서술한다.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.9
• /
• pp.67-75
• /
• 2006

A fully nonlinear model accounting for swirling effect has been applied in analyzing the dynamic response for a classical swirl injector. The current work applied highly accurate Boundary Element Methods (BEMs) in assessing its static and dynamic characteristics. On the basis of moving surface treatment method and surface instability study, which are obtained from the previous static characteristics analysis in pressure-swirl injectors, this work was expanded for analyzing the dynamics of a classical swirl injector. The dynamic response through injector components for disturbed inflow condition was investigated. The modified code was validated from comparison with the theoretical result for a typical swirl injector. Clearly the simulated result shows the interesting characteristics of swirl injectors to provide either amplification or damping of the input disturbance through each component. These results give promise in applying the current model to nonlinear dynamic characteristics of swirl injectors.

• Journal of the Korean Society for Precision Engineering
• /
• v.9 no.2
• /
• pp.108-115
• /
• 1992

이종재료의 접합면에 수직으로 존재하는 크랙에 대하여 경계요소 해석을 수행하여, 그 결과 실용가능한 수치 근사해을 얻을 수 있었다. 크랙을 정확히 모델링하기 위하여 크랙표면을 분리영역으로 하는 영역분할법을 채택하였으며, 해의 정확성을 향상시키기 위하여 등매개 2차요소로의 경계분할과 함께 크랙선단에서 표면력의 특이성을 나타내도록 하였다. 응력확대계수는 크랙표면상 절점의 상대변위를 이용하여 결정하였다. 또한 이종 재료내 크랙에 대하여 응력확대계수를 간단히 구할 수 있는 간편해석법을 제안하고 이의 적용 가능한 범위를 제 시하였다.

• Proceedings of the Korea Institute of Fire Science and Engineering Conference
• /
• 2009.04a
• /
• pp.404-409
• /
• 2009

수직벽 화재 예측의 정확성을 확인하기 위하여 화재 시뮬레이션용 전산유체역학 모델인 Fire Dynamics Simulator를 프로필렌 수직벽 화재에 적용하였다. 단위면적당 연소율 $7.0{\sim}29.29g/m^2-s$에 대한 버너 중심에서 측정한 온도분포와 비교한 결과, 최고온도가 낮게 예측되는 것 외에는 실험과 잘 일치하였다. 또 연소율의 증가에 따라 경계측의 두께가 일관되게 증가하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.6
• /
• pp.515-522
• /
• 2017

In this paper, an isogeometric analysis for multipatch problem is investigated, in which two or more geometries are connected at the interface in a conforming or non-conforming conditions. To express higher continuity at the patch interface, two approaches such as Nitsche based method and master-slave method are formulated for the linear elasticity problem and discretized using the isogeometric approach using NURBS basis functions. A short comparison between two approaches in formulations reveals the pros and cons of them with the applicability in the isogeometric multipatch problem. In addition, a NURBS based stress recovery is adopted to express a better stress continuity through the post-processing. Numerical examples indicate the effectiveness of Nitsche method in the non-conforming patch, following the exact solution well. For the stress concentration problem with the conforming patch, introduced two methodologies show comparative results, meanwhile the NURBS based stress recovery presents an improved smooth stress contour in the whole domain including the patch interface.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.1
• /
• pp.10-17
• /
• 1985

A new method is suggested for the solution of plane elasticity problems. With use of the dislocation model in the crack problems, the basic scheme of this method is to find equilibrium Burgers vectors of dislocations which are distributed along the boundary of the first fundamental boundary value problems. The stress distribution in the region can be found by superposition of the contributions of each dislocation. The method is applied to three cases with known analytical solutions, and to a V-notched specimen under uniaxial tension. The numerical results are compared with other available solutions. This method is effective and simple in its use, compared with other numerical methods. The method also provides very accurate solutions in the region except near the boundary where the discretization error is significant. The extrapolation method is suggested for the stresses in the boundary region. Extensive application are also suggested for a general estimate of the computational efficiency of the method.

• Journal of the Korean Geosynthetics Society
• /
• v.23 no.1
• /
• pp.9-16
• /
• 2024

In order to evaluate the relationship between the ground confining pressure and the shear characteristics of the pile-soil interface, this study described the comparative results of the existing experimental results and the FEA results using the strength reduction factor. The strength reduction factor was applied to simulate the shear behavior of the pile-soil interface in finite element analysis(FEA). The analysis results showed that the maximum pullout resistance decreased due to the influence of low confining pressure, as the fines content increased. This trend was similar to the previous experimental research, and this FEA model simulated with the interface strength reduction factor was evaluated as reasonable. The analysis results of the variation in the interface strength reduction factor clearly showed that the interface strength reduction factor clearly increased at a high fines content when the confining pressure was 50kPa. However, it was found that the increase rate was low when the confining pressure was 100kPa and 150kPa. Therefore, confining pressure and fines content need to be considered in FEA to evaluate the shear behavior of the pile-soil interface.

• Explosives and Blasting
• /
• v.37 no.4
• /
• pp.26-35
• /
• 2019

In designing a gravity-type anchorage of earth-anchored suspension bridge, the contact friction between a blasted rock mass and the concrete anchorage plays a key role in the stability of the entire anchorage. Therefore, it is vital to understand the shear behavior of the interface between the blasted rock mass and concrete. In this study, a portable 3D LiDAR scanner was utilized to scan the blasted bottom surfaces, and rock surface roughness was quantitatively analyzed from the scanned profiles to apply to 3D FEM modelling. In addition, based on the 3D FEM model, a three-dimensional dynamic fracture process analysis (DFPA-3D) technique was applied to study on the shear behavior of the interface between blasted rock and concrete through direct shear tests, which was analyzed under constant normal load (CNL). The effects of normal stress and the joint roughness on shear failure behavior are also analyzed.

• Journal of Korea Water Resources Association
• /
• v.48 no.12
• /
• pp.1023-1035
• /
• 2015

This study applies 3-D N-S solver based on PBM (Porous Body Model), LED-WASS-3D ver 2.0 to directly analyze non-linear interaction of seawater-freshwater-coastal aquifer in order to simulate the seawater infiltration into coastal aquifer. This numerical simulation is the first trial in Korea, as well as unusual and new numerical analysis abroad. Firstly, to validate the applied numerical model, the validity and effectiveness was verified for the numerical model by comparing and considering it with the result of laboratory experiment for seawater-freshwater interface in coastal aquifer. And then it simulated the seawater infiltration into coastal aquifer considering the changed levels of seawater and groundwater in order to analyze the distribution characteristics of flow field and seawater-freshwater interface of coastal aquifer as the level difference between seawater and groundwater and rate of seawater level (${\Delta}h/h$) increased. In addition, the characteristics of seawater infiltration were analyzed from the vertical salinity in the coastal aquifer by ${\Delta}h/h$, which cannot be obtained from existing non-diffusion numerical models. Finally, it analyzed the effect of ${\Delta}h/h$ on the seawater infiltration distance in coastal aquifer, which was indexed.