TY - JOUR

T1 - Use of post-processing to increase the order of accuracy of the trapezoidal rule at time integration of linear elastodynamics problems

AU - Idesman, A.

N1 - Funding Information:
The research has been supported in part by the Air Force Research Lab, Eglin (contract # FA8651-08-D-0108) and by Texas Tech University.

PY - 2012/4/20

Y1 - 2012/4/20

N2 - In this paper, we suggest a very simple and effective post-processing procedure to increase the order of accuracy in time for numerical results obtained by the trapezoidal rule. We first derive a new exact, closed-form, a-priori error estimator for time integration of linear elastodynamics equations by the trapezoidal rule with non-uniform time increments. Based on this error estimator, we suggest a new post-processing procedure (containing additional time integration of elastodynamics equations by the trapezoidal rule with few time increments) that systematically improves the order of accuracy of numerical results, with the increase in the number of additional time increments used for post-processing. For example, the use of just one additional time increment for post-processing after time integration with any number of uniform time increments, renders the order of accuracy of numerical results equal to 10/3. Numerical examples of the application of the new techniques to a system with a single degree of freedom and to a multi-degree system confirm the corresponding increase in the order of convergence of numerical results after post-processing. Because the same trapezoidal rule is used for basic computations and post-processing, the new technique retains all of the properties of the trapezoidal rule, requires no writing of a new computer program for its implementation, and can be easily used with any existing commercial and research codes for elastodynamics.

AB - In this paper, we suggest a very simple and effective post-processing procedure to increase the order of accuracy in time for numerical results obtained by the trapezoidal rule. We first derive a new exact, closed-form, a-priori error estimator for time integration of linear elastodynamics equations by the trapezoidal rule with non-uniform time increments. Based on this error estimator, we suggest a new post-processing procedure (containing additional time integration of elastodynamics equations by the trapezoidal rule with few time increments) that systematically improves the order of accuracy of numerical results, with the increase in the number of additional time increments used for post-processing. For example, the use of just one additional time increment for post-processing after time integration with any number of uniform time increments, renders the order of accuracy of numerical results equal to 10/3. Numerical examples of the application of the new techniques to a system with a single degree of freedom and to a multi-degree system confirm the corresponding increase in the order of convergence of numerical results after post-processing. Because the same trapezoidal rule is used for basic computations and post-processing, the new technique retains all of the properties of the trapezoidal rule, requires no writing of a new computer program for its implementation, and can be easily used with any existing commercial and research codes for elastodynamics.

KW - Elastodynamics

KW - Error estimator

KW - Finite elements

KW - Time integration

KW - Trapezoidal rule

UR - http://www.scopus.com/inward/record.url?scp=84857238046&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2011.12.036

DO - 10.1016/j.jcp.2011.12.036

M3 - Article

AN - SCOPUS:84857238046

VL - 231

SP - 3143

EP - 3165

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 8

ER -