• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.5
• /
• pp.445-455
• /
• 2010

In order to effectively predict thermo-acoustic instability within real combustors of rocket engines and gas turbines, in the present study, the Helmholtz equation in conjunction with the time lag hypothesis is discretized by the finite element method on three-dimensional hybrid unstructured mesh. Numerical nonlinearity caused by the combustion response term is linearized by an iterative method, and the large-scale eigenvalue problem is solved by the Arnoldi method available in the ARPACK. As a consequence, the final solution of complex valued eigenfrequency and acoustic pressure field can be interpreted as resonant frequency, growth rate, and modal shape for acoustic modes of interest. The predictive capabilities of the present method have been validated against two academic problems with complex impedance boundary and premixed flame, as well as an ambient acoustic test for liquid rocket combustion chamber with/without baffle.

• Structural Engineering and Mechanics
• /
• v.73 no.3
• /
• pp.271-286
• /
• 2020

Passive control technologies are commonly used in several areas to suppress structural vibrations by the addition of supplementary damping, and some modal damping may be heavy beyond critical damping even for regular structures with energy dissipation devices. The design of passive control structures is typically based on (complex) mode superposition approaches. However, the conventional mode superposition approach is predominantly applied to cases of under-critical damping. Moreover, when any modal damping ratio is equal or close to 1.0, the system becomes defective, i.e., a complete set of eigenvectors cannot be obtained such that some well-known algorithms for the quadratic eigenvalue problem are invalid. In this paper, a generalized complex mode superposition method that is suitable for under-critical, critical and over-critical damping is proposed and expressed in a unified form for structural displacement, velocity and acceleration responses. In the new method, the conventional algorithm for the eigenvalue problem is still valid, even though the system becomes defective due to critical modal damping. Based on the modal truncation error analysis, modal corrected methods for displacement and acceleration responses are developed to approximately consider the contribution of the truncated higher modes. Finally, the implementation of the proposed methods is presented through two numerical examples, and the effectiveness is investigated. The results also show that over-critically damped modes have a significant impact on structural responses. This study is a development of the original complex mode superposition method and can be applied well to dynamic analyses of non-classically damped systems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.371-383
• /
• 1999

An efficient and numerically stable eigensolution method for structures with multiple natural frequencies is presented. The proposed method is developed by improving the well-known subspace iteration method with shift. A major difficulty of the subspace iteration method with shift is that because of singularity problem, a shift close to an eigenvalue can not be used, resulting in slower convergence. In this paper, the above singularity problem has been solved by introducing side conditions without sacrifice of convergence. The proposed method is always nonsingular even if a shift is on a distinct eigenvalue or multiple ones. This is one of the significant characteristics of the proposed method. The nonsingularity is proved analytically. The convergence of the proposed method is at least equal to that of the subspace iteration method with shift, and the operation counts of above two methods are almost the same when a large number of eigenpairs are required. To show the effectiveness of the proposed method, two numerical examples are considered.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.27 no.6
• /
• pp.782-786
• /
• 1994

Integrated finite difference method (IFDM) is used to solve one dimensional free oscillation problem in the coastal area. To evaluate the solution accuracy of IFDM in free oscillation analysis, two finite difference equations based on area discretization method and point discretization method are derived from the governing equations of free oscillation, respectively. The difference equations are transformed into a generalized eigenvalue problem, respectively. A numerical example is presented, for which the analytical solution is available, for comparing IFDM to conventional finite difference equation (CFDM), qualitatively. The eigenvalue matrices are solved by sub-space iteration method. The numerical results of the two methods are in good agreement with analytical ones, however, IFDM yields better solution than CFDM in lower modes because IFDM only includes first order differential operator in finite difference equation by Green's theorem. From these results, it is concluded that IFDM is useful for the free oscillation analysis in the coastal area.

• Journal of Korean Society of Steel Construction
• /
• v.13 no.6
• /
• pp.703-713
• /
• 2001

Derivation procedures of exact dynamic stiffness matrices of thin-walled straight beams subjected to eccentrically axial forces are rigorously presented for the spatial free vibration analysis. An exact dynamic stiffness matrix is established from governing equations for a uniform beam element with nonsymmetric thin-walled cross section. First this numerical technique is accomplished via a generalized linear eigenvalue problem by introducing 14 displacement parameters and a system of linear algebraic equations with complex matrices. Thus, the displacement functions of displacement parameters are exactly derived and finally exact stiffness matrices are determined using element force-displacement relationships. The natural frequencies of nonsymmetric thin-walled straight beams are evaluated and compared with analytical solutions or results by thin-walled beam element using the cubic Hermitian polynomials and ABAQU's shell elements in order to demonstrate the validity of this study.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1995.10a
• /
• pp.225-233
• /
• 1995

This paper presents an efficient numerical method for computation of eigenpair derivatives for the real symmetric eigenvalue problem with distinct and multiple eigenvalues. The method has very simple algorithm and gives an exact solution. Furthermore, it saves computer storage and CPU time. The algorithm preserves the symmetry and band of the matrices, allowing efficient computer storage and solution techniques. Thus, the algorithm of the proposed method will be inserted easily in the commercial FEM codes. Results of the proposed method for calculating the eigenpair derivatives are compared with those of Rudisill and Chu's method and Nelson's method which is efficient one in the case of distinct natural frequencies. As an example to demonstrate the efficiency of the proposed method in the case of distinct eigenvalues, a cantilever plate is considered. The design parameter of the cantilever plate is its thickness. For the eigenvalue problem with multiple natural frequencies, the adjacent eigenvectors are used in the algebraic equation as side conditions, they lie adjacent to the m (multiplicity of multiple natural frequency) distinct eigenvalues, which appear when design parameter varies. As an example to demonstrate the efficiency of the proposed method in the case of multiple natural frequencies, a cantilever beam is considered. Results of the proposed method fDr calculating the eigenpair derivatives are compared with those of Bailey's method (an amendation of Ojalvo's work) which finds the exact eigenvector derivatives. The design parameter of the cantilever beam is its height. Data is persented showing the amount of CPU time used to compute the first ten eigenpair derivatives by each method. It is important to note that the numerical stability of the proposed method is proved.

• The KIPS Transactions:PartA
• /
• v.10A no.2
• /
• pp.137-140
• /
• 2003

Recently, CG (Conjugate Gradient) scheme for the optimization of the Rayleigh quotient has been proven a very attractive and promising technique for interior eigenvalues for the following eigenvalue problem, Ax＝λx (1) The given matrix A is assummed to be large and sparse, and symmetric. Also, the method is very amenable to parallel computations. A proper choice of the preconditioner significantly improves the convergence of the CG scheme. We compare the parallel preconditioners for the computation of the interior eigenvalues of a symmetric matrix by CG-type method. The considered preconditioners are Point-SSOR, ILU (0) in the multi-coloring order, and Multi-Color Block SSOR (Symmetric Succesive OverRelaxation). We conducted our experiments on the CRAY­T3E with 128 nodes. The MPI (Message Passing Interface) library was adopted for the interprocessor communications. The test matrices are up to $512{\times}512$ in dimensions and were created from the discretizations of the elliptic PDE. All things considered the MC-BSSOR seems to be most robust preconditioner.

• Structural Engineering and Mechanics
• /
• v.55 no.6
• /
• pp.1261-1278
• /
• 2015

Due to their low intrinsic damping, stay cables in cable-stayed bridges have often exhibited unanticipated and excessive vibrations which result in increasing maintenance frequency and disruption to normal operations of the entire bridges. Mitigation of undesired cable vibration can be achieved by attaching an external damping device near the anchorage. High Damping Rubber (HDR) dampers have many advantages such as compact size, better aesthetics, easy maintenance, temperature stability, and cost benefits; therefore, they have been widely used to increase cable damping. Although a single damper has been shown to reduce cable vibrations, it is not the most effective method due to geometric constraints. This paper proposes the use of two HDR dampers to improve effectiveness and robustness in suppressing cable vibration. Oscillation parameters of the cable-dampers system were investigated in detail by modeling the stay cable as a taut string and each HDR damper as complex-valued impedance and by using an analytical formulation of the complex eigenvalue problem. The problem of two HDR dampers arbitrarily located along a cable is solved and the solution is discussed. Asymptotic formulas to calculate the damping ratios of the cable with two HDR dampers installed near the anchorage(s) are proposed and compared with the exact solutions. Further, a design example is presented in order to justify the methodology. The results of this study show that when the two HDR dampers are installed close to each other on the same end of the cable, some interaction between the dampers leads to reduced damping ratio. When the dampers are on the opposite ends of the cable, they are effective in increasing damping ratio and can provide better vibration reduction to multiple modes.

• The Journal of the Institute of Internet, Broadcasting and Communication
• /
• v.19 no.3
• /
• pp.173-178
• /
• 2019

To analyze the diffraction properties of optical signals by periodic blazed 2D diffraction gratings, Toeplitz dielectric tensor is first defined and formulated by 2D spatial Fourier expansions associated with asymmetric blazed grating profile. The characteristic modes in each layer is then based on eigenvalue problem, and the complete solution is found rigorously in terms of modal transmission-line theory (MTLT) to address the pertinent boundary-value problems. Toeplitz matrix of symmetric and sawtooth profiles is derived from that of asymmetric blazed grating profile, and the diffraction properties for each profile are numerically simulated. The numerical results reveal that the asymmetric and symmetric profiles behave as anti-reflection GMR filter while the sawtooth profile works better as anti-transmission one rather than anti-reflection filter.

• Journal of Industrial Technology
• /
• v.17
• /
• pp.255-260
• /
• 1997

The analysis of eigenvalue and eigenvector is a crucial procedure for many electromagnetic computation problems. Although it is always the case in practice that only selected eigenpairs are needed, computation of eigenpair still seems to be a time-consuming task. In order to compute the eigenpair more quickly, there are two resorts: one is to select a good algorithm with care and another is to use parallelization technique to improve the speed of the computing. In this paper, one of the best eigensolver, the Davidson method, is parallelized on a cluster of workstations. We apply this scheme to a ridged waveguide design problem and obtain promising linear speedup and scalability.