• Journal of the Korean Geophysical Society
• /
• v.5 no.3
• /
• pp.233-245
• /
• 2002

The spectral-analysis-of-surface-waves (SASW) method is a nondestructive testing method based upon generation and detection of elastic stress waves. SASW is widely used as one of the techniques to determine stiffness profile in engineering geophysics. The essential steps involved are construction of an experimental dispersion curve from data collected in situ, and inversion of the dispersion curve to determine the stiffness profile. The main object of this study is to derive an analytical Jacobian for the inversion. If we set the subsurface to N homogeneous layer, it could save 2N times Jacobian calculation compared to numerical jacobian calculation during inversion. To reconstruct a stiffness profile, constrained damped least square method was applied for the inversion. The algorithm was tested for the numerical data and for the real asphalt and tunnel data, which were able to verify the stiffness profile. The stiffness profile reconstructed by the algorithm showed the possibility to appraise the soundness of tunnel with applications SASW.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.7
• /
• pp.1-8
• /
• 2003

Numerical characteristics of air chemistry associated with hypersonic flows are described and are compared with those of hydrogen oxygen combustion, applying the partially implicit time integration method to air chemistry. This paper reveals that the time integration of air chemistry needs a chemical Jacobian for stable calculations. However the positive real eigenvalues in air chemistry are relatively smaller than those of hydrogen combustion, and the numerical integration is less sensitive than that with combustion. lt is also found that the application of the partia1ly irnplicit method reduces the computing time without numerical instabilities.

• The Bulletin of The Korean Astronomical Society
• /
• v.35 no.2
• /
• pp.69.1-69.1
• /
• 2010

장주기 Mira 변광성의 외피층에서 발생하는 SiO 메이저에 대해 non-local한 1차원 수치계산 결과를 처음으로 제시한다. 별의 유체역학적인 맥동모델에서 얻은 시간에 따른 별 주위 가스의 속도, 온도, 밀도 분포를 사용하여 계산하였다. 임의의 속도장을 갖는 구형 분자운에서의 복사전달문제를 풀었던 이전 연구에서 개발한 수치계산코드를 이용하였고, 특정 구간에서 급격한 변화를 겪는 물리량을 잘 반영할 수 있도록 수정하였다. 또한 계산에 사용되는 거대희소행렬을 압축희소행렬로 변환하여 메모리를 절약하였고 비선형방정식의 자코비안을 해석적으로 구하여 계산속도를 향상시켰다. v=1, J=1-0, J=2-1과 v=2, J=1-0 SiO 메이저의 공간분포, 상대세기 등에 대해 이전의 LVG 모델을 이용한 연구결과와 정성적으로 비교 논의한다.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.21-28
• /
• 2004

Three-dimensional (3-D) resistivity inversion including a topographic effect can be considered theoretically to be the technique of acquiring the most accurate image in the interpretation of resistivity data, because it includes characteristic image that the actual subsurface structure is 3-D. In this study, a finite-element method was used as the numerical method in modeling, and the efficiency of Jacobian calculation has been maximized with sensitivity analysis for the destination block in inversion process. Also, during the iterative inversion, the resolution of inversion can be improved with the method of selecting the optimal value of Lagrange multiplier yielding minimum RMS(root mean square) error in the parabolic equation. In this paper, we present synthetic examples to compare the difference between the case which has the toprographic effect and the other case which has not the effect in the inversion process.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.12
• /
• pp.1433-1441
• /
• 2009

In multibody dynamic analysis, one of the most important problems is to reduce computation times for real-time simulation. This paper presents the derivation procedure of equations of motion of a 6${\times}$6 autonomous vehicle in terms of chassis local coordinates which do not require coordinates transformation matrix to enhance efficiency for real-time dynamic analysis. Also, equations of motion are derived using the VT(velocity transformation) technique and symbolic computation method coded by MATLAB. The Jacobian matrix of the equations of motion of a system is derived from symbolic operations to apply the implicit integration method. The analysis results were compared with ADAMS results to verify the accuracy and approve the feasibility of real time analysis.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.12
• /
• pp.1153-1161
• /
• 2010

An efficient solution method for harmonic balance techniques with Fourier transform is presented for periodic unsteady flow problems. The present partially-implicit harmonic balance treats the flux terms implicitly and the harmonic source term is solved explicitly. The convergence of the partially Implicit method is much faster than the explicit Runge-Kutta harmonic balance method. The method does not need to compute the additional flux Jacobian matrix from the implicit harmonic source term. Compared with fully implicit harmonic balance method, this partial approach turns out to have good convergence property. Oscillating flows over NACA0012 airfoil are considered to verify the method and to compare with results of explicit R-K(Runge-Kutta) and dual time stepping methods.

• Proceedings of the KIEE Conference
• /
• 2008.11a
• /
• pp.325-326
• /
• 2008

상정사고가 발생하거나 전력계통의 구조가 변경된 경우 빠른 시간 안에 계통의 변화 추이를 알 수 있다면 계통 운영에 많은 이점이 있다. 그러나 기존의 시간모의 방법은 모의에 많은 시간이 소요되어 충분한 검토시간을 갖지 못하는 단점이 있고, 준동적 시간모의는 빠른 시간모의가 가능하나 자코비안 특이점으로 인하여 수치해석 오류가 종종 발생하는 단점을 갖고 있다. 전력연구원에서는 이러한 단점이 극복 가능한 연속법이 적용된 QSS에 대한 연구를 진행 중에 있다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.1
• /
• pp.100-108
• /
• 2004

In this paper, the new casing oscillator, which is a construction machine and which structure is similar to that of a parallel manipulator with redundancy, is proposed. The singularity analysis of this machine is performed by two different methods. First, the singularities are found by the numerical method at configurations where the rank of the Jacobian matrix becomes deficient. The singularities are outside the workspace. To investigate the physical information on these configurations, the singularities are examined by the geometric method at configurations where the casing oscillator cannot resist the external forces and moments applied to the upper platform due to losing static equilibrium. The results of the geometric method are the same as those of the numerical method. It proves that the new casing oscillator is free from the singularity, which causes serious problems to a parallel manipulator.

• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.71-78
• /
• 1998

Inversion schemes of 2-D MT survey data generally take enormous computational time and computer memory. In addition, careful attention must be paid in handling MT data, especially in cases of TM mode, inversion results can be seriously distorted because of static effect caused by current channeling across inhomogeneous surface boundaries. There-fore inversion algorithm using the GRRI scheme for TM mode MT data was implemented. This scheme is based on a perturbation analysis with a locally 2-D analysis and local inversions were sequently performed over each divided section without additional forward modelings. The algorithm was applied to several synthetic data for the purpose of verification of its efficiency and applicability. With less computer resources than conventional schemes, it could handle static effect directly by including current channeling across inhomogeneous boundaries. Thus it is expected to be used for an useful tool such as a real-time inversion scheme in the field.