• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.2
• /
• pp.127-135
• /
• 2017

The computing environment has changed rapidly to enable large-scale finite element models to be analyzed at the PC or workstation level, such as multi-core CPU, optimal math kernel library implementing BLAS and LAPACK, and popularization of direct sparse solvers. In this paper, the design considerations on a parallel finite element code for shared memory based multi-core CPU system are proposed; (1) the use of optimized numerical libraries, (2) the use of latest direct sparse solvers, (3) parallelism using OpenMP for computing element stiffness matrices, and (4) assembly techniques using triplets, which is a type of sparse matrix storage. In addition, the parallelization effect is examined on the time-consuming works through a large scale finite element model.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.147-153
• /
• 2007

The conjugate gradient (CG) method is one of the most efficient algorithms for solving a linear system of equations. In addition to being used as a linear equation solver, it can be applied to a least-squares problem. When the CG method is applied to large-scale three-dimensional inversion of magnetotelluric data, two approaches have been pursued; one is the linear CG inversion in which each step of the Gauss-Newton iteration is incompletely solved using a truncated CG technique, and the other is referred to as the nonlinear CG inversion in which CG is directly applied to the minimization of objective functional for a nonlinear inverse problem. In each procedure we only need to compute the effect of the sensitivity matrix or its transpose multiplying an arbitrary vector, significantly reducing the computational requirements needed to do large-scale inversion.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.9
• /
• pp.775-780
• /
• 2016

In this paper, a computational algorithm of an improved and versatile structural analysis applicable for large-size flexible nonlinear structures is developed. In more detail, nonlinear finite element based on the co-rotational (CR) framework is developed. Then, a finite element tearing and interconnecting method using local Lagrange multipliers (FETI-local) is combined with the nonlinear CR finite element. The resulting computational algorithm is presented and applied for nonlinear static analyses, i.e., cantilevered beam and multibody structure. Finally, the proposed analysis is evaluated with regard to its parallel computation performance, and it is compared with those obtained by serial computation using the sparse matrix linear solver, PARDISO.

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.119-130
• /
• 2013

We developed a method for inverting magnetic data to recover the 3D susceptibility models. The major difficulty in the inversion of the potential data is the non-uniqueness and the vast computing time. The insufficient number of data compared with that of inversion blocks intensifies the non-uniqueness problem. Furthermore, there is poor depth resolution inherent in magnetic data. To overcome this non-uniqueness problem, we propose a resolution model constraint that imposes large penalty on the model parameter with good resolution; on the other hand, small penalty on the model parameter with poor resolution. Using this model constraint, the model parameter with a poor resolution can be effectively resolved. Moreover, the wavelet transform and parallel solving were introduced to save the computing time. Through the wavelet transform, a large system matrix was transformed to a sparse matrix and solved by a parallel linear equation solver. This procedure is able to enormously save the computing time for the 3D inversion of magnetic data. The developed inversion algorithm is applied to the inversion of the synthetic data for typical models of magnetic anomalies and real airborne data obtained at the Geumsan area of Korea.

• Journal of Korea Soil Environment Society
• /
• v.1 no.1
• /
• pp.89-102
• /
• 1996

An integrated model is presented to describe underground flow and mass transport, using a multicomponent multiphase approach. The comprehensive governing equation is derived considering mass and force balances of chemical species over four phases(water, oil, air, and soil) in a schematic elementary volume. Compact and systemati notations of relevant variables and equations are introduced to facilitate the inclusion of complex migration and transformation processes, and variable spatial dimensions. The resulting nonlinear system is solved by a multidimensional finite element code. The developed code with dynamic array allocation, is sufficiently flexible to work across a wide spectrum of computers, including an IBM ES 9000/900 vector facility, SP2 cluster machine, Unix workstations and PCs, for one-, two and three-dimensional problems. To reduce the computation time and storage requirements, the system equations are decoupled and solved using a banded global matrix solver, with the vector and parallel processing on the IBM 9000. To avoide the numerical oscillations of the nonlinear problems in the case of convective dominant transport, the techniques of upstream weighting, mass lumping, and elementary-wise parameter evaluation are applied. The instability and convergence criteria of the nonlinear problems are studied for the one-dimensional analogue of FEM and FDM. Modeling capacity is presented in the simulation of three dimensional composite multiphase TCE migration. Comprehesive simulation feature of the code is presented in a companion paper of this issue for the specific groundwater or flow and contamination problems.

• Composites Research
• /
• v.36 no.5
• /
• pp.303-309
• /
• 2023

In this paper, a Progressive Damage and Failure Analysis (PDFA) modeling method was developed using ABAQUS/EXPLICIT to predict in-plane damage and delamination for Open-Hole Compression (OHC) testing. The proposed PDFA model was constructed based on Hashin criteria and cohesive behavior. The strength and stiffness of OHC specimens with three types of stacking sequences [(45/-45/0₂)₃]_s , [(45/0/-45/90)₃]_s and [45/-45/0/45/-45/90/(45/-45)₂]_s were compared to comprehensively evaluate the validity of the Finite Element(FE) model of PDFA. The strength and stiffness of the OHC specimens were predicted relatively well, with less than a percentage error 10.0 %. For the numerical simulation case for each layup, the damage initiation/evolution of OHC specimens were evaluated for delamination and tension/compression matrix damage before and after failure.