Chebyshev expansion methods for the solution of the extended Graetz problem (Q1062565)

The extended Graetz problem, a classical problem described by an equation with a non-self-adjoint second-order elliptic differential operator, is solved by two numerical methods based on Chebyshev expansions: A Chebyshev-finite difference method and a Chebyshev-''finite'' element method. The latter relies entirely on Chebyshev expansions, global in the bounded (radial) coordinate direction and local in the unbounded (axial) direction and proves to be more accurate than the Chebyshev-finite difference method in resolving singularities. Both methods can be applied to general type boundary-value problems involving second-order elliptic operators, have accuracy comparable to high order finite difference schemes and are operation cost efficient.