• Title/Summary/Keyword: cubic convergence

Search Result 78, Processing Time 0.022 seconds

Dynamic response optmization using approximate search (근사 선탐색을 이용한 동적 반응 최적화)

  • Kim, Min-Soo;Choi, Dong-hoon
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.22 no.4
    • /
    • pp.811-825
    • /
    • 1998
  • An approximate line search is presented for dynamic response optimization with Augmented Lagrange Multiplier(ALM) method. This study empolys the approximate a augmented Lagrangian, which can improve the efficiency of the ALM method, while maintaining the global convergence of the ALM method. Although the approximate augmented Lagragian is composed of only the linearized cost and constraint functions, the quality of this approximation should be good since an approximate penalty term is found to have almost second-order accuracy near the optimum. Typical unconstrained optimization algorithms such as quasi-Newton and conjugate gradient methods are directly used to find exact search directions and a golden section method followed by a cubic polynomial approximation is empolyed for approximate line search since the approximate augmented Lagrangian is a nonlinear function of design variable vector. The numberical performance of the proposed approach is investigated by solving three typical dynamic response optimization problems and comparing the results with those in the literature. This comparison shows that the suggested approach is robust and efficient.

Periodic solutions of the Duffing equation

  • Tezcan, Jale;Hsiao, J. Kent
    • Structural Engineering and Mechanics
    • /
    • v.30 no.5
    • /
    • pp.593-602
    • /
    • 2008
  • This paper presents a new linearization algorithm to find the periodic solutions of the Duffing equation, under harmonic loads. Since the Duffing equation models a single degree of freedom system with a cubic nonlinear term in the restoring force, finding its periodic solutions using classical harmonic balance (HB) approach requires numerical integration. The algorithm developed in this paper replaces the integrals appearing in the classical HB method with triangular matrices that are evaluated algebraically. The computational cost of using increased number of frequency components in the matrixbased linearization approach is much smaller than its integration-based counterpart. The algorithm is computationally efficient; it only takes a few iterations within the region of convergence. An example comparing the results of the linearization algorithm with the "exact" solutions from a 4th order Runge- Kutta method are presented. The accuracy and speed of the algorithm is compared to the classical HB method, and the limitations of the algorithm are discussed.

HYBRID DIFFERENCE SCHEMES FOR A SYSTEM OF SINGULARLY PERTURBED CONVECTION-DIFFUSION EQUATIONS

  • Priyadharshini, R.Mythili;Ramanujam, N.;Tamilselvan, A.
    • Journal of applied mathematics & informatics
    • /
    • v.27 no.5_6
    • /
    • pp.1001-1015
    • /
    • 2009
  • In this paper, two hybrid difference schemes on the Shishkin mesh are constructed for solving a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a small parameter multiplying the highest derivative. We prove that the schemes are almost second order convergence in the supremum norm independent of the diffusion parameter. Error bounds for the numerical solution and its derivative are established. Numerical results are provided to illustrate the theoretical results.

  • PDF

Integral Approximate Solutions to a One-Dimensional Model for Stratified Thermal Storage Tanks (성층화된 축열조의 1차원모델에 대한 적분 근사해)

  • Chung, Jae-Dong
    • Korean Journal of Air-Conditioning and Refrigeration Engineering
    • /
    • v.22 no.7
    • /
    • pp.468-473
    • /
    • 2010
  • This paper deals with approximate integral solutions to the one-dimensional model describing the charging process of stratified thermal storage tanks. Temperature is assumed to be the form of Fermi-Dirac distribution function, which can be separated to two sets of cubic polynomials for each hot and cold side of thermal boundary layers. Proposed approximate integral solutions are compared to the previous works of the approximate analytic solutions and show reasonable agreement. The approach, however, has benefits in mathematical difficulties, complicated solution form and unstable convergence of series solution founded in the previous analytic solutions. Solutions for a semi-infinite region, which have simple closed form solutions, give close agreement to those for a finite region. Thermocline thickness is obtained in closed form and shows proportional behavior to the square root of time and inverse proportional behavior to the square root of flow rate.

Free vibration analysis of rotating cantilever plates using the p-version of the finite element method

  • Hamza-Cherif, Sidi Mohammed
    • Structural Engineering and Mechanics
    • /
    • v.22 no.2
    • /
    • pp.151-167
    • /
    • 2006
  • A p-version of the finite element method in conjunction with the modeling dynamic method using the arc-length stretch deformation is considered to determine the bending natural frequencies of a cantilever flexible plate mounted on the periphery of a rotating hub. The plate Fourier p-element is used to set up the linear equations of motion. The transverse displacements are formulated in terms of cubic polynomials functions used generally in FEM plus a variable number of trigonometric shapes functions representing the internals DOF for the plate element. Trigonometric enriched stiffness, mass and centrifugal stiffness matrices are derived using symbolic computation. The convergence properties of the rotating plate Fourier p-element proposed and the results are in good agreement with the work of other investigators. From the results of the computation, the influences of rotating speed, aspect ratio, Poisson's ratio and the hub radius on the natural frequencies are investigated.

Modal Analysis of Large Scale Multi-Machine Power System using Rayleigh Quotient and Deflation (Rayleigh Quotient와 Deflation을 이용한 대형다기(多機)전력계통의 고유치 해석)

  • Shim, Kwan-Shik;Nam, Hae-Kon
    • Proceedings of the KIEE Conference
    • /
    • 1993.07a
    • /
    • pp.76-78
    • /
    • 1993
  • This paper describes an efficient method of computing any desired number of the most unstable eigenvalues and eigenvectors of a large scale multi-machine power system. Approximate eigenvalues obtained by Hessenberg process are refined using Rayleigh quotient iteration with cubic convergence property. If further eigenvalues and eigenvectors are needed, the procedure described above are repeated with deflation. The proposed algorithm can cover all the model types of synchronous machines, exciters, speed governing system and PSS defined in AESOPS. The proposed algorithm applied to New England test system with 10 machines and 39 buses produced the results same with AESOPS in faster computation time. Also eigenvectors computed in Rayleigh quotient iteration makes it possible to make eigen-analysis for improving unstable modes.

  • PDF

Surface Relaxation of Aluminum

  • Cha, You-Yong;Han, Wone-Keun
    • Proceedings of the Korean Vacuum Society Conference
    • /
    • 2000.02a
    • /
    • pp.142-142
    • /
    • 2000
  • We performed a total energy calculation of clean alumunum surfaces of three low indices based on a density functional theory with a local density approximation, using the Ceperly-Alder exhange correlation parametrized by Perdew and Zunger. Pseudopotentials were generated for Al of which the plane wave cut-off was 15Ry. We used Gaussian broadening of a Fermi level to accelerate the convergence of our calculation with the Gaussian energy smearing parameter of 0.005Ry. First, we determine the lattice constant of the aluminum of an face-centered-cubic structure to be 3.96 which is comparable to the experimental data of 4.05 . The cohesive energy of 4.20eV/atom and the bulk modulus of 0.775$\times$1012dyne/cm2 are also comparable to the experimental values of 3.39eV/atom and 0.772$\times$1012dyne/cm2, respectively. Then we investigated the surface relaxation of (100), (110) and (111) surfaces using a 9-layer slab separated by 6-layer thick vacuum. The results are consistent with the existing experimental results.

  • PDF

Sampling-Based Sensitivity Approach to Electromagnetic Designs Utilizing Surrogate Models Combined with a Local Window

  • Choi, Nak-Sun;Kim, Dong-Wook;Choi, K.K.;Kim, Dong-Hun
    • Journal of Magnetics
    • /
    • v.18 no.1
    • /
    • pp.74-79
    • /
    • 2013
  • This paper proposes a sampling-based optimization method for electromagnetic design problems, where design sensitivities are obtained from the elaborate surrogate models based on the universal Kriging method and a local window concept. After inserting additional sequential samples to satisfy the certain convergence criterion, the elaborate surrogate model for each true performance function is generated within a relatively small area, called a hyper-cubic local window, with the center of a nominal design. From Jacobian matrices of the local models, the accurate design sensitivity values at the design point of interest are extracted, and so they make it possible to use deterministic search algorithms for fast search of an optimum in design space. The proposed method is applied to a mathematical problem and a loudspeaker design with constraint functions and is compared with the sensitivity-based optimization adopting the finite difference method.

Design of GeoPhoto Contents Markup Language for u-GIS Contents (u-GIS 콘텐츠를 위한 GeoPhoto 콘텐츠 언어의 설계)

  • Park, Jang-Yoo;Nam, Kwang-Woo;Jin, Heui-Chae
    • Journal of Korea Spatial Information System Society
    • /
    • v.11 no.1
    • /
    • pp.35-42
    • /
    • 2009
  • This paper proposes a new GeoPhoto contents markup language that can create u-GIS contents by using the spatial photos. GeoPhoto contents markup language has designed the GeoPhoto contents model and markup language for contents that can be used the spatial photos information. GeoPhoto con tents markup language is represented by the convergence of GIS information, location information, photos information respectively. GeoPhoto contents markup language to provide a variety of pictures related to the content model consists of GeoPhoto contents model and operations between the GeoPhoto contents. GeoPh oto contents model supports GeoPhoto model, CubicPhoto model, Photo model and SequenceGeoPhoto mod el. In addition, this paper propose the Annotation operation, Enlargement operation and Overlay operation for represent the GeoPhoto contents. GeoPhoto Contents Markup Language has the advantage of supportin g user custom contents model of u-GIS.

  • PDF

Structural and Electrical Properties of K(Ta,Nb)O3 Ceramics with Variation of Ag Contents for Electrocaloric Devices (전기열량소자용 Ag 첨가량에 따른 K(Ta,Nb)O3 세라믹스의 구조적·전기적 특성)

  • Lee, Min-Sung;Park, Byeong-Jun;Lim, Jeong-Eun;Lee, Sam-Haeng;Lee, Myung-Gyu;Park, Joo-Seok;Lee, Sung-Gap
    • Journal of the Korean Institute of Electrical and Electronic Material Engineers
    • /
    • v.34 no.6
    • /
    • pp.442-448
    • /
    • 2021
  • In this work, the (K1-xAgx)(Ta0.8Nb0.2)O3 (x=0.1-0.4) ceramics were fabricated using mixed-oxide method, and their structural and electrical properties were measured. All specimens represented a pseudo cubic structure with the lattice constant of 0.3989 nm. When 0.4 mol of Ag was added, second phases induced from metallic Ag and K2(Ta,Nb)6O16 phase were observed. Dielectric constant and dielectric loss of K(Ta0.8Nb0.2)O3 specimen doped with 0.3 mol of Ag were 2,737 and 0.446, respectively. The curie temperature was about -5℃, which does not change with Ag addition. The remanent polarization began to decrease sharply around 12~15℃, and the temperature at which the remanent polarization began to decrease as the applied voltage increased shifted to the high temperature side. The electrocaloric effect (ΔT) and electrocaloric efficiency (ΔT/ΔE) of the (K0.7Ag0.3)(Ta0.8Nb0.2)O3 ceramics were 0.01024℃ and 0.01825 KmV-1, respectively.