3D hybrid finite element enthalpy for anisotropic thermal conduction analysis. (June 2019)

• Record Type: Journal Article
• Authors: Erchiqui, F.; Annasabi, Z.
• Is Part Of: International journal of heat and mass transfer. Volume 136(2019)
• Page Start: 1250
• Page End: 1264
• Publication Date: 2019-06
• Subjects: Anisotropic thermal conduction -- Enthalpy -- Finite element method -- Anisotropic cylindrical -- Anisotropic spherical -- Wood thawing -- Wood freezing; Heat -- Transmission -- Periodicals; Mass transfer -- Periodicals; Chaleur -- Transmission -- Périodiques; Transfert de masse -- Périodiques; Electronic journals; 621.4022
• DOI: 10.1016/j.ijheatmasstransfer.2019.02.096
• Languages: English
• ISSNs: 0017-9310

Abstract

Highlights: 3D-Hybrid Finite Element Enthalpy (HFEM) approach is proposed for Anisotropic Thermal Conduction Analysis, with or without a phase change. HFEM using Cartesian coordinates to treat the thermal anisotropic conduction cylindrical, spherical and Cartesian. HFEM is validated with analytical testing for cylinder with convective transfer with variable circumference. HFEM is validated with analytical testing for spherical with convective transfer with variable circumference. HFEM is validated with two experimental tests related to the heating of frozen woods. Abstract: The anisotropic problem of thermal conduction in a solid is generally treated in a reference coordinate system, which adequately describes its thermal conductivity tensor (Cartesian, cylindrical or spherical). For this problem, numerical treatment is difficult, especially if the thermophysical properties are non-linear or if the anisotropic medium undergoes a phase change. In this paper, we propose an approach using a Cartesian reference system to treat the anisotropic thermal conduction of problems for which the solid medium is characterized by a set of tensors of thermal conductivity of different natures (Cartesian and/or cylindrical and/or spherical), with or without phase change. For this purpose, we the anisotropic thermal conductivity tensor, with respect to a cylindrical or spherical coordinate system, is transformed by an equivalent tensor into global Cartesian coordinates. The nonlinear heat conduction problem involving phase changes, such as wood freezing, is solved using hybrid three-dimensional volumetric specific enthalpy based on finite-element analysis. The proposed approach is validated with analytical testing for two anisotropic media and with two experimental tests related to the heating of frozen woods. As an application, we have numerically quantified, on the one hand, the minimum time required for the thaw and, on the other hand, the freezing of a log of wood, such as white pine, according to the length of its radius (7.5, 10, 15, 20 and 25 cm). The thermophysical properties are a function of temperature, moisture content and structural orientation.