scholarly journals A Weighted Residual Parabolic Acceleration Time Integration Method for Problems in Structural Dynamics

2007 ◽  
Vol 7 (3) ◽  
pp. 227-238 ◽  
Author(s):  
S.H. Razavi ◽  
A. Abolmaali ◽  
M. Ghassemieh
Keyword(s):  
Fourth Order ◽  
Initial Conditions ◽  
Time Integration ◽  
Critical Time ◽  
Linear Acceleration ◽  
Time Step ◽  
Time Integration Method ◽  
Acceleration Method ◽  
The Stability ◽  
Critical Time Step

AbstractIn the proposed method, the variation of displacement in each time step is assumed to be a fourth order polynomial in time and its five unknown coefficients are calculated based on: two initial conditions from the previous time step; satisfying the equation of motion at both ends of the time step; and the zero weighted residual within the time step. This method is non-dissipative and its dispersion is considerably less than in other popular methods. The stability of the method shows that the critical time step is more than twice of that for the linear acceleration method and its convergence is of fourth order.

Stability Analysis of Smoothed Finite Element Methods with Explicit Method for Transient Heat Transfer Problems

2019 ◽  
Vol 17 (02) ◽  
pp. 1845005 ◽  
Author(s):  
Xin Rong ◽  
Ruiping Niu ◽  
Guirong Liu
Keyword(s):  
Heat Transfer ◽  
Finite Element ◽  
Critical Time ◽  
Transient Heat Transfer ◽  
Explicit Method ◽  
Time Step ◽  
Smoothed Finite Element Method ◽  
The Stability ◽  
Transfer Problems ◽  
Critical Time Step

In this paper, transient heat transfer problems are analyzed using the smoothed finite element methods (S-FEMs) with explicit time integration. For a numerical method with spatial discretization, the computational cost per time step in the explicit method is less than that in the implicit method, but the time step is much smaller in the explicit analysis than that in the implicit analysis when the same mesh is used. This is because the stability is of essential importance. This work thus studies the stability of S-FEMs, when applied to transient heat transfer problems. Relationships are established between the critical time steps used in S-FEMs with the maximum eigenvalues of the thermal stiffness (conduction) matrix and mass matrix. It is found that the critical time step relates to the “softness” of the model. For example, node-based smoothed finite element method (NS-FEM) is softer than edge-based smoothed finite element method (ES-FEM), which leads to that the critical time step of NS-FEM is larger than that of ES-FEM. Because computing the eigenvalues and condition numbers of the stiffness matrices is very expensive but valuable for stability analysis, we proposed a concise and effective algorithm to estimate the maximum eigenvalue and condition number. Intensive numerical examples show that our scheme for computing the critical time step can work accurately and stably for the explicit method in FEM and S-FEMs.

scholarly journals Bipenalty method for time domain computational dynamics

2009 ◽  
Vol 466 (2117) ◽  
pp. 1389-1408 ◽  
Author(s):  
Harm Askes ◽  
Miguel Caramés-Saddler ◽  
Antonio Rodríguez-Ferran
Keyword(s):  
Time Integration ◽  
Critical Time ◽  
Time Step ◽  
Beam Elements ◽  
Computational Dynamics ◽  
Integration Schemes ◽  
Time Integration Schemes ◽  
The Common ◽  
Critical Time Step ◽  
Better Than

Penalty functions can be used to add constraints to systems of equations. In computational mechanics, stiffness-type penalties are the common choice. However, in dynamic applications stiffness penalties have the disadvantage that they tend to decrease the critical time step in conditionally stable time integration schemes, leading to increased CPU times for simulations. In contrast, inertia penalties increase the critical time step. In this paper, we suggest the simultaneous use of stiffness penalties and inertia penalties, which is denoted as the bipenalty method . We demonstrate that the accuracy of the bipenalty method is at least as good as (and usually better than) using either stiffness penalties or inertia penalties. Furthermore, for a number of finite elements (bar elements, beam elements and square plane stress/plane strain elements) we have derived ratios of the two penalty parameters such that their combined effect on the critical time step is neutral. The bipenalty method is thus superior to using stiffness penalties, because the decrease in critical time step can be avoided. The bipenalty method is also superior to using inertia penalties, since the constraints are realized with higher accuracy.

A Simple Explicit–Implicit Time-Marching Technique for Wave Propagation Analysis

2018 ◽  
Vol 16 (01) ◽  
pp. 1850082 ◽  
Author(s):  
Delfim Soares
Keyword(s):  
Wave Propagation ◽  
Time Integration ◽  
Critical Time ◽  
Implicit Time ◽  
Time Step ◽  
Stiffness Matrices ◽  
Time Integrators ◽  
Propagation Analysis ◽  
Time Marching ◽  
Critical Time Step

A new explicit–implicit time integration technique is proposed here for wave propagation analysis. In the present formulation, the time integrators of the model are selected at the element level, allowing each element to be considered as explicit or implicit. In the implicit elements, controllable algorithm dissipation is provided, enabling an [Formula: see text]-stable formulation. In the explicit elements, null amplitude decay is considered, enabling maximal critical time-step values. The new methodology renders a very simple and effective time-marching algorithm. Here, just displacement–velocity relations are considered, and no computation of accelerations is required. Moreover, explicit/implicit analyses can be carried out just by the tuning of local effective matrices, inputting or not stiffness matrices into their computations. At the end of the paper, numerical results are presented, illustrating the performance and potentialities of the new method.

A Precise and Stiffly Stable Time Integration Method for Vibration Equations

2001 ◽  
Author(s):  
Takeshi Fujikawa ◽  
Etsujiro Imanishi
Keyword(s):  
Integration Method ◽  
Time Integration ◽  
Low Frequency ◽  
Quadratic Function ◽  
Time Step ◽  
Integration Algorithm ◽  
Time Integration Method ◽  
Time Integration Algorithm ◽  
Time Integration Methods ◽  
Motion Problems

Abstract A method of time integration algorithm is presented for solving stiff vibration and motion problems. It is absolutely stable, numerically dissipative, and much accurate than other dissipative time integration methods. It achieves high-frequency dissipation, while minimizing unwanted low-frequency dissipation. In this method change of acceleration during time step is expressed as quadratic function including some parameters, whose appropriate values are determined through numerical investigation. Two calculation examples are demonstrated to show the usefulness of this method.

Parallel Computation of the Point Neutron Kinetic Equations Using Parallel Revisionist Integral Deferred Correction

2018 ◽  
Author(s):  
Yun Cai ◽  
Xingjie Peng ◽  
Qing Li ◽  
Zhizhu Zhang ◽  
Zhumin Jiang ◽  
...  
Keyword(s):  
Time Integration ◽  
Kinetic Equations ◽  
Accurate Method ◽  
Small Time ◽  
Time Step ◽  
Deferred Correction ◽  
Time Integration Method ◽  
Speed Up ◽  
Point Kinetics ◽  
Reactivity Insertion

The point kinetics is very important to the safety of the reactor operation. However, these equations are stiff and usually solved with very small time step. These equations are solved by Revisionist integral deferred correction (RIDC), which is a parallel time integration method. RIDC is a highly accurate method, and it reduces the error by iteration. Based on C++ and MPI, a four-core fourth-order RIDC is implemented and tested by several cases, such as step, ramp, and sinusoidal reactivity insertion. Compared with other methods, the time step of RIDC in the step reactivity insertion case is smaller, but it’s larger in the case of the sinusoidal reactivity insertion. RIDC can keep high accuracy while the time step is appropriately large. The numerical results also show that the speed-up ratio can achieve 2 when 4 processors are used.

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.

Implicit-Explicit Finite Elements in Transient Analysis: Stability Theory

1978 ◽  
Vol 45 (2) ◽  
pp. 371-374 ◽  
Author(s):  
T. J. R. Hughes ◽  
W. K. Liu
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Transient Analysis ◽  
Critical Time ◽  
Unconditional Stability ◽  
Time Step ◽  
Explicit Finite Element ◽  
New Family ◽  
Critical Time Step ◽  
Analysis Stability

A stability analysis is carried out for a new family of implicit-explicit finite-element algorithms. The analysis shows that unconditional stability may be achieved for the implicit finite elements and that the critical time step of the explicit elements governs for the system.

scholarly journals Mathematical modeling of temperature changes impact on artificial ice islands

2021 ◽  
Vol 13 (1) ◽  
pp. 79-86
Author(s):  
Maxim V. Muratov ◽  
◽  
Vladimir A. Biryukov ◽  
Denis S. Konov ◽  
Igor B. Petrov ◽  
...  
Keyword(s):  
Numerical Solution ◽  
Stefan Problem ◽  
Initial Conditions ◽  
Economic Feasibility ◽  
The Arctic ◽  
Seasonal Temperature ◽  
Time Step ◽  
Temperature Changes ◽  
The Stability ◽  
The Impact

The article is devoted to the numerical solution of the Stefan problem for thermal effects on an artificial ice island. For modern tasks of the development of the Arctic, associated with the exploration and production of minerals, it is important to create artificial ice islands in the Arctic shelf, due to the speed of their construction, economic feasibility and other factors. The most important task for the exploitation of such islands is their stability, including against melting. This paper discusses the issue of the stability of ice islands to melting. For this, the Stefan problem on the change in the phase state of matter is formulated. An enthalpy solution method is constructed, and the applicability of this method is considered. For the numerical solution, the Peasman-Reckford scheme is used, which is unconditionally spectrally stable in the two-dimensional case, which allows to freely choose the time step. In addition, the developed approach takes into account the flow of water and the flow of melted water, which is important in the task at hand. The developed computational algorithms are parallelized for use on modern multiprocessor computing systems An approach is implemented for modeling thermal processes in the thickness of an arbitrary mass of substances, taking into account arbitrary initial conditions, environmental conditions, tidal currents of water, and solar radiation. This approach was used to calculate the temperature distribution in the thickness of the ice island, as well as to study the impact of seasonal temperature changes on the stability of the island.

Sign in / Sign up

Export Citation Format

Share Document