The Effect of Time-Stepping on the Accuracy of Green Element Formulation of Unsteady Convective Transport

Despite the significant number of boundary element method (BEM) solutions of time-dependent problems, certain concerns still need to be addressed. Foremost among these is the impact of different time discretization schemes on the accuracy of BEM modeling. Although very accurate for steady-state problems, the boundary element methods more often than not are computationally challenged when applied to transient problems. For the work reported herein, we investigate the level of accuracy achieved with different time-discretization schemes for the Green element method (GEM) solution of the unsteady convective transport equation. The Green element method (a modified BEM formulation) solves the boundary integral theory (A Fredholm integral equation of the second kind) on a generic element of the problem domain in a way that is typical of the finite element method (FEM). In this integration process a new system of discrete equations is produced which is banded and hence amenable to matrix manipulations. This is subsequently deployed to investigate the proper resolution in both space and time for the chosen transient 1D transport problems especially those involving shock wave propagation and different types of boundary conditions. It is found that for three out of the four numerical models developed in this study, the new system of discrete element equations generated for both space and temporal domains exhibits accurate characteristics even for cases involving advection-dominant transport. And for all the cases considered, the overall performance relies heavily on the temporal discretization scheme adopted.

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

Onyejekwe, O. (2014) The Effect of Time-Stepping on the Accuracy of Green Element Formulation of Unsteady Convective Transport. Journal of Applied Mathematics and Physics, 2, 621-633. doi: 10.4236/jamp.2014.27068.

 [1] Clearly, R.W. and Ungs, M.J. (1978) Analytical Models for Groundwater Pollution and Hydrology. Report 78-WR-15, Water Resources Program, Princeton University, Princeton. [2] van Genuchten, M.Th. and Alves, W.J. (1982) Analytical Solutions of the One-Dimensional Convective-Dispersive Solute Transport Equation. Technical Bulletin 1661, US Department of Agriculture, Washington DC. [3] Peaceman, D.W. and Rachford Jr., H.H. (1962) Numerical Calculation of Multi-Dimensional Miscible Displacement. Society of Petroleum Engineers Journal, 2, 327-339 [4] Price, H.S., Cavendish, J.C. and Varga, R.S. (1968) Numerical Methods of Higher Order Accuracy for ConvectionDiffusion Equations. Society of Petroleum Engineers Journal, 16, 293-300. [5] Morton, K.W. (1996) Numerical Solution of Convection-Diffusion Problems. Chapman and Hall, London. [6] Li, X.-G. and Chan, C.K. (2004) The Finite Element Method with Weighted Basis Function for Singularly Perturbed Convection Diffusion Problems. Journal of Computational Physics, 195, 773-789. http://dx.doi.org/10.1016/j.jcp.2003.10.028 [7] Solin, P. (2005) Partial Differential Equations and the Finite Element Methods. John Wiley and Sons. http://dx.doi.org/10.1002/0471764108 [8] John, V., Mitkhova, T., Roland, M., Sundmacher, K., Tobiska, L. and Voigt, A. (2009) Simulations of Population Balance Systems with One Internal Coordinate Using Finite Element Methods. Chemical Engineering Science, 64, 733-741. [9] Fu, H.F. and Rui, H.X. (2012) A Mass-Conservative Characteristic Finite Element Scheme for Optimal Control Governed by Convection-Diffusion Equations. Compt. Methods. Applied Mechanical Engineering, 241, 82. [10] Taigbenu, A.E. and Onyejekwe, O.O. (1997) Transient 1D Transport Equation Simulated by a Mixed Green Element Formulation. International Journal for Numerical Methods in Fluids, 25, 437-454.http://dx.doi.org/10.1002/(SICI)1097-0363(19970830)25:4<437::AID-FLD570>3.0.CO;2-J [11] Onyejekwe, O.O. (2005) Green Element Method for 2D Helmholtz and Convection-Diffusion Problems with Variable Coefficients. Numerical Methods Partial Differential Equations, 21, 229-241. http://dx.doi.org/10.1002/num.20034 [12] Banerjee, P.K. (1994) The Boundary Element Methods. McGraw-Hill, London, New York. [13] Wrobel, L.C. and Brebbia, C.A. (1987) Dual Reciprocity Boundary Element Formulation for Nonlinear Diffusion Problems. Computer Methods in Applied Mechanics and Engineering, 65, 147-164. http://dx.doi.org/10.1016/0045-7825(87)90010-7 [14] Power, H. and Patridge, P.W. (1994) “Use of Stokes” Fundamental Solution for the Boundary Only Element Formulation of the Three-Dimensional Navier-Stokes Equations for Moderate Reynolds Numbers. International Journal for Numerical Methods in Engineering, 37, 1825-1840. [15] Bulgakov, V., Sarler, B. and Kuhn, G. (1998) Iterative Solution of Systems of Equations in the Dual Reciprocity Boundary Element Method for the Diffusion Equation. International Journal for Numerical Methods in Engineering, 43, 713-732. http://dx.doi.org/10.1002/(SICI)1097-0207(19981030)43:4<713::AID-NME445>3.0.CO;2-8 [16] Grigoriev, M.M. and Dargush, G.F. (2003) Boundary Element Methods for Transient Convective Diffusion, Part 1: General Formulation and 1D Implementation. Computer Methods in Applied Mechanics and Engineering, 192, 4281-4298. [17] Augustin, M., Caiazzo, A., Fiebach, A., Fuhrmann, J., John, V., Linke, A. and Ulma, R. (2011) An Assessment of Discretizations for Convection-Dominated Convection-Diffusion Equations. Computer Methods in Applied Mechanics and Engineering, 200, 3395-3409. http://dx.doi.org/10.1016/j.cma.2011.08.012 [18] Feng, X.F. and Tian, Z.F. (2006) Alternating Group Explicit Method with Exponential-Type for the Convection-Diffusion Equation. International Journal of Computer Mathematics, 83, 765-775. [19] Brebbia, C.A. and Skerget, P. (1984) Diffusion-Convection Problems Using Boundary Elements. In: Liable, J.P., et al., Eds., Proceedings: Finite elements in Water Resources, Springer-Verlag, New York, 747-768. http://dx.doi.org/10.1007/978-3-662-11744-6_63 [20] Taigbenu, A.E. and Liggett, J.A. (1986) An Integral Solution for the Diffusion-Advection Equation. Water Resources Research, 22, 1237-1246. http://dx.doi.org/10.1029/WR022i008p01237 [21] Grigoriev, M.M. (2000) Higher-Order Boundary Element Methods for Unsteady Convective Transport. Ph.D. Thesis, Faculty of the Graduate School of the State University of New York, New York. [22] Peratta, A. and Popov, V. (2003) Numerical Stability of the BEM for Advection-Diffusion Problems. http://dx.doi.org/10.1002/num.20009 [23] Taigbenu, A.E. (1995) The Green Element Method. International Journal for Numerical Methods in Engineering, 38, 2241-2263. http://dx.doi.org/10.1002/nme.1620381307 [24] Onyejekwe, O.O. (1998) A Boundary-Element-Finite-Element Equations Solutions for Flow in Heterogeneous Porous Media. Transport in Porous Media, 31, 293-312. http://dx.doi.org/10.1023/A:1006529122626 [25] Onyejekwe, O.O. (2012) Combined Effects of Shear and Buoyancy for Mixed Convection in an Enclosure. Advances in Engineering Software, 47, 188-193. http://dx.doi.org/10.1016/j.advengsoft.2011.11.002 [26] Tian, Z.F. and Yu, P.X. (2011) A High-Order Exponential Scheme for Solving 1D Unsteady Convection-Diffusion Equation. Journal of Computational Mathematics, 235, 2477-2491. http://dx.doi.org/10.1016/j.cam.2010.11.001