• 한국지반환경공학회 논문집
• /
• 제25권4호
• /
• pp.21-28
• /
• 2024

본 연구에서는 SPH 해석을 위한 다면체영역분할 기법이 소개된다. SPH 기법은 유체 유동 모사를 위한 수치해석기법으로 무요소기법(meshless method) 중 하나이다. 유동성 지반 또는 고체-유체 상호작용 해석 등에 유용하게 쓰일 수 있다. SPH는 입자기반 해석이기 때문에 입자가 많을수록 결과의 정확도는 높아지지만 수치적 효율성은 떨어진다. 일반적으로 해석의 효율성을 높이기 위해 병렬 프로세싱 알고리즘과 함께 쓰이는데 직교좌표계 기반의 영역분할 기법이 대표적이다. 그러나 복잡한 기하학적 형태나 동적 경계조건에서 유동 모사 등을 병렬 해석하기 위해서는 직교좌표계 영역분할 방법이 적합하지 않다. 소개하는 다면체영역분할 기법은 이와 같은 문제에서 병렬효율성을 높일 수 있는 장점을 갖는다. 다양한 형태의 3차원 다면체 요소로 분할하여 문제에 적합하게 모델링할 수 있다. SPH 입자들의 물리적 값들은 smoothing 길이 이내의 주위 입자들 정보를 이용하여 계산된다. 영역분할 시 물리적으로 분리될 수 있는 입자정보들을 코어간 공유할 수 있는 방법과 병렬효율성이 떨어질 수 있는 cross-point에서의 정보공유 방법이 소개된다. 수치해석 예제를 통하여 제안된 방법의 병렬효율성은 12코어까지 95%에 근접하였다. 이후 코어가 증가할수록 코어간 공유되는 정보량이 많아져 병렬효율성이 떨어지는 문제가 발생되기도 하였다.

• Structural Engineering and Mechanics
• /
• 제68권5호
• /
• pp.603-619
• /
• 2018

In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

• 한국가시화정보학회지
• /
• 제19권3호
• /
• pp.92-98
• /
• 2021

Large Eddy Simulation (LES) has been popularly applied and used in the last several decades to simulate turbulent boundary layer in the numerical domain. A fully developed turbulent boundary layer has also been applied to predict the complicated wake flow behind bluff bodies. In this study we aimed to generate an artificial turbulent boundary layer, which is based on an exponential correlation function, and generates a series of realistic three-dimensional velocity data in two-dimensional inlet section which are correlated both in space and in time. The results suggest its excellent capability for high Reynolds number flows. To make an effective generation, a hexahedral mesh has been used and Cholesky decomposition was applied to possess suitable turbulent statistics such as the randomness and correlation of turbulent flow. As a result, the flow characteristics in the domain and fluctuating pressure near the wall are very close to those of fully developed turbulent boundary layers.

• Structural Engineering and Mechanics
• /
• 제67권4호
• /
• pp.317-326
• /
• 2018

The iterative decomposition coupling formulation of the precise integration finite element method (FEM) and the time domain boundary element method (TD-BEM) is presented for elstodynamic problems. In the formulation, the FEM node and the BEM node are not required to be coincident on the common interface between FEM and BEM sub-domains, therefore, the FEM and BEM are independently discretized. The force and displacement converting matrices are used to transfer data between FEM and BEM nodes on the common interface between the FEM and BEM sub-domains, to renew the nodal variables in the process of the iterations for the un-coincident FEM node and BEM node. The iterative coupling formulation for elastodynamics in current paper is of high modeling accuracy, due to the semi-analytical solution incorporated in the precise integration finite element method. The decomposition coupling formulation for elastodynamics is verified by examples of a cantilever bar under a Heaviside-type force and a harmonic load.

• Structural Engineering and Mechanics
• /
• 제20권2호
• /
• pp.225-239
• /
• 2005

Applications of the Trefftz boundary element method (BEM) to anti-plane electroelastic problems are presented in this paper. Both direct and indirect methods with domain decomposition are discussed in details. Each crack is treated as semi-infinite thin slit defined in a subregion, for which a particular solution of the anti-plane problem, satisfying exactly the crack-face condition, is derived. The stress intensity factors defined at each crack tip can be directly computed from the coefficients of the particular solution. The performance of the proposed formulation is assessed by two examples and comparison is made with results obtained by other approaches. The Trefftz boundary element approach is demonstrated to be suitable for the analysis of the anti-plane problem of piezoelectric materials.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1995년도 추계 학술대회논문집
• /
• pp.82-89
• /
• 1995

In the present paper, a method of nearly orthogonal grid generation in an arbitrary simply-connected 3D domain will be presented. The method is a new direct and non-iterative scheme based on the concept of the decomposition of the global orthogonal transformation into consecutive mapping of a conformal mapping and an auxiliary orthogonal mapping, which was suggested by King and Leal [4]. In our numerical scheme. Kang and Leal's method is extended from 2D problems to 3D problems while the advantage of the non-iterative algorithm is maintained. The essence of the present mapping method is that an iterative scheme can be avoided by introducing a preliminary step. This preliminary step corresponds to a conformal map and is based on the boundary element method(BEM). This scheme is applied to generate several nearly-orthogonal grid systems which are orthogonal at boundaries.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2006년도 정기 학술대회 논문집
• /
• pp.291-298
• /
• 2006

A parallelized topology design optimization method is developed on a distributed memory system. The parallelization is based on a domain decomposition method and a boundary communication scheme. For the finite element analysis of structural responses and design sensitivities, the PCG method based on a Krylov iterative scheme is employed. Also a parallelized optimization method of optimality criteria is used to solve large-scale topology optimization problems. Through several numerical examples, the developed method shows efficient and acceptable topology optimization results for the large-scale problems.

• Journal of Mechanical Science and Technology
• /
• 제20권10호
• /
• pp.1773-1778
• /
• 2006

Steady-state natural convection heat transfer characteristics from cylinders in a multiply-connected bounded region are clarified. A spectral finite difference scheme (spectral decomposition of the system of partial differential equations, semi-implicit time integration) is applied in numerical analysis, with a boundary-fitted conformal coordinate system through a Jacobian elliptic function with a successive transformation to formulate a system of governing equations in terms of a stream function, vorticity and temperature. Multiplicity of the domain is expressed explicitly.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제7권6호
• /
• pp.1034-1043
• /
• 2015

This is the continuation of a numerical study on vortex shedding from a blunt trailing-edge of a hydrofoil. In our previous work (Lee et al., 2015), numerical schemes for efficient computations were successfully implemented; i.e. multiple domains, the approximation of domain boundary conditions using cubic spline functions, and particle-based domain decomposition for better load balancing. In this study, numerical results through a hybrid particle-mesh method which adopts the Vortex-In-Cell (VIC) method and the Brinkman penalization model are further rigorously validated through comparison to experimental data at the Reynolds number of $2{\times}10^6$. The effects of changes in numerical parameters are also explored herein. We find that the present numerical method enables us to reasonably simulate vortex shedding phenomenon, as well as turbulent wakes of a hydrofoil.

• 한국산학기술학회논문지
• /
• 제4권3호
• /
• pp.263-269
• /
• 2003

이 논문은 3차원 정수 웨이브릿 변환을 이용한 손실 의료 영상 압축에 대한 방법을 보여준다. 의료영상에 3차원 웨이브릿 분할법과 3차원 spatial dependence tree를 이용한 Set Partitioning in Hierarchical Trees(SPIHT) algorithm을 제공한다. 또한 3차원 웨이브릿 분할법에서 정수 웨이브릿 필터들을 이용한 lifting방법을 이용하고, unitary 변환을 만들기 위하여 정확한 scaling을 이용한다. 압축률이 증가하면 할수록 인접한 coding unit사이에선 boundary effect가 생긴다. Video와 같지 않아서 인접한 coding unit사이에서의 boundary artifact는 보여서는 안 된다. 이러한 현상을 제거하기 위해서 인접한 coding unit사이에 axial domain으로 overlapping방법을 사용한다. 또한 코딩 할 때 여러 종류의 정수 필터들을 사용한다. 결과로 어떤 특정한 필터를 사용할 때 좋은 결과를 얻었고 overlapping방법을 사용했기 때문에 ringing artifacts는 거의 찾아 볼 수가 없게 되었다. 그리고 어두운 면을 코딩 할 때도 좋은 결과를 얻었다.