Transient 1D transport equation simulated by a mixed Green element formulation

✍ Scribed by Akpofure Efemena Taigbenu; Okey Oseloka Onyejekwe

• Publisher: John Wiley and Sons
• Year: 1997
• Tongue: English
• Weight: 286 KB
• Volume: 25
• Category: Article
• ISSN: 0271-2091
• DOI: 10.1002/(sici)1097-0363(19970830)25:4<437::aid-fld570>3.0.co;2-j

✦ Synopsis

New discrete element equations or coef®cients 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 coef®cient matrix which is amenable to ef®cient matrix solvers. It is herein derived for the transient 1D transport equation with uniform and non-uniform ambient ¯ow conditions and in which ®rst-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 ®rst 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 modi®ed fully implicit scheme with a difference weighting value of two give superior solutions in all simulated examples.