Complex Three-Dimensional Mathematical Model of the Ignition of a Coniferous Tree via a Cloud-to-Ground Lightning Discharge: Electrophysical, Thermophysical and Physico-Chemical Processes

: Thunderstorms are the main natural source of forest fires. The ignition mechanism of trees begins with the impact of cloud-to-ground lightning discharge. A common drawback of all predicting systems is that they ignore the physical mechanism of forest fire as a result of thunder-storm activity. The purpose of this article is to develop a physically based mathematical model for the ignition of a coniferous tree via cloud-to-ground lightning discharge, taking into account thermophysical, electrophysical, and physicochemical processes. The novelty of the article is explained by the development of an improved mathematical model for the ignition of coniferous trees via cloud-to-ground lightning discharge, taking into account the processes of soot formation caused by the thermal decomposition phase of dry organic matter. Mathematically, the process of tree ignition is described by a system of non-stationary nonlinear differential equations of heat conduction and diffusion. In this research, a locally one-dimensional method is used to solve three-dimensional partial differential equations. The finite difference method is used to solve one-dimensional heat conduction and diffusion equations. Difference analogues of the equations are solved using the marching method. To resolve nonlinearity, a simple iteration method is used. Temperature distributions in a structurally inhomogeneous trunk of a coniferous tree, as well as distributions of volume fractions of phases and concentrations of gas mixture components, are obtained. The conditions for tree trunk ignition under conditions of thunderstorm activity are determined. As a result, a complex three-dimensional mathematical model is developed, which makes it possible to identify the conditions for the ignition of a coniferous tree trunk via cloud-to-ground lightning discharge.


Introduction
Predicting forest fires can help to reduce damage [1].There are various methods for predicting forest fires, ranging from empirical to complex deterministic-probabilistic ones [2,3].As a rule, a specific information-computing system is based on some specific method of predicting forest fire danger.For example, in the Russian Federation, there is an information system for remote the monitoring of forest fires known as ISDM Rosleskhoz [4].The predictive component of this system is based on the Nesterov index [5].
It should be noted that in recent years, a deterministic-probabilistic method for predicting forest fire danger has been intensively developed [6].This method is based on deterministic mathematical models of the drying and ignition of forest fuels and probabilistic criteria [7].It is known [8] that forest fires from thunderstorms often occur in remote areas.Therefore, it is important to develop mathematical models for the ignition of Citation: Baranovskiy trees in conditions of lightning activity.The ignition of a tree by lightning is a complex phenomenon that includes various electrophysical, thermophysical and physicochemical processes.These processes occur both in the tree trunk and in the near-wall region of the trunk in the gas phase [9][10][11][12][13][14][15].
In the forestry structures of various states in the world community, instrumental systems are used to obtain information on lightning discharge.First of all, forest services are interested in information on the distribution of cloud-to-ground lightning discharges in time and space.There are examples of lightning direction-finding systems that provide information online: - On every continent, there are regions with a higher density of lightning discharges that lead to forest fires.For example, in Austria there is the Austrian Lightning Detection and Information System (ALDIS).This system can detect lightning strikes and provide access to information on lightning current, polarity, and complexity.In addition [16], information is analyzed using the time and place of fire in the forest, the type of vegetation, and the area covered by the fire.Using this information, the following [17] were calculated: FFMC-Fine Fuel Moisture Code; BUI-Build Up Index; and FWI-Fire Weather Index.An analysis of the application of ALDIS showed that this system makes it possible to recognize more than 90% of cloud-to-ground lightning discharges over the territory of Austria [18,19].A comparative analysis showed that this is one of the most accurate systems for the direction finding of such lightning discharges [20][21][22][23][24][25].
In turn, a lightning direction-finding system operates in Asia using SAFIR 300 sensors, which work according to the IMPACT method [26][27][28].In addition, the LINET method [29][30][31][32][33] is used, which makes it possible to determine almost half of cloud-to-ground lightning discharges.It should be noted that similar studies were carried out in Malaysia and Colombia [27,29].Lightning direction-finding systems can use two methods to determine the point of application of lightning discharge [16]: Magnetic Direction Finding (MDF) and Time of Arrival (TOA).These methods allow the system to determine the following discharge characteristics [21]: geographic coordinates, polarity, current, and number of strokes in a complex discharge.In addition, similar studies using lightning direction-finding systems were carried out in the territory of the Iberian Peninsula [34][35][36] and in the mountainous regions of the Alps [37].It should be noted that there is also a European Cooperation for Lightning Detection network (EUCLID) project in Europe [38].
The assessment of the impact of thunderstorm activity on the occurrence of forest fires can be carried out using maps of the density of lightning discharges passed from cloud to ground over a controlled forested area [39].It should be noted that precipitation is of great importance for the occurrence of forest fires under conditions of thunderstorm activity [40,41].On the one hand, a positive relationship has been established between precipitation and thunderstorm activity [42][43][44][45][46].This reduces the probability of forest fires occurring in rainy conditions.On the other hand, there are so-called dry thunderstorms, when the probability of ignition of forest fuels via cloud-to-ground lightning discharge is high [47].There are also more complex approaches to modeling thunderstorm activity in order to predict forest fires [48].For example, in Portugal a system has been developed using the CELLS atmospheric electrical model and the Meso-NH non-hydrostatic atmospheric model [49,50] to provide meteorological parameters based on the SURFEX model [51] and CELLS [52].
In addition, studies on the analysis of the characteristics of lightning discharges were carried out in the USA [53][54][55][56] and China [57].For example, an analysis of the characteristics of cloud-to-ground lightning discharges was carried out using the YNLDN [57,58].It was found that the maximum current is >70 kA and the minimum is <30 kA.For this analysis, a dataset of 400,000 lightning discharges was used [58].
This review provides information on lightning observation and direction-finding systems in various countries of the world community.Existing systems, such as WWLLN, can identify the duration and location of a lightning strike, as well as the strength of the discharge.There are mathematical expressions that allow researchers to obtain the value of peak lightning current.The duration and discharge current are the input data for the proposed mathematical model.Currently, a scenario-modeling approach is used when these parameters change within certain ranges.However, in future studies, when organizing interactions with GIS systems to predict forest fire danger, there will be a geographic reference to the locations of registered cloud-to-ground lightning discharges in a controlled forested area.Therefore, the review is directly related to the research topic.
Shortcomings or open questions in previously published work should be noted.Despite the fact that there has been intensive development in recent decades in the field of predicting forest fire danger, almost all existing systems do not have a physical basis.The forest fire danger prediction systems used in various countries do not have software modules that would allow researchers to take into account physical and chemical processes when a lightning discharge impacts trees.Typically, forest fire statistics and lightning direction-finding data are used.Therefore, it is necessary to develop new systems for predicting forest fire danger based on the real physical and chemical processes that occur during the impact of a cloud-to-ground lightning discharge with a coniferous tree.As a consequence, this is impossible without the development of complex mathematical models of the ignition of coniferous trees by lightning discharges.
The novelty of the article is explained by the development of an improved mathematical model for the ignition of coniferous trees via cloud-to-ground lightning discharge, taking into account the processes of soot formation caused by the thermal decomposition phase of dry organic matter.The development of this mathematical model allows us to get closer to the creation of hardware and software systems for the visual observation of the presence of soot particles in the observation zone using LiDAR atmospheric sensing technologies.The formation of soot emissions can be monitored, and the occurrence of a surface forest fire can be predicted.The mathematical model can be used to develop a new information and computing system for the scenario modeling of forest fire danger in conditions of thunderstorm activity.
The purpose of this article is to develop a physically based mathematical model for the ignition of coniferous trees via cloud-to-ground lightning discharge, taking into account thermophysical, electrophysical, and physicochemical processes.
Research objectives: (1) The formulation of a physical model of the coniferous tree trunk lightning ignition, taking into account the heating of wood due to release according to the Joule-Lenz law, thermal decomposition of dry organic matter, chemical reaction in the gas phase, and the formation of soot particles.(2) The development of a spatial mathematical model, taking into account the indicated physical effects and its software implementation.(3) The scenario modeling of heat and mass transfer in a tree trunk under the influence of a cloud-to-ground lightning discharge.(4) The formulation of conclusions and proposals for the practical application of the developed mathematical model in order to predict forest fires from the observation of thunderstorms.

Physical Statement
System of basic assumptions: (1) The subcrustal zone contains moisture and is a conductor of electric current, like a resistor [11,12]; (2) In the structure of the trunk, there is heterogeneity in the lower part of the branch in the form of reactive wood [11]; (3) In the structure of the cortical layer, there is a heterogeneity in the form of a crack that runs along the entire vertical of the trunk (Figure 1); (4) Only heat transfer in that part of the branch that grows out of the trunk is taken into account; (5) The main gaseous combustible component of the thermal decomposition products of dry organic matter is carbon monoxide; (6) The oxidation reaction of carbon monoxide into carbon dioxide is the leading reaction in the gas phase; (7) Heating of wood occurs as a result of heat release according to the Joule-Lenz law [12]; (8) Pyrolysis products enter the gas phase instantly; (9) Thermophysical characteristics do not depend on temperature.
The ignition of a coniferous tree trunk occurs when a critical temperature is reached in the zone of localization of the leading reaction and when the heat flux from the chemical reaction exceeds the heat flux from the subcrustal zone.

Mathematical Statement
The process of tree lightning ignition is described via a system of three-dimensional non-stationary nonlinear equations of heat conduction and diffusion (1)- (19).A locally one-dimensional finite-difference method was used for numerical implementation [59].Difference analogs of one-dimensional heat conduction equations are solved by using the marching method in combination with the simple iteration method [59].exp exp exp exp exp Initial and boundary conditions for the system of Equations ( 1)- (19).Boundary conditions outside the branch area: Boundary conditions on the inner side of the branch: Conditions on the boundary of the right side of the branch and the crack: Boundary conditions on the outer cut of the branch: Boundary conditions on the lower edge of the branch: Boundary conditions on the left side of the branch: Boundary conditions on the upper face of the branch: Conditions at the crack boundary, excluding the right side of the branch: ) Г25, e g T T = , (62) Initial conditions: Equation ( 1) describes the process of heat transfer in space and time in the core zone of a coniferous tree trunk, taking into account the heat drain due to the reaction of thermal decomposition of dry organic matter of wood.
Equation (2) describes heat transfer in time and space in the subcrustal zone of a coniferous tree trunk, taking into account the heat release according to the Joule-Lenz law and thermal decomposition of the dry organic matter of wood.
Equation (3) describes heat transfer in time and space in the region of the bark layer of the trunk, taking into account the thermal decomposition of dry organic matter of wood.
Equation (4) describes heat transfer in time and space in the zone of the upper part of the branch (grown into the core), taking into account the thermal decomposition of the dry organic matter of the branch.
Equation ( 5) describes heat transfer in space and time in the area of the upper part of the branch (grown into the subcrustal zone of the trunk), taking into account heat release according to the Joule-Lenz law and thermal decomposition of the dry organic matter of the branch.
Equation (6) describes heat transfer in time and space in the upper part of the branch (into the bark layer), taking into account the thermal decomposition of the dry organic matter of the branch wood.
Equation (7) describes heat transfer in time and space in the lower part of the branch (grown into the core), taking into account the thermal decomposition of the dry organic matter of the branch.
Equation (8) describes heat transfer in time and space in the lower part of the branch (grown into the subcrustal zone), taking into account the heat release according to the Joule-Lenz law and the thermal decomposition of the dry organic matter of the branch.Equation (9) describes heat transfer in time and space in the lower part of the branch (grown into the bark), taking into account the thermal decomposition of the dry organic matter of the branch.
Equation (10) describes heat transfer in time and space in the zone of the gas mixture in a crack in the bark layer, taking into account the heat release due to the oxidation reaction of carbon monoxide to carbon dioxide.
Equation (11) describes the diffusion transfer of a substance (air oxygen), taking into account the chemical reaction in the gas phase of carbon monoxide with air oxygen.
Equation ( 12) describes the diffusion transfer of a substance (carbon monoxide), taking into account the chemical reaction in the gas phase of carbon monoxide with atmospheric oxygen.
Relation (13) allows us to determine the concentration of inert components based on the balance of concentrations of the gas mixture.
Relation (14) allows us to determine the volume fraction of the gas phase based on the balance of the volume fractions of the phases.
Equation ( 15) allows us to determine the volume fraction of dry organic matter and describes the kinetics of thermal decomposition of this phase.
Equation ( 16) allows us to determine the rate of the chemical reaction of the oxidation of carbon monoxide into carbon dioxide using atmospheric oxygen.
Relationship ( 17) is an auxiliary expression to simplify the expression for the rate of the chemical reaction of the oxidation of carbon monoxide into carbon dioxide using atmospheric oxygen.
Equation ( 18) is the equation of state for a gas mixture.Equation ( 19) allows us to determine the volume fraction of soot particles and describes the kinetics of soot formation, which depends on the thermal decomposition of dry organic matter.
Designations of variables and indexes are presented in Appendix A. Roadmap of mathematical model design presented in Figure 3.The initial stage of the computational procedure involves establishing initial data for use in this program.Then, it is necessary to construct a three-dimensional mesh in a cylindrical system of coordinates r, φ, and z.Three-dimensional arrays of floating variables were also constructed for temperature, concentrations, and volume fractions within this stage.At the moment of time equal to zero, the initial data are arranged into arrays of variables.The next multistage step is to circle until the time is equal to the final time of tree exposure to cloud-to-ground lightning discharge.
It is necessary to compute temperature, concentrations and volume fractions for the next temporal layer.First, a set of heat conduction equations was solved using the finite difference method and locally one-dimensional method.Difference analogues of mathematical equations were solved using the marching method and the simple iteration method.Equations for core, subcrustal zone, bark, and air gap subsequently were solved at this stage.Second, a set of diffusion equations was solved using the finite difference method and locally one-dimensional method.As at the previous stage, difference analogues of mathematical equations were solved by the marching method in conjunction with the simple iteration method.The equations for carbon monoxide and oxygen were subsequently solved.Then, the concentration of inert components was calculated using the balance expression for concentration.Third, a set of ordinary differential equations were solved in order to compute the volume fractions of gas phase, dry organic matter and soot particles.Then, in the case of time being equal to exposure time, the resultant arrays were saved to files.Otherwise, the next iteration of circling was started.
Development environment: Embarcadero RAD Studio 10 Architect Edition (Delphi package) (Embarcadero, Austin, TX, USA) with Object Pascal programming language.Compilation was performed using standard Delphi tools.The visualization of the results was performed in the Origin Pro 2022 package (OriginLab Corporation, Northampton, MA, USA).Description of workstation is presented in Appendix B.

Results and Discussion
The distribution of temperature along the radius and height of the tree trunk at different points in time is shown in Figure 5: (a) t = 0.01 s; (b) 0.1 s; (c) 0.3 s; and (d) 0.5 s.The top figure corresponds to the section of the trunk of a coniferous tree in which the branch is missing.The lower figure corresponds to the section of the trunk of a coniferous tree in which a branch is present.An analysis of temperature distributions shows that, the subcrustal zone of a tree trunk warms up over time due to heat release according to the Joule-Lenz law.An exposure interval of lightning discharge up to 500 ms was considered.Not all cloud-to-ground lightning discharges are characterized by such an exposure time, nor do all such discharges lead to the ignition of a tree and a subsequent forest fire.The figures show that when the exposure time of lightning discharge is about 0.01 s, there is insignificant heating of the subcrustal zone of a coniferous tree trunk and practically no difference between areas with and without branches.At the time instant of 0.1 s, the temperature field at its maximum corresponds to approximately 400 K. Differences in the temperature of the regions with and without branches are already noticeable.In addition, nearby areas of the core and bark layer begin to warm up.At a time of 0.3 s, there is a noticeable difference in the temperatures of the regions with and without branches.This is due to the fact that the reactive wood of the branch has different thermophysical characteristics and this leads to a change.Namely, there is a decrease in temperature in the reactive wood zone of the coniferous branch.The maximum temperatures in the subcrustal zone and the maximum differences in branch temperatures are observed at an exposure time of 0.5 s.Also, temperature distributions along the horizontal cut of a tree trunk at different points in time are obtained.There is also an increase in temperature in the subcrustal zone over time.The maximum differences in the temperature of the areas with and without branches are typical for an exposure time of lightning discharge of 0.5 s.Moreover, temperature distributions along the radial cut of a tree trunk at different points in time were obtained.This section was isolated in the subcrustal zone of a coniferous tree trunk.It can be seen that the reactive wood of the branch has a lower temperature.That is, a full spatial analysis of the temperature field shows that the part of the branch with reactive wood is characterized by a field of lower temperature.All areas around it have a higher temperature due to the difference in the thermophysical characteristics of various fragments of the wood of the trunk of a coniferous tree, especially the reactive wood of the branch.Based on this fact, it can be assumed that more branches there are, the less likely that the trunk of a coniferous tree will ignite when exposed to lightning discharge.In principle, this fact can be used in the analysis of descriptions of forest areas and to identify the most fire-prone trees in the controlled forested area.This can be used as the basis for a geoinformation system for assessing the fire-related danger to forests under conditions of thunderstorm activity.In addition, in relation to the most fire dangerous areas, certain measures may be introduced to care for forest areas by forest fire services.It may make sense to forcibly plant species of coniferous trees with spreading crowns and thus a large number of branches during reforestation or around industrial infrastructure.This could be a deterrent to the passage of a thunderstorm front in a controlled forested area.
The radial distribution of temperature at the moment of ignition presented in Figure 6: 1-section outside the crack in the bark; 2-section passing through a crack in the bark.An analysis of the temperature curves in Figure 6 shows that, in the area of the subcrustal zone in the absence of a crack in the bark, an inert field of elevated temperature is formed.Under these conditions, the ignition of the tree trunk does not occur since there is no access for the oxidizer to the heated subcrustal zone.In the case of a crack in the bark layer, the situation is diametrically opposite.Gaseous pyrolysis products diffuse into the crack area and, mixing with the oxidizing agent, enter a chemical reaction at an elevated temperature.It is in the crack that the place of localization of ignition is located.
The fields of volume fractions of dry organic matter and the gas phase in the region of the subcrustal zone of a coniferous tree trunk and a crack in the bark layer are presented in Figure 7.The most intense thermal decomposition processes occur in the subcrustal zone of a coniferous tree trunk, releasing gaseous combustible pyrolysis products.It should be noted here that there are different approaches to modeling the composition of dry organic matter [60].For example, the following components can be suggested: lignin, cellulose and hemicellulose.Each of these components is characterized by a certain set of products that are formed during the thermal decomposition of a certain component.
In this work, we used the gross substance approach, during which the thermal decomposition of this substance produces gaseous and condensed pyrolysis products.It should be noted that the formation of liquid pyrolysis products is not taken into account in this work.It is believed that the gaseous product of pyrolysis is carbon monoxide, which reacts with oxygen and oxidizes into carbon dioxide.Other gaseous combustible components are neglected since their proportion is noticeably smaller than the formation of carbon monoxide during the thermal decomposition of dry organic matter.Accordingly, the condensed phase is represented by residues of dry organic matter and dispersed soot particles.Moreover, the part of soot particles that is formed during the period of intense thermal decomposition is considered.The soot formation process is considered to follow the thermal decomposition process.
It is known that there are video surveillance systems for smoke plumes in forest areas near wildland-urban interface areas.Perhaps the development of software based on the proposed mathematical model could improve the technical characteristics of such video surveillance systems for forest areas.In this article, a console application has been developed using the high-level programming language Object Pascal in RAD Studio Delphi 10 Architect.There are no obstacles to making a software module for such a system based on the software implementation of the developed mathematical model for the ignition of coniferous trees via cloud-to-ground lightning discharge.
The dependences of the volume fractions of dry organic matter and the gas phase in the subcrustal zone of a coniferous tree trunk at the location of a crack in the bark is presented in Figure 8.It can be seen that the most significant processes of thermal decomposition of dry organic matter occur at exposure times above 0.3 s.That is, short-term discharge will not lead to the ignition of a coniferous tree trunk.This is due not only to insufficient heating of the wood in the zone of the electric current flow of a cloud-to-ground lightning discharge, but also to an insufficient concentration of combustible gaseous pyrolysis products of dry organic matter.
However, even prolonged exposure to lightning discharge does not lead to complete thermal decomposition of dry organic matter in the subcrustal zone of a coniferous tree trunk.In fact, an unstable porous structure of the material is formed in this zone, which can be subject to destruction.Previously, field observations of thunderstorm activity in the Timiryazevsky forestry of the Tomsk region were conducted [61].As a result of numerous observations, a strike by such discharge into the trunk of a coniferous tree was recorded.A subsequent inspection of the adjacent territory showed the presence of charred wood fragments in the ground cover, which could potentially lead to a forest fire in dry thunderstorm conditions.This confirms the possibility of the destruction of the tree trunk in the subcrustal zone.Inspection of the tree trunk also confirmed the breakdown of the bark layer and the presence of traces of wood destruction in the subcrustal zone.
The dependences of the concentrations of the gas phase components at the moment of ignition presented in Figure 9: 1-section in the crack; 2-section outside the crack.An analysis of the results of numerical simulation showed that, over time, the formation of combustible gaseous pyrolysis products occurs in the subcrustal zone.It is believed that such products instantly appear in the near-wall region in the bark layer.When considering a section of the bark with no cracks, it is found that the combustible gaseous products of pyrolysis cannot reach the required critical ignition temperature.In addition, there is no oxidizing agent in this zone, which also does not contribute to ignition.In the area with a crack in the bark layer, a certain concentration of reacting components is also achieved, but this occurs when already in the presence of an oxidizing agent and when the critical temperature of the gas mixture is reached.It should be noted that the option of having a bark barrier is considered when the crack is partially filled with bark.It is found that, after reaching a thickness of the bark barrier of more than 1.5 mm, the gas mixture may not be heated to a critical temperature.As a result, the process of ignition of a coniferous tree trunk does not occur.The dependence of the volume fractions of phases are presented in Figure 10, taking into account the formation of soot particles in the subcrustal zone of a coniferous tree trunk.These are typical simulation results.The analysis shows that the greatest soot formation processes begin when the exposure time of lightning discharge is more than 0.3 s.In addition, scenario calculations are carried out to take into account the intensity of the soot formation process.This factor is taken into account by varying the dispersion coefficient, which leads to a change in the volume fraction of soot formed when a coniferous tree trunk is ignited via cloud-to-ground lightning discharge.Tables 1 and 2 show the ignition delays of a coniferous tree depending on the parameters of lightning discharge.The conditions for the ignition of a coniferous tree in the context of the current-voltage characteristics of a cloud-to-ground lightning discharge are established.It is known [62] that the ignition conditions for large wood fragments can be established by comparing the ignition surface temperature and the heat flux.Submodules are included in the computer program to calculate these parameters along with the calculation of the ignition delay time.Analysis of the data presented in Tables 1 and 2 shows good agreement between both the experimental conditions for the ignition of wood and the conditions theoretically determined as a result of solving the presented mathematical model.Therefore, the developed mathematical model and its software implementation can be used to calculate the conditions and parameters of the ignition of a tree by an electric current of a lightning discharge passing through the tree trunk.The reliability of the numerical results is shown both when varying the electric current and the voltage of the lightning discharge.
The following indicators for assessing the results obtained have been introduced: -Possibility of practical use; -Reliability, -Completeness of information.
First of all, such an indicator as the completeness of information obtained using the proposed mathematical model must be discussed.
Numerical implementation of the mathematical model and scenario modeling make it possible to determine the following parameters as a result of modeling: -Ignition delay.This parameter can be used to parametrically assess the probability of a forest fire.-Temperature field in the tree trunk structure.This parameter allows us to determine the location of elevated temperatures, which can also be used to determine the probability of a forest fire.However, it must be noted that, for scenarios of modeling the formation of heated wood particles that fall on the ground layer of forest fuel material, it is also necessary to model the stress-strain state of the tree trunk, and such work will be carried out in future studies.It will then be possible to simulate the ignition of the ground cover, taking into account the size and temperature of the heated particles formed as a result of the destruction of the trunk wood.Future research will determine thermal stress in the trunk of a coniferous tree.-Volume fraction of soot particles.Information about this parameter can be used when setting up information and computing systems within the framework of hardware and software systems for the detection of soot and smoke particles in a controlled area.That is, information about the location of the emission of soot particles can be used to detect potential forest fires.
It is also necessary to discuss such a parameter as the possibility of the practical use of the proposed mathematical model of ignition of coniferous trees via cloud-to-ground lightning discharge.Below comparative information is presented on the execution time of program code developed on the basis of one-dimensional, two-dimensional and three-dimensional mathematical models.The following modeling scenario was used, reflecting typical situations under conditions of thunderstorm activity: Table 3 presents the results of a comparative analysis of the execution time of various software implementations.Thus, the use of a three-dimensional mathematical model is possible only with parallel implementation on supercomputers with multiprocessor and multicore computer system architecture.It should be noted that parallel implementation of the proposed mathematical model can be performed using massive parallelism.In this case, iterations of cyclic constructions are spread across the different computing cores of the computer system.Since the three-dimensional equations of thermal conductivity and diffusion are solved using a locally one-dimensional method, three blocks are distinguished for which parallel calculations can be carried out: these blocks are divided along the coordinates r, φ, and z.Taking into account the given data on the spatial grid covering the modeling area, it is necessary to use at least 10 computing nodes of a multiprocessor system.If necessary, you can make a flexible algorithm that will use 100 processor cores when modeling a block at the r coordinate, 10 processor cores when modeling at the φ coordinate, and 20 processor cores when modeling at the z coordinate.However, even in the worst-case scenario, the acceleration will approach the theoretical value of 10 times.The difference from theoretical estimates is due to the fact that part of the time will be spent on sending data at the beginning of the parallel program, collecting data at the end of calculations, and exchanging data between the blocks of the parallel program.Now, as for the reliability of the results obtained.As has already been shown in the tables, there is qualitative agreement with the results of experimental data published in [62].The volume fraction of soot particles also qualitatively corresponds to the scenarios for modeling soot formation under the influence of a forest fire sustained by forest fuels [13] and the removal of heated particles from the forest fire front [14].It is clear that there is a difference in the time coordinate, since the process of the impact of a lightning discharge on a tree trunk is significantly shorter than the processes indicated above.However, qualitative agreement between the results can be established.This indicates the physical adequacy and satisfactory reliability of the proposed mathematical model.However, in the future it will be necessary to conduct another series of field observations of cloud-to-ground lightning discharges in a controlled forested area.zone of a coniferous tree trunk.Thirdly, the proportion of soot formed during the combustion of dry organic matter is not considered, although a slight increase in the total amount of soot particles can be expected.
In addition, future research should be discussed.This paper presents a fairly complex mathematical model for the ignition of a coniferous tree trunk in conditions of lightning activity, taking into account electrophysical, thermophysical, and physicochemical processes.However, even this mathematical model can be improved.This can be performed, first of all, by taking into account the more plausible geometry of the study area, that is, the structure of the trunk of a coniferous tree.Secondly, for practical purposes, it is necessary to develop an application with a graphical user interface capable of organizing interactions with a lightning discharge direction-finding system, such as WWLLN [64].

Conclusions
In the present work, a generalized three-dimensional mathematical model for the ignition of coniferous trees via cloud-to-ground lightning discharge was developed, taking into account a complex of electrophysical, thermophysical, and physicochemical processes.It was found that the time spent on calculations is significant and the use of single-processor computing technology does not have great prospects in the context of predicting forest fire danger.It would be more logical in the future to develop a parallel software implementation of the proposed mathematical model.
Key findings and conclusions: (1) The presence in the branch structure of heterogeneities like reactive wood causes the formation of a low-temperature field.It is possible that trees with spreading crowns are less likely to be ignited by lightning discharge.(2) The electric current passes in the subcrustal zone and leads to the formation of a zone of elevated temperature and the thermal decomposition of dry organic matter, with the formation of gaseous combustible pyrolysis products.(3) Gaseous products of thermal decomposition enter the gas phase and ignite, which leads to the ignition of the coniferous tree trunk in the crack.(4) The presence of a crustal barrier in a crack makes it impossible to heat the reactants of the gas mixture in the crack and to ignite them in the gas phase.(5) During the thermal decomposition of dry organic matter, a certain proportion of soot particles is formed.It is possible that the remote sensing of the surface layer of the atmosphere for the presence of soot particles can be used to monitor the occurrence of forest fires.

Figure 1 .Figure 2 .
Figure 1.The structure of the bark layer of a coniferous tree (Republic of Altay, Russian Federation).The solution area is shown in Figure2aand the boundaries of the regions are indicated in Figure2b.

Figure 3 .
Figure 3. Roadmap of mathematical model design.The workflow diagram of the presented mathematical model is shown in Figure 4.

Figure 5 .
Figure 5. Temperature distribution along the radius and height of the tree trunk at different points in time: (a) t = 0.01 s; (b) 0.1 s; (c) 0.3 s; (d) 0.5 s.

Figure 6 .
Figure 6.Radial temperature distribution at the moment of ignition: 1-section outside the crack in the bark; 2-section passing through a crack in the bark.

Figure 7 .
Figure 7. Distribution of volume fractions of phases: (a) dry organic matter; (b) gas.

Figure 8 .
Figure 8. Dependences of volume fractions of phases on time: 1-volume fraction of dry organic matter; 2-volume fraction of gas phase.

Figure 9 .
Figure 9. Dependences of the concentrations of the gas phase components at the moment of ignition: 1-section in the crack; 2-section outside the crack.

Figure 10 .
Figure 10.Typical dependences taking into account soot formation: 1-volume fraction of dry organic matter; 2-volume fraction of gas phase; 3-volume fraction of soot particles.

Table 1 .
Wood ignition delays depending on discharge voltage at current J = 23.5 kA.

Table 2 .
Ignition delays of a tree trunk depending on the current at a voltage of U = 100 kV.

Table 3 .
Execution time of various software implementations.