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

A domain/boundary decomposition technique is applied to carry out efficient finite element analyses of 3-D contact problems. Appropriate penalty functions are selected for connecting an interface and contact interfaces with neighboring subdomains that satisfy continuity constraints. As a consequence, all the effective stiffness matrices have positive definiteness, and computational efficiency can be improved to a considerable degree. If necessary, any complex-shaped 3-D domain can be divided into several simple-shaped subdomains without considering the conformity of meshes along the interface. With a set of numerical examples, the basic characteristics of computational efficiency are investigated carefully.

• Nuclear Engineering and Technology
• /
• 제52권11호
• /
• pp.2667-2677
• /
• 2020

A domain decomposition (DD) scheme for GPU-based Monte Carlo (MC) calculation which is essential for whole-core depletion is introduced within the framework of the modified history-based tracking algorithm. Since GPU-offloaded MC calculations suffer from limited memory capacity, employing DDMC is inevitable for the simulation of depleted cores which require large storage to save hundreds of newly generated isotopes. First, an automated domain decomposition algorithm named wheel clustering is devised such that each subdomain contains nearly the same number of fuel assemblies. Second, an innerouter iteration algorithm allowing overlapped computation and communication is introduced which enables boundary neutron transactions during the tracking of interior neutrons. Third, a bank update scheme which is to include the boundary sources in a way to be adequate to the peculiar data structures of the GPU-based neutron tracking algorithm is presented. The verification and demonstration of the DDMC method are done for 3D full-core problems: APR1400 fresh core and a mock-up depleted core. It is confirmed that the DDMC method performs comparably with the standard MC method, and that the domain decomposition scheme is essential to carry out full 3D MC depletion calculations with limited GPU memory capacities.

• 한국전산구조공학회논문집
• /
• 제22권1호
• /
• pp.35-44
• /
• 2009

본 논문에서는 대규모 3차원 구조해석에 필요한 병렬 유한요소해석을 위한 영역분할법의 적용에 대해 묘사하였다. 영역분할법을 사용한 병렬 유한요소법 시스템을 개발하였다. 절점 생성시, 절점들간의 거리가 특정절점에서의 공간함수와 같아지면 절점이 생성되어 진다. 이 절점공간함수는 퍼지지식처리에 의해 조절되어 진다. 기본적인 요소생성은 데로우니 삼각화 기법을 적용하였다. 자동요소생성 시스템을 이용한 영역분할법은 3차원 해석에 큰 도움이 된다. 공간함수와 유사하게 절점들간의 유한요소해석을 위한 병렬 수치 알고리즘으로서 영역분할법을 전체의 해석영역을 완전히 여러 개의 작은 영역으로 겹치지 않게 나누는 공역구배인 반복적 솔버와 결합시켰다. 개발된 시스템의 효용성에 대한 성능을 몇 가지 예를 통해 제시하였다.

• 한국전산구조공학회논문집
• /
• 제20권4호
• /
• pp.469-476
• /
• 2007

많은 계산량이 요구되는 삼차원 접촉 문제의 효율적인 유한요소 해석을 위하여 영역/경계 분할 기법을 적용하였다. 접촉 경계면의 부등식 적합 조건과 부영역, 공유면, 접촉 공유면의 등식 적합 조건을 모두 벌칙 함수로 처리하였다. 이에 따라 모든 유효 강성 행렬이 양 정치화되므로, 역행렬과 같은 각종 행렬 연산이 매우 간편해진다. 또한 전체 영역의 형상이 복잡하더라도, 임의의 부영역, 공유면, 접촉 공유면 단위로 쉽게 유한요소 모델링할 수 있다. 즉, 관련 지배 방정식은 물론 경계 조건도 독립적으로 이산화할 수 있으므로, 국부적인 비선형 접촉 조건에 대한 효율적인 해석이 가능하다. 간단한 수치 예제를 통하여 삼차원 접촉 해석의 효율성에 관한 기본적인 경향을 검토하였다.

• 대한수학회지
• /
• 제61권4호
• /
• pp.659-681
• /
• 2024

We propose V-cycle multigrid methods for vector field problems arising from the lowest order hexahedral Nédélec finite element. Since the conventional scalar smoothing techniques do not work well for the problems, a new type of smoothing method is necessary. We introduce new smoothers based on substructuring with nonoverlapping domain decomposition methods. We provide the convergence analysis and numerical experiments that support our theory.

• Nuclear Engineering and Technology
• /
• 제48권3호
• /
• pp.635-641
• /
• 2016

Because of prohibitive data storage requirements in large-scale simulations, the memory problem is an obstacle for Monte Carlo (MC) codes in accomplishing pin-wise three-dimensional (3D) full-core calculations, particularly for whole-core depletion analyses. Various kinds of data are evaluated and quantificational total memory requirements are analyzed based on the Reactor Monte Carlo (RMC) code, showing that tally data, material data, and isotope densities in depletion are three major parts of memory storage. The domain decomposition method is investigated as a means of saving memory, by dividing spatial geometry into domains that are simulated separately by parallel processors. For the validity of particle tracking during transport simulations, particles need to be communicated between domains. In consideration of efficiency, an asynchronous particle communication algorithm is designed and implemented. Furthermore, we couple the domain decomposition method with MC burnup process, under a strategy of utilizing consistent domain partition in both transport and depletion modules. A numerical test of 3D full-core burnup calculations is carried out, indicating that the RMC code, with the domain decomposition method, is capable of pin-wise full-core burnup calculations with millions of depletion regions.

• 한국전자파학회논문지
• /
• 제25권3호
• /
• pp.357-364
• /
• 2014

본 연구에서는 빠르게 전자파 해석을 수행할 수 있는 병렬 유한차분 시간영역(Finite-Difference Time-Domain: FDTD) 알고리즘을 구현하기 위하여 CPU 클러스터를 구축하였다. 병렬 FDTD 알고리즘은 단일 프로세서를 이용한 FDTD 알고리즘에 비해 해석 시간을 크게 줄일 수 있으며, 전기적으로 매우 큰 구조물에 대한 전자파 해석도 가능하다. 본 연구팀에서는 CPU 클러스터 기반의 병렬 FDTD 알고리즘에서 요구되는 프로세스 간의 통신을 위해 MPI(Message Passing Interface) 라이브러리를 이용하였으며, 3차원 공간분할을 적용하여 프로세스 간의 통신 시간을 최소화하였다. 단일 프로세서를 이용한 FDTD 알고리즘 대비 CPU 클러스터 기반의 병렬 FDTD 알고리즘의 계산속도 향상도를 기본 모드와 하이퍼 모드에서 분석하였으며, 전기적으로 매우 큰 콘크리트 구조물의 전자파 해석을 하였다.

• Wind and Structures
• /
• 제5권2_3_4호
• /
• pp.379-392
• /
• 2002

Numerical flow computations around an aeroelastic 3D square cylinder immersed in the turbulent boundary layer are shown. Present computational code can be characterized by three numerical aspects which are 1) the method of artificial compressibility is adopted for the incompressible flow computations, 2) the domain decomposition technique is used to get better grid point distributions, and 3) to achieve the conservation law both in time and space when the flow is computed a with moving and transformed grid, the time derivatives of metrics are evaluated using the time-and-space volume. To provide time-dependant inflow boundary conditions satisfying prescribed time-averaged velocity profiles, a convenient way for generating inflow turbulence is proposed. The square cylinder is modeled as a 4-lumped-mass system and it vibrates with two-degree of freedom of heaving motion. Those blocks which surround the cylinder are deformed according to the cylinder's motion. Vigorous oscillations occur as the vortex shedding frequency approaches cylinder's natural frequencies.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집A
• /
• pp.894-899
• /
• 2000

An automatic mesh generation system with unstructured quadrilateral elements on trimmed NURBS surfaces has been developed.. In this paper, NURBS surface geometries in the IGES format have been used to represent model shape. NURBS surface is represented as parametric surface. So each surface could be mapped to a 2D parametric plane through the parametric domain. And then meshes with quadrilateral elements are constructed in this plane. Finally, the constructed meshes are mapped back to the original 3D surface through the parametric domain. In this paper, projection plane, quasi-expanded plane and parametric Plane are used as 2D mesh generation plane. For mapping 3D surface to parametric domain, Newton-Rhapson Method is employed. For unstructured mesh generation with quadrilateral elements on 2D plane, a domain decomposition algorithm using loop operators has been employed. Sample meshes are represented to demonstrate the effectiveness of the proposed algorithm.

• 소성∙가공
• /
• 제9권4호
• /
• pp.362-371
• /
• 2000

In order to reduce weight of a high-speed railroad vehicle, the main body has been manufactured by hollow section extrusion using aluminum alloys. A porthole die has utilized for the hollow section extrusion process, which causes complicated die geometry and flow characteristics. Design of porthole die is very difficult due to such a complexity. The three-dimensional finite element analysis for hollow section is also an arduous job from the viewpoint of appropriate mesh construction and tremendous computation time. In the present work, mismatching refinement, an efficient domain decomposition method with different mesh density for each subdomain, is implemented for the analysis of the hollow section extrusion process. In addition, a modified grid-based approach with the surface element layer is utilized lot three-dimensional mesh generation of a complicated shape with hexahedral elements. The effects of porthole design are discussed through the simulation for extrusion of an underframe part of a railroad vehicle. An experiment has also been carried out for the comparison. Comparing the velocity distribution at the outlet with the thickness variation of the extruded part, it is concluded that the analysis results can provide reliable measures whether the die design is acceptable to obtain uniform part thickness. The analysis results are then successfully reflected on the industrial porthole die design.