scholarly journals Damping Acoustic Modes in Compressible Horizontally Explicit Vertically Implicit (HEVI) and Split-Explicit Time Integration Schemes

2018 ◽  
Vol 146 (6) ◽  
pp. 1911-1923 ◽  
Author(s):  
Joseph B. Klemp ◽  
William C. Skamarock ◽  
Soyoung Ha

Although the equations of motion for a compressible atmosphere accommodate acoustic waves, these modes typically play an insignificant role in atmospheric processes of physical interest. In numerically integrating the compressible equations, it is often beneficial to filter these acoustic modes to control acoustic noise and prevent its artificial growth. Here, a new technique is proposed for filtering the 3D divergence that may damp acoustic modes more effectively than filters previously implemented in numerical modes using horizontally explicit vertically implicit (HEVI) and split-explicit time integration schemes. With this approach, a divergence damping term is added as a final adjustment to the horizontal velocity at the new time level after completing the vertically implicit portion of the time step. In this manner, the divergence used in the filter term has exactly the same numerical form as that used in the discrete pressure equation. Analysis of the dispersion equation for this form of the filter documents its stability characteristics and confirms that it effectively damps acoustic modes with little artificial influence on the amplitude or propagation of the gravity wave modes that are of physical interest. Some specific aspects of the implementation of the filter in the Model for Prediction Across Scales (MPAS) are discussed, and results are presented to illustrate some of the beneficial aspects of suppressing acoustic noise.

1986 ◽  
Vol 65 (2) ◽  
pp. 253-272 ◽  
Author(s):  
L. Garcia ◽  
H.R. Hicks ◽  
B.A. Carreras ◽  
L.A. Charlton ◽  
J.A. Holmes

1981 ◽  
Vol 48 (4) ◽  
pp. 909-914 ◽  
Author(s):  
A. Libai

The field equations of motion and compatibility for the nonlinear dynamics of doubly curved shells are recast in an intrinsic form, in terms of the metric and curvature functions of their reference surfaces. For appropriate input, the motion of the shell is described without the need for an external reference coordinate system or the use of vector quantities such as position, velocity, and acceleration. The equations are shown to be readily applicable to time integration schemes. Such cases, as the (spatially) constant load problem and the inextensional dynamics problem, are also considered. The need for further work in the area is emphasized.


Author(s):  
J. M. Rodriguez ◽  
S. Larsson ◽  
J. M. Carbonell ◽  
P. Jonsén

AbstractThis work presents the development of an explicit/implicit particle finite element method (PFEM) for the 2D modeling of metal cutting processes. The purpose is to study the efficiency of implicit and explicit time integration schemes in terms of precision, accuracy and computing time. The formulation for implicit and explicit time marching schemes is developed, and a detailed study on the explicit solution steps is presented. The PFEM remeshing procedures for insertion and removal of particles have been improved to model the multiple scales of time and/or space of the solution. The detection and treatment of the rigid tool contact are presented for both, implicit and explicit schemes. The performance of explicit/implicit integration is studied with a set of different two-dimensional orthogonal cutting tests of AISI 4340 steel at cutting speeds ranging from 1 m/s up to 30 m/s. It was shown that if the correct selection of the time integration scheme is made, the computing time can decrease up to 40 times. It allows us to affirm that the computing time of the PFEM simulations can be excessive due to the used time marching scheme independently of the meshing process. As a practical result, a set of recommendations to select the time integration schemes for a given cutting speed are given. This is intended to minimize one of the negative constraints pointed out by the industry when using metal cutting simulators.


2019 ◽  
Vol 54 (2) ◽  
pp. 116-129 ◽  
Author(s):  
Roberto Ortega ◽  
Geraldine Farías ◽  
Marcela Cruchaga ◽  
Matías Rivero ◽  
Mariano Vázquez ◽  
...  

The focus of this work is on the computational modeling of a pendulum made of a hyperelastic material and the corresponding experimental validation with the aim of contributing to the study of a material commonly used in seismic absorber devices. From the proposed dynamics experiment, the motion of the pendulum is recorded using a high-speed camera. The evolution of the pendulum’s positions is recovered using a capturing motion technique by tracking markers. The simulation of the problem is developed in the framework of a parallel multi-physics code. Particular emphasis is placed on the analysis of the Newmark integration scheme and the use of Rayleigh damping model. In particular, the time step size effect is analyzed. A strong time step size dependency is obtained for dissipative time integration schemes, while the Rayleigh damping formulation without time integration dissipation shows time step–independent results when convergence is achieved.


Author(s):  
Olivier A. Bauchau ◽  
Alexander Epple ◽  
Carlo L. Bottasso

This paper addresses practical issues associated with the numerical enforcement of constraints in flexible multibody systems, which are characterized by index-3 differential algebraic equations (DAEs). The need to scale the equations of motion is emphasized; in the proposed approach, they are scaled based on simple physical arguments, and an augmented Lagrangian term is added to the formulation. Time discretization followed by a linearization of the resulting equations leads to a Jacobian matrix that is independent of the time step size, h; hence, the condition number of the Jacobian and error propagation are both O(h0): the numerical solution of index-3 DAEs behaves as in the case of regular ordinary differential equations (ODEs). Since the scaling factor depends on the physical properties of the system, the proposed scaling decreases the dependency of this Jacobian on physical properties, further improving the numerical conditioning of the resulting linearized equations. Because the scaling of the equations is performed before the time and space discretizations, its benefits are reaped for all time integration schemes. The augmented Lagrangian term is shown to be indispensable if the solution of the linearized system of equations is to be performed without pivoting, a requirement for the efficient solution of the sparse system of linear equations. Finally, a number of numerical examples demonstrate the efficiency of the proposed approach to scaling.


2014 ◽  
Vol 11 (05) ◽  
pp. 1350104 ◽  
Author(s):  
SACHIN S. GAUTAM ◽  
ROGER A. SAUER

Simulation of dynamic adhesive peeling problems at small scales has attracted little attention so far. These problems are characterized by a highly nonlinear response. Accurate and stable time integration schemes are required for simulation of dynamic peeling problems. In the present work, a composite time integration scheme is proposed for the simulation of dynamic adhesive peeling problems. It is shown through numerical examples that the proposed scheme remains stable and also has some gain in accuracy. The performance of the scheme is compared with two collocation-based schemes, i.e., Newmark scheme and Bathe composite scheme. It is shown that the proposed scheme and Bathe composite scheme perform equally. However, the proposed scheme adds very little to the computational cost of Newmark scheme. Through a numerical simulation of the peeling of a gecko spatula from a rigid substrate it is shown that the proposed scheme and the Bathe composite scheme are able to simulate the complete peeling process for given time step whereas the Newmark scheme diverges. It is also shown that the maximum pull-off force is within the range reported in the literature.


Sign in / Sign up

Export Citation Format

Share Document