Premium
Transient 1D transport equation simulated by a mixed Green element formulation
Author(s) -
Taigbenu Akpofure Efemena,
Onyejekwe Okey Oseloka
Publication year - 1997
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/(sici)1097-0363(19970830)25:4<437::aid-fld570>3.0.co;2-j
Subject(s) - mathematics , mathematical analysis , boundary element method , finite element method , matrix (chemical analysis) , boundary value problem , partial differential equation , convection–diffusion equation , integral equation , diffusion equation , physics , materials science , composite material , thermodynamics , economy , economics , service (business)
New discrete element equations or coefficients are derived for the transient 1D diffusion–advection or transport equation based on the Green element replication of the differential equation using linear elements. The Green element method (GEM), which solves the singular boundary integral theory (a Fredholm integral equation of the second kind) on a typical element, gives rise to a banded global coefficient matrix which is amenable to efficient matrix solvers. It is herein derived for the transient 1D transport equation with uniform and non‐uniform ambient flow conditions and in which first‐order decay of the containment is allowed to take place. Because the GEM implements the singular boundary integral theory within each element at a time, the integrations are carried out in exact fashion, thereby making the application of the boundary integral theory more utilitarian. This system of discrete equations, presented herein for the first time, using linear interpolating functions in the spatial dimensions shows promising stable characteristics for advection‐dominant transport. Three numerical examples are used to demonstrate the capabilities of the method. The second‐order‐correct Crank–Nicolson scheme and the modified fully implicit scheme with a difference weighting value of two give superior solutions in all simulated examples. © 1997 John Wiley & Sons, Ltd.