• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.2
• /
• pp.117-124
• /
• 2019

Parallel sparse solvers are essential for solving large-scale finite element models. This paper introduces the combination of iterative and direct solver that can be applied efficiently to problems that require continuous solution for a subtly changing sequence of systems of equations. The iterative-direct sparse solver combination technique, proposed and implemented in the parallel sparse solver package, PARDISO, means that iterative sparse solver is applied for the newly updated linear system, but it uses the direct sparse solver's factorization of previous system matrix as a preconditioner. If the solution does not converge until the preset iterations, the solution will be sought by the direct sparse solver, and the last factorization results will be used as a preconditioner for subsequent updated system of equations. In this study, an improved method that sets the maximum number of iterations dynamically at the first Krylov iteration step is proposed and verified thereby enhancing calculation efficiency by the frequency domain analysis.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.3 s.246
• /
• pp.310-317
• /
• 2006

Currently most practical vibration and structural problems in automotive suspensions require the use of the finite element method to obtain their structural responses. When the finite element model has a very large number of degrees of freedom the harmonic and dynamic analyses are computationally too expensive to repeat within a feasible design process time. To alleviate the computational difficulty, this paper presents a moment-matching based model order reduction (MOR) which reduces the number of degrees of freedom of the original finite element model and speeds up the necessary simulations with the reduced-size models. The moment-matching model reduction via the Arnoldi process is performed directly to ANSYS finite element models by software mor4ansys. Among automotive suspension components, a knuckle is taken as an example to demonstrate the advantages of this approach for vibration simulation. The frequency and transient dynamic responses by the MOR are compared with those by the mode superposition method.

• Journal of the Korean Institute of Navigation
• /
• v.11 no.2
• /
• pp.61-72
• /
• 1987

When a ship is running in following seas, the encounter frequency is reduced to a very low one. In that case broaching, surfiding and capsizing phenomena are most likely to occur due to wave exciting forces acting on a ship in following seas. In this paper, the emphasis is mainly laid upon transverse stability of ships following seas, which is related to capszing phenomenon. The authors take the case that ship speed is equal to the wave celerity, i.e., the encounter frequency is zero. Hydrostatic force and moment due to Froude-Krylov hypothesis are calculated by line intergral method. Transverse stability is evaluated from hydrostatic force and moment. Through the application of present calculation method to box-shaped vessel, it is confirmed that the transversestability of a vessel can be reduced to critical level at wave crest.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2001.10a
• /
• pp.164-168
• /
• 2001

The motion of a ship advancing in regular waves is analyzed in the time-domain using the convolution integral of the radiation forces. The memory effect functions and infinite frequency added masses are obtained from the solution of the three dimensional improved Green integral equation in the frequency domain by making use of the Fourier transformation. The ship motions in regular waves have been calculated by both the time and frequency domain methods. It has been shown that they agree very well with each other. The present time-domain method can be used to predict the time histories of unsteady motions in irregular waves. It can also be used to calculate the hydrostatic and Froude-Krylov forces over the instantaneous wetted surface of the ship hull to predict large ship motions, in a practical sense, advancing in large amplitude waves.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1181-1188
• /
• 2017

This paper proposes the adoption of Monte Carlo perturbation theory to approximate the Jacobian matrix of coupled neutronics/thermal-hydraulics problems. The projected Jacobian is obtained from the eigenvalue decomposition of the fission matrix, and it is adopted to solve the coupled problem via the Newton method. This avoids numerical differentiations commonly adopted in Jacobian-free Newton-Krylov methods that tend to become expensive and inaccurate in the presence of Monte Carlo statistical errors in the residual. The proposed approach is presented and preliminarily demonstrated for a simple two-dimensional pressurized water reactor case study.

• The Transactions of the Korean Institute of Electrical Engineers A
• /
• v.49 no.4
• /
• pp.168-176
• /
• 2000

This paper presents the Hessenberg method, a new sparsity-based small signal stability analysis program for large interconnected power systems. The Hessenberg method as well as the Arnoldi method computes the partial eigen-solution of large systems. However, the Hessenberg method with pivoting is numerically very stable comparable to the Householder method and thus re-orthogonalization of the krylov vectors is not required. The fractional transformation with a complex shift is used to compute the modes around the shift point. If only the dominant electromechanical oscillation modes are of concern, the modes can be computed fast with the shift point determined by Fourier transforming the time simulation results for transient stability analysis, if available. The program has been successfully tested on the New England 10-machine 39-bus system and Korea Electric Power Co. (KEPCO) system in the year of 2000, which is comprised of 791-bus, 1575-branch, and 215-machines. The method is so efficient that CPU time for computing five eigenvalues of the KEPCO system is 3.4 sec by a PC with 400 MHz Pentium IIprocessor.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.3
• /
• pp.177-183
• /
• 2018

This study simulated ice load and the motion response of a moored semi-submersible rig in pack-ice conditions using a finite element method. Ice flows of random size and shape were modeled, and interactions for ice-sea, ice-structure, ice-ice were simulated using a simplified method. Parameters for the simplified method such as drag force coefficient and the pressure-penetration relation were obtained based on the result of detailed analysis using the coupled Eulerian-Lagrangian method. The mooring lines were modeled by spring elements based on their stiffness. As a result of the simulation over 1,400 seconds, the force and motion response of the rig were obtained and validated using discrete elements and compared with the results found by the Krylov State Research Centre.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.1
• /
• pp.53-64
• /
• 2011

The present study considers multi-level approach for the analysis of parametric roll phenomena. Three kinds of computation method, GM variation, impulse response function (IRF), and Rankine panel method, are applied for the multi-level approach. IRF and Rankine panel method are based on the weakly nonlinear formulation which includes nonlinear Froude-Krylov and restoring forces. In the computation result of parametric roll occurrence test in regular waves, IRF and Rankine panel method show similar tendency. Although the GM variation approach predicts the occurrence of parametric roll at twice roll natural frequency, its frequency criteria shows a little difference. Nonlinear roll motion in bichromatic wave is also considered in this study. To prove the unstable roll motion in bichromatic waves, theoretical and numerical approaches are applied. The occurrence of parametric roll is theoretically examined by introducing the quasi-periodic Mathieu equation. Instability criteria are well predicted from stability analysis in theoretical approach. From the Fourier analysis, it has been verified that difference-frequency effects create the unstable roll motion. The occurrence of unstable roll motion in bichromatic wave is also observed in the experiment.

• Journal of applied mathematics & informatics
• /
• v.17 no.1_2_3
• /
• pp.457-473
• /
• 2005

In each step of the quasi-minimal residual (QMR) method which uses a look-ahead variant of the nonsymmetric Lanczos process to generate basis vectors for the Krylov subspaces induced by A, it is necessary to decide whether to construct the Lanczos vectors $v_{n+l}\;and\;w{n+l}$ as regular or inner vectors. For a regular step it is necessary that $D_k\;=\;W^{T}_{k}V_{k}$ is nonsingular. Therefore, in the floating-point arithmetic, the smallest singular value of matrix $D_k$, ${\sigma}_min(D_k)$, is computed and an inner step is performed if $\sigma_{min}(D_k)<{\epsilon}$, where $\epsilon$ is a suitably chosen tolerance. In practice it is absolutely impossible to choose correctly the value of the tolerance $\epsilon$. The subject of this paper is to show how discrete stochastic arithmetic remedies the problem of this tolerance, as well as the problem of the other tolerances which are needed in the other checks of the QMR method with the estimation of the accuracy of some intermediate results. Numerical examples are used to show the good numerical properties.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.148-154
• /
• 2004

This article reviews the development of three-dimensional (3-D) magnetotelluric (MT) modeling. The 3-D modeling of electromagnetic fields is essential in understanding the physics of MT soundings, and in implementing an inversion method to reconstruct a 3-D resistivity image. Although various numerical schemes have been developed over the last two decades, practical methods have been quite limited. However, the recent rapid improvement in computer speed and memory, as well as the advance in iterative solution algorithms for a large system of equations, makes it possible to model the MT responses of complex 3-D structures, which have been very difficult to simulate before. The use of staggered grids in finite difference method has become popular, conserving a magnetic flux and an electric current and allowing for realistic discontinuous fields. The convergence of numerical solutions has been greatly accelerated by adopting Krylov subspace methods, proper preconditioning techniques, and static divergence corrections. The vector finite-element method using edge elements is also free from the discontinuity problem, and seems a natural choice for modeling complex structures including irregular topography because its flexibility allows one to capture full geometric complexity.