Kemerovo State University(Russia) Mathematical Modeling of Large Forest Fires

1. Kemerovo State University(Russia)Mathematical Modeling of Large Forest Fires Valeriy A. Perminov pva@belovo.kemsu.ru, p_valer@mail.ru

3. Mathematical modeling using the deterministic approach contains the following relatively independent stages: • Physical analysis of the phenomenon of forest fire spread; • definition of the medium type (biogeocenosis), • and creation of a physical model of the phenomenon. 2. Determination of the reaction and thermophysical properties of the medium, the transfer coefficients and structural parameters of the medium, and deduction of the basic system of equations with corresponding additional (boundary and initial) conditions. 3. Selection of a method of numerical solution of the problem, and derivation of differentialequations approximating the basic system of equations.

4. 4. Programming, test check of program; evaluation of the accuracy of the difference scheme; numerical solution of the system of differential equations. 5. Testing to see how well the derived results comply with the real phenomenon, their physical interpretation; development of new technical suggestions for ways of fighting forest fires.

5. Mathematical models of forest fires 1) models for initiation of forest fire; 2) models for predicting the rate of spread of fire and the contours of forest fires; 3) models for predicting the characteristics of flow, heat and mass transfer at the forest fires due to the natural and technogeneous catastrophes .

6. The basic assumptions adopted during the deduction of equations, and boundary and initial conditions: 1) the forest represents a multi-phase, multistoried, spatially heterogeneous medium; 2) in the fire zone the forest is a porous-dispersed, multi-phase, two-temperature, single-velocity, reactive medium; 3) the forest canopy is supposed to be non - deformed medium, affects only the magnitude of the force of resistance in the equation of conservation of momentum in the gas phase; 4) the energetics and physico-chemical processes at the forest fire front are took into account; 5) the length of the free path of a photon during a forest fire is much smaller than the characteristic dimensions of the forest biogeocenosis l* and , where l* is the horizontal dimension of the forest massif, and is the average height of the forest fuel layer;

7. 6) let there be a so-called “ventilated” forest massif, in which the volume of fractions of condensed forest fuel phases, consisting of dry organic matter, water in liquid state, solid pyrolysis products, and ash, can be neglected compared to the volume fraction of gas phase (components of air and gaseous pyrolysis products); 7) the flow has a developed turbulent nature and molecular transfer is neglected; 8) gaseous phase density doesn’t depend on the pressure because of the low velocities of the flow in comparison with the velocity of the sound; 9) diffusion approximation is applied to describe energy transfer by radiation.

8. Mathematical model of forest fires

12. Mathematical modeling of forest fire initiation (axisymmetrical formulation)

15. The initial and boundary conditions

16. The local – equilibrium model of turbulence.

17. Fig.1 a,b,c. Relationships of dimensionless temperatures , concentrations and volume fractions in the lower boundary of the forest canopy at r=0 :

18. Fig.2 a,b. Relationships between and on the axis of symmetry and dimensionless height .

19. The distribution of temperatures and velocities for axisymmetrical case.

20. I II I II I II Figure 4. Isotherms of gas (I) and condensed (II) phases at different instants of time t=2.0 s (a), t=2.2 s (b), t=2.4 s (c). Arabic numerals 1,2,3 corresponded to the isotherms 2., 3. and 5 respectively.

21. I II Figure 5. Lines of equal relative concentrations of oxygen and carbon monoxide at the same moments as in Figure 2. Arabic numerals correspond to -I: , II: .

27. Vectorial fields of velocity and distributions of temperatures

28. Mass concentrations of oxygenand combustible products of pyrolysis

31. Mathematical Modeling of Forest Fire Spread Initiation Because of the horizontal sizes of forest massif more than height of forest – h - average value of .

34. Figure 2. Isotherms of solid phase at t=5 s and 10 s: Figure1. Isotherms of the forest fire for t=5 s and 10 s:

35. Figure The distribution of oxygen at t=5 and 10 s: Figure The distribution of at t=5 and 10 s:

37. Mathematical Modeling of Large Forest Fire Initiation

40. 1 – x, y = 0; 2 – x = -10 km, y=0; 3 – x= - 15 km, y = 0. 1 – x, y = 0; 2 – x = -10 km, y=0; 3 – x= - 15 km, y = 0; .

41. 1 – x, y = 0; 2 – x = -10 km, y=0; 3 – x= - 15 km, y = 0.

42. 1 – t=7.0 sec, 2 – t=5 sec, 3 t=4.3 sec. 1 – t=7.0 sec, 2 – t=5 sec, 3 t=4.3 sec. 1 – t=7.0 sec, 2 – t=5 sec, 3 t=4.3 sec.

44. The distribution of basic functions in the region of forest ignition.