This paper presents a characteristic Galerkin finite element method with an
implicit algorithm for solving multidimensional, time-dependent convection
-diffusion equations. The method is formulated on the basis of the combinat
ion of both the precise and the implicit numerical integration procedures a
iming to reference particles. The precise integration procedure with a 2(N)
algorithm is taken as a tool to determine the material (Lagrangian) deriva
tive of the convective function in the operator splitting procedure. The st
ability analysis of the algorithm and numerical results illustrate good per
formance of the present method in stability and accuracy. Copyright (C) 200
0 John Wiley & Sons, Ltd.