• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.660-665
• /
• 2004

본 연구에서는 비선형적인 메커니즘을 갖는 수문현상의 불확실성을 해석하고자 하는 목적으로 새로운 개념의 지배방정식이 유도된다. 제안된 모형의 불확실성은 토양 특성치의 공간적 변동성에 기인하고 있는 것으로 가정하여, 유도된 방정식은 Fokker-Planckl 방정식의 형태를 가지고 있다. 실제 유역단위에서 토양 내 수분 흐름의 연직방향 흐름을 모의하기 위해 미소단위에서 유도된 Richards 방정식은 토양의 공간적 변동성으로 말미암아 불확실한 매개변수를 갖는 비선형 추계학적 편미분방정식의 형태를 갖게 된다. 이는 먼지 수직 방향적분을 통하여 단순화된 비선형 추계학적 상미분방정식으로 전환되고, 이렇게 전환된 비선형 추계학적 상미분방정식은 다시 추계학적 Liouville 방정식을 이용하여 선형 추계학적 편미분방정식으로 전환되어진다. 최종적으로 cumulant 급수방법을 이용하여 상기 방정식을 선형 결정론적 편미분방정식으로 전환시킴으로써, 강우 시 토양 내 수분 침투현상을 모형화할 경우 유역단위에서 토양의 공간적 변동성을 설명할 수 있는 지배방정식을 유도할 수 있다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.414-414
• /
• 2020

고전 천수방정식을 적용한 2차원 도시홍수 모델링에서는 지형의 정확한 반영을 위해 고해상도 격자가 요구되며 이는 많은 계산시간과 노력을 필요로 한다. 최근에는 다공성 천수방정식을 적용한 도시홍수해석으로 많은 계산 노력이 요구되는 도시홍수모델링의 한계를 극복한 연구가 많이 이루어지지고 있다. 이러한 연구는 도시 홍수에서 흐름이 공간적으로 변화할 때 불균일 공극이 존재하므로 격자의 크기를 다르게 하여 이러한 불균일성을 해결하고자 하는 등방성 천수모형의 적용에서 시작되었다. 하지만, 등방성 공극을 고려한 도시홍수 해석모형은 대표요소체적(REV)보다 더 큰 격자의 적용을 해야 하는 제한성을 가진다. 반면, 비등방성 공극은 대표요소체적의 적용이 필요하지 않아 불균일 공극의 크기에 관계없이 이론상으로는 동일한 해상도의 격자가 사용가능하긴 하지만, 실제 도시홍수 해석에서 중요하면서도 도전적인 연구이다. 본 연구에서는 도시홍수의 효율적 계산을 위해 비등방성 공극을 고려한 적분형 다공성 천수방정식을 기반으로 하는 2차원 도시홍수 해석모형을 개발하였다. 모형의 개발을 위해, 적용 격자내에서 도시지역의 건물이 차지하는 길이 및 면적을 산정하고 그 값을 2차원 천수방정식에 적용 가능하도록 체적공극(𝜙_j)와 면적공극(𝜓_k)을 2차원 고전 천수방정식에 추가하였다. 개발모형은 고전 천수방정식, 등방성 공극 고려(미분형 다공성) 천수방정식 및 비등방성 공극 고려(적분형 다공성) 천수방정식의 적용이 가능하여, 각 모형에 적합한 2차원 격자 생성, 각 모형의 매개변수를 보정 그리고 정확성, 효율성, 적용성이 비교 가능하다. 각 모형의 정확성과 효율성 비교를 위해 3가지의 오차 비교 (구조적 오차, 격자크기 오차, 공극 모형 오차), 계산시간 비교, 공간 변동성 검증을 위한 수심 종단형상 비교하였다.

• Proceedings of the Korean Institute of Communication Sciences Conference
• /
• 1990.06a
• /
• pp.403-406
• /
• 1990

• Proceedings of the Korean Institute of Communication Sciences Conference
• /
• 1987.10a
• /
• pp.258-261
• /
• 1987

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.3
• /
• pp.8-16
• /
• 2002

The implicit formulation of rotor dynamics for helicopter flight simulation has been derived and and presented. The generalized vector kinematics regarding the relative motion between coordinates were expressed as a unified matrix operation and applied to get the inertial velocities and accelerations at arbitaty rotor blade span position. Based on these results the rotor aeromechanic equations for flapping dynamics, lead-lag dynamics and torque dynamics were formulated as an implicit form. Spatial integration methods of rotor dynamic equations along blade span and the expanded applicability of the present implicit formulations for arbitrary hings geometry and hinge sequences have been investigated. Time integration methods for present DAE(Differential Algebraic Equation) to calculate dynamic response calculation are recommenaded as future works.

• Computational Structural Engineering
• /
• v.7 no.3
• /
• pp.83-93
• /
• 1994

To design and study the dynamic performance of safety structures, crash tests are needed. Method to get the angular accelerations at the time of impact by integating the Euler equations are introduced. Numerically stable 9-array system contains several 7 and 8-array sub-systems in it. Numerical stability of those latent sub-systems are studied using test files. All of the 8-array subsystems were found to be numerically stable. Six of the 7-array sub-systems were stable and other six of the 7-array sub-systems were unstable. Using this findings fail-safe measurement system can be developed.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1429-1437
• /
• 1992

A new finite element formulation based on the relaxation type hereditary integral is presented for a time-domain analysis of isotropic, linear viscoelastic problems. The semi-discrete variational approximation and elastic-viscoelastic correspondence principle are used in the theoretical development of the proposed method. In a time-stepping procedure of final, linear algebraic system equations, only a small additional computation for past history is required since the equivalent stiffness matrix is constant. The viscoelasticity matrices are derived and the stress computation algorithm is given in matrix form. The effect of time increment and Gauss point numbers on the numerical accuracy is examined. Two dimensional numerical examples of plane strain and plane stress are solved and compared with the analytical solutions to demonstrate the versatility and accuracy of the present method.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• 2009.10a
• /
• pp.821-824
• /
• 2009

Conformal mapping is useful to solve problems in physics, engineering and so on. This paper is to discuss the numerical conformal mapping from the unit disk onto Jordan region, which can be solved by Theodorsen equation. Wegmann's method has been known as the most efficient one for the Theodorsen equation. However, we found divergence through numerical experiments by the iterative method of Wegmann. The divergence occurs especially when some degree of difficulty is high. We analyze the cause of divergence and propose an improved method by applying a low frequency pass filter to Wegmann's method. By this proposed method we can get a stable convergence for all the problems which was unstable with the Wegmann's method.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.3
• /
• pp.41-48
• /
• 2014

Experimental study of sedimentation and self-weight consolidation has been primary research area in dredged soil. However, good quality of the dredged soil and minimum water pollution caused by the pumping of reclaimed soil require intensive study of the flow characteristics of dredged material due to dumping. In this study, continuity and the equilibrium equations for mass flow assuming single phase was derived to simulate mass flow in dredged containment area. To optimize computation and modeling time for three dimensional geometry and boundary conditions, depth integration is applied to governing equations to consider three dimensional topography of the site. Petrov-Galerkin formulation is applied in spatial discretization of governing equations. Generalized trapezoidal rule is used for time integration, and Newton iteration process approximated the solution. DG and CDG technique were used for weighting matrix in discontinuous test function in dredged flow analysis, and numerical stability was evaluated by performed a square slump simulation. A comparative analysis for numerical methods showed that DG method applied to SU / PG formulation gives minimal pseudo oscillation and reliable numerical results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.