Received 24.01.2024, Revised 18.04.2024, Accepted 30.05.2024

Calculation of three-dimensional stationary potential thermal fields by the finite element method

Volodymyr Karashetsky, Volodymyr Yarkun

The article describes a clear algorithm for forming a system of nonlinear algebraic equations from the condition of the minimum of the functional to calculate the temperature distribution inside the three-dimensional domain filled with hysteresis-free nonlinear anisotropic environments. The boundary value problem of calculating the three-dimensional stationary potential temperature field is formulated. The finite element model for calculating the temperature distribution inside a three-dimensional domain filled with hysteresis-free nonlinear anisotropic environments is built. Basic formulas of the finite element method for the boundary value problem of calculating three-dimensional stationary potential thermal fields in domains filled with hysteresis-free nonlinear anisotropic environments with using Lagrangian tetrahedrons 1– 4 orders as finite elements and cubature formulas of numerical integration over the volume of the Lagrangian tetrahedron are derived. The algorithm determining the contribution of each finite element to the vector of residuals and matrix Jacobi  nonlinear system of equations solved by Newton’s method using the elements of the tensor of the differential thermal conductivity of the environment was considered. The cubature formula of numerical integration over the volume of the Lagrangian tetrahedron based on the interpolation complete polynomial for the Lagrangian finite element of the first order was applied. This algorithm is suitable when using Lagrangian finite elements of the second, third, and fourth orders with the use of the corresponding cubature formulas for numerical integration over the volume of the Lagrangian tetrahedron. Formulas for calculating tensors of differential thermal conductivity for non-linear isotropic and  linear environments are given. The boundary conditions of Dirichlet (first kind), Neumann (second kind), third and fourth kind and their consideration are described

potential thermal field, Lagrangian tetrahedron, finite element method, cubature formula, boundary conditions, Newton’s method
139-145
Karashetsky, V., & Yarkun, V. (2024). Calculation of three-dimensional stationary potential thermal fields by the finite element method. Information Technologies and Computer Engineering, 21(1), 139-145. https://doi.org/10.31649/1999-9941-2024-59-1-139-145

References

[1] Khalesidoost, S., Faiz, J., & Mazaheri, E. (2022). An overview of thermal modelling techniques for permanent magnet machines. IET Science, Measurement & Technology, 16(4), 219-241. doi: 10.1049/smt2.12099.

[2] Prajwal, K.T., & Bhat, P. (2021). Thermal analysis of a Thermoelectric Generator (TEG) using FEM technique. In First international conference on electrical energy and power engineering (ICEEPE 2020). IOP conference series: Materials science and engineering (Vol. 1045, article number 012018). Penang: IOP. doi 10.1088/1757-899X/1045/1/012018.

[3] Petresevics, F., & Nagy, B. (2022). FEM-based evaluation of the point thermal transmittance of various types of ventilated façade cladding fastening systems. Buildings, 12(8), article number 1153. doi: 10.3390/buildings12081153.

[4] Jindra, D., Hradil, P., Kala, J., & Salajka, V. (2020). Nonlinear FEM analysis οf composite concrete slab exposed τo extreme thermal load. In International conference of numerical analysis and applied mathematics ICNAAM 2019 (Vol. 2293, article number 130018). doi: 10.1063/5.0026880.

[5] Ponsati, T.L., Bahman, A.S., & Iannuzzo, F. (2021). Thermal modeling of large electrolytic capacitors using FEM and considering the internal geometry. In IEEE Journal of Emerging and Selected Topics in Power Electronics, 9(5), 6315-6328. doi: 10.1109/JESTPE.2021.3089899.

[6] Boukounacha, A.Y., Zegnini, B., Belkacem, Y., & Tahar, S. (2023). The effect of temperature on the thermal conductivity of transformer oils using the finite element method. In 1st international conference on materials sciences and applications (ICMSA2023) (article number 444816). Khenchela: ICMSA.

[7] Silvester, P.P., & Ferrari, R.L. (1996). Finite elements for electrical engineers. Cambridge: Cambridge University Press. doi: 10.1017/CBO9781139170611.

[8] Gallagher, R.H. (1975). Finite element analysis. Fundamentals (1st ed.). Victoria: Pearson College Div.

[9] Segerlind, L.J. (1984). Applied finite element analysis. Hoboken: J. Wiley & Sons

[10] Silvester, P., Cabayan, H.S., & Browne, B.T. (1973). Efficient techniques for finite element analysis of electric machines. – IEEE Transactions on Power Apparatus and Systems, PAS-92(4), 1274-1281. doi: 10.1109/TPAS.1973.293811.

[11] Dyshovyi, R.V. (1983). Calculation of the static magnetic field in implicit-pole electric machines by the differential grid method. (Abstract of Doctoral dissertation, Lviv Polytechnic National University, Lviv, Ukraine).

[12] Karashetsky, V.P. (2007). Cubature formulas for numerical integration over the volume of a tetrahedron based on interpolating full polynomials. Scientific Bulletin of the National Technical University of Ukraine, 17(6), 258-264.

[13] Yushko, S.V, Borshch, O.E., & Yushko, M.A. (2011). Stationary thermal conductivity. Kharkiv: NTU “KhPI”.