• Journal of Korean Association for Spatial Structures
• /
• v.2 no.3 s.5
• /
• pp.115-122
• /
• 2002

It is often hard to obtain analytical solutions of boundary value problems of shells. Introducing some approximations into the governing equations may allow us to get analytical solutions of boundary value problems. Instead of an analytical procedure, we can apply a numerical method to the governing equations. Since the governing equations of shells of revolution under symmetric load are expressed by ordinary differential equations, a numerical solution of ordinary differential equations is applicable to solve the equations. In this paper, the governing equations of orthotropic spherical shells under symmetric load are derived from the classical theory based on differential geometry, and the analysis is numerically carried out by computer program of Runge-Kutta methods. The numerical results are compared to the solutions of a commercial analysis program, SAP2000, and show good agreement.

• The Korean Journal of Rheology
• /
• v.8 no.3_4
• /
• pp.187-198
• /
• 1996

최근에 본 연구자에 의해서 단축 스크류 공정에서 카오스 스크류라고 명명되어진 카오스 혼합장치가 성공적으로 개발되었다. 기하학적 조건이나 공정조건에 대한 설계변수로 카오스 스크류를 설계하기 위하여 체류시간, 포인카레 단면 그리고 혼합패턴등에 대한 계산 과 해석이 이루어져야 하는데 이를 단지 Runge-Kutta 방법에 의해 속도장을 적분한다면 상당한 계산시간이 소비된다. 이러한 수치문제를 극복하기 위하여 본논문에서는 새로운 사 상법을 제안한다. 이 방법으　 사용하면 벽면 근처의 특이점 영역에서도 수치문제가 해결된 다. 본 논문에서 제안하는 수치사상법은 Runge-Kutta 방법에 비하여 수치계산의 효율성과 정확도 면에서, 특히 유안요소법으로 얻은 속도장에 대하여 우수함이 밝혀졌다. 이러한 사상 법은 공간주기 유동장뿐만 아니라 시간주기 유동장에서 적용할수 있다.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.87-99
• /
• 1994

Reynolds Averaged Navier-Stokes equations are solved numerically for the computation of turbulent flow around a Wigley double model. A second order finite difference method is applied for the spatial discretization on the nonstaggered grid system and 4-stage Runge-Kutta scheme for the numerical integration in time. In order to increase the time step, residual averaging scheme of Jameson is adopted. Pressure field is obtained by solving the pressure-Poisson equation with the appropriate Neumann boundary condition. For the turbulence closure, 0-equation turbulence model of Baldwin-Lomax is used. Numerical computation is carried out for the Reynolds number of 4.5 million. Comparisons of the computed results with the available experimental data show good agreements for the velocity and pressure distributions.

• 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
• /
• 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법과 제안한 방법을 적용시켜서 제안한 방법의 정화함을 보여준다.

• 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.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.28 no.6
• /
• pp.573-578
• /
• 2010

Precise determination of satellite positions is necessary to improve positioning accuracy in GNSS. In this study, GLONASS orbits were predicted from broadcast ephemeris using the 4th-order Runge-Kutta numerical integration method and their accuracy dependence on the integration step and the integration time was analyzed. The 3D RMS (Root Mean Square) differences between the results from I-second integration step and 300-second integration step was about 3 cm, but the processing time was one hundred times less for the I-second integration time case. For trials of different integration times, the 3D RMS errors were 8.3 m, 187.3 m, and 661.5 m for 30-, 150-, and 300-minutes of integration time, respectively. Though this integration-time analysis, we concluded that the accuracy gets higher with a shorter integration time. Thus we suggest forward and backward integration methods to improve GLONASS positioning accuracy, and with this method we can achieve a 5-meter level of 3-D orbit accuracy.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.11 no.2
• /
• pp.107-115
• /
• 1999

In order to develop a convenient method for the estimation of transport distance of settling stones in quiescent water or flowing water, introduced was the simple but relatively accurate equation of drag coefficient. The equation of drag coefficient introduced was confirmed to give relatively accurate evaluation for the drag force of smooth-surface sphere, and the effects of surface roughness and shape can be considered by adjusting empirical parameters. A theoretical equation has been developed for the settling velocity or settling distance of smooth-surface sphere in quiescent fluid, and the computation results have been obtained by adjusting the empirical parameter for the settling distance of stone in quiescent water. The 2nd order ordinary differential equation has been developed for the case of settling stones in flowing fluid, and a numerical model has been developed by using Runge-Kutta method for its solution. A number of cases have been tested by adjusting the empirical parameter.

• Journal of Korea Water Resources Association
• /
• v.39 no.7 s.168
• /
• pp.593-603
• /
• 2006

The basin response to storm is regarded as nonlinearity inherently. In addition, the consistent nonlinearity of hydrologic system response to rainfall has been very tough and cumbersome to be treated analytically. The thing is that such nonlinear models have been avoided because of computational difficulties in identifying the model parameters from recorded data. The parameters of nonlinear system considered as dynamic effects in the conceptual model are optimized as the sum of errors between the observed and computed runoff is minimized. For obtaining the optimal parameters of functions, the historical data for the Bocheong watershed in the Geum river basin were tested by applying the numerical methods, such as quasi-linearization technique, Runge-Kutta procedure, and pattern-search method. The estimated runoff carried through from the storage function with dynamic effects was compared with the one of 1st-order differential equation model expressing just nonlinearity, and also done with Nash model. It was found that the 2nd-order model yields a better prediction of the hydrograph from each storm than the 1st-order model. However, the 2nd-order model was shown to be equivalent to Nash model when it comes to results. As a result, the parameters of nonlinear 2nd-order differential equation model performed from the present study provided not only a considerable physical meaning but also a applicability to Korean watersheds.

• 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.