Abstract
We present three time discretization schemes for the Green element solution of the linear conduction equation. Numerical results obtained from the three methods are assessed by their convergence and stability-related properties, namely, numerical amplitude and amplification factor through Fourier analysis. The main emphasis is the ability of any of the schemes to handle conduction problems typified by those properties that are known to be taxing to most numerical schemes. This was found in some cases to be related to how accurately the harmonics of the Fourier waves are propagated.