An enriched finite element model with q-refinement for radiative boundary layers in glass cooling
Abstract
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 P{sub 1} 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 ismore »
- Authors:
-
- Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom)
- School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom)
- Publication Date:
- OSTI Identifier:
- 22230861
- Resource Type:
- Journal Article
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Volume: 258; Other Information: Copyright (c) 2013 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICAL METHODS AND COMPUTING; APPROXIMATIONS; BOUNDARY LAYERS; COMPARATIVE EVALUATIONS; EQUATIONS; FINITE DIFFERENCE METHOD; FINITE ELEMENT METHOD; GLASS; HEAT TRANSFER; RADIATIVE COOLING; REDUCTION
Citation Formats
Mohamed, M. Shadi, Seaid, Mohammed, Trevelyan, Jon, and Laghrouche, Omar. An enriched finite element model with q-refinement for radiative boundary layers in glass cooling. United States: N. p., 2014.
Web. doi:10.1016/J.JCP.2013.11.005.
Mohamed, M. Shadi, Seaid, Mohammed, Trevelyan, Jon, & Laghrouche, Omar. An enriched finite element model with q-refinement for radiative boundary layers in glass cooling. United States. https://doi.org/10.1016/J.JCP.2013.11.005
Mohamed, M. Shadi, Seaid, Mohammed, Trevelyan, Jon, and Laghrouche, Omar. 2014.
"An enriched finite element model with q-refinement for radiative boundary layers in glass cooling". United States. https://doi.org/10.1016/J.JCP.2013.11.005.
@article{osti_22230861,
title = {An enriched finite element model with q-refinement for radiative boundary layers in glass cooling},
author = {Mohamed, M. Shadi and Seaid, Mohammed and Trevelyan, Jon and Laghrouche, Omar},
abstractNote = {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 P{sub 1} 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 P{sub 1} 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.},
doi = {10.1016/J.JCP.2013.11.005},
url = {https://www.osti.gov/biblio/22230861},
journal = {Journal of Computational Physics},
issn = {0021-9991},
number = ,
volume = 258,
place = {United States},
year = {Sat Feb 01 00:00:00 EST 2014},
month = {Sat Feb 01 00:00:00 EST 2014}
}