UNSTEADY FLOW COMPUTATIONS OVER MOVING BODY USING DYNAMIC MESHES

2004 ◽  
Vol 01 (03) ◽  
pp. 507-518
Author(s):  
J. C. MANDAL ◽  
J. BALLMANN
Keyword(s):  
Compressible Flows ◽  
Time Integration ◽  
Time Discretization ◽  
Second Order ◽  
Integration Scheme ◽  
Arbitrary Lagrangian Eulerian ◽  
Time Step ◽  
Naca0012 Airfoil ◽  
Moving Body ◽  
Time Integration Scheme

An efficient implicit unstructured grid algorithm for solving unsteady inviscid compressible flows over moving body employing an Arbitrary Lagrangian Eulerian formulation is presented. In the present formulation, the time discretization is performed using a second-order accurate 3-point time integration scheme and the upwind-biased space discretization using second-order accurate finite volume formulation with Venkatakrishnan limiter. The face-velocities of the control volumes are computed using Geometric Conservation Laws. The nonlinear system arising from the implicit formulation is solved using an ILU preconditioned Newton–Krylov iteration at every time step. The computed results for two test cases involving harmonically oscillating NACA0012 airfoil are presented in order to demonstrate the efficacy of the present solver.

A New Unconditionally Stable Time Integration Method for Analysis of Nonlinear Structural Dynamics

2013 ◽  
Vol 80 (2) ◽  
Author(s):  
Ali Akbar Gholampour ◽  
Mehdi Ghassemieh ◽  
Mahdi Karimi-Rad
Keyword(s):  
Time Integration ◽  
Second Order ◽  
Numerical Dispersion ◽  
Computational Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Similar Time ◽  
Unconditionally Stable ◽  
Average Acceleration

A new time integration scheme is presented for solving the differential equation of motion with nonlinear stiffness. In this new implicit method, it is assumed that the acceleration varies quadratically within each time step. By increasing the order of acceleration, more terms of the Taylor series are used, which are expected to have responses with better accuracy than the classical methods. By considering this assumption and employing two parameters δ and α, a new family of unconditionally stable schemes is obtained. The order of accuracy, numerical dissipation, and numerical dispersion are used to measure the accuracy of the proposed method. Second order accuracy is achieved for all values of δ and α. The proposed method presents less dissipation at the lower modes in comparison with Newmark's average acceleration, Wilson-θ, and generalized-α methods. Moreover, this second order accurate method can control numerical damping in the higher modes. The numerical dispersion of the proposed method is compared with three unconditionally stable methods, namely, Newmark's average acceleration, Wilson-θ, and generalized-α methods. Furthermore, the overshooting effect of the proposed method is compared with these methods. By evaluating the computational time for analysis with similar time step duration, the proposed method is shown to be faster in comparison with the other methods.

scholarly journals Simulation of Wire Bonding Process Using Explicit Fem with Ale Remeshing Technology

2019 ◽  
Vol 36 (1) ◽  
pp. 47-54
Author(s):  
C. C. Yang ◽  
Y. F. Su ◽  
Steven Y. Liang ◽  
K. N. Chiang
Keyword(s):  
Time Integration ◽  
Wire Bonding ◽  
Bonding Process ◽  
Integration Scheme ◽  
Arbitrary Lagrangian Eulerian ◽  
Cracking Behavior ◽  
Free Air ◽  
Time Integration Scheme ◽  
Time Marching ◽  
The Wire

ABSTRACTThermosonic wire bonding is a common fabrication process for connecting devices in electronic packaging. However, when the free air ball (FAB) is compressed onto the I/O pad of the chip during bonding procedure, chip cracking may occur if the contact pressure is too large. This study proposes an effective simulation technique that can predict the wire ball geometry after bonding in an accurate range. The contact force obtained in the simulation can be used for possible die cracking behavior evaluation. The simulation in this study used the explicit time integration scheme to deal with the time marching problem, and the second-order precision arbitrary Lagrangian-Eulerian (ALE) algorithm was used to deal with the large deformation of the wire ball during the bonding process. In addition, the equilibrium smoothing algorithm in LS-DYNA can make the contact behavior and geometry of the bonding wire almost the same as the experiment, which can also significantly reduce the distortion of the mesh geometry after remeshing.

Newmark-Beta Method in Discrete Elastic Rods Algorithm to Avoid Energy Dissipation

2019 ◽  
Vol 86 (8) ◽  
Author(s):  
Weicheng Huang ◽  
Mohammad Khalid Jawed
Keyword(s):  
Energy Dissipation ◽  
Time Integration ◽  
Integration Scheme ◽  
Integration Step ◽  
Elastic Rods ◽  
Computationally Efficient ◽  
Step Size ◽  
Time Step ◽  
Time Step Size ◽  
Time Integration Scheme

Discrete elastic rods (DER) algorithm presents a computationally efficient means of simulating the geometrically nonlinear dynamics of elastic rods. However, it can suffer from artificial energy loss during the time integration step. Our approach extends the existing DER technique by using a different time integration scheme—we consider a second-order, implicit Newmark-beta method to avoid energy dissipation. This treatment shows better convergence with time step size, specially when the damping forces are negligible and the structure undergoes vibratory motion. Two demonstrations—a cantilever beam and a helical rod hanging under gravity—are used to show the effectiveness of the modified discrete elastic rods simulator.

scholarly journals Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation

2012 ◽  
Vol 2012 ◽  
pp. 1-12 ◽  
Author(s):  
Zhongdi Cen ◽  
Anbo Le ◽  
Aimin Xu
Keyword(s):  
Difference Scheme ◽  
Time Integration ◽  
Second Order ◽  
Integration Scheme ◽  
Uniform Mesh ◽  
Central Difference ◽  
Exponential Time ◽  
Black Scholes ◽  
Time Integration Scheme ◽  
Exponential Time Integration

We apply an exponential time integration scheme combined with a central difference scheme on a piecewise uniform mesh with respect to the spatial variable to evaluate a generalized Black-Scholes equation. We show that the scheme is second-order convergent for both time and spatial variables. It is proved that the scheme is unconditionally stable. Numerical results support the theoretical results.

Studies on the Accuracy Stability and Efficiency of a New Time Integration Scheme for Ballast Modeling Using the Discrete Element Method

2014 ◽  
Author(s):  
G. F. Mathews ◽  
R. L. Mullen ◽  
D. C. Rizos
Keyword(s):  
Particle Interaction ◽  
Time Integration ◽  
Discrete Element ◽  
Critical Discussion ◽  
Implicit Time ◽  
Integration Scheme ◽  
Computationally Efficient ◽  
Time Step ◽  
Central Difference ◽  
Time Integration Scheme

This paper presents the development of a semi-implicit time integration scheme, originally developed for structural dynamics in the 1970’s, and its implementation for use in Discrete Element Methods (DEM) for rigid particle interaction, and interaction of elastic bodies that are modeled as a cluster of rigid interconnected particles. The method is developed in view of ballast modeling that accounts for the flexibility of aggregates and the arbitrary shape and size of granules. The proposed scheme does not require any matrix inversions and is expressed in an incremental form making it appropriate for non-linear problems. The proposed method focuses on improving the efficiency, stability and accuracy of the solutions, as compared to current practice. A critical discussion of the findings of the studies is presented. Extended verification and assessment studies demonstrate that the proposed algorithm is unconditionally stable and accurate even for large time step sizes. It is demonstrated that the proposed method is at least as computationally efficient as the Central Difference Method. Guidelines for the implementation of the method to ballast modeling are discussed.

A COMPOSITE TIME INTEGRATION SCHEME FOR DYNAMIC ADHESION AND ITS APPLICATION TO GECKO SPATULA PEELING

2014 ◽  
Vol 11 (05) ◽  
pp. 1350104 ◽  
Author(s):  
SACHIN S. GAUTAM ◽  
ROGER A. SAUER
Keyword(s):  
Nonlinear Response ◽  
Time Integration ◽  
Computational Cost ◽  
Integration Scheme ◽  
Time Step ◽  
Integration Schemes ◽  
Newmark Scheme ◽  
Highly Nonlinear ◽  
Time Integration Schemes ◽  
Time Integration Scheme

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.

A Novel Sub-Stepping Method with Numerical Dissipation Control for Time Integration of Highly Flexible Structures

2018 ◽  
Vol 10 (10) ◽  
pp. 1850106 ◽  
Author(s):  
Saeed Mohammadzadeh ◽  
Mehdi Ghassemieh
Keyword(s):  
Equilibrium Points ◽  
Flexible Structures ◽  
Time Integration ◽  
Single Step ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Step Method ◽  
Time Step ◽  
Time Increment ◽  
Time Integration Scheme

Sub-stepping time integration methods attempt to march each time step with multiple sub-steps. Generally, for the first sub-step, a single-step method is applied and the following sub-steps are solved using a method that utilizes the data obtained from two or three previous equilibrium points. Despite the robust stability in problems, control of numerical dissipation in sub-stepping schemes is a tough task due to applying different algorithms on a time increment. In order to overcome this insufficiency, a new sub-stepping time integration scheme, which uses two sub-steps in each time increment, is proposed. Newmark and quadratic acceleration methods are applied on the first and second sub-steps, respectively. Both methods utilize constant parameters that enable the control of numerical dissipation in the analysis. For the proposed method, the stability analysis revealed the unconditional stability region for the pertinent parameters. Additionally, the precision investigation disclosed an advantage of the proposed method with the presence of minor period elongation error. Finally, the application of the proposed method is illuminated via several numerical examples. In addition to numerical dissipation control, the proposed method proved to have an outstanding advantage over other methods in solving highly flexible structures more efficiently and more accurately.

Numerical Method

2013 ◽  
Author(s):  
Gary A. Glatzmaier
Keyword(s):  
Finite Difference ◽  
Finite Difference Method ◽  
Numerical Method ◽  
Fluid Velocity ◽  
Time Integration ◽  
Spectral Space ◽  
Integration Scheme ◽  
Time Step ◽  
Time Integration Scheme ◽  
Solving Equations

This chapter describes a numerical method for solving equations of thermal convection on a computer. It begins by introducing the vorticity-streamfunction formulation as a means of conserving mass. The approach involves updating for the vorticity first and then solving for the fluid velocity each time step. The chapter continues with a discussion of two very different spatial discretizations, whereby the vertical derivatives are approximated with a finite-difference method and the horizontal derivatives with a spectral method. The nonlinear terms are computed in spectral space. The chapter also considers the Adams-Bashforth time integration scheme and explains how the Poisson equation can be solved at each time step for the updated streamfunction given the updated vorticity.

Sign in / Sign up

Export Citation Format

Share Document