• Atmosphere
• /
• v.33 no.4
• /
• pp.331-341
• /
• 2023

A numerical forecasting models usually predict future states by performing time integration considering fixed static time-steps. A time-step that is too long can cause model instability and failure of forecast simulation, and a time-step that is too short can cause unnecessary time integration calculations. Thus, in numerical models, the time-step size can be determined by the CFL (Courant-Friedrichs-Lewy)-condition, and this condition acts as a necessary condition for finding a numerical solution. A static time-step is defined as using the same fixed time-step for time integration. On the other hand, applying a different time-step for each integration while guaranteeing the stability of the solution in time advancement is called an adaptive time-step. The adaptive time-step algorithm is a method of presenting the maximum usable time-step suitable for each integration based on the CFL-condition for the adaptive time-step. In this paper, the adaptive time-step algorithm is applied for the Korean Integrated Model (KIM) to determine suitable parameters used for the adaptive time-step algorithm through the monthly verifications of 10-day simulations (during January and July 2017) at about 12 km resolution. By comparing the numerical results obtained by applying the 25 second static time-step to KIM in Supercomputer 5 (Nurion), it shows similar results in terms of forecast quality, presents the maximum available time-step for each integration, and improves the calculation efficiency by reducing the number of total time integrations by 19%.

• Earthquakes and Structures
• /
• v.13 no.3
• /
• pp.279-288
• /
• 2017

Having the ability to keep on yielding stable solutions in problems involving high potential of instability, composite time integration methods have become very popular among scientists. These methods try to split a time step into multiple sub-steps so that each sub-step can be solved using different time integration methods with different behaviors. This paper proposes a new composite time integration in which a time step is divided into two sub-steps; the first sub-step is solved using the well-known Newmark method and the second sub-step is solved using Simpson's Rule of integration. An unconditional stability region is determined for the constant parameters to be chosen from. Also accuracy analysis is perform on the proposed method and proved that minor period elongation as well as a reasonable amount of numerical dissipation is produced in the responses obtained by the proposed method. Finally, in order to provide a practical assessment of the method, several benchmark problems are solved using the proposed method.

• Smart Structures and Systems
• /
• v.14 no.6
• /
• pp.1269-1289
• /
• 2014

With the sub-stepping technique, the numerical analysis in real-time dynamic hybrid testing is split into the response analysis and signal generation tasks. Two target computers that operate in real-time may be assigned to implement these two tasks, respectively, for fully extending the simulation scale of the numerical substructure. In this case, the integration time-step of solving the dynamic response of the numerical substructure can be dozens of times bigger than the sampling time-step of the controller. The time delay between the real and desired feedback forces becomes more striking, which challenges the well-developed delay compensation methods in real-time dynamic hybrid testing. This paper focuses on displacement prediction and force correction for delay compensation in the real-time dynamic hybrid testing with a large integration time-step. A new displacement prediction scheme is proposed based on recently-developed explicit integration algorithms and compared with several commonly-used prediction procedures. The evaluation of its prediction accuracy is carried out theoretically, numerically and experimentally. Results indicate that the accuracy and effectiveness of the proposed prediction method are of significance.

• Earthquakes and Structures
• /
• v.8 no.6
• /
• pp.1325-1348
• /
• 2015

Real-time hybrid testing (RTHT) involves virtual splitting of the structure into two parts: physical substructure that contains the key region of interest which is tested in a laboratory and numerical substructure that contains the remaining part of the structure in the form of a numerical model. This paper numerically assesses four step-by-step integration methods (Central difference method (CDM), Operator splitting method (OSM), Rosenbrock based method (RBM) and CR-integration method (CR)) which are widely used in RTHT. The methods have been assessed in terms of stability and accuracy for various realistic damping ratios of the physical substructure. The stability is assessed in terms of the spectral radii of the amplification matrix while the accuracy in terms of numerical damping and period distortion. In order to evaluate the performance of the methods, five carefully chosen examples have been studied - undamped SDOF, damped SDOF, instantaneous softening, instantaneous hardening and hysteretic system. The performance of the methods is measured in terms of a non-dimensional error index for displacement and velocity. Based on the error indices, it is observed that OSM and RBM are robust and performs fairly well in all the cases. CDM performed well for undamped SDOF system. CR method can be used for the system showing softening behaviour. The error indices indicate that accuracy of OSM is more than other method in case of hysteretic system. The accuracy of the results obtained through time integration methods for different damping ratios of the physical substructure is addressed in the present study. In the presence of a number of integration methods, it is preferable to have criteria for the selection of the time integration scheme. As such criteria are not available presently, this paper attempts to fill this gap by numerically assessing the four commonly used step-by-step methods.

• Journal of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.29 no.10
• /
• pp.89-96
• /
• 2015

Explicit numerical integration methods for power system transient stability simulation require very small time steps to avoid numerical instability. The EXST1 exciter model is a primary source of fast dynamics in power system transients. In case of the EXST1, the required small integration time step for entire system simulation increases the computational demands in terms of running time and storage. This paper presents a practical exciter model reduction approach which allows the increase of the required step size and thus the method can decrease the computational demands. The fast dynamics in the original EXST1 are eliminated in the reduced exciter model. The use of a larger time step improves the computational efficiency. This paper describes the way to eliminate the fast dynamics from the original exciter model based on linear system theory. In order to validate the performance of the proposed method, case studies with the GSO-37 bus system are provided. Comparisons between the original and reduced models are made in simulation accuracy and critical clearing time.

• 한국기계연구소 소보
• /
• s.18
• /
• pp.55-66
• /
• 1988

Many direct integration methods are used for numerical analyses of dynamic motion. In these methods, the governing equations of a dynamic system are integrated successively using a step-by-step numerical integration procedure. Time derivatives in the equations are generally approximated using difference formulas involving one or more increments of the time. Time increment has closely relationship with the accuracy of the motion analysis. In this paper, a 4th order Houbolt direct integration method is derived. For a spring-mass system, the motion of the system are analyzed from the 3rd order Houbolt and the 4th order Houbolt approaches respectively. Finally the paper proposes the optimal time-increment based on the accuracy of numerical analyses.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.4
• /
• pp.190-200
• /
• 1996

Step-by-step time integration methods are widely used for solving structural dynamics problem. One difficult yet critical choice an analyst must make is to decide an appropriate time step size. The choice of time step size has a significant effect on solution accuracy and computational expense. The objective of this research is to derive error estimate for newly developed time integration method and develop automatic time step size control algorithm for structural dynamics. A formula for computing error tolerance is derived based on desired period resolution. An automatic time step size control strategy is proposed based on a normalized local error estimate for the generalized-α method. Numerical examples demonstrate the developed strategy satisfies general design criteria for time step size control algorithm for dynamic problem.

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

• Smart Structures and Systems
• /
• v.14 no.6
• /
• pp.1131-1150
• /
• 2014

One of the issues in extending the range of applicable problems of real-time hybrid simulation is the computation speed of the simulator when large-scale computational models with a large number of DOF are used. In this study, functionality of real-time dynamic simulation of MDOF systems is achieved by creating a logic circuit that performs the step-by-step numerical time integration of the equations of motion of the system. The designed logic circuit can be implemented to an FPGA-based system; FPGA (Field Programmable Gate Array) allows large-scale parallel computing by implementing a number of arithmetic operators within the device. The operator splitting method is used as the numerical time integration scheme. The logic circuit consists of blocks of circuits that perform numerical arithmetic operations that appear in the integration scheme, including addition and multiplication of floating-point numbers, registers to store the intermediate data, and data busses connecting these elements to transmit various information including the floating-point numerical data among them. Case study on several types of linear and nonlinear MDOF system models shows that use of resource sharing in logic synthesis is crucial for effective application of FPGA to real-time dynamic simulation of structural response with time step interval of 1 ms.

• Earthquakes and Structures
• /
• v.7 no.2
• /
• pp.159-178
• /
• 2014

Although the family methods with unconditional stability and numerical dissipation have been developed for structural dynamics they all are implicit methods and thus an iterative procedure is generally involved for each time step. In this work, a new family method is proposed. It involves no nonlinear iterations in addition to unconditional stability and favorable numerical dissipation, which can be continuously controlled. In particular, it can have a zero damping ratio. The most important improvement of this family method is that it involves no nonlinear iterations for each time step and thus it can save many computationally efforts when compared to the currently available dissipative implicit integration methods.