• Title/Summary/Keyword: Runge-Kutta법

Search Result 4, Processing Time 0.017 seconds

A New Numerical Method for Solving Differential Equation by Quadratic Approximation (포물선 근사법에 의한 상태방정식의 새로운 수치해석적 접근법에 관한 연구)

  • Lee, Jong-Gi;Kwon, Yong-Jun;Choi, Byoung-Kon;Moon, Young-Hyun
    • Proceedings of the KIEE Conference
    • /
    • 2003.07a
    • /
    • pp.107-109
    • /
    • 2003
  • 전력계통의 과도 안정도 해석의 접근방법에는 SI(Simultaneous Implicit)법과 PE(Partitioned Explicit)법 두 가지방법을 사용해오고 있다. SI법에는 Trapezoidal법 등이 있고, PE법에는 Runge-Kutta법, Euler법등이 사용되고 있다. SI법인 Trapezoidal법은 PE법의 Runge-Kutta법 또는 Euler법에 비해 시간간격을 크게 해서 계산속도를 줄일 수 있다는 장점이 있지만, 정화도면에서는 신뢰한 수 없는 단점이 있다. 이 논문에서는 포물선 사법을 이용하여 Trapezoidal법의 정확도를 개선학 수 있는 방법을 제시하고 명확한 수학적 증명을 통해 타당성을 보여준다. 연속함수와 불연속함수에 대해서 Runge-Kutta법과 Trapezoidal법과 제안한 방법을 적용시켜서 제안한 방법의 정화함을 보여준다.

  • PDF

Development of a Three-Dimensional Euler Solver for Analysis of Basic Contraction Flow (수축부 기초 유동 해석을 위한 삼차원 Euler 방정식 풀개 개발)

  • Kim J.;Kim H. T.
    • Journal of computational fluids engineering
    • /
    • v.2 no.1
    • /
    • pp.8-12
    • /
    • 1997
  • The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

  • PDF

Buckling Loads and Post-Buckling Behaviors of Shear Deformable Columns with Regular Cross-Section (전단변형을 고려한 정다각형 단면 기둥의 좌굴하중 및 후좌굴 거동)

  • Lee, Byeoung Koo;Lee, Tae Eun;Kwon, Yun Sil;Kim, Sun Gi
    • Journal of Korean Society of Steel Construction
    • /
    • v.13 no.6
    • /
    • pp.683-691
    • /
    • 2001
  • Numerical methods are developed for solving the elastica and buckling load of tapered columns with shear deformation, subjected to a compressive end load. The linear, parabolic and sinusoidal tapers with the regular polygon cross-sections are considered, whose material volume and span length are always held constant. The differential equations governing the elastica of buckled column are derived. The Runge-Kutta method is used to integrate the differential equations, and the Regula-Falsi method is used to determine the rotation at left end and the buckling load, respectively. The numerical methods developed herein for computing the elastica and the buckling loads of the columns are found to be efficient and reliable.

  • PDF

Strongest Simple Beams with Constant Volume (일정체적 단순지지 최강보)

  • Lee, Byoung Koo;Lee, Tae Eun;Kim, Young Il
    • KSCE Journal of Civil and Environmental Engineering Research
    • /
    • v.29 no.2A
    • /
    • pp.155-162
    • /
    • 2009
  • This paper deals with the strongest beams with the solid regular polygon cross-section, whose volumes are always held constant. The differential equation of the elastic deflection curve of such beam subjected to the concentrated and trapezoidal distributed loads are derived and solved numerically. The Runge-Kutta method and shooting method are used to integrate the differential equation and to determine the unknown initial boundary condition of the given beam. In the numerical examples, the simple beams are considered as the end constraint and also, the linear, parabolic and sinusoidal tapers are considered as the shape function of cross sectional depth. As the numerical results, the configurations, i.e. section ratios, of the strongest beams are determined by reading the section ratios from the numerical data related with the static behaviors, under which static maximum behaviors become to be minimum.