Title: An Analysis of Time Discretization in the Finite Element Solution of Hyperbolic Problems
Authors: SELMIN V.
Citation: Journal of Computational Physics vol. 70
Publication Year: 1987
JRC N°: JRC4206
URI: http://publications.jrc.ec.europa.eu/repository/handle/JRC4206
Type: Articles in Journals
Abstract: 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.
JRC Institute:Joint Research Centre Historical Collection

Files in This Item:
There are no files associated with this item.


Items in repository are protected by copyright, with all rights reserved, unless otherwise indicated.