scholarly journals The Selection of Time Step in Runge Kutta Fourth Order for Determine Deviation in the Weapon Arm Vehicle

2015 ◽  
Vol 68 ◽  
pp. 363-369 ◽  
Author(s):  
Aditya Sukma Nugraha
Keyword(s):  
Fourth Order ◽  
Time Step ◽  
Runge Kutta ◽  
Selection Of

Long-time symplectic integration: the example of four-vortex motion

1991 ◽  
Vol 432 (1886) ◽  
pp. 481-494 ◽  
Keyword(s):  
Fourth Order ◽  
Initial Conditions ◽  
Hamiltonian Formulation ◽  
Superior Performance ◽  
Integration Time ◽  
Characteristic Time Scale ◽  
Symplectic Integration ◽  
Time Step ◽  
Long Time ◽  
Runge Kutta

Detailed comparisons are made between long-time numerical integration of the motion of four identical point vortices obtained using both a fourth-order symplectic integration method of the implicit Runge-Kutta type and a standard fourth-order explicit Runge-Kutta scheme. We utilize the reduced hamiltonian formulation of the four-vortex problem due to Aref & Pomphrey. Initial conditions which give both fully chaotic and also quasi-periodic motions are considered over integration times of order 10 6 -10 7 times the characteristic time scale of the evolution. The convergence, as the integration time step is decreased, of the Poincaré section is investigated. When smoothness of the section compared to the converged image, and the fractional change in the hamiltonian H are used as diagnostic indicators, it is found that the symplectic scheme gives substantially superior performance over the explicit scheme, and exhibits only an apparent qualitative degrading in results up to integration time steps of order the minimum timescale of the evolution. It is concluded that this performance derives from the symplectic rather than the implicit character of the method.

On the accuracy of some explicit and implicit methods for the inviscid GRLW equation subject to initial Gaussian conditions

2016 ◽  
Vol 26 (3/4) ◽  
pp. 698-721 ◽  
Author(s):  
J I Ramos
Keyword(s):  
Finite Difference ◽  
Fourth Order ◽  
Finite Difference Methods ◽  
Second Order ◽  
Time Step ◽  
Content Type ◽  
Difference Methods ◽  
Temporal And Spatial ◽  
Runge Kutta ◽  
Grlw Equation

Purpose – The purpose of this paper is to both determine the effects of the nonlinearity on the wave dynamics and assess the temporal and spatial accuracy of five finite difference methods for the solution of the inviscid generalized regularized long-wave (GRLW) equation subject to initial Gaussian conditions. Design/methodology/approach – Two implicit second- and fourth-order accurate finite difference methods and three Runge-Kutta procedures are introduced. The methods employ a new dependent variable which contains the wave amplitude and its second-order spatial derivative. Numerical experiments are reported for several temporal and spatial step sizes in order to assess their accuracy and the preservation of the first two invariants of the inviscid GRLW equation as functions of the spatial and temporal orders of accuracy, and thus determine the conditions under which grid-independent results are obtained. Findings – It has been found that the steepening of the wave increase as the nonlinearity exponent is increased and that the accuracy of the fourth-order Runge-Kutta method is comparable to that of a second-order implicit procedure for time steps smaller than 100th, and that only the fourth-order compact method is almost grid-independent if the time step is on the order of 1,000th and more than 5,000 grid points are used, because of the initial steepening of the initial profile, wave breakup and solitary wave propagation. Originality/value – This is the first study where an accuracy assessment of wave breakup of the inviscid GRLW equation subject to initial Gaussian conditions is reported.

Conservative Transport Schemes for Spherical Geodesic Grids: High-Order Flux Operators for ODE-Based Time Integration

2011 ◽  
Vol 139 (9) ◽  
pp. 2962-2975 ◽  
Author(s):  
William C. Skamarock ◽  
Almut Gassmann
Keyword(s):  
Fourth Order ◽  
Time Integration ◽  
Higher Order ◽  
Second Order ◽  
Time Step ◽  
Diffusion Term ◽  
The Third ◽  
Integration Schemes ◽  
Order Scheme ◽  
Runge Kutta

Higher-order finite-volume flux operators for transport algorithms used within Runge–Kutta time integration schemes on irregular Voronoi (hexagonal) meshes are proposed and tested. These operators are generalizations of third- and fourth-order operators currently used in atmospheric models employing regular, orthogonal rectangular meshes. Two-dimensional least squares fit polynomials are used to evaluate the higher-order spatial derivatives needed to cancel the leading-order truncation error terms of the standard second-order centered formulation. Positive definite or monotonic behavior is achieved by applying an appropriate limiter during the final Runge–Kutta stage within a given time step. The third- and fourth-order formulations are evaluated using standard transport tests on the sphere. The new schemes are more accurate and significantly more efficient than the standard second-order scheme and other schemes in the literature examined by the authors. The third-order formulation is equivalent to the fourth-order formulation plus an additional diffusion term that is proportional to the Courant number. An optimal value for the coefficient scaling this diffusion term is chosen based on qualitative evaluation of the test results. Improvements using the higher-order scheme in place of the traditional second-order centered approach are illustrated within 3D unstable baroclinic wave simulations produced using two global nonhydrostatic models employing spherical Voronoi meshes.

A Two-Stage Fourth-Order Multimoment Global Shallow-Water Model on the Cubed Sphere

2020 ◽  
Vol 148 (10) ◽  
pp. 4267-4279
Author(s):  
Yuzhang Che ◽  
Chungang Chen ◽  
Feng Xiao ◽  
Xingliang Li ◽  
Xueshun Shen
Keyword(s):  
Shallow Water ◽  
Fourth Order ◽  
Time Integration ◽  
Water Model ◽  
Shallow Water Model ◽  
Time Step ◽  
Two Stage ◽  
Order Scheme ◽  
Runge Kutta ◽  
Cubed Sphere

AbstractA new multimoment global shallow-water model on the cubed sphere is proposed by adopting a two-stage fourth-order Runge–Kutta time integration. Through calculating the values of predicted variables at half time step t = tn + (1/2)Δt by a second-order formulation, a fourth-order scheme can be derived using only two stages within one time step. This time integration method is implemented in our multimoment global shallow-water model to build and validate a new and more efficient numerical integration framework for dynamical cores. As the key task, the numerical formulation for evaluating the derivatives in time has been developed through the Cauchy–Kowalewski procedure and the spatial discretization of the multimoment finite-volume method, which ensures fourth-order accuracy in both time and space. Several major benchmark tests are used to verify the proposed numerical framework in comparison with the existing four-stage fourth-order Runge–Kutta method, which is based on the method of lines framework. The two-stage fourth-order scheme saves about 30% of the computational cost in comparison with the four-stage Runge–Kutta scheme for global advection and shallow-water models. The proposed two-stage fourth-order framework offers a new option to develop high-performance time marching strategy of practical significance in dynamical cores for atmospheric and oceanic models.

Is Hydraulics Over-Used or Under-Used in Storm Drainage?

1994 ◽  
Vol 29 (1-2) ◽  
pp. 53-61
Author(s):  
Ben Chie Yen
Keyword(s):  
Unsteady Flow ◽  
Computational Time ◽  
Time Step ◽  
Flow Routing ◽  
Kinematic Wave Equation ◽  
Storm Drainage ◽  
Saint Venant Equations ◽  
Unsteady Flow Routing ◽  
Selection Of

Urban drainage models utilize hydraulics of different levels. Developing or selecting a model appropriate to a particular project is not an easy task. Not knowing the hydraulic principles and numerical techniques used in an existing model, users often misuse and abuse the model. Hydraulically, the use of the Saint-Venant equations is not always necessary. In many cases the kinematic wave equation is inadequate because of the backwater effect, whereas in designing sewers, often Manning's formula is adequate. The flow travel time provides a guide in selecting the computational time step At, which in turn, together with flow unsteadiness, helps in the selection of steady or unsteady flow routing. Often the noninertia model is the appropriate model for unsteady flow routing, whereas delivery curves are very useful for stepwise steady nonuniform flow routing and for determination of channel capacity.

scholarly journals Krylov SSP Integrating Factor Runge–Kutta WENO Methods

Mathematics ◽  
2021 ◽  
Vol 9 (13) ◽  
pp. 1483
Author(s):  
Shanqin Chen
Keyword(s):  
Large Time ◽  
Hyperbolic Equations ◽  
Strong Stability ◽  
Matrix Exponential ◽  
Linear Component ◽  
Nonlinear Hyperbolic Equations ◽  
Time Step ◽  
Integrating Factor ◽  
The Matrix ◽  
Runge Kutta

Weighted essentially non-oscillatory (WENO) methods are especially efficient for numerically solving nonlinear hyperbolic equations. In order to achieve strong stability and large time-steps, strong stability preserving (SSP) integrating factor (IF) methods were designed in the literature, but the methods there were only for one-dimensional (1D) problems that have a stiff linear component and a non-stiff nonlinear component. In this paper, we extend WENO methods with large time-stepping SSP integrating factor Runge–Kutta time discretization to solve general nonlinear two-dimensional (2D) problems by a splitting method. How to evaluate the matrix exponential operator efficiently is a tremendous challenge when we apply IF temporal discretization for PDEs on high spatial dimensions. In this work, the matrix exponential computation is approximated through the Krylov subspace projection method. Numerical examples are shown to demonstrate the accuracy and large time-step size of the present method.

Application of the piecewise constant method in nonlinear dynamics of drill string

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Jialin Tian ◽  
Jie Wang ◽  
Yi Zhou ◽  
Lin Yang ◽  
Changyue Fan ◽  
...  
Keyword(s):  
Nonlinear Dynamics ◽  
Dynamic Model ◽  
Dynamic Characteristics ◽  
Drill String ◽  
Vibration Velocity ◽  
Kutta Method ◽  
Time Step ◽  
Piecewise Constant ◽  
Runge Kutta Method ◽  
Runge Kutta

Abstract Aiming at the current development of drilling technology and the deepening of oil and gas exploration, we focus on better studying the nonlinear dynamic characteristics of the drill string under complex working conditions and knowing the real movement of the drill string during drilling. This paper firstly combines the actual situation of the well to establish the dynamic model of the horizontal drill string, and analyzes the dynamic characteristics, giving the expression of the force of each part of the model. Secondly, it introduces the piecewise constant method (simply known as PT method), and gives the solution equation. Then according to the basic parameters, the axial vibration displacement and vibration velocity at the test points are solved by the PT method and the Runge–Kutta method, respectively, and the phase diagram, the Poincare map, and the spectrogram are obtained. The results obtained by the two methods are compared and analyzed. Finally, the relevant experimental tests are carried out. It shows that the results of the dynamic model of the horizontal drill string are basically consistent with the results obtained by the actual test, which verifies the validity of the dynamic model and the correctness of the calculated results. When solving the drill string nonlinear dynamics, the results of the PT method is closer to the theoretical solution than that of the Runge–Kutta method with the same order and time step. And the PT method is better than the Runge–Kutta method with the same order in smoothness and continuity in solving the drill string nonlinear dynamics.

scholarly journals Two-Way Contact Network Modeling for Identifying the Route of COVID-19 Community Transmission

Informatics ◽  
2021 ◽  
Vol 8 (2) ◽  
pp. 22
Author(s):  
Sung Jin Lee ◽  
Sang Eun Lee ◽  
Ji-On Kim ◽  
Gi Bum Kim
Keyword(s):  
Simulated Data ◽  
Contact Network ◽  
Transmission Route ◽  
Time Step ◽  
Source Of Infection ◽  
Health Authorities ◽  
Control And Prevention ◽  
Community Transmission ◽  
Centrality Analysis ◽  
Selection Of

In this study, we address the problem originated from the fact that “The Corona 19 Epidemiological Research Support System,” developed by the Korea Centers for Disease Control and Prevention, is limited to analyzing the Global Positioning System (GPS) information of the confirmed COVID-19 cases alone. Consequently, we study a method that the authority predicts the transmission route of COVID-19 between visitors in the community from a spatiotemporal perspective. This method models a contact network around the first confirmed case, allowing the health authorities to conduct tests on visitors after an outbreak of COVID-19 in the community. After securing the GPS data of community visitors, it traces back to the past from the time the first confirmed case occurred and creates contact clusters at each time step. This is different from other researches that focus on identifying the movement paths of confirmed patients by forward tracing. The proposed method creates the contact network by assigning weights to each contact cluster based on the degree of proximity between contacts. Identifying the source of infection in the contact network can make us predict the transmission route between the first confirmed case and the source of infection and classify the contacts on the transmission route. In this experiment, we used 64,073 simulated data for 100 people and extracted the transmission route and a top 10 list for centrality analysis. The contacts on the route path can be quickly designated as a priority for COVID-19 testing. In addition, it is possible for the authority to extract the subjects with high influence from the centrality theory and use them for additional COVID-19 epidemic investigation that requires urgency. This model is expected to be used in the epidemic investigation requiring the quick selection of close contacts.

Sign in / Sign up

Export Citation Format

Share Document