An Analysis of Time Discretization in the Finite Element Solution of Hyperbolic Problems

THE PROBLEM OF THE TIME DISCRETIZATION OF HYPERBOLIC EQUATIONS WHEN FINITE ELEMENTS ARE USED TO REPRESENT THE SPATIAL DEPENDENCE IS CRITICALLY EXAMINED. A MODIFIED EQUATION ANALYSIS REVEALS THAT THE CLASSICAL, SECOND-ORDER ACCU- RATE, TIME-STEPPING ALGORITHMS, I.E., THE LAX-WENDROFF, LEAP-FROG, AND CRANK- NICOLSON METHODS, PROPERLY COMBINE WITH PIECEWISE LINEAR FINITE ELEMENTS IN ADVECTION PROBLEMS ONLY FOR SMALL VALUES OF THE TIME STEP. ON THE CONTRARY, AS THE COURANT NUMBER INCREASES, THE NUMERICAL PHASE ERROR DOES NOT DECREASE UNIFORMLY AT ALL WAVELENGTHS SO THAT THE OPTIMAL STABILITY LIMIT AND THE UNIT CFL PROPERTY ARE NOT ACHIEVED. THESE FUNDAMENTAL NUMERICAL PROPERTIES CAN HOWEVER BE RECOVERED, WHILE STILL REMAINING IN THE STANDARD GALERKIN FINITE ELEMENT SETTING, BY INCREASING THE ORDER OF ACCURACY OF THE TIME DISCRETI- ZATION. THIS IS ACCOMPLISHED BY EXPLOITING THE TAYLOR SERIES EXPANSION IN THE TIME INCREMENT UP TO THE THIRD ORDER BEFORE PERFORMING THE GALERKIN SPATIAL DISCRETIZATION USING PIECEWISE LINEAR LOCAL INTERPOLATIONS. AS A RESULT, IT APPEARS THAT THE PROPER FINITE ELEMENT EQUIVALENT OF SECOND-ORDER FINITE DIFFERENCE SCHEMES ARE IMPLICIT METHODS OF INCREMENTAL TYPE HAVING THIRD- AND FOURTH-ORDER GLOBAL ACCURACY (TAYLOR-GALERKIN METHODS). NUMERICAL RESULTS FOR SEVERAL LINEAR EXAMPLES ARE PRESENTED TO ILLUSTRATE THE PROPERTIES OF THE TAYLOR-GALERKIN SCHEMES IN ONE- AND TWO-DIMENSIONAL CALCULATIONS.

SELMIN V.; 

1995-03-15

JRC4206

https://publications.jrc.ec.europa.eu/repository/handle/JRC4206,