An improved finite-element flux-corrected transport algorithm

An improved finite-element flux-corrected transport (FE-FCT) method for the numerical solution of hydrodynamic conservation equations is described, ba sed on the method developed by Lohner and his collaborators to solve conser vation equations in fluid mechanics, and its application is extended to gas discharge problems. The high- and low-order schemes used are the ones prop osed by Lohner who adds diffusion to the high-order scheme by subtracting t he lumped-mass matrix from the consistent-mass matrix to give the low-order scheme; the diffusion coefficient is adjusted globally. A variable diffusi on coefficient is introduced; it is assumed to be constant in each element and is shown to transform the high-order solution to a scheme equivalent to an upwind scheme which has minimal diffusion but ensures positive results. This avoids the complexity of upwinding in FE, especially in two dimension s. It is also shown that the correct amount of "real" diffusion may be easi ly added to the algorithm when required, for example, for electrons. Result s are presented which show that the high-order scheme reduces to the upwind difference scheme when the new diffusion is used. The proposed FCT scheme is shown to give similar results, in comparison with a fourth-order FD-FCT algorithm, Finally, the new method is applied to a streamer propagation pro blem in one dimension, and the results obtained are shown to agree well wit h previously published results. (C) 1999 Academic Press.