• Title/Summary/Keyword: truncation error

Search Result 94, Processing Time 0.021 seconds

A substructure formulation for the earthquake -induced nonlinear structural pounding problem

  • Shi, Jianye;Bamer, Franz;Markert, Bernd
    • Earthquakes and Structures
    • /
    • v.17 no.1
    • /
    • pp.101-113
    • /
    • 2019
  • Earthquake-induced pounding is one of the major reasons for structural failure in earthquake prone cities. An accurate description of the pounding phenomenon of two buildings requires the consideration of systems with a large number of degrees of freedom including adequate contact impact formulations. In this paper, firstly, a node to surface formulation for the realization of state-of-the-art pounding models for structural beam elements is presented. Secondly, a hierarchical substructure technique is introduced, which is adapted to the structural pounding problem. The numerical accuracy and efficiency of the method, especially for the contact forces, are verified on an academic example, applying four different impact elements. Error estimations are carried out and compared with the classical modal truncation method. It is demonstrated that the hierarchical substructure method is indeed able to significantly speed up the numeric integration procedure by preserving a required level of accuracy.

Numerical Solution of Nonlinear Diffusion in One Dimensional Porous Medium Using Hybrid SOR Method

  • Jackel Vui Lung, Chew;Elayaraja, Aruchunan;Andang, Sunarto;Jumat, Sulaiman
    • Kyungpook Mathematical Journal
    • /
    • v.62 no.4
    • /
    • pp.699-713
    • /
    • 2022
  • This paper proposes a hybrid successive over-relaxation iterative method for the numerical solution of a nonlinear diffusion in a one-dimensional porous medium. The considered mathematical model is discretized using a computational complexity reduction scheme called half-sweep finite differences. The local truncation error and the analysis of the stability of the scheme are discussed. The proposed iterative method, which uses explicit group technique and modified successive over-relaxation, is formulated systematically. This method improves the efficiency of obtaining the solution in terms of total iterations and program elapsed time. The accuracy of the proposed method, which is measured using the magnitude of absolute errors, is promising. Numerical convergence tests of the proposed method are also provided. Some numerical experiments are delivered using initial-boundary value problems to show the superiority of the proposed method against some existing numerical methods.

Using the Monte Carlo method to solve the half-space and slab albedo problems with Inönü and Anlı-Güngör strongly anisotropic scattering functions

  • Bahram R. Maleki
    • Nuclear Engineering and Technology
    • /
    • v.55 no.1
    • /
    • pp.324-329
    • /
    • 2023
  • Different types of deterministic solution methods were used to solve neutron transport equations corresponding to half-space and slab albedo problems. In these types of solution methods, in addition to the error of the numerical solutions, the obtained results contain truncation and discretization errors. In the present work, a non-analog Monte Carlo method is provided to simulate the half-space and slab albedo problems with Inönü and Anlı-Güngör strongly anisotropic scattering functions. For each scattering function, the sampling method of the direction of the scattered neutrons is presented. The effects of different beams with different angular dependencies and the effects of different scattering parameters on the reflection probability are investigated using the developed Monte Carlo method. The validity of the Monte Carlo method is also confirmed through the comparison with the published data.

A Study on the Periodic Transient Response Characteristics in Annular Fin with Uniform Thickness (均一두께의 環狀흰에서 週期的 過渡應答 特性에 관한 硏究)

  • 김광수
    • Transactions of the Korean Society of Mechanical Engineers
    • /
    • v.12 no.2
    • /
    • pp.338-348
    • /
    • 1988
  • This study presents an analysis of periodic heat diffusion in an annular fin with uniform thickness. When the temperature of the fin base is changed in the form of a sinusoidal function, the exact temperature solution can be obtained by Laplace transformation in terms of the dimensionless parameters in the infinite series. Local heat flux and average heat flux, local fin efficiency and average fin efficiency were obtained. Particularly, the table of eigenvalues that are the indispensable condition in solving the heat transfer problem of annular fin in a transient state with convection phenomena at the fin edge is provided. The tables of heat fluxes and average heat fluxes, fin efficiencies and average fin efficiencies are also provided from the computed results. Also, substituting the variations of dimensionless parameters into the these exact solutions, the characteristics of these response are investigated.

A Case Study on the Importance of Residual Compensation in FRF-based Substructuring (주파수 응답함수를 이용한 부분구조 합성에서 모드자름 오차 보정에 관한 수치적 연구)

  • 박윤식;김경호
    • Transactions of the Korean Society for Noise and Vibration Engineering
    • /
    • v.12 no.4
    • /
    • pp.302-309
    • /
    • 2002
  • A FRF-based substructuring method attempts to predict the dynamic characteristics of a complex structure from predetermined FRFs of the comprising uncoupled substructures. Although this method has the advantage of being able to incorporate experimental component FRFs directly, it is prone to errors : measurement errors, coordinate incompleteness, modal incompleteness, etc. Among the various sources of errors, this paper deals with the problem of modal incompleteness (or residual problem) of which importance is underestimated compared to others. It is a well-known rule of thumb that such a problem can be overcome by including modes up to 2 or 3 times the upper frequency of interest. Using a simulated case study, it is demonstrated that even including modes up to 20 times the upper frequency of interest does not guarantee a satisfactory result. A method to compensate the residual errors is introduced. This method requires the whole FRF matrices of substructures which is practically impossible for a complex structure. An applicable alternative is suggested and applied successfully to the case study. Finally, the effects of measurement errors on the residual compensation are also discussed.

A Study on Transient Analysis of Linear Induction Motor with Ununiform Airgap for Shallow-depth Underground Train (저심도철도용 선형유도전동기의 공극 불균일 과도특성 분석 연구)

  • Lee, Hyung-Woo;Park, Chan-Bae;Won, Sunghong
    • The Transactions of The Korean Institute of Electrical Engineers
    • /
    • v.62 no.5
    • /
    • pp.723-729
    • /
    • 2013
  • This paper presents an parallel type Linear Induction Motor with ununiform airgap for a shallow-depth underground train with 100‰ grade and 15 m curvature. This parallel type LIM has enough gradability but has inherently ununiform airgap between center and end parts. Consequently, performance when the train passes curved section should be considered with transient analysis. Moreover, general parallel operation, 1C2M which is usually used for train operation, deteriorates LIM performance because of different line velocity between inner and outer LIMs. Transient analysis has many problems such as huge model, lots of meshes, very long calculation time, truncation error and so on. This paper has presented a novel technique using equivalent linear rotating model in order to solve these problems and has analyzed parallel type LIM by using the proposed technique. Finally, LIM performance according to independent operating control has been investigated.

(Frequency Weighted Reduction Using Iterative Approach of BMI) (BMI의 반복적 해법을 이용한 주파수하중 차수축소)

  • Kim, Yong-Tae;O, Do-Chang;Park, Hong-Bae
    • Journal of the Institute of Electronics Engineers of Korea SC
    • /
    • v.39 no.1
    • /
    • pp.33-41
    • /
    • 2002
  • In this paper, we present a frequency weighted model reduction using LMIs for minimizing the H$\infty$ weighted model error compared with the methods of frequency weighted balanced truncation and frequency weighted Hankel norm approximation. The proposed algorithm, its form is equal to the sufficient condition of performance preserving controller approximation, is based on an iterative two-step LMI scheme induced from bound real lemma. So, it can be applied to the problem of the performance preserving controller approximation. The controller reduction is useful in a practical control design and ensures its easy implementation and high reliability The validity of the proposed algorithm is shown through numerical examples. Additionaly, we extend the proposed algorithm to performance preserving controller approximation by applying to the HIMAT(highly maneuverable aircraft technology) system.

Optimized Compact Finite Difference Schemes for Computational Aeroacoustics (수치적 공력음향해석을 위한 집적유한차분법의 개발)

  • Lee D. J.;Kim J. W.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 1995.10a
    • /
    • pp.197-202
    • /
    • 1995
  • 비정상(unsteady) 압축성(compressible) 유동에 의한 공력음향(aeroacoustics)을 모사하여 공력소음원을 해석하기 위해서는 고차(high order)의 정확도와 높은 해상도(resolution)를 가지며, 상대적으로 계산시간을 많이 필요로 하지 않는 외재적(explicit) 유한차분법이 필수적으로 요구된다. 이것은 주어진 차분방식과 격자계로써 공간과 시간상에 존재하는 미소크기의 파동성분들을 충분히 구현하여야 만족할 만한 수치해를 얻을 수 있기 때문이다. 본 연구에서는, 그러한 유한차분법 중 최근에 관심의 대상이 되고있는 삼각(tridiagonal)또는 오각(pentadiagonal) 집적유한차분법(compact finite difference scheme)이 최대의 해상도를 갖도록 하는 수학적인 방법을 개발하고, 이 방법으로써 새롭게 집적유한차분법을 최적화하였다. 개발된 최적화 방법은, 푸리에 해석법(Fourier analysis)을 통하여 파동수(wavenumber) 영역에서 수학적으로 계산된 위상오차(phase error)를 최소화하는 것이며, 이러한 개념과 방법은 본 연구에서 처음으로 집적유한차분법에 적용되었다. 여러가지 절단정확도(truncation order)에 대해서 최적화 된 집적유한차분법들이 실제 공간과 시간상에서 보여주는 정확도와 오차특성을 알아보기 위하여, 이 방법들을 1차원 선형파동방정식에 적용하였고, 이 결과를 통하여 가장 정확하고 효과적인 절단정확도의 집적유한차분법을 선별하였다. 특히, 오각(pentadiagonal)법에 비해 더욱 효율적인 6차 삼각(tridiagonal)법을 1차원 Euler방정식에 적용하여, 비선형 파동에 대한 모사를 수행할 수 있었다.

  • PDF

RESEARCH ON THE WAVELET METHOD FOR THE IMPROVEMENT OF COMPUTATIONAL EFFICIENCY OF TWO DIMENSIONAL FLOW PROBLEMS (2차원 비정상 유동 해석 효율 향상을 위한 Wavelet 기법 응용 연구)

  • Kang, H.M.;Hong, S.W.;Jeong, J.H.;Kim, K.H.;Lee, D.H.;Lee, D.H.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2008.03a
    • /
    • pp.42-49
    • /
    • 2008
  • A wavelet method is presented in order to improve the computational efficiency of two dimensional unsteady flow problems while maintaining the order of accuracy of conventional CFD schemes. First, by using the interpolating wavelet transformation including decomposition and thresholding, an adaptive dataset to a solution is constructed. Then, inviscid and viscous fluxes are calculated only at the points within an adaptive dataset, which enhances the computational efficiency. Second, thresholding step is modified to maintain the spatial and temporal accuracy of conventional CFD schemes automatically by selecting the threshold value between user-defined value and the magnitude of spatial or temporal truncation error. The wavelet method suggested in this study is successfully applied to various unsteady flow problems and it is shown that the computational efficiency is enhanced with maintaining the computational accuracy of CFD schemes.

  • PDF

RESEARCH ON THE WAVELET METHOD FOR THE IMPROVEMENT OF COMPUTATIONAL EFFICIENCY OF TWO DIMENSIONAL FLOW PROBLEMS (2차원 비정상 유동 해석 효율 향상을 위한 Wavelet 기법 응용 연구)

  • Kang, H.M.;Hong, S.W.;Jeong, J.H.;Kim, K.H.;Lee, D.H.;Lee, D.H.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 2008.10a
    • /
    • pp.42-49
    • /
    • 2008
  • A wavelet method is presented in order to improve the computational efficiency of two dimensional unsteady flow problems while maintaining the order of accuracy of conventional CFD schemes. First, by using the interpolating wavelet transformation including decomposition and thresholding, an adaptive dataset to a solution is constructed. Then, inviscid and viscous fluxes are calculated only at the points within an adaptive dataset, which enhances the computational efficiency. Second, thresholding step is modified to maintain the spatial and temporal accuracy of conventional CFD schemes automatically by selecting the threshold value between user-defined value and the magnitude of spatial or temporal truncation error. The wavelet method suggested in this study is successfully applied to various unsteady flow problems and it is shown that the computational efficiency is enhanced with maintaining the computational accuracy of CFD schemes.

  • PDF