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

• Bulletin of the Korean Space Science Society
• /
• 2008.10a
• /
• pp.26.4-27
• /
• 2008

The completion ('initiation' de facto) of the KASI Orbit Propagator and Estimator (KASIOPEA) has been delayed for several reasons unfortunately. Due to the lack of working staffs and the Division priority rearrangement, the initial plan was dismantled and ignored for many years. However, fundamental researches regarding the essential parts of KASIOPEA has been done by author. The numerical integration module of the KASIOPEA is the most sensitive part in the precision of the final output in general. There is no silver bullet in the numerical integration in an orbit propagation as a non-stiff ODE case. Many numerical integration method like single-step methods, multi-step method, and extrapolation methods have been used in overly populated orbit propagator or estimator. In this study, several popular methods from single-step, multi-step, and extrapolation methods have been tested in numerical accuracy and stability.

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

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

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.17-22
• /
• 2002

This paper describes numerical research on transverse jet behind rearward-facing step in turbulent supersonic flowfields without chemical reaction. The purpose of transverse jet behind rearward-facing step is to improve mixing of the fuel in the combustor. Two-dimensional unsteady flowfields generated by slot injection into supersonic flow are numerically simulated by integration of Navier-Stokes equation. Final-scale turbulence effects are modeled with two-equation $\kappa-\epsilon$ model. Numerical methods are modeled high-order upwind TVDschemes. A total of 4 cases are computed, comprising slot momentum flux ratios at four step heights downstream of the step. These numerical results are represented periodic phenomenon in unsteady flowfields.

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

• Structural Engineering and Mechanics
• /
• v.22 no.6
• /
• pp.701-717
• /
• 2006

In the existing reports regarding free transverse vibrations of the Euler-Bernoulli beams, most of them studied a uniform beam carrying various concentrated elements (such as point masses, rotary inertias, linear springs, rotational springs, spring-mass systems, ${\ldots}$, etc.) or a stepped beam with one to three step changes in cross-sections but without any attachments. The purpose of this paper is to utilize the numerical assembly method (NAM) to determine the exact natural frequencies and mode shapes of the multiple-step Euler-Bernoulli beams carrying a number of lumped masses and rotary inertias. First, the coefficient matrices for an intermediate lumped mass (and rotary inertia), left-end support and right-end support of a multiple-step beam are derived. Next, the overall coefficient matrix for the whole vibrating system is obtained using the numerical assembly technique of the conventional finite element method (FEM). Finally, the exact natural frequencies and the associated mode shapes of the vibrating system are determined by equating the determinant of the last overall coefficient matrix to zero and substituting the corresponding values of integration constants into the associated eigenfunctions, respectively. The effects of distribution of lumped masses and rotary inertias on the dynamic characteristics of the multiple-step beam are also studied.

• 한국기계연구소 소보
• /
• 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.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.453-458
• /
• 2004

This study proposes a new two-level condensation scheme for the construction of a reduced system. In the first step, the candidate area is selected for the construction of the reduced system by energy estimation in element-level. In the second step, primary degrees of freedom are selected by sequential elimination from the candidate degrees of freedom linked to the selected elements. Numerical examples demonstrate that the proposed method saves the computational cost effectively and provides a reduced system which predicts the eigenvalues accurately. Moreover, the well-constructed reduced system can present the reliable behavior of the structure under arbitrary dynamic loads comparing to that of global system. Time integration in a reduced system can save the computing time remarkably. Through a few numerical examples, the efficiency and reliability of the proposed scheme are verified.