• Geosystem Engineering
• /
• 제21권6호
• /
• pp.309-317
• /
• 2018

We employ a three-dimensional indirect boundary element method (BEM) to simulate temperature change around an underground liquefied natural gas storage cavern. The indirect BEM (IBEM) uses fictitious heat source strength on boundary elements as basic variables which are solved from equations of boundary conditions and then used to compute the temperature change at other points in the considered problem domain. The IBEM requires evaluation of singular integration for temperature change due to heat conduction from a constant heat source on a planar (triangular) region. The singularity can be eliminated by a semi-analytical integration scheme. However, it is found that the semi-analytical integration scheme yields sharp temperature gradient for points close to vertices of triangle. This affects the accuracy of heat flux, if they are evaluated by finite difference method at these points. This difficulty can be overcome by a combination of using a direct numerical integration for these points and the semi-analytical scheme for other points distance away from the vertices. The IBEM and the hybrid integration scheme have been verified with an analytic solution and then used to the application of the underground storage.

• Journal of Astronomy and Space Sciences
• /
• 제28권1호
• /
• pp.37-54
• /
• 2011

Two semi-analytic solutions for a perturbed two-body problem known as Lagrange planetary equations (LPE) were compared to a numerical integration of the equation of motion with same perturbation force. To avoid the critical conditions inherited from the configuration of LPE, non-singular orbital elements (EOE) had been introduced. In this study, two types of orbital elements, classical Keplerian orbital elements (COE) and EOE were used for the solution of the LPE. The effectiveness of EOE and the discrepancy between EOE and COE were investigated by using several near critical conditions. The near one revolution, one day, and seven days evolutions of each orbital element described in LPE with COE and EOE were analyzed by comparing it with the directly converted orbital elements from the numerically integrated state vector in Cartesian coordinate. As a result, LPE with EOE has an advantage in long term calculation over LPE with COE in case of relatively small eccentricity.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.278-285
• /
• 2011

In the present study, the numerical prediction of the oil amount leaked from the hole of a damaged tank is investigated using the improved MPS (Moving Particle Semi-implicit) method, which was originally proposed by Koshizuka and Oka (1996) for incompressible flow. The governing equations, which consist of the continuity and Navier-Stokes equations, are solved by Lagrangian moving particles, and all terms expressed by differential operators should be replaced by the particle interaction models based on a Kernel function. The simulation results are validated though the comparison with the analytic solution based on Torricelli's equilibrium relation. Furthermore, a series of numerical simulations under the various conditions are performed in order to estimate more accurately the initial amount of leaked oil.

• 한국항공우주학회지
• /
• 제50권8호
• /
• pp.531-540
• /
• 2022

본 논문에서는 IMM 필터 기반으로 장사정포의 탄종을 식별하고 탄착점을 신속하게 예측하는 알고리즘을 제시한다. 탄도궤적 방정식을 시스템 모델로 사용하고, 각각 다른 탄도계수 값을 갖는 3가지 모델을 IMM 필터에 적용한다. 가속도를 중력, 공기저항, 양력에 의한 3가지 성분으로 나누고 양력가속도를 새로운 상태변수로 추가하여 추정한다. 속도벡터와 양력가속도가 수직이라는 운동학 조건을 유사 측정값으로 추가한 측정방정식을 다룬다. IMM 필터를 통해 추정된 상태변수와 모드 확률이 가장 높은 모델의 탄도계수를 기반으로 탄착점을 예측한다. 탄착점 예측을 위해 일반적으로 사용되는 룽게-쿠타 수치적분 대신, 준해석적인 방법을 사용하여 적은 계산량으로 탄착점을 예측할 수 있음을 설명한다. 마지막으로 최소제곱법을 이용한 상태변수 초기화 방법에 대해 제안하고 성능을 확인하였다. 탄종식별, 탄착점 예측 및 초기화를 포함한 통합 알고리즘을 제시하고 시뮬레이션을 통해 제안한 방법의 타당성을 검증하였다.