scholarly journals A novel explicit three-sub-step time integration method for wave propagation problems

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


2019 ◽  
Vol 26 (3-4) ◽  
pp. 161-174
Author(s):  
Taufeeq Ur Rehman Abbasi ◽  
Hui Zheng

Engineering systems for different levels of energy dissipation use internal variable models, which may lead to tremendous problems in accurate analysis. This article aims to provide an alternative direct integration method for the analysis of systems involving an anelastic displacement field model. A new state-space formulation built on an augmented set of anelastic variables for asymmetric systems is developed. Then, a precise time integration method based on state-space matrix formulation is proposed by introducing a Legendre–Gauss quadrature. The new integration method in terms of numerical stability and its implementation is discussed. The effect of sensitivity of the selection of the time-step and computational time on the performance of the new method is investigated by using a multi-degree-of-freedom system. The performance of the new method is also evaluated in terms of both computational accuracy and efficiency at higher degrees of freedom by using a continuum system. It is demonstrated that the computational accuracy and efficiency of the new method on large-scale problems are higher than that of the direct integration linear displacement–velocity method.


2013 ◽  
Vol 29 (3) ◽  
pp. 551-558
Author(s):  
J.-Y. Chang ◽  
C.-C. Tsai ◽  
T.-W. Hsu

AbstractIn this study, the fictitious time integration method (FTIM) is applied to investigate wave propagation over an arbitrary bathymetry with measured uncertainty. The FTIM is used to convert the higher-order elliptic mild-slope equation (EMSE) into a FTIM like EMSE (FTIMEMSE). It has the advantage to describe wave transformation from deep water to shallow water region in a large coastal area with numerical efficiency. The validity of the noise resistance for the measured uncertainty of the bathymetry is also studied. In addition, typical examples for waves propagating over an elliptic shoal rest on a horizontal and sloping bottom is presented. It is concluded that the FTIM is robust in the numerical stability and capable of against the noise of the measurement.


1999 ◽  
Author(s):  
Bertrand Tchamwa ◽  
Ted Conway ◽  
Christian Wielgosz

Abstract The purpose of this paper is to introduce a new simple explicit single step time integration method with controllable high-frequency dissipation. As opposed to the methods generally used in structural dynamics, with a consistency experimentally chosen of second order, the new method is only first-order-consistent but yields smaller numerical errors in low frequencies and is therefore very efficient for structural dynamic analysis. The new method remains explicit for any structural dynamics problem, even when a non-diagonal damping matrix is used in linear structural dynamics problem or when the non-linear internal force vector is a function of velocities. Convergence and spectral properties of the new algorithm are discussed and compared to those of some well-known algorithms. Furthermore, the validity and efficiency of the new algorithm are shown in a non-linear dynamic example by comparison of phase portraits.


Sign in / Sign up

Export Citation Format

Share Document