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

A weakly nonlinear seakeeping methodology for predicting motions and loads is presented in this paper. This methodology assumes linear radiation and diffraction forces, calculated in the frequency domain, and fully nonlinear Froude-Krylov and hydrostatic forces, evaluated in the time domain. The particular approach employed here allows to overcome numerical problems connected to the determination of the impulse response functions. The procedure is divided into three consecutive steps: evaluation of dynamic sinkage and trim in calm water that can significantly influence the final results, a linear seakeeping analysis in the frequency domain and a weakly nonlinear simulation. The first two steps are performed employing a three-dimensional Rankine panel method. Nonlinear Froude-Krylov and hydrostatic forces are computed in the time domain by pressure integration on the actual wetted surface at each time step. Although nonlinear forces are evaluated into the time domain, the equations of motion are solved in the frequency domain iteratively passing from the frequency to the time domain until convergence. The containership S175 is employed as a test case for evaluating the capability of this methodology to correctly predict the nonlinear behavior related to wave induced motions and loads in head seas; numerical results are compared with experimental data provided in literature.

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.38 no.12
• /
• pp.56-65
• /
• 2001

This paper proposes an algorithm for the capacitance extraction of general 3-dimensional conductors in an ideal uniform dielectric that uses a high-order quadrature approximation method combined with the typical first-order collocation method to enhance the accuracy and adopts an efficient matrix-vector product algorithm for the model-order reduction to achieve efficiency. The proposed method enhances the accuracy using the quadrature method for interconnects containing corners and vias that concentrate the charge density. It also achieves the efficiency by reducing the model order using the fact that large parts of system matrices are of numerically low rank. This technique combines an SVD-based algorithm for the compression of rank-deficient matrices and Gram-Schmidt algorithm of a Krylov-subspace iterative technique for the rapid multiplication of matrices. It is shown through the performance evaluation procedure that the combination of these two techniques leads to a more efficient algorithm than Gaussian elimination or other standard iterative schemes within a given error tolerance.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.9
• /
• pp.933-941
• /
• 2014

This paper provides a comparison between the Krylov subspace method (KSM) and modal truncation method (MTM), which are typical projection-based model order reduction methods. The frequency responses are compared to determine the numerical accuracies and efficiencies. In order to compare the numerical accuracies of the KSM and MTM, the frequency responses and relative errors according to the order of the reduced model and frequency of interest are studied. Subsequently, a numerical examination shows whether a reduced order can be determined automatically with the help of an error convergence indicator. As for the numerical efficiency, the computation time needed to generate the projection matrix and the solution time to perform a frequency response analysis are compared according to the reduced order. A finite element model for a car suspension is considered as an application example of the numerical comparison.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.16 no.3
• /
• pp.533-540
• /
• 2021

To evaluate the performance of a system, one-dimensional 3-neighborhood cellular automata(CA) based pseudo-random generators are widely used in many fields. Although two-dimensional CA and one-dimensional 5-neighborhood CA have been applied for more effective key sequence generation, designing symmetric one-dimensional 5-neighborhood CA corresponding to a given primitive polynomial is a very challenging problem. To solve this problem, studies on one-dimensional 5-neighborhood CA synthesis, such as synthesis method using recurrence relation of characteristic polynomials and synthesis method using Krylov matrix, were conducted. However, there was still a problem with solving nonlinear equations. To solve this problem, a symmetric one-dimensional 5-neighborhood CA synthesis method using a transition matrix of 90/150 CA and a block matrix has recently been proposed. In this paper, we detail the theoretical process of the proposed algorithm and use it to obtain symmetric one-dimensional 5-neighborhood CA corresponding to high-order primitive polynomials.

• Journal of applied mathematics & informatics
• /
• v.8 no.2
• /
• pp.305-326
• /
• 2001

We introduce a class of multilevel recursive incomplete LU preconditioning techniques (RILUM) for solving general sparse matrices. This techniques is based on a recursive two by two block incomplete LU factorization on the coefficient martix. The coarse level system is constructed as an (approximate) Schur complement. A dynamic preconditioner is obtained by solving the Schur complement matrix approximately. The novelty of the proposed techniques is to solve the Schur complement matrix by a preconditioned Krylov subspace method. Such a reduction process is repeated to yield a multilevel recursive preconditioner.

• Communications of the Korean Mathematical Society
• /
• v.26 no.4
• /
• pp.695-708
• /
• 2011

Based on the multiple-time-scale (MTS) method, a general formula has been presented for solving an n-th, n = 2, 3, ${\ldots}$, order ordinary differential equation with strong linear damping forces. Like the solution of the unified Krylov-Bogoliubov-Mitropolskii (KBM) method or the general Struble's method, the new solution covers the un-damped, under-damped and over-damped cases. The solutions are identical to those obtained by the unified KBM method and the general Struble's method. The technique is a new form of the classical MTS method. The formulation as well as the determination of the solution from the derived formula is very simple. The method is illustrated by several examples. The general MTS solution reduces to its classical form when the real parts of eigen-values of the unperturbed equation vanish.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.283-290
• /
• 1998

A solution method is presented to solve the eigenproblem arising in tile dynamic analysis of non-proportional damping systems with symmetric matrices. The method is based on tile use of Lanczos method to generate a Krylov subspace of trial vectors, witch is then used to reduce a large eigenvalue problem to a much smaller one. The method retains the η order quadratic eigenproblem, without the need to the method of matrix augmentation traditionally used to cast the problem as a linear eigenproblem of order 2n. In the process, the method preserves tile sparseness and symmetry of the system matrices and does not invoke complex arithmetics, therefore, making it very economical for use in solving large problems. Numerical results are presented to demonstrate the efficiency and accuracy of the method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.223-231
• /
• 2013

This paper presents a time-domain Gauss-Newton full-waveform inversion method for the material profile reconstruction in heterogeneous semi-infinite solid media. To implement the inverse problem in a finite computational domain, perfectly-matchedlayers( PMLs) are introduced as wave-absorbing boundaries within which the domain's wave velocity profile is to be reconstructed. The inverse problem is formulated in a partial-differential-equations(PDE)-constrained optimization framework, where a least-squares misfit between measured and calculated surface responses is minimized under the constraint of PML-endowed wave equations. A Gauss-Newton-Krylov optimization algorithm is utilized to iteratively update the unknown wave velocity profile with the aid of a specialized regularization scheme. Through a series of one-dimensional examples, the solution of the Gauss-Newton inversion was close enough to the target profile, and showed superior convergence behavior with reduced wall-clock time of implementation compared to a conventional inversion using Fletcher-Reeves optimization algorithm.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.11
• /
• pp.1852-1861
• /
• 1997

An enhanced subspace iteration algorithm has been developed to solve eigenvalue problems reliably and efficiently. Basic subspace iteration algorithm has been improved by eliminating recalculation of converged eigenvectors, using Krylov sequence as initial vectors and incorporating with shifting techniques. The number of iterations and computational time have been considerably reduced when compared with the original one, and reliability for catching copies of the multiple roots has been retained successfully. Further research would be required for mathematical justification of the present method.

• Journal of applied mathematics & informatics
• /
• v.36 no.5_6
• /
• pp.459-474
• /
• 2018

Cao et al. in (Numer. Linear. Algebra Appl. 18 (2011) 875-895) proposed a splitting method for saddle point problems which unconditionally converges to the solution of the system. It was shown that a Krylov subspace method like GMRES in conjunction with the induced preconditioner is very effective for the saddle point problems. In this paper we first modify the iterative method, discuss its convergence properties and apply the induced preconditioner to the problem. Numerical experiments of the corresponding preconditioner are compared to the primitive one to show the superiority of our method.