• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.2
• /
• pp.153-163
• /
• 2010

In this paper a frequency response analysis using Krylov subspace-based model reduction and its design sensitivity analysis with respect to design variables are presented. Since the frequency response and its design sensitivity information are necessary for a gradient-based optimization, problems of high computational cost and resource may occur in the case that frequency response of a large sized finite element model is involved in the optimization iterations. In the suggested method model order reduction of finite element models are used to calculate both frequency response and frequency response sensitivity, therefore one can maximize the speed of numerical computation for the frequency response and its design sensitivity. As numerical examples, a semi-monocoque shell and an array-type $4{\times}4$ MEMS resonator are adopted to show the accuracy and efficiency of the suggested approach in calculating the FRF and its design sensitivity. The frequency response sensitivity through the model reduction shows a great time reduction in numerical computation and a good agreement with that from the initial full finite element model.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.9
• /
• pp.878-885
• /
• 2009

One of important factors in designing MEMS resonators for RF filters is obtaining a desired frequency response function (FRF) within a specific frequency range of interest. Because various array-type MEMS resonators have been recently introduced to improve the filter characteristics such as bandwidth, pass-band, and shape factor, the degrees of freedom (DOF) of finite elements for their FRF calculation dramatically increases and therefore raises computational difficulties. In this paper the Krylov subspace-based model order reduction using moment-matching with non-zero expansion points is represented as a numerical solution to perform the frequency response analyses of those array-type MEMS resonators in an efficient way. By matching moments at a frequency around the specific operation range of the array-type resonators, the required FRF can be efficiently calculated regardless of their operating frequency from significantly reduced systems. In addition, because of the characteristics of the moment-matching method, a minimal order of reduced system with a prearranged accuracy can be determined through an error indicator using successive reduced models, which is very useful to automate the order reduction process and FRF calculation for structural optimization iterations. We also found out that the presented method could obtain the FRF of a $6\times6$ array-type resonator within a seventieth of the computational time necessary for the direct method and in addition FRF calculation by the mode superposition method could not even be completed because of a data overflow with a half after calculation of 9,722 eigenmodes.

• Coupled systems mechanics
• /
• v.1 no.2
• /
• pp.183-204
• /
• 2012

A method for the analysis of wave scattering in 3-D elastic full space is developed by means of the coupled boundary-volume integral equation, which takes into account the effects of both the boundary of inclusions and the uctuation of the wave field. The wavenumber domain formulation is used to construct the Krylov subspace by means of FFT. In order to achieve the wavenumber domain formulation, the boundary-volume integral equation is transformed into the volume integral equation. The formulation is also focused on this transform and its numerical implementation. Several numerical results clarify the accuracy and effectiveness of the present method for scattering analysis.

• Journal of applied mathematics & informatics
• /
• v.32 no.3_4
• /
• pp.389-403
• /
• 2014

In this paper, we first provide comparison results of preconditioned AOR methods with direct preconditioners $I+{\beta}L$, $I+{\beta}U$ and $I+{\beta}(L+U)$ for solving a linear system whose coefficient matrix is a large sparse irreducible L-matrix, where ${\beta}$ > 0. Next we propose how to find a near optimal parameter ${\beta}$ for which Krylov subspace method with these direct preconditioners performs nearly best. Lastly numerical experiments are provided to compare the performance of preconditioned iterative methods and to illustrate the theoretical results.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.2
• /
• pp.45.2-45.2
• /
• 2019

The main-belt asteroid (5247) Krylov is known as a Non-Principal Axis (NPA) rotator. However, the shape model and spin state of this asteroid were not revealed. The physical model of an asteroid including spin state and shape is regarded to be important to understand its physical properties and dynamical evolution. Thus, in order to reconstruct the physical model of Kryolv, we applied the light curve inversion method using not only the optical light curves observed with ground-based telescopes in three apparitions during 2006, 2016, and 2017, but also the infrared light curves obtained with the Wide-field Infrared Survey Explorer (WISE) in 2010. We found that it is rotating in Short Axis Mode (SAM) with the rotation and precession periods of 368.71 hr and 67.277 hr, respectively. The orientation of the angular momentum vector is (298°, -58°) in the ecliptic coordinate system. The ratio of moments of inertia of the longest axis to the shortest axis is I_a/I_c = 0.36; the ratio of moments of inertia of the intermediate axis to the shortest axis is I_b/I_c = 0.96. Finally, the excitation level of this asteroid is found to be rather low with a ratio of the rotational kinetic energy to the basic spin state energy as E/E₀ ≃ 1.024. We will briefly discuss the possible evolutionary process of Krylov in this presentation.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.258-270
• /
• 2020

A new task of using the Jacobian-Free-Newton-Krylov (JFNK) method for the PWR core transient simulations involving neutronics, thermal hydraulics and mechanics is conducted. For the transient scenario of PWR, normally the Picard iteration of the coupled coarse-mesh nodal equations and parallel channel TH equations is performed to get the transient solution. In order to solve the coupled equations faster and more stable, the Newton Krylov (NK) method based on the explicit matrix was studied. However, the NK method is hard to be extended to the cases with more physics phenomenon coupled, thus the JFNK based iteration scheme is developed for the nodal method and parallel-channel TH method. The local gap conductance is sensitive to the gap width and will influence the temperature distribution in the fuel rod significantly. To further consider the local gap conductance during the transient scenario, a 1D mechanics model is coupled into the JFNK scheme to account for the fuel thermal expansion effect. To improve the efficiency, the physics-based precondition and scaling technique are developed for the JFNK iteration. Numerical tests show good convergence behavior of the iterations and demonstrate the influence of the fuel thermal expansion effect during the rod ejection problems.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.2
• /
• pp.127-136
• /
• 2016

Springing is the resonance phenomenon of a ship hull girder with incoming waves having the same natural frequency of the ship. In this study, a simple and reliable calculation method was developed based on quadratic strip theory using the Timoshenko beam approach as an elastic hull girder. Second-order hydrodynamic forces and Froude-Krylov forces were applied as the external force. To improve the accuracy of the strip method, the variation in the added mass along the ship hull longitudinal direction, so called tip-effect, was considered. The J-factor was also employed to compensate for the effect of three-dimensional fluid motion on the two-node vibration of the ship. Using the developed method, the first- and second-order vertical bending moments of the Flokstra ship were compared. A comparative study was also carried out for a uniform barge ship and a 10,000 TEU container ship with the respective methods including the J-factor and tip-effect.

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

In the proposed paper numerical calculations are carried out using two versions of a three-dimensional, timedomain panel method developed by the group of Prof. P. Sclavounos at MIT, i.e. the linear code SWAN2, enabling optionally the use of the instantaneous non-linear Froude-Krylov and hydrostatic forces and the fully non-linear SWAN4. The analytical results are compared with experimental results for three hull forms with increasing geometrical complexity, the Series 60, a reefer vessel with stern bulb and a modern fast ROPAX hull form with hollow bottom in the stern region. The details of the geometrical modeling of the hull forms are discussed. In addition, since SWAN4 does not support transom sterns, only the two versions of SWAN2 were evaluated over experimental results for the parent hull form of the NTUA double-chine, wide-transom, high-speed monohull series. The effect of speed on the numerical predictions was investigated. It is concluded that both versions of SWAN2 the linear and the one with the non-linear Froude-Krylov and hydrostatic forces provide a more robust tool for prediction of the dynamic response of the vessels than the non-linear SWAN4 code. In general, their results are close to what was expected on the basis of experience. Furthermore, the use of the option of non-linear Froude-Krylov and hydrostatic forces is beneficial for the accuracy of the predictions. The content of the paper is based on the Diploma thesis of the second author, supervised by the first one and further refined by the third one.

• Journal of applied mathematics & informatics
• /
• v.28 no.1_2
• /
• pp.351-361
• /
• 2010

The solution of large sparse linear systems is one of the most important problems in large scale scientific computing. Among the many methods developed, the preconditioned Krylov subspace methods are considered the preferred methods. Selecting a suitable preconditioner with appropriate parameters for a specific sparse linear system presents a challenging task for many application scientists and engineers who have little knowledge of preconditioned iterative methods. The prediction of ILU type preconditioners was considered in [27] where support vector machine(SVM), as a data mining technique, is used to classify large sparse linear systems and predict best preconditioners. In this paper, we apply the data mining approach to the sparse approximate inverse(SAI) type preconditioners to find some parameters with which the preconditioned Krylov subspace method on the linear systems shows best performance.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.5
• /
• pp.541-549
• /
• 2007

In the case of control for mechanical systems, it is highly useful to be able to provide a compact model of the mechanical system to control engineers using the smallest number of state variables, while still providing an accurate model. The reduced mechanical model can then be inserted into the complete system models and used for extended system-level dynamic simulation. In this paper, moment-matching based model order reductions (MOR) using Krylov subspaces, which reduce the number of degrees of freedom of an original finite element model via the Arnoldi process, are presented to study the eigenvalue and frequency response problems of a HDD actuator and suspension system.