Skip to main content

Research Repository

Advanced Search

An enriched finite element model with q-refinement for radiative boundary layers in glass cooling

Mohamed, M.S.; Seaid, M.; Trevelyan, J.; Laghrouche, O.

An enriched finite element model with q-refinement for radiative boundary layers in glass cooling Thumbnail


M.S. Mohamed

O. Laghrouche


Radiative cooling in glass manufacturing is simulated using the partition of unity finite element method. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary simplified P1 approximation for the radiation in non-grey semitransparent media. To integrate the coupled equations in time we consider a linearly implicit scheme in the finite element framework. A class of hyperbolic enrichment functions is proposed to resolve boundary layers near the enclosure walls. Using an industrial electromagnetic spectrum, the proposed method shows an immense reduction in the number of degrees of freedom required to achieve a certain accuracy compared to the conventional h -version finite element method. Furthermore the method shows a stable behaviour in treating the boundary layers which is shown by studying the solution close to the domain boundaries. The time integration choice is essential to implement a q -refinement procedure introduced in the current study. The enrichment is refined with respect to the steepness of the solution gradient near the domain boundary in the first few time steps and is shown to lead to a further significant reduction on top of what is already achieved with the enrichment. The performance of the proposed method is analysed for glass annealing in two enclosures where the simplified P1 approximation solution with the partition of unity method, the conventional finite element method and the finite difference method are compared to each other and to the full radiative heat transfer as well as the canonical Rosseland model.


Mohamed, M., Seaid, M., Trevelyan, J., & Laghrouche, O. (2014). An enriched finite element model with q-refinement for radiative boundary layers in glass cooling. Journal of Computational Physics, 258, 718-737.

Journal Article Type Article
Acceptance Date Nov 6, 2013
Online Publication Date Nov 12, 2013
Publication Date Feb 1, 2014
Deposit Date Jan 22, 2014
Publicly Available Date Nov 12, 2014
Journal Journal of Computational Physics
Print ISSN 0021-9991
Publisher Elsevier
Peer Reviewed Peer Reviewed
Volume 258
Pages 718-737
Keywords Finite element method, Partition of unity method, Radiative heat transfer, Simplified P1 approximation, Glass cooling, Boundary layers.


Accepted Journal Article (4.6 Mb)

Copyright Statement
NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Computational Physics. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Computational Physics, 258, 2014, 10.1016/

You might also like

Downloadable Citations