• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1634-1638
• /
• 2008

하천의 2차원 흐름 해석, 유사이동 해석, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있다. 유체에 대한 수치해석 기법으로 전통적으로 가장 많이 사용되고 있는 방법은 유한차분법이지만, 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 다양한 형태의 하천 해석에는 더욱 적합할 것이다. 본 연구에서는 비구조적 요소망을 advanced front method를 이용하여 생성해 보았다. Advanced front method는 해석하고자 하는 영역에 적절한 절점들을 생성한 후 삼각 요소망을 구성하는 grid based advanced front method와 절점들을 생성하지 않고 해석 영역에 삼각 요소를 바로 구성하는 element based advanced front method로 나눌 수 있다. Grid based advanced front method에서 해석 영역에 적절한 절점을 생성하는 방법으로는 일반적인 격자 구조의 절점 생성 방법을 적용하였으며 경계와의 거리가 가까운 절점은 생성되지 않으며, 삼각 요소를 구성할 때에는 두 개의 인접 절점을 비교하여 최적의 삼각 요소를 구성하게 된다. 단 두 개의 인접 절점만을 비교함으로서 비교적 빠른 시간 안에 최적의 삼각 요소망을 구성할 수 있다. 삼각 요소망을 생성한 후에는 Laplacian smoothing을 이용하여 삼각 요소망의 형질을 개선하였다. Element based advanced front method는 외부 경계에서부터 시작된 Front가 내부 영역으로 확대되어지며 각 Front에서 적절한 절점을 직접 생성하여 바로 삼각 요소를 구성하게 된다. Front에서 생성된 절점은 인접 절점들이 있는지 검색하여 인접 절점이 있다면 생성된 절점은 삭제되어지며 인접 절점이 삼각 요소를 위한 나머지 한 점으로 채택되어진다. Front는 외부 경계와 교차되어지지 않아야 하며 또한 연속된 Front를 효율적으로 관리하기 위해 list 자료 구조를 활용하였다.

• Proceedings of the KIEE Conference
• /
• 2004.10a
• /
• pp.55-59
• /
• 2004

유한요소법을 이용하여 전자장을 해석할 경우 전류원이 전 영역에 비해 극히 작은 영역이면, 요소분할 과정에서 소스부분을 세분하여야 하므로 결국 미지수의 증가를 가져오게 된다. 또한, 선전류 문제의 경우 2차원 유한 요소 해석이 용이하지 않다. 이를 보안하기 위해 본 논문에서는 소스가 선전류이고 관심 영역이 선전류원으로부터 떨어져 있는 경우, 소스 영역은 해석해를 적용하여 유한요소법과 결합하는 방법을 제시하였다. 해석적인 해는 원통좌표계에서 반정에 대한 멱함수와 회전각도에 대한 삼각함수의 곱의 형태로 표현된다. 이때 두 종류의 적분 상수가 있는데, 이는 경계상의 포텐셜값과 유한요소법의 경계 적분항을 푸리에급수로 전개한 계수로 표현된다. 제안한 알고리즘의 검증을 위하여 해석해가 존재하는 모델을 설정하여 해석적인 방법, 기존의 유한요소 법 및 결합 방법에 의한 해를 비교 검증하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.2
• /
• pp.199-208
• /
• 2008

In this paper, to analyze the initial valued non-homogeneous elastic half space by the scaled boundary analysis, the infinite element approach was introduced. The free surface of the initial valued non-homogeneous elastic half space was modeled as a circumferential direction of boundary scaled boundary coordinate. The infinite element was used to represent the infinite length of the free surface. The initial value of material property(elastic modulus) was considered by the combination of the position of the scaling center and the power function of the radial direction. By use of the mapping type infinite element, the consistent elements formulation could be available. The performance and the feasibility of proposed approach are examined by two numerical examples.

• Journal of Ocean Engineering and Technology
• /
• v.23 no.4
• /
• pp.6-11
• /
• 2009

In this study, the internal waves in two-density layered fluids were analyzed using the Numerical Wave Tank (NWT) technique in the frequency domain. The NWT is based on a two-domain Boundary Element Method with the potential fluids using the whole-domain matrix scheme. From the mathematical solution of the two-domain boundary integral equation, two different wave modes could be classified: a surface wave mode and an internal wave mode, and each mode were shown to have a wave number determined by a respective dispersion relation. The magnitudes of the internal waves against surface waves were investigated for various fluid densities and water depths. The calculated results are compared with available theoretical data.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.93-99
• /
• 1996

In the particular situations where the crack is terminated at an interface of two materials, the order of stress singularity depends on the elastic constants which specify the properties of two materials. A multidomain boundary element technique is used to solve a crack normal to bimaterial interface. A correct order of shape function is used for displacement by using the isoparametric elements by shifting adequately the side nodes adjacent to this crack tip. A shape function containing the same order of singularity as that in the interface crack is also used for the interpolation of traction. Numerical testing of a binaterial with a crack normal to the interface is carried out with three-node elements. The results obtained are compared with the previous solutions.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.12 no.3
• /
• pp.139-148
• /
• 2000

In this paper, the interaction of oblique incident waves with a bottom-mounted and fluid-filled flexible membrane structure is investigated in the frame of linear hydro-elastic theory. The static shape of a membrane structure containing the fluid of a specific density is initially unknown and must be calculated before the hydrodynamic analysis. To solve hydrodynamic problem, the fluid domain is divided into the inner and outer region. The inner solution based on discrete membrane dynamic model and simple-source distribution over the entire fluid boundaries is matched to the outer solution ba~ed on an eigenfunction expansion method. The numerical results were compared to a series of Ohyama's experimental results. The measured reflection and tran¬smission coefficients reasonably follow the trend of predicted values. Using the computer program developed, the performance of a bottom-mounted and fluid-filled flexible membrane strocture is tested with various system parameters (membrane shape, internal pressure, density ratio) and wave characteristics (wave frequencies, incident wave angle). It is found that a bottom-mounted and fluid-filled flexible membrane structure can be an effel;tive wave barrier if properly designed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.499-512
• /
• 2016

The boundary reaction method(BRM) is a substructure time domain method, it removes global iterations between frequency and time domain analyses commonly required in the hybrid approaches, so that it operates as a two-step uncoupled method. The BRM offers a two-step method as follows: (1) the calculation of boundary reaction forces in the frequency domain on an interface of linear and nonlinear regions, (2) solving the wave radiation problem subjected to the boundary reaction forces in the time domain. In the time domain analysis, the near-field soil is modeled to simulate the wave radiation problem. This paper evaluates the performance of the BRM according to modeling extent of near-field soil for the nonlinear SSI analysis of base-isolated NPP structure. For this purpose, parametric studies are performed using equivalent linear SSI problems. The accuracy of the BRM solution is evaluated by comparing the BRM solution with that of conventional SSI seismic technique. The numerical results show that the soil condition affects the modeling range of near-field soil for the BRM analysis as well as the size of the basemat. Finally, the BRM is applied for the nonlinear SSI analysis of a base-isolated NPP structure to demonstrate the accuracy and effectiveness of the method.

• The Journal of the Acoustical Society of Korea
• /
• v.20 no.6
• /
• pp.94-103
• /
• 2001

The goal of using numerical methods in this study is two-fold: to replicate a set of measured, individualized HRTFs by a computer simulation, and also to visualise the resultant sound field around the head. Two methods can be wed: the Boundary Element Method (BEM) and the Infinite-Finite Element Method (IFEM). This paper presents the results of a preliminary study carried out on a KEMAR dummy-head, the geometry of which was captured with a high accuracy 3-D laser scanner and digitiser. The scanned computer model was converted to a few valid BEM and IFEM meshes with different polygon resolutions, enabling us to optimise the simulation for different frequency ranges. The results show a good agreement between simulations and measurements of the sound pressure at the blocked ear-canal of the dummy-head. The principle of reciprocity provides an effect method to simulate HRTF database. The BEM was also used to investigate the total sound field around the head, providing a tool to visualise the sound field for different arrangements of virtual acoustic imaging systems.

• Composites Research
• /
• v.23 no.3
• /
• pp.19-30
• /
• 2010

The present study describes a method for analytical solution to elastic field in structural elements of general symmetric laminated composite materials. The two dimensional plane stress elasticity problems under mixed boundary conditions are reduced to the solution of a single fourth order partial differential equation, expressed in terms of a single unknown function, called displacement potential function. In addition, all the components of stress and displacement are expressed in terms of the same displacement potential function, which makes the method suitable for any boundary conditions. The method is applied to obtain analytical solutions to two particular problems of structural elements consisting of an angle-ply laminate and a cross-ply laminate, respectively. Some numerical results are presented for both the problems with reference to the glass/epoxy composite. The results are highly accurate and reliable as all the boundary conditions including those in the critical regions of supports and loads are satisfied exactly. This verifies the method as a simple and reliable one as well as capable to obtain exact analytical solution to elastic field in structural elements of composite materials under mixed and any other boundary conditions.

• Journal of Navigation and Port Research
• /
• v.29 no.10 s.106
• /
• pp.853-858
• /
• 2005

Recently, various novel numerical models based on Navier-Stokes equation have been developed for calculating wave motions in the sea with coastal or ocean structures. Among those models, Volume Of Fluid (VOF) method might be the most popular one, and it has been used for numerical simulations of wave motions including complicated phenomena of wave breakings. VOF method, however, needs enormous computation time and large computational storage memories in general, thus it is practically difficult to use this method for calculations in the case of random waves because long and stable computation (e.g for more than 100 significant wave periods) is required to obtain statistically meaningful results. On the other hand if the wave motion is potential motion, Boundary Element Method (BEM), which is a much faster and more accurate method than VOF method, can be effectively used. The aim of this study is to develop a new efficient model applicable to calculations of wave motion and/or wave-structure interactions under random waves. To achieve this, a strictly combined BEM-VOF model has been developed by making the best use of both methods' merits; VOF method is used in a restricted fluid domain around a structure where complicated phenomena of wave breakings may exist, and BEM is used in the other domains far from the disturbance where the wave motion may be assumed to be potential. The verification of the model was performed with numerical results for Stokes' 5th order wave propagation and a random wave propagation.