• Structural Engineering and Mechanics
• /
• 제76권5호
• /
• pp.643-651
• /
• 2020

Generally, it is necessary to perform transient structural analysis in order to verify and improve the seismic performance of high-rise buildings and bridges against earthquake loads. In this paper, we propose the model order reduction (MOR) method using the Krylov vectors to perform seismic analysis for linear and elastic systems in an efficient way. We then compared the proposed method with the mode superposition method (MSM) by using the limited numbers of modal vectors (or eigenvectors) calculated from the modal analysis. In the calculation, the data of the El Centro earthquake in 1940 were adopted for the seismic loading in the transient analysis. The numerical accuracy and efficiency of the two methods were compared in detail in the case of a simplified high-rise building.

• Journal of applied mathematics & informatics
• /
• 제25권1_2호
• /
• pp.293-304
• /
• 2007

In this paper, we present a Krylov subspace based projection method for reduced-order modeling of large scale bilinear multi-input multi-output (MIMO) systems. The reduced-order bilinear system is constructed in such a way that it can match a desired number of moments of multi-variable transfer functions corresponding to the kernels of Volterra series representation of the original system. Numerical examples report the effectiveness of this method.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.

• 전자공학회논문지
• /
• 제50권11호
• /
• pp.28-35
• /
• 2013

Krylov-Schur 반복법을 활용하여 2-차원 사각 도파관에서 나타나는 고유특성을 밝혔다. 고유 행렬 방정식은 삼각형 그물 요소의 접선을 기저벡터로 사용한 FEM(유한요소법)으로 구성하였다. 우선 Arnoldi 분해법을 이용하여 이 방정식에 대한 상위 Hessenberg 행렬을 구하였다. 그리고 QR 알골리즘을 통하여 이것을 삼각형 대각 행렬인 Shur 형태로 변형하였다. 수렴 조건에 부합된 몇몇 고유 값들이 삼각형 대각 행렬의 대각 요소에 나타났다. 이들에 대응하는 고유 모드들을 역-반복법으로 구하였다. 수렴조건에 부합되는 고유 값들은 Shur 행렬의 대각선 선두 부분으로 재배열시켰다. 이들은 나머지 고유값 및 고유모드의 쌍을 구하는 반복 과정에서 변형되지 않도록 배제되었다. 이 과정이 연속하여 서너 번 반복되었는데, 그 결과 충분한 신뢰도를 갖는 주요한 몇 개의 TM 및 TE 고유 쌍들이 구하여졌다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제3권1호
• /
• pp.111-115
• /
• 2011

When formulating a general, non-linear mathematical model of ship dynamics in waves the hydrostatic forces and moments along with the Froude-Krylov part of wave load are usually concerned. Normally radiation and the diffraction forces are regarded as linear ones. The paper discusses briefly few approaches, which can be used in this respect. The concerned models attempt to model the non-linearities of the surface waves; both regular and the irregular ones, and the nonlinearities of the restoring forces and moments. The approach selected in the Laidyn method, which is meant for the evaluation of large amplitude motions in the 6 degrees-of-freedom, is presented in a bigger detail. The workability of the method is illustrated with the simulation of ship motions in irregular stern quartering waves.

• Structural Engineering and Mechanics
• /
• 제81권4호
• /
• pp.411-428
• /
• 2022

This paper presents the Campbell diagram analysis of the rotordynamic system using the full order model (FOM) and the reduced order model (ROM) techniques to determine the critical speeds, identify the stability and reduce the computational time. Due to the spin-speed-dependent matrices (e.g., centrifugal stiffening matrix), several model order reduction (MOR) techniques may be considered, such as the modal superposition (MS) method and the Krylov subspace-based MOR techniques (e.g., Ritz vector (RV), quasi-static Ritz vector (QSRV), multifrequency quasi-static Ritz vector (MQSRV), multifrequency/ multi-spin-speed quasi-static Ritz vector (MMQSRV) and the combined Ritz vector & modal superposition (RV+MS) methods). The proposed MMQSRV method in this study is extended from the MQSRV method by incorporating the rotational-speed-dependent stiffness matrices into the Krylov subspace during the MOR process. Thus, the objective of this note is to respond to the question of whether to use the MS method or the Krylov subspace-based MOR technique in establishing the Campbell diagram of the shaft-disc-blade assembly systems in three-dimensional (3D) finite element analysis (FEA). The Campbell diagrams produced by the FOM and various MOR methods are presented and discussed thoroughly by computing the norm of relative errors (ER). It is found that the RV and the MS methods are dominant at low and high rotating speeds, respectively. More precisely, as the spinning velocity becomes large, the calculated ER produced by the RV method is significantly increased; in contrast, the ER produced by the MS method is smaller and more consistent. From a computational point of view, the MORs have substantially reduced the time computing considerably compared to the FOM. Additionally, the verification of the 3D FE rotordynamic model is also provided and found to be in close agreement with the existing solutions.

• 한국전자통신학회논문지
• /
• 제15권6호
• /
• pp.1105-1112
• /
• 2020

1차원 3-이웃 셀룰라 오토마타(Cellular Automata, 이하 CA) 기반의 의사난수 생성기는 시스템의 성능을 평가하기 위한 테스트 패턴 생성과 암호 시스템의 키수열 생성기 등에 많이 응용되고 있다. 본 논문에서는 더 복잡하고 혼돈스러운 수열을 생성할 수 있는 CA기반의 키 수열 생성기를 설계하기 위해 각 셀의 상태전이에 영향을 주는 이웃을 5개로 확장한 1차원 대칭 5-이웃 CA에 대해 연구한다. 특히 대칭 5-이웃 CA를 합성하기 위해 Krylov 행렬을 이용하는 대수적인 방법과 Cho et al.의 알고리즘을 기반으로 한 1차원 n셀 대칭 5-이웃 CA 합성 알고리즘을 제안한다.

• Nuclear Engineering and Technology
• /
• 제54권11호
• /
• pp.4280-4295
• /
• 2022

In this paper, a new multi-group neutron-gamma transport calculation code system STRAUM-MATXST for complicated geometrical problems is introduced and its development status including numerical tests is presented. In this code system, the MATXST (MATXS-based Cross Section Processor for S_N Transport) code generates multi-group neutron and gamma cross sections by processing MATXS format libraries generated using NJOY and the STRAUM (S_N Transport for Radiation Analysis with Unstructured Meshes) code performs multi-group neutron-gamma coupled transport calculation using tetrahedral meshes. In particular, this work presents the recent implementation and its test results of the Krylov subspace methods (i.e., Bi-CGSTAB and GMRES(m)) with preconditioners using DSA (Diffusion Synthetic Acceleration) and TSA (Transport Synthetic Acceleration). In addition, the Krylov subspace methods for accelerating the energy-group coupling iteration through thermal up-scatterings are implemented with new multi-group block DSA and TSA preconditioners in STRAUM.

• 대한수학회논문집
• /
• 제10권4호
• /
• pp.981-995
• /
• 1995

This paper first establishes some conditions for preconditioner under which PGCR does not break down. Next, VPGCR algorithm whose preconditioners can be easily obtained is introduced and then its breakdown and convergence properties are discussed. Lastly, implementation details of VPGCR are described and then numerical results of VPGCR with a certain criterion guaranteeing no breakdown are compared with those of restarted GMRES.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 추계 학술대회논문집
• /
• pp.153-159
• /
• 1998

The Newton-Krylov method on the unstructured grid flow solver using the cell-centered spatial discretization oi compressible Euler equations is presented. This flow solver uses the reconstructed primitive variables to get the higher order solutions. To get the quadratic convergence of Newton method with this solver, the careful linearization of face flux is performed with the reconstructed flow variables. The GMRES method is used to solve large sparse matrix and to improve the performance ILU preconditioner is adopted and vectorized with level scheduling algorithm. To get the quadratic convergence with the higher order schemes and to reduce the memory storage. the matrix-free implementation and Barth's matrix-vector method are implemented and compared with the traditional matrix-vector method. The convergence and computing times are compared with each other.