An unconditionally stable time integration method with controllable dissipation for second-order nonlinear dynamics

Author(s):  
Yi Ji ◽  
Yufeng Xing ◽  
Marian Wiercigroch
2021 ◽  
Vol 2090 (1) ◽  
pp. 012145
Author(s):  
Ryuma Honda ◽  
Hiroki Suzuki ◽  
Shinsuke Mochizuki

Abstract This study presents the impact of the difference between the implicit and explicit time integration methods on a steady turbulent flow field. In contrast to the explicit time integration method, the implicit time integration method may produce significant kinetic energy conservation error because the widely used spatial difference method for discretizing the governing equations is explicit with respect to time. In this study, the second-order Crank-Nicolson method is used as the implicit time integration method, and the fourth-order Runge-Kutta, second-order Runge-Kutta and second-order Adams-Bashforth methods are used as explicit time integration methods. In the present study, both isotropic and anisotropic steady turbulent fields are analyzed with two values of the Reynolds number. The turbulent kinetic energy in the steady turbulent field is hardly affected by the kinetic energy conservation error. The rms values of static pressure fluctuation are significantly sensitive to the kinetic energy conservation error. These results are examined by varying the time increment value. These results are also discussed by visualizing the large scale turbulent vortex structure.


1990 ◽  
Vol 57 (3) ◽  
pp. 553-561 ◽  
Author(s):  
I. Nishiguchi ◽  
T.-L. Sham ◽  
E. Krempl

A one-step time integration method is developed for the finite deformation theory of viscoplasticity based on overstress (FVBO) described in Part I. This time integration method is based on a forward gradient approximation and it leads to explicit expressions of the tangent operators suitable for finite element implementation. Numerical experiments and closed-form solutions for a hypoelastic material in homogeneous deformation states are presented. The FVBO is applied to the modeling of second-order effects in torsion. The numerical results show that a modification of the Jaumann rate and second-order terms of the inelastic rate of deformation are necessary to model the observed effects.


Author(s):  
Huimin Zhang ◽  
Runsen Zhang ◽  
Andrea Zanoni ◽  
Yufeng Xing ◽  
Pierangelo Masarati

AbstractA novel explicit three-sub-step time integration method is proposed. From linear analysis, it is designed to have at least second-order accuracy, tunable stability interval, tunable algorithmic dissipation and no overshooting behaviour. A distinctive feature is that the size of its stability interval can be adjusted to control the properties of the method. With the largest stability interval, the new method has better amplitude accuracy and smaller dispersion error for wave propagation problems, compared with some existing second-order explicit methods, and as the stability interval narrows, it shows improved period accuracy and stronger algorithmic dissipation. By selecting an appropriate stability interval, the proposed method can achieve properties better than or close to existing second-order methods, and by increasing or reducing the stability interval, it can be used with higher efficiency or stronger dissipation. The new method is applied to solve some illustrative wave propagation examples, and its numerical performance is compared with those of several widely used explicit methods.


2021 ◽  
Vol 11 (4) ◽  
pp. 1932
Author(s):  
Weixuan Wang ◽  
Qinyan Xing ◽  
Qinghao Yang

Based on the newly proposed generalized Galerkin weak form (GGW) method, a two-step time integration method with controllable numerical dissipation is presented. In the first sub-step, the GGW method is used, and in the second sub-step, a new parameter is introduced by using the idea of a trapezoidal integral. According to the numerical analysis, it can be concluded that this method is unconditionally stable and its numerical damping is controllable with the change in introduced parameters. Compared with the GGW method, this two-step scheme avoids the fast numerical dissipation in a low-frequency range. To highlight the performance of the proposed method, some numerical problems are presented and illustrated which show that this method possesses superior accuracy, stability and efficiency compared with conventional trapezoidal rule, the Wilson method, and the Bathe method. High accuracy in a low-frequency range and controllable numerical dissipation in a high-frequency range are both the merits of the method.


Sign in / Sign up

Export Citation Format

Share Document