Русская версия English version

Study of thermal explosion conditions in a plate with variable properties and temperature-dependent heat source

E.V. Kotova, V.A. Kudinov, K.V. Kolotilkina, E.V. Stefanyuk, I.V. Kudinov

Vestnik IGEU, 2026 issue 5, pp. 36—45

Download PDF

Abstract in English: 

Background. The interest to obtain analytical solutions to heat conduction problems with variable physical properties and a temperature-dependent internal heat source, is relevant due to the need to study high-temperature processes of thermal ignition and explosion that occur in blow-up regimes. In these regimes heat becomes localized in a bounded region of space and leads to an unbounded increase in temperature (thermal explosion) within a short time interval. The study of such problems is mainly carried out by numerical methods, since their solution using known analytical methods is not feasible. In this regard, the development of analytical methods to get solutions to such problems is of current importance.

Materials and methods. To solve the specified heat conduction problems, an additional function (AF) and additional boundary conditions (ABCs) are introduced. These conditions consist of coordinate relations that ensure the heat conduction equation is satisfied at the boundaries. The AF represents the temporal variation of temperature at a certain point of the spatial variable. Its application makes it possible to reduce the order of the original equation, transforming it from a partial differential equation into an ordinary differential equation. The application of ABCs allows one to eliminate the process of integrating the original differential equation with respect to the spatial variable, as the equation is satisfied at the boundaries.

Results. A proof is presented for a theorem stating that if the equation holds on the boundary, it holds throughout the entire domain of the spatial variable. A solution has been found for the heat conduction problem involving variable physical properties and a temperature-dependent internal heat source. The obtained solution makes it possible to separate stationary regimes from regimes of unbounded temperature growth, which occur when the first eigenvalue is positive. It has been found that the temperature rise occurs within a narrow spatial region over a short time interval, that is, in a blow-up regime. A method is employed that makes it possible to eliminate integration of the equation with respect to the coordinate when a large number of additional boundary conditions are used. This is confirmed by the proof of the theorem. The obtained solution has been used to determine the heat source power at which the temperature gradient of the methane pyrolysis reactor wall does not exceed the permissible (specified) value.

Conclusions. By employing additional boundary conditions and an auxiliary unknown function, an analytical solution has been obtained for a highly complex nonlinear heat conduction problem involving a nonlinear heat source. The method described enables the analysis of combustion and explosion processes characterized by unbounded temperature rise, as well as other thermal processes in structures with internal heat sources, such as the induction heating of metals.

References in English: 

1. Lykov, A.V. Teoriya teploprovodnosti [Theory of Heat Conductivity]. Moscow: Vysshaya shkola, 1967. 601 p.

2. Grigor′ev, L.Ya., Man′kovskiy, O.N. Inzhenernye zadachi nestatsionarnogo teploobmena [Engineering Problems of Non-Stationary Heat Transfer]. Leningrad: Energiya, 1968. 84 p.

3. Kudinov, V.A., Averin, B.V., Stefanyuk, E.V. Teploprovodnost′ i termouprugost′ v mnogosloynykh konstruktsiyakh [Thermal Conductivity and Thermoelasticity in Multilayer Structures]. Moscow: Vysshaya shkola, 2008. 391 p.

4. Kotova, E.V., Eremin, A.V., Kudinov, V.A., Tkachev, V.K., Kuznetsova, A.E. Metod dopolnitel′nykh iskomykh funktsiy v zadachakh teploprovodnosti s peremennymi fizicheskimi svoystvami sredy [Method of additional unknown functions in heat onductivity problems with variable physical properties of the medium]. Vestnik IGEU, 2019, issue 2, pp. 59–70.

5. Kudinov, I.V., Kotova, E.V., Kudinov, V.A. Metod polucheniya analiticheskikh resheniy kraevykh zadach na osnove opredeleniya dopolnitel′nykh granichnykh usloviy i dopolnitel′nykh iskomykh funktsiy [A method for obtaining analytical solutions to boundary value problems based on determining additional boundary conditions and additional sought-for functions]. Sibirskiy zhurnal vychislitel′noy matematiki, 2019, vol. 22, no. 2, pp. 153–165.

6. Fedorov, F.M. Granichnyy metod resheniya prikladnykh zadach matematicheskoy fiziki [Boundary value method for solving applied problems of mathematical physics]. Novosibirsk: Nauka, 2000. 220 p.

7. Eremin, A., Krasnova, N. One method of heat transfer process mathematical modeling in solids. 2019 International Multi-Conference on Industrial Engineering and modern Technologies (FarEastCon). IEEE, 2019, pp. 1–5.

8. Pashchenko, D., Eremin, A. Heat flow inside a catalyst particle for steam methane reforming: CFD modeling and analytical solution. International Journal of Heat and Mass Transfer, 2021, vol. 165, p. 120617.

9. Vasyutin, M.A., Danilova, E.V., Kuz′michev, N.D. Otsenka maksimal′nogo nagreva plenki NbN tokovym impul′som pri nizkikh temperaturakh na osnove uravneniya teploprovodnosti [Estimation of the maximum heating of an NbN film by a current pulse at low temperatures based on the heat conduction equation]. Teplofizika i aeromekhanika, 2022, vol. 29, no. 3, pp. 473–480.

10. Karelin, V.A., Salomatov, V.V. Matematicheskaya model′ perenosa tepla v sloe dielektrika pri SVCH-obluchenii [Mathematical model of heat transfer in a dielectric layer under microwave irradiation]. Teplofizika i aeromekhanika, 2023, vol. 30, no. 2, pp. 375–380.

11. Karelin, V.A., Salomatov, V.V. Postroenie priblizhennykh analiticheskikh resheniy dlya modeli perenosa tepla v sloe l′da pri SVCH-obluchenii [Construction of approximate analytical solutions for the model of heat transfer in an ice layer under microwave irradiation]. Teplofizika i aeromekhanika, 2023, vol. 30, no. 6, pp. 1189–1196.

12. Dutta, D. Matematicheskaya model′ temperaturnogo polya, voznikayushchego pri dugovoy svarke pod flyusom [Mathematical model of the temperature field arising during submerged arc welding]. Prikladnaya mekhanika i tekhnicheskaya fizika, 2024, vol. 65, no. 3(385), pp. 74–82.

13. Kurdyumov, S.P., Malinetskiy, G.G., Potapov, A.B., Samarskiy, A.A. Struktury v nelineynykh sredakh [Structures in nonlinear media]. Komp′yutery i nelineynye yavleniya: Informatika i sovremennoe estestvoznanie [Computers and nonlinear phenomena: Information technology and modern natural science]. Moscow: Nauka, 1988. 192 p.

Key words in Russian: 
теплопроводность, пластина, переменные свойства, температурозависящий источник, дополнительные граничные условия, дополнительная функция, интеграл теплового баланса
Key words in English: 
thermal conductivity, plate, variable properties, temperature-dependent source, additional boundary conditions, additional function, heat balance integral
The DOI index: 
10.17588/2072-2672.2026.5.036-045
Downloads count: 
4