Моделювання процесу розповсюдження лісової пожежі
Розглянуто вирішення важливої задачі визначення контуру низинної лісової пожежі з огляду на здійснення засобів боротьби із лісовими пожежами. Розглянуто числово-аналітичний метод визначення температурного поля у зоні пожежі. Запропоновано числово-аналітичний метод розв’язання крайової задачі із...
Збережено в:
| Опубліковано в: : | Збірник наукових праць Інституту проблем моделювання в енергетиці ім.Г.Є.Пухова НАН України |
|---|---|
| Дата: | 2009 |
| Автор: | |
| Формат: | Стаття |
| Мова: | Українська |
| Опубліковано: |
Інститут проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України
2009
|
| Онлайн доступ: | https://nasplib.isofts.kiev.ua/handle/123456789/27768 |
| Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Цитувати: | Моделювання процесу розповсюдження лісової пожежі / В.О. Лщина // Збірник наукових праць Інституту проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України. — К.: ІПМЕ ім. Г.Є. Пухова НАН України, 2009. — Вип. 50. — С. 166-170. — Бібліогр.: 2 назв. — укр. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1859720767045894144 |
|---|---|
| author | Лщина, В.О. |
| author_facet | Лщина, В.О. |
| citation_txt | Моделювання процесу розповсюдження лісової пожежі / В.О. Лщина // Збірник наукових праць Інституту проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України. — К.: ІПМЕ ім. Г.Є. Пухова НАН України, 2009. — Вип. 50. — С. 166-170. — Бібліогр.: 2 назв. — укр. |
| collection | DSpace DC |
| container_title | Збірник наукових праць Інституту проблем моделювання в енергетиці ім.Г.Є.Пухова НАН України |
| description | Розглянуто вирішення важливої задачі визначення контуру низинної
лісової пожежі з огляду на здійснення засобів боротьби із лісовими пожежами.
Розглянуто числово-аналітичний метод визначення температурного поля у зоні
пожежі. Запропоновано числово-аналітичний метод розв’язання крайової задачі
із рухомими межами для визначення фронту низинної пожежі.
The decision of important task of determination the contour of low-laying area forest fire is considered, taking into account realization of facilities of fight against forest fires. It is considered numerical-analytical method of determination of the temperature field in the area of fire. It is offered numerical-analytical method of decision of regional task with mobile limits for determination of front of low-laying area fire.
|
| first_indexed | 2025-12-01T09:44:31Z |
| format | Article |
| fulltext |
УДК 621.397.3
В.О. Ліщина, Луцький інститут розвитку людини ВМУРоЛ "Україна",
м. Луцьк
МОДЕЛЮВАННЯ ПРОЦЕСУ РОЗПОВСЮДЖЕННЯ ЛІСОВОЇ
ПОЖЕЖІ
Розглянуто вирішення важливої задачі визначення контуру низинної
лісової пожежі з огляду на здійснення засобів боротьби із лісовими пожежами.
Розглянуто числово-аналітичний метод визначення температурного поля у зоні
пожежі. Запропоновано числово-аналітичний метод розв'язання крайової задачі
із рухомими межами для визначення фронту низинної пожежі.
The decision of important task of determination the contour of low-laying area
forest fire is considered, taking into account realization of facilities of fight against
forest fires. It is considered numerical-analytical method of determination of the
temperature field in the area of fire. It is offered numerical-analytical method of
decision of regional task with mobile limits for determination of front of low-laying
area fire.
Вступ
Дослідження процесу розповсюдження лісової пожежі ставить за мету
визначення контуру лісової пожежі з метою організації боротьби з лісовими
пожежами. Контур лісової пожежі можна визначити шляхом розв'язання
загальної системи рівнянь теплового балансу у зоні пожежі [1]. Замість цієї
складної системи, з якої визначаються тривимірні поля швидкостей,
температур та концентрацій, можна використовувати простішу систему, яку
отримують осередненням за висотою шару рослинності.
Постановка задачі
Контур пожежі у кожний момент часу можна розглядати як лінію рівня
на площині (х, y) або як поверхню у просторі (x, y, t), яку можна задати у
вигляді (р(х, y, t) = 0 або у явному вигляді y = f (x, t). Якщо визначити
кромку пожежі як ізотерму, що відповідає температурі горіння шару лісового
займистого матеріалу (ЛЗМ), то для опису лінії контуру досить рівняння
теплового балансу на площині x, y :
cpmpmpm { j t f + Vmgrad Tm j = div ( ( (1 - Y)grad Tm ) + Q , (1)
де Qm - тепловий потік, що виділяється при горінні ЛЗМ, cm, pm ,Хт -
теплофізичні характеристики ЛЗМ.
Основна частина
Для пожеж, розмір яких у плані набагато ширше за ширину горіння,
величиною теплопровідності твердого палива À.m можна знехтувати. Тоді
166 © В.О. Ліщина
рівняння (1) набуває вигляду гіперболічного рівняння
дТ
Срт Рвт + ^^гай Тт j = ат . (2)
Початкові умови для цього рівняння мають задаватися у вигляді
т< т Т ПРи Х,У Є Во,
т
( Х , У , 0 ) = І Т „ В
[Тр при Х,У <£ Во.
де Тр с - температура горіння шару; В0 с - задана область на площині
(х, у) , межа якої С являє собою кромку пожежі у початковий момент часу.
Якщо р(х, у, ґ) = 0 є рівняння контуру, функція р має задовольняти
умову Гамільтона—Якобі для рівняння (2) :
+ vgraй р = 0 . (3)
дґ
Це рівняння із початковою умовою р( х, у,0) = р0( х, у) розв'язується за
методом характеристик, що і визначає шуканий контур розповсюдження
лісової пожежі. У разі явного завдання контуру у = /(х, ґ),
р(х, у, ґ) = /(х, ґ) - у рівняння контуру набуває вигляду
— + V.. 1 + № I2 = 0
дґ \ {дх)
з початковою умовою у( х,0) = / (х,0) = / 0 ( х).
Виконано моделювання цих рівнянь з метою визначення контурів для
випадків, коли початкові умови мають вигляд кола і еліпса.
У роботі також розглядається загальніший підхід, коли контур має
складний вигляд, тобто (р0 - початковий контур осередку пожежі Г 0 . Цей
контур розбивається на елементарні відрізки в околі точок М і 0 . Для
визначення контуру лісової пожежі запишемо систему рівнянь відносно
температури середовища перед фронтом і за фронтом пожежі:
С1 = ^ ^ + ат( - Т ) х, < 1 (4) дґ 1 дх, 4 дх2 к
дТ2, дТ2 і „ д 2Т2, а
+ V, ^ = + к ( - Т 0 ) , х > х* 1 . (5) - 2 1 дґ 1 дх1 2 1 дх2 к
У цій системі рівнянь параметри с1}, с2}, У1 } , У2} ,Я1/ ,Л2і вважаються
відомими (у загальному випадку вони визначаються як розв'язання рівнянь
гідродинаміки).
Початкова умова: Т (х і , 0) = Т п . Межові умови задаються у вигляді
167
Т(0,і) = Тс; т\х
К
дТ,
дх. + а-
йх„ у
йі
•=К
= Т.;
)
дТ
дх..
(6)
(7)
У рівняннях (4)-(6) Т п = т на контурі Г 0 ; - координати фронту
лісової пожежі, що вираховуються від точок Мо по нормалі до контуру Г 0 .
Задача (4)-(7) - задача із рухомим межами, відома як задача типу
Стефана. Розв'язання цієї задачі дозволяє визначити х„Ді) у будь-який
момент часу і, отже, визначити контур Г шляхом інтерполяції.
З метою спрощення запису системи рівнянь та межових умов опустимо
індекс ] і запровадимо заміну змінних:
Т1(х,і) = е^/(2Л)хМі6>(х,і); Т2(х,і) = ^ /(2К)х+^6>(х,і); ёк = V2/(4\ск) +ак /к . (8)
Отримуємо крайову задачу із змінними межами у вигляді
1 ді
х < 4 ;
д61 3 дГВХ а
= К " к Т 0 ,
д̂2 з д2В2 «2 Т <4
дґ 1 д г
2 к
Початкові умови:
6>(х,0) = е"/(2^)хТт; 0 2 ^ , 0 ) = е ^ / ( 2^ гТ т ;
Межові умови:
ВХ (0, ґ) = е"Л1% : В2 (0, ґ) = е"Лг'Тс;
6>(4,ґ) = с>1,Тт : 02(4,ґ) = е " ^ 2 ^ 4 - ^ Т т .
Умови сполучення на межі розділу фаз набувають вигляд
"є V /(2К)4 д6> у_
дх 2Л1
'=4
= а—— + Ке
йі 2
• / (2К)4 дв2 V — 2
ді 23
2=4
Власні функції для крайових задач мають вигляд:
2 2
X ( х ) = 4 І п Р п х ; Zn
(г) = - в і п р п і .
_ 4 _ 4
Ч.Рп,ґ) = х,ґ)Хп(х)йх; в2(Рп,ґ) = \в2(2,ґ)Іп(2)й2 .
0 0
Рівняння (9), (10) набувають вигляд:
^ ^ + Д & = Я ( 4 ) - О Н ; к = 1,2; " йі кК
Як (4) = е^к і [е±v /(2Кк)4Т0 - Тт ] ;
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
168
— — — 2
й (в ,0) = Тс; Тс = Г-8ІП рпхв± /(2Ак )х*х. (17)
0
І
Із урахуванням виразів (16),(17) розв'язання рівняння (15) запишеться як
ґ л
й в , ') = е
- а к п' — а
Т +——Т
С і с
На,_
А
На,
Тс + 1 ( - е-'к')/(2-*Т0 - Тт ] . (18)
*к - ак,п
А У цьому виразі акп =-к- вп.
' ск
Розв'язання крайової задачі (9)-(13) має вигляд:
й ( х , ' ) = Х|5ІП
рпхв1(вп,'); вг(2,') = х | 5 і п Рп2в2(Рп,'). (19)
п 7 п 7
Ці розв'язання мають задовольняти умову сполучення (14). Маємо:
дй . 2 дй - , 2
- 1 = Х в І Рпхй(Рп,'); --і = Х в п І Р Й 2 , ' ) .
дх п І & п І
Підстановка цих виразів в умову сполучення (14) призводить до
диференціального рівняння відносно межі розділу І ( ' ) фази горіння Т1(І,') і
фази гасіння пожежі Т2 (І , ' ) .
^ =1 -Це/(2-)ІУ -
* Я [ п І
- 1 А
Я I и д
в „ « и в п і - З К 5 І П в І
РП c o s в п І - в п І
й(вп , ' ) -
й2(вп , ' ) . (20)
Це рівняння перепишемо у такій формі.
* ( І 2 / 2 )
= 1 АХ
Я п
в п е - / ( 2 А ) І С ° 8 в п І - 2 - е " / ( 2 А ) І 8 І П в І
- 1 А Х
Я п
РПЄ
' / (2-2 ) І V Л З Ю І С І
2 -
®ІП в п І
й(вп, І) -
Й(вп , ' ) = ^ ( І , ' ) . (21)
У виразі (17)
Т ( І ) = ^^^^ [~1 - ( -1 ) п - 1 е ± У / ( 2 - ) І ] .
І 2 ( V / 2 А к ) 2 + п 2 П 1 ]
(22)
таке:
З урахуванням того, що вп = пп/І, замість рівняння (21) отримаємо
^ ^ ^ = 1 Х в п (-1)п-1 [ - е " " / ( 2 А ) І Й ( Р П , ' ) - - 2 е" / ( 2 -2 ) Ій2(вп , ' ) ] .
Я п = 1
Інтегрування цього рівняння дозволяє визначити закон зміни у часі межі
169
розділу фаз.
2 Ь = 1 } X П ^ Г 1 [ Л ^ Ч в , 0 - Л е ^ Ч в , /) ] л . (23)
2 Ч 0 и = 1 ь
Оскільки пряме інтегрування рівняння (21) утруднюється за наявності
функції Ь у знаменнику виразу (22), доцільно це рівняння апроксимувати
різницевим рівнянням із кроком т :
Ь+1 = Ь + т ^ ь , ^).
Розв'язання (19) та (23) отримують для кожної із ділянок контуру Г 0 .
Далі можна виконати інтерполяцію отриманих частинних контурів.
1. Доррер Г.А. Математические модели динамики лесных пожаров.—М.: Лесн. пром-
сть, 1979.—161 с.
2. Зеленський К.Х., Ігнатенко В.М., Коц О.П. Комп'ютерні методи прикладної
математики. - К.: Академперіодика, 2002.—480 с.
Поступила 19.01.2009р.
УДК 621.3
В.О. Пелішок , к.т.н., доцент, НУ «Львівська Політехніка» каф. ТК
ВИЗНАЧЕННЯ КОЕФІЦІЄНТА НАПРАВЛЕНОЇ ДІЇ
НЕПЕРЕРВНИХ АНТЕННИХ РЕШІТОК
The purpose of this article is to show that existing methods of antennas
characteristic calculation provide some inaccuracies. Thus, was designed new method,
which provides more precise calculation results. This method and example of
calculation are shown in this article.
В даній роботі показано, що загальноприйняті способи визначення КНД
неперервних лінійних АР забезпечують занижені, або завищені значення
реального КНД множника АР. Тому запропоновано спосіб визначення реального
значення КНД множника досліджуваної неперервної АР заданої довжини Ln.
Приведено приклади визначення реального значення КНД.
1. Відмінність між КНД простих антен та антенних решіток
В табл. 1 приведені значення коефіцієнта направленої дії (КНД) деяких
простих випромінювачів та антенн.
Варто зауважити, що на практиці часто виникає потреба в значно
більших значеннях КНД (10 та більше). Але для простих антен значення
КНД, приведені в табл.1, являються максимальними. Зовсім інша ситуація
виникає при застосуванні антенних решіток (АР). В них немає принципових
обмежень для забезпечення будь-якого значення КНД. Справа в тому, що для
170 © В.О. Пелішок
|
| id | nasplib_isofts_kiev_ua-123456789-27768 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | XXXX-0067 |
| language | Ukrainian |
| last_indexed | 2025-12-01T09:44:31Z |
| publishDate | 2009 |
| publisher | Інститут проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України |
| record_format | dspace |
| spelling | Лщина, В.О. 2011-10-17T15:07:38Z 2011-10-17T15:07:38Z 2009 Моделювання процесу розповсюдження лісової пожежі / В.О. Лщина // Збірник наукових праць Інституту проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України. — К.: ІПМЕ ім. Г.Є. Пухова НАН України, 2009. — Вип. 50. — С. 166-170. — Бібліогр.: 2 назв. — укр. XXXX-0067 https://nasplib.isofts.kiev.ua/handle/123456789/27768 621.397.3 Розглянуто вирішення важливої задачі визначення контуру низинної лісової пожежі з огляду на здійснення засобів боротьби із лісовими пожежами. Розглянуто числово-аналітичний метод визначення температурного поля у зоні пожежі. Запропоновано числово-аналітичний метод розв’язання крайової задачі із рухомими межами для визначення фронту низинної пожежі. The decision of important task of determination the contour of low-laying area forest fire is considered, taking into account realization of facilities of fight against forest fires. It is considered numerical-analytical method of determination of the temperature field in the area of fire. It is offered numerical-analytical method of decision of regional task with mobile limits for determination of front of low-laying area fire. uk Інститут проблем моделювання в енергетиці ім. Г.Є. Пухова НАН України Збірник наукових праць Інституту проблем моделювання в енергетиці ім.Г.Є.Пухова НАН України Моделювання процесу розповсюдження лісової пожежі Article published earlier |
| spellingShingle | Моделювання процесу розповсюдження лісової пожежі Лщина, В.О. |
| title | Моделювання процесу розповсюдження лісової пожежі |
| title_full | Моделювання процесу розповсюдження лісової пожежі |
| title_fullStr | Моделювання процесу розповсюдження лісової пожежі |
| title_full_unstemmed | Моделювання процесу розповсюдження лісової пожежі |
| title_short | Моделювання процесу розповсюдження лісової пожежі |
| title_sort | моделювання процесу розповсюдження лісової пожежі |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/27768 |
| work_keys_str_mv | AT lŝinavo modelûvannâprocesurozpovsûdžennâlísovoípožeží |