Publication Cover
Numerical Heat Transfer, Part B: Fundamentals
An International Journal of Computation and Methodology
Volume 70, 2016 - Issue 1
116
Views
4
CrossRef citations to date
0
Altmetric
Original Articles

New unconditionally stable scheme for solving the convection–diffusion equation based on the Associated Hermite orthogonal functions

&
Pages 80-90 | Received 05 Dec 2015, Accepted 20 Feb 2016, Published online: 23 Jun 2016
 

ABSTRACT

In this article, a new unconditionally stable scheme, based on the Associated Hermite orthogonal functions combined with first-order upwind scheme (AH-FUS), is proposed for solving the convection–diffusion equation. To eliminate the time variable from the computations, the time derivatives are expanded by Hermite functions, and a Galerkin's temporal testing procedure is introduced to the expanded equation. A set of implicit difference equations are derived in AH domain under no convergent conditions, and the numerical results can be obtained by solving the expanded coefficients recursively. Two numerical examples were considered to verify the accuracy and the efficiency of the proposed scheme.

Nomenclature

a=

diffusion coefficient

A=

banded matrix [Eq. (16)]

c=

convection coefficient

H=

Hermite polynomial

I=

unit matrix

S=

heat flux

t=

time

=

transformed time variable

T=

time interval

Tf=

time translating parameter

u=

solution of equation

α=

tridiagonal matrix defined in Eq. (15)

β=

unit matrix defined in Eq. (15)

λ=

matrix coefficient [Eq. (15)]

σ=

time-scaling parameter

ϕ=

Associated Hermite basis function

Ω=

computational domain

Subscripts and Superscripts=
i, j=

value at grid spatial points

k, m, n, p=

order of expansion coefficient

Nx=

number of grid points in x direction

Ny=

number of grid points in y direction

Nomenclature

a=

diffusion coefficient

A=

banded matrix [Eq. (16)]

c=

convection coefficient

H=

Hermite polynomial

I=

unit matrix

S=

heat flux

t=

time

=

transformed time variable

T=

time interval

Tf=

time translating parameter

u=

solution of equation

α=

tridiagonal matrix defined in Eq. (15)

β=

unit matrix defined in Eq. (15)

λ=

matrix coefficient [Eq. (15)]

σ=

time-scaling parameter

ϕ=

Associated Hermite basis function

Ω=

computational domain

Subscripts and Superscripts=
i, j=

value at grid spatial points

k, m, n, p=

order of expansion coefficient

Nx=

number of grid points in x direction

Ny=

number of grid points in y direction

Acknowledgment

This work was supported by PLA University of Science and Technology.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.