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

본 논문에서는 유동체를 해석할 수 있는 수치해석기법 중 하나인 SPH(Smoothed Particle Hydrodynamics)의 병렬해석 알고리즘이 소개된다. 무요소법(meshless method)의 SPH는 연속체 거동을 입자기반으로 표현하기 때문에 컴퓨팅하는데 높은 자원을 요구한다. 그래서 병렬해석 알고리즘은 SPH 시뮬레이션에서 필수적으로 고려되어야 한다. 계산영역을 일정한 간격으로 분할시켜 독립적으로 해석하는 영역분할 알고리즘은 병렬해석 알고리즘 중에 가장 대표적인 방법이다. 그리고 그 중 Cartesian 좌표계의 영역분할 방법은 입자들의 좌표를 빠르고 편리하게 검색할 수 있는 장점이 있어, DEM(Discrete Element Method)이나 MD(Molecular Dynamics)에서 대중적으로 사용되고 있다. 그러나 SPH의 경우 입자들이 smoothing 길이 이내의 주위 입자 정보가 필요하기 때문에 분할 영역 간의 입자정보 공유가 중요하다. 그리고 이에 따른 CPU의 로드밸런스가 중요하다. 본 연구에서는 직교 영역분할의 크기를 동적으로 미소화 시켜 잉여 CPU가 발생하지 않도록 하는 높은 병렬효율성의 알고리즘이 제안되었다. 그리고 수치해석 모델을 통하여 효율성을 검증하였다. 유동체 모델에 대해 총 30 CPU까지 제안된 방법의 병렬효율성을 검토하였고, 28개의 물리적 코어 수까지 90%의 병렬효율성을 얻을 수 있었다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회B
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Korea-Australia Rheology Journal
• /
• 제21권1호
• /
• pp.1-12
• /
• 2009

This paper reports the suitability of a domain decomposition technique for the hybrid simulation of dilute polymer solution flows using Eulerian Brownian dynamics and Radial Basis Function Networks (RBFN) based methods. The Brownian Configuration Fields (BCF) and RBFN method incorporates the features of the BCF scheme (which render both closed form constitutive equations and a particle tracking process unnecessary) and a mesh-less method (which eliminates element-based discretisation of domains). However, when dealing with large scale problems, there appear several difficulties: the high computational time associated with the Stochastic Simulation Technique (SST), and the ill-condition of the system matrix associated with the RBFN. One way to overcome these disadvantages is to use parallel domain decomposition (DD) techniques. This approach makes the BCF-RBFN method more suitable for large scale problems.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.55-58
• /
• 2002

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.

• 한국전산구조공학회논문집
• /
• 제24권4호
• /
• pp.365-373
• /
• 2011

병렬과 영역분할을 이용한 폭발하중을 받는 철근콘크리트패널의 손상을 분석하였다. 폭풍파는 극도로 짧은시간 동안에 발생되기 때문에 수치해석을 통한 결과값은 폭풍파의 메쉬크기에 영향을 받는다. 그러므로 폭풍파 메쉬크기의 영향을 분석하기 위해 explicit 유한요소해석 프로그램인 AUTODYN을 이용하여 기존 실험결과와 메쉬크기에 따른 해석결과가 비교되었다. 폭발해석결과 메쉬크기가 작을수록 정확도가 높았으나 수행시간이 증가하여 효율성이 떨어졌다. 추가로 수치해석의 효율성을 높이기 위해 영역별 Euler와 Lagrange 기법을 달리하는 병렬해석이 수행되었다. 결과로, 폭풍파영역에서는 영역분할된 Euler 메쉬를 사용하고 구조물영역에서는 영역 분할된 Lagrange 메쉬를 사용하는 것이 구조물영역에서 영역 분할된 Lagrange 메쉬만을 사용한 것보다 수치효율성이 가장 높았다.

• International Journal of Aeronautical and Space Sciences
• /
• 제16권2호
• /
• pp.177-189
• /
• 2015

This paper presents an all-direct domain decomposition approach for large-scale structural analysis. The proposed approach achieves computational robustness and efficiency by enforcing the compatibility of the displacement field across the sub-domain boundaries via local Lagrange multipliers and augmented Lagrangian formulation (ALF). The proposed domain decomposition approach was compared to the existing FETI approach in terms of the computational time and memory usage. The parallel implementation of the proposed algorithm was described in detail. Finally, a preliminary validation was attempted for the proposed approach, and the numerical results of two- and three-dimensional problems were compared to those obtained through a dual-primal FETI approach. The results indicate an improvement in the performance as a result of the implementing the proposed approach.

• Journal of applied mathematics & informatics
• /
• 제8권1호
• /
• pp.27-44
• /
• 2001

A fast Poisson solver on irregular domains, based on bound-ary methods, is presented. The harmonic polynomial approximation of the solution of the associated homogeneous problem provides a good practical boundary method which allows a trivial parallel processing for solution evaluation or straightfoward computations of the interface values for domain decomposition/embedding. AMS Mathematics Subject Classification : 65N35, 65N55, 65Y05.

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1999년도 유체기계 연구개발 발표회 논문집
• /
• pp.200-207
• /
• 1999

An efficient inverse design technique based on the MGM (Modified Garabedian-McFadden) method has been developed. The 2-D Navier-Stokes equations are solved for obtaining the surface pressure distributions and coupled with the MGM method to perform the inverse design. The solver is parallelized by using the domain decomposition method and the standard MPI library for communications between the processors. The MGM method is a residual-correction technique, in which the residuals are the difference between the desired and the computed pressure distribution. The developed code was applied to several airfoil shapes and the axial blade. It has been found that they are well converged to their target pressure distribution.

• Structural Engineering and Mechanics
• /
• 제38권4호
• /
• pp.405-416
• /
• 2011

An efficient two-level domain decomposition parallel algorithm is suggested to solve large-DOF structural problems with nonlinear material models generating unsymmetric tangent matrices, such as a group of plastic-damage material models. The parallel version of the stabilized bi-conjugate gradient method is developed to solve unsymmetric coarse problems iteratively. In the present approach the coarse DOF system is solved parallelly on each processor rather than the whole system equation to minimize the data communication between processors, which is appropriate to maintain the computing performance on a non-supercomputer level cluster system. The performance test results show that the suggested algorithm provides scalability on computing performance and an efficient approach to solve large-DOF nonlinear structural problems on a cluster system.

• 한국전산구조공학회논문집
• /
• 제23권5호
• /
• pp.475-482
• /
• 2010

본 논문에서는 대규모 자유도 시스템의 병렬처리를 위하여 2단계로 이루어진 영역분할법(Domain Decomposition Method) 기반의 병렬 알고리즘을 제안하였다. 분할된 영역의 내부 및 외부 경계를 상위영역문제로 정의하고 국부영역문제는 변위 경계조건이 모두 주어지는 분할영역에서의 Dirichlet 문제로 구성한다. 상위영역에서는 전체 상위영역에 대한 강성 행렬의 어셈블이 필요없는 반복법을 통하여 변위를 구하고, 이를 바탕으로 국부영역에서 Multi-Frontal Sparse Solver (MFSS)를 이용하여 변위를 계산한다. 상위영역문제의 연산에서 프로세서 간의 데이터 교환을 최소화하여 계산효율을 유지하며, 동시에 해석 가능한 자유도를 증대시키는 병렬 PCG(Preconditioned Conjugate Gradient)법 기반의 알고리즘을 개발하였다. 제안된 알고리즘을 적용하여 수치해석을 수행한 결과, 프로세서 수가 증가할수록 계산성능의 손실없이 해석 가능한 자유도가 비례하여 증가하는 선형 확장성을 관찰할 수 있었으며, 대규모 자유도 문제에 효과적으로 사용 가능함을 확인하였다.