• Journal of Positioning, Navigation, and Timing
• /
• v.8 no.4
• /
• pp.201-208
• /
• 2019

Numerical integration is necessary for satellite orbit determination and its prediction. The numerical integration algorithm can be divided into single-step and multi-step method. There are lots of single-step and multi-step methods. However, the Runge-Kutta method in single-step and the Adams method in multi-step are generally used in global navigation satellite system (GNSS) satellite orbit. In this study, 4th and 8th order Runge-Kutta methods and various order of Adams-Bashforth-Moulton methods were used for GLObal NAvigation Satellite System (GLONASS) orbit integration using its broadcast ephemeris and these methods were compared with international GNSS service (IGS) final products for 7days. As a result, the RMSE of Runge-Kutta methods were 3.13m and 4th and 8th order Runge-Kutta results were very close and also 3rd to 9th order Adams-Bashforth-Moulton results. About result of computation time, this study showed that 4th order Runge-Kutta was the fastest. However, in case of 8th order Runge-Kutta, it was faster than 14th order Adams-Bashforth-Moulton but slower than 13th order Adams-Bashforth-Moulton in this study.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.4
• /
• pp.238-266
• /
• 2013

The Interaction Picture (IP) method is a valuable alternative to Split-step methods for solving certain types of partial differential equations such as the nonlinear Schr$\ddot{o}$dinger equation or the Gross-Pitaevskii equation. Although very similar to the Symmetric Split-step (SS) method in its inner computational structure, the IP method results from a change of unknown and therefore do not involve approximation such as the one resulting from the use of a splitting formula. In its standard form the IP method such as the SS method is used in conjunction with the classical 4th order Runge-Kutta (RK) scheme. However it appears to be relevant to look for RK scheme of higher order so as to improve the accuracy of the IP method. In this paper we investigate 5th order Embedded Runge-Kutta schemes suited to be used in conjunction with the IP method and designed to deliver a local error estimation for adaptive step size control.

• Steel and Composite Structures
• /
• v.15 no.4
• /
• pp.439-449
• /
• 2013

In this study simply supported nonlinear Euler-Bernoulli beams resting on linear elastic foundation and subjected to the axial loads is investigated. A new kind of analytical technique for a non-linear problem called He's Energy Balance Method (EBM) is used to obtain the analytical solution for non-linear vibration behavior of the problem. Analytical expressions for geometrically non-linear vibration of Euler-Bernoulli beams resting on linear elastic foundation and subjected to the axial loads are provided. The effect of vibration amplitude on the non-linear frequency and buckling load is discussed. The variation of different parameter to the nonlinear frequency is considered completely in this study. The nonlinear vibration equation is analyzed numerically using Runge-Kutta $4^{th}$ technique. Comparison of Energy Balance Method (EBM) with Runge-Kutta $4^{th}$ leads to highly accurate solutions.

• Journal of Ocean Engineering and Technology
• /
• v.4 no.2
• /
• pp.73-84
• /
• 1990

A circular cylinder is oscillated in th otherwise quiescent viscous fluid. Numerical analysis performed for this problem by using the fourth-order Runge-kutta method for the unsteady Navier-stokes equations. For K(Kelegan-Carpenter's No.)=5, the flow developed symmetrically, while for K=10, it revealed random patterns. The coefficient of the rms force is overestimated by 20-30% compared with the experimental result.

• 제어로봇시스템학회:학술대회논문집
• /
• 1991.10a
• /
• pp.383-388
• /
• 1991

This paper analizes two numerical integration methods, both based on the Runge Kutta 4-th order formula for deterministic systems, for digital simulation of a differential equation driven by white noise. It is shown that a "standard' Runge Kutta method for stochasitic systems yields solutions of Stratonovich differential equations, while Riggs and Phillips' method results in solutions of Ito differential equations. Therefore the white noise differential equation must be converted into the equivalent Ito equation before the latter method is used. Digital simulation results for a simple differential equation are also presented.nted.

• 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 computational fluids engineering
• /
• v.15 no.2
• /
• pp.55-61
• /
• 2010

A dependence of weighting parameter in preconditioning method for solving low Mach number flow with incompressible flow nature is investigated. The present preconditioning method employs a finite-difference method applied Roe‘s flux difference splitting approximation with the MUSCL-TVD scheme and 4th-order Runge-Kutta method in curvilinear coordinates. From the computational results of benchmark flows through a 2-D backward-facing step duct it is confirmed that there exists a suitable value of the weighting parameter for accurate and stable computation. A useful method to determine the weighting parameter is introduced. With this method, high accuracy and stable computational results were obtained for the flow with low Mach number in the range of Mach number less than 0.3.

• Journal of Ocean Engineering and Technology
• /
• v.4 no.2
• /
• pp.223-223
• /
• 1990

A circular cylinder is oscillated in th otherwise quiescent viscous fluid. Numerical analysis performed for this problem by using the fourth-order Runge-kutta method for the unsteady Navier-stokes equations. For K(Kelegan-Carpenter''s No.)=5, the flow developed symmetrically, while for K=10, it revealed random patterns. The coefficient of the rms force is overestimated by 20-30% compared with the experimental result.

• Geomechanics and Engineering
• /
• v.16 no.4
• /
• pp.355-361
• /
• 2018

In this paper, nonlinear vibration of Euler-Bernoulli beams resting on linear elastic foundation is studied. It has been tried to prepare a semi-analytical solution for whole domain of vibration. Only one iteration lead us to high accurate solution. The effects of linear elastic foundation on the response of the beam vibration are considered and studied. The effects of important parameters on the ratio of nonlinear to linear frequency of the system are studied. The results are compared with numerical solution using Runge-Kutta $4^{th}$ technique. It has been shown that the Max-Min approach can be easily extended in nonlinear partial differential equations.

• Journal of Ocean Engineering and Technology
• /
• v.9 no.2
• /
• pp.53-60
• /
• 1995

A multi-leg spread mooring system for floating offshore structures is important, but the multi-leg static analysis is complicated due to the nonlinear behavior of each line and the effect of current which affects each line differently. The pretensioned position of the multi-leg mooring system obtained from the static equilibrium condition changes into a different position due to external loads and current. In this paper, the new position and the static tension at each line are caculated. The relation between the initial static equilibrium position and the new position due to the external loads is expressed in terms of the Taylor's series expansion. The Runge-Kutta $4^{th}$ method is employed in analyzing the 3-dimensional static cable nonlinear equations.