TY - GEN
T1 - Optimal reduction of numerical dispersion for wave propagation problems. Application to isogeometric elements
AU - Idesman, A.
N1 - Funding Information:
The research has been supported in part by the Air Force Office of Scientific Research (contract FA9550-16-1-0177) and by Texas Tech University.
PY - 2017
Y1 - 2017
N2 - Based on the optimal coefficients of the stencil equation, a numerical technique for the reduction of the numerical dispersion error has been suggested. New isogeometric elements with the reduced numerical dispersion error for wave propagation problems have been developed with the suggested approach. By the minimization of the order of the dispersion error of the stencil equation, the order of the dispersion error is improved from order 2p (the conventional isogeometric elements) to order 4p (the isogeometric elements with reduced dispersion) where p is the order of the polynomial approximations. Because all coefficients of the stencil equation are obtained from the minimization procedure, the obtained accuracy is maximum possible. The corresponding elemental mass and stiffness matrices of the isogeometric elements with reduced dispersion are calculated with help of the optimal coefficients of the stencil equation. The analysis of the dispersion error of the isogeometric elements with the lumped mass matrix has also shown that independent of the procedures for the calculation of the lumped mass matrix, the second order of the dispersion error cannot be improved with the conventional stiffness matrix. However, the dispersion error with the lumped mass matrix can be improved from the second order to order 2p by the modification of the stiffness matrix. The numerical examples confirm the computational efficiency of the isogeometric elements with reduced dispersion that significantly reduce the computation time at a given accuracy. The numerical results obtained by the new and conventional isogeometric elements may include spurious oscillations due to the dispersion error. These oscillations can be quantified and filtered by the two-stage time-integration technique developed recently. The approach developed can be directly applied to other space-discretization techniques with similar stencil equations.
AB - Based on the optimal coefficients of the stencil equation, a numerical technique for the reduction of the numerical dispersion error has been suggested. New isogeometric elements with the reduced numerical dispersion error for wave propagation problems have been developed with the suggested approach. By the minimization of the order of the dispersion error of the stencil equation, the order of the dispersion error is improved from order 2p (the conventional isogeometric elements) to order 4p (the isogeometric elements with reduced dispersion) where p is the order of the polynomial approximations. Because all coefficients of the stencil equation are obtained from the minimization procedure, the obtained accuracy is maximum possible. The corresponding elemental mass and stiffness matrices of the isogeometric elements with reduced dispersion are calculated with help of the optimal coefficients of the stencil equation. The analysis of the dispersion error of the isogeometric elements with the lumped mass matrix has also shown that independent of the procedures for the calculation of the lumped mass matrix, the second order of the dispersion error cannot be improved with the conventional stiffness matrix. However, the dispersion error with the lumped mass matrix can be improved from the second order to order 2p by the modification of the stiffness matrix. The numerical examples confirm the computational efficiency of the isogeometric elements with reduced dispersion that significantly reduce the computation time at a given accuracy. The numerical results obtained by the new and conventional isogeometric elements may include spurious oscillations due to the dispersion error. These oscillations can be quantified and filtered by the two-stage time-integration technique developed recently. The approach developed can be directly applied to other space-discretization techniques with similar stencil equations.
KW - High-order Elements
KW - Numerical Dispersion
KW - Wave Propagation
UR - http://www.scopus.com/inward/record.url?scp=85042508343&partnerID=8YFLogxK
U2 - 10.7712/120117.5468.16785
DO - 10.7712/120117.5468.16785
M3 - Conference contribution
AN - SCOPUS:85042508343
T3 - COMPDYN 2017 - Proceedings of the 6th International Conference on Computational Methods in Structural Dynamics and Earthquake Engineering
SP - 946
EP - 954
BT - COMPDYN 2017 - Proceedings of the 6th International Conference on Computational Methods in Structural Dynamics and Earthquake Engineering
A2 - Papadrakakis, M.
A2 - Fragiadakis, Michalis
PB - National Technical University of Athens
T2 - 6th International Conference on Computational Methods in Structural Dynamics and Earthquake Engineering, COMPDYN 2017
Y2 - 15 June 2017 through 17 June 2017
ER -