
книги / Методы математического моделирования рудничных аэрологических процессов и их численная реализация в аналитическом комплексе Аэросеть
..pdf
становится возможным определить мощность теплообмена, приходящуюся на длину одного участка
Wpart |
Lpart |
W , |
Lend Lstart Lrib |
где Lpart – длина одного участка в выработке, W – заданная
мощность источника.
Также есть возможность вычислить значение Lportion по границам Lpart _start и Lpart _end текущего участка:
|
1) если |
источник |
находится |
правее |
( Lpart _end Lstart ), |
то |
||||||||
Lportion 0; |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2) если |
источник |
находится |
|
левее |
(Lpart _start Lend ), |
то |
|||||||
Lportion 0; |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
3) если |
участок помещается |
внутри |
источника |
целиком |
|||||||||
( Lstart Lpart _ start и Lpart _end |
Lend ), тогда Lportion 1; |
|
|
|||||||||||
|
4) если |
источник |
выходит за |
границу участка |
справа |
|||||||||
(L |
L |
part _end |
), то L |
portion |
|
Lpart _end |
Lstart |
|
; |
|
|
|||
|
|
|
|
|
||||||||||
end |
|
|
|
Lpart _end |
Lpart _start |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
5) если |
источник |
выходит за границу участка слева, |
то |
||||||||||
Lportion |
Lend |
Lpart _start |
. |
|
|
|
|
|
|
|
||||
Lpart _end Lpart _ start |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
Приведенный точный способ учета пересечений участков с областями теплообмена в выработках позволяет обеспечить получение достоверных результатов даже при условии выделения очень малого числа участков.
341

При моделировании пожара задается его мощность и протяженность, после чего используется та же самая зависимость,
что и для калориферной установки, только в Lsource полагается равной не 2 м, а заданной протяженности пожара.
Для двигателя внутреннего сгорания мощность тепловыделения рассчитывается по формуле
W Pconsumption Ccombustion 1000,
3600
где Pconsumption – средний расход топлива в кг в час, Ccombustion –
удельная теплота сгорания топлива в МВт/кг. Причем, средний расход топлива рассчитывается на основе справочных данных
Pconsumption Wengine pconsumption firregularity ,
где Wengine – номинальная мощность двигателя внутреннего сго-
рания в кВт, pconsumption – удельный расход топлива в кг на кВт-ч, firregularity – безразмерный коэффициент, учитывающий неравно-
мерность работы двигателя. Протяженность источника теплообмена для двигателя внутреннего сгорания принимается равной 2 м.
Двигатель внутреннего сгорания фигурирует в модельной схеме газораспределения в качестве отрицательного источника газовыделения, если речь идет о содержании кислорода в воздухе.
342

По умолчанию двигатель никак не влияет на концентрацию газа в выработке, однако при включении соответствующей опции, интенсивность потребления кислорода будет рассчитываться по формуле
E |
|
Pconsumption moxygen |
, |
emission |
|
air 3600 |
|
|
|
где Eemission – интенсивность потребления кислорода двигателем внутреннего сгорания в м3/с, Pconsumption – средний расход топлива в кг/ч, air – плотность воздуха, moxygen – масса кислорода, необ-
ходимая для сжигания одного кг топлива (принимается равной 14,42 кг кислорода на 1 кг топлива).
Применительно к конвейеру вся потребляемая им мощность на приводе расходуется в основном в результате тепловыделения на разных участках. Так, по умолчанию считается, что сам привод тратит на нагрев 15 % своей номинальной мощности, а лента конвейера – 85 %. Скорректировать эти соотношения можно вручную в свойствах соответствующих объектов.
343

В остальном же тепловой расчет конвейера ничем не отличается от расчета простого источника тепловыделения.
Для источников газовыделения задается их протяженность и интенсивность. На каждом шаге по времени производится пересчет распределения концентрации газа по длине выработки
|
|
Cn V |
part |
L |
portion |
V |
partgas |
|
Cn 1 |
|
i |
|
|
, |
|||
i |
|
|
|
Vpart |
|
|
|
|
|
|
|
|
|
|
|
где Cin 1 – концентрация газа (от 0 до 1) на i-м участке после газообмена, Cin – концентрация газа (от 0 до 1) на i-м участке до газообмена, Vpart – объем i-го элементарного объема, Lportion – доля (от 0 до 1), которую занимает область газообмена на текущем участке, Vpart _ gas – объем выделившегося газа за время га-
зообмена.
Величина Lportion рассчитывается аналогично тому, как это
делается для источников тепловыделения, а объем выделившихся газов Vpart _ gas вычисляется по формуле
Vpart _ gas |
Eemission t Lpart |
, |
Lend Lstart Lrib |
где Eemission – общая интенсивность газовыделения, t – время газообмена, Lpart – длина одного участка, Lstart – относительное
расположение (от 0 до 1) начала области газовыделения в выработке, Lend – относительное расположение (от 0 до 1) окончания
области газовыделения в выработке, Lrib – длина выработки, рассчитывается так же, как и для источников выде-

По умолчанию все источники тепло- и газовыделения, расположенные на схеме, начинают свою работу в начальный момент времени и действуют на протяжении всего периода моделирования, о чем сообщает соответствующая надпись на вкладке «Расписания» в свойствах теплофизической модели.
При необходимости можно задать расписание работы каждого источника по отдельности, и в модели будут учитываться только те источники, у которых оно задано. После выбора источника на схеме расписание его работы добавляется в список.
345

По умолчанию поля с временем начала и окончания действия источника пусты, что означает активность источника на протяжении всего времени моделирования. Если указать только время начала, то источник будет активен с указанного момента модельного времени до окончания моделирования. Если есть необходимость, чтобы источник работал какое-то время, потом прекращал свою работу, а затем снова ее начинал и т.д., то допускается добавлять несколько расписаний для одного и того же источника. Возможность гибкой настройки времени действия источников тепло- и газовыделения позволяет моделировать, например, разные сценарии развития пожаров в руднике.
8.3.5.Источники тяги переменного напора
Впрограмме предусмотрена возможность задания переменных напоров источников тяги. Например, этот функционал востребован при моделировании реверсирования или остановки вентилятора главного проветривания. Функционал активируется кнопкой «Добавить расписание» на вкладке «Вентиляторы»
уформы параметров теплофизической модели с выбором нужного вентилятора на схеме.
346

Далее заполняется список напоров вентилятора с привязкой к моментам модельного времени. При этом для всех промежуточных моментов времени напор вентилятора интерполируется между указанными значениями. Схема изменения напора вентилятора отображается графически. Чтобы добавить новую точку с привязкой напора вентилятора к моменту модельного времени, нужно ввести значения в пустые поля внизу списка.
8.3.6. Выбор модели теплообмена
Если в теплофизической модели включена опция «Тепло-
обмен со стенками выработок», то температура воздуха на участках будет пересчитываться с учетом взаимодействия с окружающей выработку породой. Предусмотрены два основных способа расчета теплообмена: упрощенный с использованием коэффициента теплоотдачи стенок или более точный с использованием коэффициента нестационарного теплообмена. В первом способе на основе справочных значений задается коэффициент теплоотдачи и указывается температура стенок выработки. Во втором способе задается температура удаленной породы, а коэффициент нестационарного теплообмена рассчитывается.
347

8.3.6.1. Упрощенный способ моделирования теплообмена на основе использования коэффициента теплоотдачи
По умолчанию теплообмен со стенками рассчитывается на основе проектных данных с использованием коэффициента теплоотдачи стенок выработки, величина которого зависит от расхода воздуха по выработке
|
|
Q |
0.8 |
||
0 |
|
|
, |
||
|
|||||
|
Q0 |
|
|
где – коэффициент теплоотдачи у стенок выработки при расходе воздуха, равном Q , 0 – указанный пользователем коэффициент теплоотдачи при расходе, равном Q0 .
Температура воздуха на участках рассчитывается следующим образом:
T n |
|
|
air |
c |
V |
T n 1 |
S t T |
|
|
|
air |
air |
air |
wall |
, |
||
air |
|
|
|
air cair Vair |
S t |
|||
|
|
|
|
где Tairn – неизвестная температура воздуха на текущем участке на текущем шаге модельного времени, Tairn 1 – известная температура воздуха на текущем участке на предыдущем шаге модельного времени, air – плотность воздуха, cair – удельная теп-
348

лоемкость воздуха, Vair – объем воздуха на текущем участке, Twall – температура стенок выработки в том месте, где находится
текущий участок, S – площадь соприкосновения текущего элементарного объема со стенками выработки, t – время теплообмена со стенками выработки.
Температура стенок по умолчанию рассчитывается на основе температуры, заданной в узлах выработок («Вид» -> «Тепло-
физика» -> «Исходные данные» -> «Конечные вершины (индикаторы)» -> «Фактическая температура стенок выработок»).
Twall i |
Twall |
|
n i Twall |
i 1 |
|
|
from |
to |
, |
||
|
|
|
|
||
|
|
|
n 1 |
|
|
где Twall i – температура стенок выработки рядом с i-м участком, Twall from – температура стенок в начальном узле выработки, Twallto – температура стенок в конечном узле выработки, n – общее число участков в выработке.
8.3.6.2.Модель нестационарного теплообмена
сучетом распространения тепла в породном массиве
При теплообмене в течение длительных промежутков времени температура породного массива меняется, поэтому моделирование на основе задания постоянной температуры стенок выработок является некорректным. Проблема решается применением более детализированной модели нестационарного теплообмена, учитывающей распространение тепла как в воздухе,
349

так и в массиве. Выбор модели осуществляется в окне «Тепло-
обмен со стенками»:
Результат моделирования теплообмена воздуха с породой можно просмотреть при помощи соответствующей кнопки внизу окна.
По умолчанию считается, что вычисленная или заданная температура стенок при дальнейшем теплофизическом моделировании в сетевой постановке не меняется. Это допустимо для небольших интервалов времени порядка нескольких часов, однако при увеличении таких интервалов до нескольких суток и месяцев такое приближение становится грубым. Для увеличения точности необходимо перейти на динамический тип модели:
350