TRANSIENT 1D TRANSPORT EQUATION SIMULATED BY A MIXED GREEN ELEMENT FORMULATION

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.

[1]  Green element simulations of the transient nonlinear unsaturated flow equation , 1995 .

[2]  戸田 保幸.  The Finite Analytic Method , 1992 .

[3]  D. Spalding.  A novel finite difference formulation for differential expressions involving both first and second derivatives , 1972 .

[4]  H. H. Rachford,et al.  Numerical Calculation of Multidimensional Miscible Displacement , 1962 .

[5]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[6]  O. Zienkiewicz,et al.  The coupling of the finite element method and boundary solution procedures , 1977 .

[7]  W. Yeh,et al.  A proposed upstream weight numerical method for simulating pollutant transport in groundwater , 1983 .

[8]  W. Gray,et al.  An analysis of the numerical solution of the transport equation , 1976 .

[9]  O. Onyejekwe,et al.  A MIXED GREEN ELEMENT FORMULATION FOR THE TRANSIENT BURGERS EQUATION , 1997 .

[10]  R. Banks,et al.  A solution of the differential equation of longitudinal dispersion in porous media , 1961 .

[11]  Mary P. Anderson,et al.  Using models to simulate the movement of contaminants through groundwater flow systems , 1979 .

[12]  Richard S. Varga,et al.  Numerical methods of higher-order accuracy for diffusion-convection equations , 1968 .

[13]  Okey Oseloka Onyejekwe,et al.  GREEN ELEMENT DESCRIPTION OF MASS TRANSFER IN REACTING SYSTEMS , 1996 .

[14]  An Integral Solution for the Diffusion-Advection Equation , 1986 .

[15]  Derek Elsworth,et al.  A boundary element‐finite element procedure for porous and fractured media flow , 1987 .

[16]  J. Z. Zhu,et al.  The finite element method , 1977 .

[17]  I. Raju,et al.  Coupling finite and boundary element methods for two‐dimensional potential problems , 1993 .

[18]  E. Dick.  Accurate Petrov‐Galerkin methods for transient convective diffusion problems , 1983 .

[19]  A. Taigbenu.  The green element method , 1999 .

[20]  O. C. Zienkiewicz,et al.  An ‘upwind’ finite element scheme for two‐dimensional convective transport equation , 1977 .