• Title/Summary/Keyword: Numerical inversion method

Search Result 148, Processing Time 0.027 seconds

Analysis of Linear Consolidation Problems by the Boundary Element Method (경계요소법에 의한 선형 압밀문제의 해석)

  • 서일교
    • Computational Structural Engineering
    • /
    • v.8 no.4
    • /
    • pp.129-136
    • /
    • 1995
  • This paper presents a boundary element method for obtaining approximate solutions of 2-dimensional consolidation problems based on the Biot's linear theory. Laplace transform is applied to differential equation system in order to eliminate the time dependency. The boundary integral equations in transformed space are formulated and the fundamental solutions are shown in a closed form. In order to convert the transformed solutions to the ones in real space, the Hosono's numerical Laplace transform inversion method is applied. As a numerical example, a half-space consolidation problem subjected to a strip local load is selected and the applicability of the method is demonstrated through the comparison with the exact solutions.

  • PDF

Axisymmetric analysis of multi-layered transversely isotropic elastic media with general interlayer and support conditions

  • Lee, J.S.;Jiang, L.Z.
    • Structural Engineering and Mechanics
    • /
    • v.2 no.1
    • /
    • pp.49-62
    • /
    • 1994
  • Based on the transfer matrix approach and integral transforms, a solution method is developed for the stress analysis of axisymmetrically loaded transversely isotropic elastic media with generalized interlayer and support conditions. Transfer functions (Green's functions in the transformed domain) are obtained in explicit integral form. For several problems of practical interest with different loading and support conditions, solutions are worked out in detail. For the inversion operation, an efficient technique is introduced to remedy the slow convergence of numerical integrals involving oscillating functions. Several illustrative examples are considered and numerical results are presented.

Application of Effective Regularization to Gradient-based Seismic Full Waveform Inversion using Selective Smoothing Coefficients (선택적 평활화 계수를 이용한 그래디언트기반 탄성파 완전파형역산의 효과적인 정규화 기법 적용)

  • Park, Yunhui;Pyun, Sukjoon
    • Geophysics and Geophysical Exploration
    • /
    • v.16 no.4
    • /
    • pp.211-216
    • /
    • 2013
  • In general, smoothing filters regularize functions by reducing differences between adjacent values. The smoothing filters, therefore, can regularize inverse solutions and produce more accurate subsurface structure when we apply it to full waveform inversion. If we apply a smoothing filter with a constant coefficient to subsurface image or velocity model, it will make layer interfaces and fault structures vague because it does not consider any information of geologic structures and variations of velocity. In this study, we develop a selective smoothing regularization technique, which adapts smoothing coefficients according to inversion iteration, to solve the weakness of smoothing regularization with a constant coefficient. First, we determine appropriate frequencies and analyze the corresponding wavenumber coverage. Then, we define effective maximum wavenumber as 99 percentile of wavenumber spectrum in order to choose smoothing coefficients which can effectively limit the wavenumber coverage. By adapting the chosen smoothing coefficients according to the iteration, we can implement multi-scale full waveform inversion while inverting multi-frequency components simultaneously. Through the successful inversion example on a salt model with high-contrast velocity structures, we can note that our method effectively regularizes the inverse solution. We also verify that our scheme is applicable to field data through the numerical example to the synthetic data containing random noise.

Inverse Scattering of Two-Dimensional Objects Using Linear Sampling Method and Adjoint Sensitivity Analysis

  • Eskandari, Ahmadreza;Eskandari, Mohammad Reza
    • Journal of Electrical Engineering and Technology
    • /
    • v.10 no.1
    • /
    • pp.308-313
    • /
    • 2015
  • This paper describes a technique for complete identification of a two-dimensional scattering object and multiple objects immersed in air using microwaves where the scatterers are assumed to be a homogenous dielectric medium. The employed technique consists of initially retrieving the shape and position of the scattering object using a linear sampling method and then determining the electric permittivity and conductivity of the scatterer using adjoint sensitivity analysis. Incident waves are assumed to be TM (Transverse Magnetic) plane waves. This inversion algorithm results in high computational speed and efficiency, and it can be generalized for any scatterer structure. Also, this method is robust with respect to noise. The numerical results clearly show that this hybrid approach provides accurate reconstructions of various objects.

Experimental Study on a Monte Carlo-based Recursive Least Square Method for System Identification (몬테카를로 기반 재귀최소자승법에 의한 시스템 인식 실험 연구)

  • Lee, Sang-Deok;Jung, Seul
    • The Transactions of The Korean Institute of Electrical Engineers
    • /
    • v.67 no.2
    • /
    • pp.248-254
    • /
    • 2018
  • In this paper, a Monte Carlo-based Recursive Least Square(MC-RLS) method is presented to directly identify the inverse model of the dynamical system. Although a RLS method has been used for the identification based on the deterministic data in the closed loop controlled form, it would be better for RLS to identify the model with random data. In addition, the inverse model obtained by inverting the identified forward model may not work properly. Therefore, MC-RLS can be used for the inverse model identification without proceeding a numerical inversion of an identified forward model. The performance of the proposed method is verified through experimental studies on a control moment gyroscope.

Indirect Input Identification by Modal Filter Technique (모드필터방법에 의한 간접적 입력규명)

  • 김영렬;김광준
    • Journal of KSNVE
    • /
    • v.9 no.2
    • /
    • pp.377-386
    • /
    • 1999
  • This paper is a study on model method for estimating system inputs from vibration responses, which is one of indirect input identification methods in frequency domain. The method has advantages over direct inverse method especially when points of operational inputs are inaccessible so that artificial excitation forces cannot be applied to obtain frequency response functions of the complete system. Procedures of extended modal model method are proposed and checked by numerical experiment. Mechanisms of error propagation, i.e., how errors in modal parameters such as poles nad mode shape vectors affect estimation of the input forces, are illustrated. Then, in order to counteract the error propagation, discrete modal filter approach is taken in this paper to compute the inversion of modal matrix in which the most serious errors seem to be generated. Further, a Reduced form of Modified Reciprocal Modal Vector(RMRMV) is proposed for estimating multiple inputs. It is shown to have smaller orthogonality error than MRMV.

  • PDF

Solving Time-dependent Schrödinger Equation Using Gaussian Wave Packet Dynamics

  • Lee, Min-Ho;Byun, Chang Woo;Choi, Nark Nyul;Kim, Dae-Soung
    • Journal of the Korean Physical Society
    • /
    • v.73 no.9
    • /
    • pp.1269-1278
    • /
    • 2018
  • Using the thawed Gaussian wave packets [E. J. Heller, J. Chem. Phys. 62, 1544 (1975)] and the adaptive reinitialization technique employing the frame operator [L. M. Andersson et al., J. Phys. A: Math. Gen. 35, 7787 (2002)], a trajectory-based Gaussian wave packet method is introduced that can be applied to scattering and time-dependent problems. This method does not require either the numerical multidimensional integrals for potential operators or the inversion of nearly-singular matrices representing the overlap of overcomplete Gaussian basis functions. We demonstrate a possibility that the method can be a promising candidate for the time-dependent $Schr{\ddot{o}}dinger$ equation solver by applying to tunneling, high-order harmonic generation, and above-threshold ionization problems in one-dimensional model systems. Although the efficiency of the method is confirmed in one-dimensional systems, it can be easily extended to higher dimensional systems.

An efficient method for computation of receptances of structural systems with sparse, non-proportional damping matrix (성긴 일반 감쇠행렬을 포함하는 구조물에 대한 효율적인 주파수 응답 계산 방법)

  • Park, Jong-Heuck;Hong, Seong-Wook
    • Journal of the Korean Society for Precision Engineering
    • /
    • v.12 no.7
    • /
    • pp.99-106
    • /
    • 1995
  • Frequency response functions are of great use in dynamic analysis of structural systems. The present paper proposes an efficient method for computation of the frequency rewponse functions of linear structural dynamic models with a sparse, non-proportional damping matrix. An exact condensation procedure is proposed which enables the present method to condense the matrices without resulting in any errors. Also, an iterative scheme is proposed to be able to avoid matrix inversion in computing frequency response matrix. The proposed method is illustrated through a numerical example.

  • PDF

Simulation for nanoimprint lithography process using temperature controlled nonequilibrium molecular dynamics (온도 제어 비평형 분자동역학 방법을 이용한 나노임프린트 리소그라피 공정의 전산모사)

  • Kwon, Sung-Jin;Lee, Young-Min;Im, Se-Young
    • Proceedings of the KSME Conference
    • /
    • 2007.05a
    • /
    • pp.332-336
    • /
    • 2007
  • Temperature is an essential process variable in nanoimprint lithography(NIL) where the temperature varies between room temperature and above the glass transition temperature. To simulate NIL process, we employ both the Nose-Poincare method for temperature controlled molecular dynamics(MD) and force field for polymer material i.e. polymethyl methacrylate(PMMA), which is most widely selected as NIL resist. Nose-Poincare method, which convinces the conservation of Hamiltonian structure and time-reversal symmetry, overcomes the drawbacks inherent in the conventional methods such as Nose thermostat and Nose-Hoover thermostat. Thus, this method exhibits enhanced numerical stability even when the temperature fluctuation is large. To describe PMMA, we adopt the force field which account for bond stretch, bending, torsion, inversion, partial charge, and van der Waals energy.

  • PDF

On the Numerical Inversion of the Laplace Transform by the Use of an Optimized Legendre Polynomial

  • Al-Shuaibi, Abdulaziz
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • v.4 no.1
    • /
    • pp.49-65
    • /
    • 2000
  • A method for inverting the Laplace transform is presented, using a finite series of the classical Legendre polynomials. The method recovers a real-valued function f(t) in a finite interval of the positive real axis when f(t) belongs to a certain class ${\mathcal{W}}_{\beta}$ and requires the knowledge of its Laplace transform F(s) only at a finite number of discrete points on the real axis s > 0. The choice of these points will be carefully considered so as to improve the approximation error as well as to minimize the number of steps needed in the evaluations. The method is tested on few examples, with particular emphasis on the estimation of the error bounds involved.

  • PDF