V.I. Havrysh

Èlektron. model. 2023, 45(5):03-19


Linear and non-linear mathematical models for the determination of the temperature field, and subsequently for the analysis of temperature regimes in isotropic media with semi-transparent foreign inclusions that are also subject to external heat load, have been developed. To solve the nonlinear boundary-value problem, a linearizing function was introduced, using which the original nonlinear heat conduction equation and nonlinear boundary conditions were linearized, and as a result, a partially linearized second-order differential equation with partial derivatives and discontinuous and singular coefficients relative to the linearizing function and partially linearized boundary conditions was obtained. For the final linearization of the partially linearized differential equation and boundary conditions, the approximation of the temperature according to one of the spatial coordinates on the boundary surfaces of the inclusion was performed by piecewise constant functions. To solve the obtained linear boundary value problem, the Hankel integral transformation method was used, as a result of which an analytical solution was obtained, which determines the introduced linearizing function. For the numerical analysis of temperature behavior and heat exchange processes caused by external heat load, software tools have been developed, which are used to create a geometric image of temperature distribution depending on spatial coordinates. The obtained results indicate the correspondence of the developed mathematical models of the analysis of heat exchange processes in heterogeneous media with external heating to a real physical process. The obtained results indicate the correspondence of the developed mathematical models of the analysis of heat exchange processes in heterogeneous media with external heating to a real physical process. The developed models make it possible to analyze environments with foreign elements under external thermal loads regarding their thermal resistance. As a result, it becomes possible to increase it and protect it from overheating, which can cause the destruction of not only individual elements, but also the entire structure.


temperature field; isotropic heterogeneous medium; thermal conductivity; convective heat exchange; perfect thermal contact; external heating; heat flow; thermosensitivity, half-through foreign inclusion.


  1. Car­pin­te­ri, A., & Pag­gi, M. (2008). Ther­mo­elas­tic mis­match in non­ho­mo­ge­ne­ous be­ams. Jo­ur­nal of En­gi­ne­ering Mat­he­ma­tics, 61(2-4), 371-384
  2. No­da, N. (1991). Ther­mal stres­ses in ma­te­ri­als with tem­pe­ra­tu­re-de­pen­dent pro­per­ti­es. Appli­ed Mec­ha­nics Re­vi­ews, 44, 383-397
  3. Otao, Y., Ta­ni­ga­wa, O., & Is­hi­ma­ru, O. (2000). Op­ti­mi­za­ti­on of ma­te­ri­al com­po­si­ti­on of functi­ona­lity gra­ded pla­te for ther­mal stress re­la­xa­ti­on using a ge­ne­tic al­go­rithm. Jo­ur­nal of Ther­mal Stres­ses, 23, 257 - 271
  4. Ta­ni­ga­wa, Y., Akai, T., & Ka­wa­mu­ra, R. (1996). Tran­si­ent he­at con­duc­ti­on and ther­mal stress prob­lems of a non­ho­mo­ge­ne­ous pla­te with tem­pe­ra­tu­re-de­pen­dent ma­te­ri­al pro­per­ti­es. Jo­ur­nal of Ther­mal Stres­ses, 19(1), 77-102
  5. Ta­ni­ga­wa, Y., & Otao, Y. (2002). Tran­si­ent ther­mo­elas­tic analysis of functi­onally gra­ded pla­te with tem­pe­ra­tu­re-de­pen­dent ma­te­ri­al pro­per­ti­es ta­king in­to ac­co­unt the ther­mal ra­di­ati­on. Ni­hon Ki­kai Gak­kai Nen­ji Ta­ikai Ko­en Ron­bunshu, 2,133-134.­mec­jo.2002.2.0_133
  6. Yan­gi­an, Xu. (2009). Analysis of ste­ady ther­mal stress in a ZrO2/FGM/Ti-6Al-4V compo­site ECBF plate with temperature-dependent material properties by NFEM. & Daihui, Tu (Ред.), 2009-WA­SE Int. Conf. on In­for­ma. Eng, Vol. 2, (433-436).­TA.2009.842
  7. Dovbnia, K.M., & Dun­dar, O.D. (2016). Stat­si­onarnyi tep­lo­ob­min tonkykh po­lohykh izot­ropnykh obo­lo­nok, ya­ki znak­ho­di­at­sia pid diieiu dzhe­rel tep­la, zo­se­redzhenykh po dvovy­mir­nii ob­las­ti. Visnyk Donetckogo ­natcionalnogo univercutety. Seriya A: Pryrodnychi na­u­ky, 1-2, 107-112.
  8. Aza­ren­kov, V.I. (2012). Issle­do­va­nie i raz­ra­bot­ka tep­lo­voi mo­de­li i me­to­dov ana­li­za tem­pe­ra­tur­nikh po­lei konstruktcii ra­dioelektron­noi ap­pa­ra­tu­ri. Techno­logy au­dit and pro­duc­ti­on re­ser­ves, 3/1(5), 39-
  9. Havrysh, V.I., & Grys­juk, Y.I. (2022). Temperature fields in heterogeneous enviroments­ with consideration of thermal sensitivity. Lviv: Pub­lis­hing hou­se of Lviv Po­li­technic Na­ti­o­­nal Uni­ver­sity.
  10. Vasyl Havrysh. (2022). Modelling temperature field in spatial thermosensitive elements of electronic devices with local thermal heating. У & Orest Kochan (Ред.), Computer Scien­ce and Information Technologies. Materials of the XVIIth International Scientific and Tech­nical Conference, (CSIT-2022), Lviv, Ukraine 10-12 November (547-550).
  11. Vasyl Havrysh, Lubov Kolyasa, & Svitlana Vozna (2021). Temperature field in a layered plate with local heating. International scientific journal “Mathematical modeling”. Vol. 5, Issue 3, 90-94.

Full text: PDF