ePrints@IIScePrints@IISc Home | About | Browse | Latest Additions | Advanced Search | Contact | Help

An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence

Ganesan, Sashikumaar (2012) An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence. In: ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE, 46 (6). pp. 1447-1465.

[img] PDF
esaim_math_mod_num_ana_46_1447-1465_2012.pdf - Published Version
Restricted to Registered users only

Download (299kB) | Request a copy
Official URL: http://dx.doi.org/10.1051/m2an/2012012

Abstract

We present a heterogeneous finite element method for the solution of a high-dimensional population balance equation, which depends both the physical and the internal property coordinates. The proposed scheme tackles the two main difficulties in the finite element solution of population balance equation: (i) spatial discretization with the standard finite elements, when the dimension of the equation is more than three, (ii) spurious oscillations in the solution induced by standard Galerkin approximation due to pure advection in the internal property coordinates. The key idea is to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In the proposed splitting scheme, the shape of the physical domain can be arbitrary, and different discretizations can be applied to the low-dimensional equations. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. The stability and error estimates of the Galerkin/SUPG finite element discretization of the population balance equation are derived. It is shown that a slightly more regularity, i.e. the mixed partial derivatives of the solution has to be bounded, is necessary for the optimal order of convergence. Numerical results are presented to support the analysis.

Item Type: Journal Article
Publication: ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE
Publisher: CAMBRIDGE UNIV PRESS
Additional Information: Copy right for this article belongs to CAMBRIDGE UNIV PRESS
Keywords: Population balance equations;operator-splitting method;error analysis;streamline Upwind Petrov Galerkin finite element methods;backward Euler scheme
Department/Centre: Division of Interdisciplinary Sciences > Supercomputer Education & Research Centre
Date Deposited: 21 Jul 2012 06:48
Last Modified: 21 Jul 2012 06:51
URI: http://eprints.iisc.ac.in/id/eprint/44815

Actions (login required)

View Item View Item