Dynamic Relaxation Using Continuous Kinetic Damping—Part I: Basic Algorithm

Author(s):  
Samuel Jung ◽  
Tae-Yun Kim ◽  
Wan-Suk Yoo

Dynamic relaxation (DR) is the most widely used approach for static equilibrium analyses. Specifically, DR compels dynamic systems to converge to a static equilibrium through the addition of fictitious damping. DR methods are classified by the method in which fictitious damping is applied. Conventional DR methods use a fictitious mass matrix to increase the fictitious damping while maintaining numerical stability. There are many calculation methods for the fictitious mass matrix; however, it is difficult to select the appropriate method. In addition, these methods require a stiffness matrix of a model, which makes it difficult to apply nonlinear models. To resolve these problems, a new DR method that uses continuous kinetic damping (CKDR) is proposed in this study. The proposed method does not require the fictitious mass matrix and any tuning coefficients, and it possesses a second-order convergence rate. The aforementioned advantages are unique and significant when compared to those of conventional methods. The stability and convergence rate were analyzed by using an eigenvalue analysis and demonstrated by simulating nonlinear models of a pendulum and cable. Simple but representative models were used to clearly demonstrate the features of the proposed DR method and to enable the reproducibility of the verification results.

2018 ◽  
Vol 2018 ◽  
pp. 1-9 ◽  
Author(s):  
Samuel Jung ◽  
Tae-Yun Kim ◽  
Wan-Suk Yoo

Dynamic relaxation (DR) is a widely used numerical method to determine the static equilibrium of a dynamic system. However, it is difficult to apply conventional DR methods to nonlinear models because they require estimation of a stiffness matrix of the model. To resolve the forementioned problem, a new dynamic relaxation method using continuous kinetic damping (CKDR) was proposed in previous research. The CKDR method does not require any model parameters including the stiffness matrix, and it possesses absolute stability and a second-order convergence rate. However, the convergence rate is proportional to square of the step size, and this may result in a low convergence rate if the selected step size is excessively small. This problem leads to difficulties in the practical use of CKDR. Thus, an adaptive step-size method is proposed in this paper to control the convergence rate of CKDR. The proposed method estimates natural frequency of the model and determines adaptive step size. Static equilibrium simulations were performed for three different models to verify the method. The results revealed that the computational cost of CKDR with a variable step size was very efficient when compared to fixed step sizes and that the convergence rate was also controlled as intended. In addition, the lowest natural frequencies of models in static equilibrium were accurately estimated.


2010 ◽  
Vol 145 ◽  
pp. 128-133
Author(s):  
Chun Jiang Zhao ◽  
Lian Yun Jiang ◽  
Jin Zhi Zhang ◽  
Qing Xue Huang ◽  
Xiao Kai Yu

Based on the theory of mathematical analysis, I find the rolling disturbance can be measured. Then the mathematical model of dynamic setting AGC is gotten by recursive methods. By the mathematical model I find out the influence of model parameters on the stability and convergence rate of the control system. When the system is stable, an influence of model parameters and parameter of the control system on steel strap thickness have been obtained, which will be helpful for us to choose suitable parameters in the end.


2019 ◽  
Vol 1 (1) ◽  
pp. 49-60
Author(s):  
Simon Heru Prassetyo ◽  
Ganda Marihot Simangunsong ◽  
Ridho Kresna Wattimena ◽  
Made Astawa Rai ◽  
Irwandy Arif ◽  
...  

This paper focuses on the stability analysis of the Nanjung Water Diversion Twin Tunnels using convergence measurement. The Nanjung Tunnel is horseshoe-shaped in cross-section, 10.2 m x 9.2 m in dimension, and 230 m in length. The location of the tunnel is in Curug Jompong, Margaasih Subdistrict, Bandung. Convergence monitoring was done for 144 days between February 18 and July 11, 2019. The results of the convergence measurement were recorded and plotted into the curves of convergence vs. day and convergence vs. distance from tunnel face. From these plots, the continuity of the convergence and the convergence rate in the tunnel roof and wall were then analyzed. The convergence rates from each tunnel were also compared to empirical values to determine the level of tunnel stability. In general, the trend of convergence rate shows that the Nanjung Tunnel is stable without any indication of instability. Although there was a spike in the convergence rate at several STA in the measured span, that spike was not replicated by the convergence rate in the other measured spans and it was not continuous. The stability of the Nanjung Tunnel is also confirmed from the critical strain analysis, in which most of the STA measured have strain magnitudes located below the critical strain line and are less than 1%.


2021 ◽  
Vol 3 (2) ◽  
Author(s):  
Abdullah Al-Mamun ◽  
S. M. Arifuzzaman ◽  
Sk. Reza-E-Rabbi ◽  
Umme Sara Alam ◽  
Saiful Islam ◽  
...  

AbstractThe perspective of this paper is to characterize a Casson type of Non-Newtonian fluid flow through heat as well as mass conduction towards a stretching surface with thermophoresis and radiation absorption impacts in association with periodic hydromagnetic effect. Here heat absorption is also integrated with the heat absorbing parameter. A time dependent fundamental set of equations, i.e. momentum, energy and concentration have been established to discuss the fluid flow system. Explicit finite difference technique is occupied here by executing a procedure in Compaq Visual Fortran 6.6a to elucidate the mathematical model of liquid flow. The stability and convergence inspection has been accomplished. It has observed that the present work converged at, Pr ≥ 0.447 indicates the value of Prandtl number and Le ≥ 0.163 indicates the value of Lewis number. Impact of useful physical parameters has been illustrated graphically on various flow fields. It has inspected that the periodic magnetic field has helped to increase the interaction of the nanoparticles in the velocity field significantly. The field has been depicted in a vibrating form which is also done newly in this work. Subsequently, the Lorentz force has also represented a great impact in the updated visualization (streamlines and isotherms) of the flow field. The respective fields appeared with more wave for the larger values of magnetic parameter. These results help to visualize a theoretical idea of the effect of modern electromagnetic induction use in industry instead of traditional energy sources. Moreover, it has a great application in lung and prostate cancer therapy.


2009 ◽  
Vol 131 (3) ◽  
Author(s):  
Mario Durán ◽  
Jean-Claude Nédélec ◽  
Sebastián Ossandón

An efficient numerical method, using integral equations, is developed to calculate precisely the acoustic eigenfrequencies and their associated eigenvectors, located in a given high frequency interval. It is currently known that the real symmetric matrices are well adapted to numerical treatment. However, we show that this is not the case when using integral representations to determine with high accuracy the spectrum of elliptic, and other related operators. Functions are evaluated only in the boundary of the domain, so very fine discretizations may be chosen to obtain high eigenfrequencies. We discuss the stability and convergence of the proposed method. Finally we show some examples.


Author(s):  
David Roos Launchbury ◽  
Luca Mangani ◽  
Ernesto Casartelli ◽  
Francesco Del Citto

Abstract In the industrial simulation of flow phenomena, turbulence modeling is of prime importance. Due to their low computational cost, Reynolds-averaged methods (RANS) are predominantly used for this purpose. However, eddy viscosity RANS models are often unable to adequately capture important flow physics, specifically when strongly anisotropic turbulence and vortex structures are present. In such cases the more costly 7-equation Reynolds stress models often lead to significantly better results. Unfortunately, these models are not widely used in the industry. The reason for this is not mainly the increased computational cost, but the stability and convergence issues such models usually exhibit. In this paper we present a robust implementation of a Reynolds stress model that is solved in a coupled manner, increasing stability and convergence speed significantly compared to segregated implementations. In addition, the decoupling of the velocity and Reynolds stress fields is addressed for the coupled equation formulation. A special wall function is presented that conserves the anisotropic properties of the model near the walls on coarser meshes. The presented Reynolds stress model is validated on a series of semi-academic test cases and then applied to two industrially relevant situations, namely the tip vortex of a NACA0012 profile and the Aachen Radiver radial compressor case.


Sign in / Sign up

Export Citation Format

Share Document