• 제목/요약/키워드: Eigenvalue problems

검색결과 210건 처리시간 0.024초

A MIXED METHOD OF SUBSPACE ITERATION FOR DIRICHLET EIGENVALUE PROBLEMS

  • Lee, Gyou-Bong;Ha, Sung-Nam;Hong, Bum-Il
    • Journal of applied mathematics & informatics
    • /
    • 제4권1호
    • /
    • pp.243-248
    • /
    • 1997
  • A full multigrid scheme was used in computing some eigenvalues of the Laplace eigenvalues problem with the Dirichlet bound-ary condition. We get a system of algebraic equations with an aid of finite difference method and apply subspace iteration method to the system to compute first some eigenvalues. The result shows that this is very effective in calculating some eigenvalues of this problem.

THE CONVERGENCE OF A DUAL ALGORITHM FOR NONLINEAR PROGRAMMING

  • Zhang, Li-Wei;He, Su-Xiang
    • Journal of applied mathematics & informatics
    • /
    • 제7권3호
    • /
    • pp.719-738
    • /
    • 2000
  • A dual algorithm based on the smooth function proposed by Polyak (1988) is constructed for solving nonlinear programming problems with inequality constraints. It generates a sequence of points converging locally to a Kuhn-Tucker point by solving an unconstrained minimizer of a smooth potential function with a parameter. We study the relationship between eigenvalues of the Hessian of this smooth potential function and the parameter, which is useful for analyzing the effectiveness of the dual algorithm.

수직펌프의 진동해석 (Vibration Analysis of Vertical Pumps)

  • 홍진선
    • 소음진동
    • /
    • 제9권5호
    • /
    • pp.928-935
    • /
    • 1999
  • Vertical pumps are widely used owing to the fact that they occupy small floor space. In this type of pumps, however, the vibrational problems are very important, since, in many cases, they have less stiffness in comparison with lateral pumps. This study presents a simple solution method for calculating the natural frequencies and modes of vertical pumps. In this study, a mode of a vertical pump was developed and the nondimensional parameters for the vibrational characteristics of it were determined. Added mass was calculated for the effects of water and the transfer matrix method was used.

  • PDF

고유치 솔버 Davidson Method 의 병렬화 (A Parallel Algorithm of Davidson Method for Eigenproblems)

  • 김형중
    • 대한전기학회:학술대회논문집
    • /
    • 대한전기학회 1997년도 하계학술대회 논문집 A
    • /
    • pp.12-14
    • /
    • 1997
  • The analysis of eigenvalue and eigenvector is a crucial procedure for many electromagnetic computation problems. However, eigenpair computation is timing-consuming task. Thus, its parallelization is required for designing large-scale and precision three-dimensional electromagnetic machines. In this paper, the Davidson method is parallelized on a cluster of workstations. Performance of the parallelization scheme is reported. This scheme is applied to a ridged waveguide design problem.

  • PDF

대형 시스템에서의 다단계 부분구조 기법을 이용한 시스템 축소기법에 관한 연구 (Study on the Structural System Condensation Using Multi-level Sub-structuring Scheme in Large-scale Problems)

  • 백승민;조맹효;김현기
    • 한국전산구조공학회논문집
    • /
    • 제21권3호
    • /
    • pp.281-285
    • /
    • 2008
  • 축소시스템 기법은 전체 구조의 거동을 나타내는 저차 고유모드를 근사화한다. 지난 연구에서 축소시스템을 구축하기 위한 2단계 축소기법을 제안하였다. 또, 기존의 2단계 축소기법을 반복적 IRS기법을 통해 중간 주파수 대역의 고유모드에 대한 해의 정확도를 높이는 방안에 대해 연구가 제안되었다. 본 연구에서는 기존의 향상된 2단계 축소기법에 다단계 부구조화 기법을 적용하는 기법을 제안한다. 첫 단계에서는 전체 시스템을 그래프 분할을 통해 계층적으로 부구조로 분할되고, 두 번째 단계에서는 각각의 부구조를 개선된 2단계 축소기법을 이용하여 축소한다. 각각의 축소된 분절화된 고유치문제의 조합을 총해 최종적 축소시스템을 구축하고 이렇게 구한 축소된 고유치 문제를 란초스 기법(ARPACK)을 통해 해석한다. 최종적으로 제안된 기법의 성능을 수치 예제를 통해 검증한다.

Stabilizing Solutions of Algebraic Matrix riccati Equations in TEX>$H_\infty$ Control Problems

  • Kano, Hiroyuki;Nishimura, Toshimitsu
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 제어로봇시스템학회 1994년도 Proceedings of the Korea Automatic Control Conference, 9th (KACC) ; Taejeon, Korea; 17-20 Oct. 1994
    • /
    • pp.364-368
    • /
    • 1994
  • Algebraic matrix Riccati equations of the form, FP+PF$^{T}$ -PRP+Q=0. are analyzed with reference to the stability of closed-loop system F-PR. Here F, R and Q are n * n real matrices with R=R$^{T}$ and Q=Q$^{T}$ .geq.0 (nonnegative-definite). Such equations have been playing key roles in optimal control and filtering problems with R .geq. 0. and also in the solutions of in H$_{\infty}$ control problems with R taking the form R=H$_{1}$$^{T}$ H$_{1}$-H$_{2}$$^{T}$ H$_{2}$. In both cases an existence of stabilizing solution, i.e. the solution yielding asymptotically stable closed-loop system, is an important problem. First, we briefly review the typical results when R is of definite form, namely either R .geq. 0 as in LQG problems or R .leq. 0. They constitute two extrence cases of Riccati to the cases H$_{2}$=0 and H$_{1}$=0. Necessary and sufficient conditions are shown for the existence of nonnegative-definite or positive-definite stabilizing solution. Secondly, we focus our attention on more general case where R is only assumed to be symmetric, which obviously includes the case for H$_{\infty}$ control problems. Here, necessary conditions are established for the existence of nonnegative-definite or positive-definite stabilizing solutions. The results are established by employing consistently the so-called algebraic method based on an eigenvalue problem of a Hamiltonian matrix.x.ix.x.

  • PDF

Airy 함수와 WKB 해를 이용한 언덕형 굴절율 분포를 갖는 광섬유 해석 (Analysis of Optical Fibers with Graded-Index Profile By a Combination of Modified Airy Functions and WKB Solutions)

  • 정민섭;김영문;김창민
    • 대한전자공학회논문지SD
    • /
    • 제37권2호
    • /
    • pp.28-37
    • /
    • 2000
  • 수정된 Airy 함수들과 WKB 시도해를 조합하여 언덕형 굴절율 분포를 가지는 광섬유에 대한 매우 정확한 고유방정식을 수학적으로 유도하였다. 적절한 경계조건을 적용하여 일반적인 WKB 방법이 가지는 고유한 오차 문제를 개선하도록 위상 천이 보정항 δ를 이끌어냈다. 모의 전산을 통하여, 유도된 고유방정식의 결과가 유한차분법에 의한 결과와 잘 일치함을 보였다.

  • PDF

비부합 결합을 이용한 구조물 변경법 (Structural dynamics modification using non-matching substructure synthesis.)

  • 정의일;박윤식
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 2002년도 춘계학술대회논문집
    • /
    • pp.666-671
    • /
    • 2002
  • For a large structure, substructure based SDM(structural dynamics modification) method is very effective to raise its dynamic characteristics. Dividing into smaller substructures has a major advantage in the aspect of computation especially for getting sensitivities, which are in the core of SDM process. But quite often, non-matching nodes problem occurs in the process of synthesizing substructures. The reason is that, in general, each substructure is modelled separately, then later combined together to form a entire structure model under interface constraint conditions. Without solving the non-matching nodes problem, the substructure based SDM can not be processed. In this work, virtual node concept is introduced. Lagrange multipliers are used to enforce the interface compatibility constraint. The governing equation of whole structure is derived using hybrid variational principle. The eigenvalues of whole structure are calculated using determinant search method. The number of degrees of freedom of the eigenvalue problem can be drastically reduced to just the number of interface degree of freedom. Thus, the eigenvalue sensitivities can be easily calculated, and further SDM can be efficiently performed. Some numerical problems are tested to show the effectiveness of handling non-matching nodes.

  • PDF

Convergence study of traditional 2D/1D coupling method for k-eigenvalue neutron transport problems with Fourier analysis

  • Boran Kong ;Kaijie Zhu ;Han Zhang ;Chen Hao ;Jiong Guo ;Fu Li
    • Nuclear Engineering and Technology
    • /
    • 제55권4호
    • /
    • pp.1350-1364
    • /
    • 2023
  • 2D/1D coupling method is an important neutron transport calculation method due to its high accuracy and relatively low computation cost. However, 2D/1D coupling method may diverge especially in small axial mesh size. To analyze the convergence behavior of 2D/1D coupling method, a Fourier analysis for k-eigenvalue neutron transport problems is implemented. The analysis results present the divergence problem of 2D/1D coupling method in small axial mesh size. Several common attempts are made to solve the divergence problem, which are to increase the number of inner iterations of the 2D or 1D calculation, and two times 1D calculations per outer iteration. However, these attempts only could improve the convergence rate but cannot deal with the divergence problem of 2D/1D coupling method thoroughly. Moreover, the choice of axial solvers, such as DGFEM SN and traditional SN, and its effect on the convergence behavior are also discussed. The results show that the choice of axial solver is a key point for the convergence of 2D/1D method. The DGFEM SN based 2D/1D method could converge within a wide range of optical thickness region, which is superior to that of traditional SN method.