• Coupled systems mechanics
• /
• v.1 no.1
• /
• pp.19-37
• /
• 2012

In this work, the iterative coupling of finite element and boundary element methods for the investigation of coupled fluid-fluid, solid-solid and fluid-solid wave propagation models is reviewed. In order to perform the coupling of the two numerical methods, a successive renewal of the variables on the common interface between the two sub-domains is performed through an iterative procedure until convergence is achieved. In the case of local nonlinearities within the finite element sub-domain, it is straightforward to perform the iterative coupling together with the iterations needed to solve the nonlinear system. In particular, a more efficient and stable performance of the coupling procedure is achieved by a special formulation that allows to use different time steps in each sub-domain. Optimized relaxation parameters are also considered in the analyses, in order to speed up and/or to ensure the convergence of the iterative process.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• 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.

• Steel and Composite Structures
• /
• v.10 no.6
• /
• pp.541-561
• /
• 2010

This paper presents finite element analyses, experimental measurements and finite element model updating of an arch type steel laboratory bridge model using semi-rigid connections. The laboratory bridge model is a single span and fixed base structure with a length of 6.1 m and width of 1.1m. The height of the bridge column is 0.85 m and the maximum arch height is 0.95 m. Firstly, a finite element model of the bridge is created in SAP2000 program and analytical dynamic characteristics such as natural frequencies and mode shapes are determined. Then, experimental measurements using ambient vibration tests are performed and dynamic characteristics (natural frequencies, mode shapes and damping ratios) are obtained. Ambient vibration tests are performed under natural excitations such as wind and small impact effects. The Enhanced Frequency Domain Decomposition method in the frequency domain and the Stochastic Subspace Identification method in the time domain are used to extract the dynamic characteristics. Then the finite element model of the bridge is updated using linear elastic rotational springs in the supports and structural element connections to minimize the differences between analytically and experimentally estimated dynamic characteristics. At the end of the study, maximum differences in the natural frequencies are reduced on average from 47% to 2.6%. It is seen that there is a good agreement between analytical and experimental results after finite element model updating. Also, connection percentages of the all structural elements to joints are determined depending on the rotational spring stiffness.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.33-51
• /
• 2008

This paper investigates the modelling of coupled soil-structure interaction problems by domain decomposition techniques. It is assumed that the soil-structure system is physically partitioned into soil and structure subdomains, which are independently modelled. Coupling of the separately modelled partitioned subdomains is undertaken with various algorithms based on the sequential iterative Dirichlet-Neumann sub-structuring method, which ensures compatibility and equilibrium at the interface boundaries of the subdomains. A number of mathematical and computational characteristics of the coupling algorithms, including the convergence conditions and choice of algorithmic parameters leading to enhanced convergence of the iterative method, are discussed. Based on the presented coupling algorithms a simulation environment, utilizing discipline-oriented solvers for nonlinear structural and geotechnical analysis, is developed which is used here to demonstrate the performance characteristics and benefits of various algorithms. Finally, the developed tool is used in a case study involving nonlinear soil-structure interaction analysis between a plane frame and soil subjected to ground excavation. This study highlights the relative performance of the various considered coupling algorithms in modelling real soil-structure interaction problems, in which nonlinearity arises in both the structure and the soil, and leads to important conclusions regarding their adequacy for such problems as well as the prospects for further enhancements.

• Nuclear Engineering and Technology
• /
• v.46 no.5
• /
• pp.655-666
• /
• 2014

The CUPID code has been developed at KAERI for a transient, three-dimensional analysis of a two-phase flow in light water nuclear reactor components. It can provide both a component-scale and a CFD-scale simulation by using a porous media or an open media model for a two-phase flow. In this paper, recent advances in the CUPID code are presented in three sections. First, the domain decomposition parallel method implemented in the CUPID code is described with the parallel efficiency test for multiple processors. Then, the coupling of CUPID-MARS via heat structure is introduced, where CUPID has been coupled with a system-scale thermal-hydraulics code, MARS, through the heat structure. The coupled code has been applied to a multi-scale thermal-hydraulic analysis of a pool mixing test. Finally, CUPID-SG is developed for analyzing two-phase flows in PWR steam generators. Physical models and validation results of CUPID-SG are discussed.

• Current Applied Physics
• /
• v.18 no.11
• /
• pp.1327-1337
• /
• 2018

This paper investigates the vibration behavior of pristine and defected triangular graphene sheets; which has recently attracted the attention of researchers and compare these two types in natural frequencies and sensitivity. Here, the molecular dynamics method has been employed to establish a virtual laboratory for this purpose. After measuring the different parameters obtained by the molecular dynamics approach, these data have been analyzed by using the frequency domain decomposition (FDD) method, and the dominant frequencies and mode shapes of the system have been extracted. By analyzing the vibration behaviors of pristine triangular graphene sheets in four cases (right angle of 45-90-45 configuration, right angle of 60-90-30 configuration, equilateral triangle and isosceles triangle), it has been demonstrated that the natural frequencies of these sheets are higher than the natural frequency of a square sheet, with the same number of atoms, by a minimum of 7.6% and maximum of 26.6%. Therefore, for increasing the resonance range of sensors based on 2D materials, nonrectangular structures, and especially the triangular structure, can be considered as viable candidates. Although the pristine and defective equilateral triangular sheets have the highest values of resonance, the sensitivity of defective (45,90,45) triangular sheet is more than other configurations and then, defective (45,90,45) sheet is the worst choice for sensor applications.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.49-54
• /
• 1998

3-dimensional Euler solver is parallelized. The spatial discretization method is the 2nd order TVD scheme and DADI method with multigrid is used as a time integration. In order to parallelize this solver, the domain decomposition method with overlapped grid and message passing techniques are used. The informations on the each inter-processor bound-aries are communicated with MPI library. Finally, the parallel performance repsented by calculating the ONERA M6 wing at transonic flow condition using CRAY T3E and C90.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• 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.

• Geomechanics and Engineering
• /
• v.7 no.2
• /
• pp.149-164
• /
• 2014

One of the most advanced classes of techniques for ground response analysis is based on the use of Transfer Functions. They represent the ratio of Fourier spectrum of amplitude motion at the free surface to the corresponding spectrum of the bedrock motion and they are applied in frequency domain usually by FFT method. However, Fourier spectrum only shows the dominant frequency in each time step and is unable to represent all frequency contents in every time step and this drawback leads to inaccurate results. In this research, this process is optimized by decomposing the input motion into different frequency sub-bands using Wavelet Multi-level Decomposition. Each component is then processed with transfer Function relating to the corresponding component frequency. Taking inverse FFT from all components, the ground motion can be recovered by summing up the results. The nonlinear behavior is approximated using an iterative procedure with nonlinear soil properties. The results of this procedure show better accuracy with respect to field observations than does the Conventional method. The proposed method can also be applied to other engineering disciplines with similar procedure.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.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.