Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы математического моделирования рудничных аэрологических процессов и их численная реализация в аналитическом комплексе Аэросеть

..pdf
Скачиваний:
4
Добавлен:
19.11.2023
Размер:
53.4 Mб
Скачать

Расстояния до центров ближайших участков:

Sold _left li _center

jold _left

loldcenter

,

Sold _right jold _right

lold

center

li

.

 

 

center

 

Значение теплофизического параметра Tnew в следующий

момент времени для участка с индексом i вычисляется по следующей формуле

 

T

j

 

 

 

T

j

 

 

old oldleft

 

 

old oldright

 

Tnew i

 

Soldleft

 

 

 

Soldright

 

.

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Soldleft

 

Soldright

 

 

 

 

 

 

 

 

 

 

 

Подобный метод интерполяции массива теплофизических параметров позволяет добиться непрерывности изменения этих параметров при резком увеличении числа участков, однако несколько замедляет расчетный процесс. Поэтому если расходы в выработках меняются незначительно, то интерполяцию данных на участках целесообразно отключать.

В процессе расчета теплофизические параметры воздуха меняются под действием различных факторов, и на каждом шаге по времени участки сдвигаются на одну позицию по направлению течения воздуха. В местах слияния воздушных потоков теплофизические параметры на участках пересчитываются с учетом перемешивания воздуха.

331

На приведенной схеме корректировка параметров с учетом перемешивания воздуха будет выполняться на участке, закрашенном синим цветом, по формуле

T

Qinput T

i

i ,

mixed

Qinput

 

 

i

 

где Tmixed – значение температуры, абсолютной влажности воздуха или концентрации газа после перемешивания воздушных потоков, Qiinput – расход воздуха во входящем в узел i-м потоке, Ti – значение соответствующего параметра во входящем в узел

i-м потоке.

По умолчанию на каждом шаге по времени рассчитывается стационарное воздухораспределение. Однако этот способ моделирования является нерациональным в случае отсутствия причин для изменения расходов в выработках. После отключения этой опции расчет стационарного воздухораспределения будет проведен только один раз в начале, а далее на всех остальных шагах по времени будут использоваться вычисленные на первом шаге модельные расходы.

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

332

Вместе с пересчетом воздухораспределения деактивируют-

ся также и опции «Учитывать естественную тягу» и «Инерци-

онность воздуха», так как теряют свой смысл в случае постоянных расходов.

Учет естественной тяги инициирует корректировку расходов в выработках с учетом меняющихся тепловых депрессий. В отличие от стационарного расчета воздухораспределения, в теплофизическом расчете температура воздуха задается не в узлах выработки и равномерно по ней распределяется, а вычисляется в каждом элементарном объеме воздуха внутри выработки.

Учет инерционности воздуха означает, что при резком изменении перепада давлений расход изменяется не мгновенно по причине наличия у воздуха массы. Например, при реверсировании вентилятора к напору в каждой выработке добавляется инерционная поправка

P

 

L Qcurrent Qprevious

,

 

inertia

 

g S t

 

 

 

 

где – плотность воздуха, L – длина выработки, Qcurrent – расход воздуха в выработке на текущем шаге моделирования, Qprevious – расход воздуха в выработке на предыдущем шаге мо-

делирования, g – ускорение свободного падения, S – площадь

поперечного сечения выработки, t – шаг модельного времени. В результате расходы приходят к своим модельным значениям с некоторым запаздыванием.

333

Для каждого элементарного объема воздуха в выработках отдельно выбираются те теплофизические параметры, которые будут отслеживаться.

Отключение любой из этих опций высвобождает память и процессорное время.

При включении учета температуры воздуха в выработках становятся доступны дополнительные опции. Первая из них включает или отключает изменение температуры воздуха в выработках в результате теплообмена с источниками тепловыделения, установленными на схеме. Данная опция используется, например, для учета влияния пожаров, двигателей внутреннего сгорания, кондиционеров и пр. Вторая указывает, нужно ли моделировать теплообмен со стенками выработок, третья задает необходимость учета нагревания воздуха вследствие его гидростатического сжатия. При этом температура корректируется аналогично тому, как это делается в модели нестационарного теплообмена с породой за исключением того, что высотная отметка задается на панели свойств узла на схеме.

334

Что касается учета источников влаговыделения, то при их инициации влажность воздуха будет меняться, в том числе в результате влагообмена со стенками выработок. Учет же источников газовыделения включает или отключает не только газообмен с соответствующими источниками, но и собственно учет концентрации газа в выработках.

Предусмотрены два способа инициализации теплофизических параметров в выработках в начале моделирования. Первый способ связан с использованием параметров, указанных на схеме.

335

В частности, температура воздуха в элементарных объемах задается исходя из температуры воздуха в конечных узлах выработки.

При этом используется следующая зависимость

T 0

 

Tfrom n i Tto i 1

 

i

 

n 1

 

 

где Ti0 – начальная температура воздуха в i-м участке выработки, n – число участков, выделенных в выработке, Tfrom – темпе-

ратура воздуха в начальном узле выработки, Tto – температура

воздуха в конечном узле выработки.

Влажность воздуха инициализируется на основе значения, указанного в качестве влажности атмосферного воздуха. Концентрация газа в начальный момент времени считается равной нулю по умолчанию, но это значение может быть отредактировано в параметрах атмосферы.

336

Второй способ инициализации теплофизических параметров связан с использованием уже полученных результатов теплофизического моделирования. В этом случае в качестве начальных берутся параметры участков в конечный момент времени для выбранной модели-предшественницы. Причем в качестве такой модели может выступать и текущая модель, если для нее уже был выполнен теплофизический расчет.

8.3.3. Параметры атмосферного воздуха

Если в свойствах узла выработки указано, что он имеет связь с атмосферой, то в теплофизической модели на участке, граничащем с этим узлом в сторону направления движения воздуха, используются теплофизические параметры атмосферы. Концентрация газа в атмосфере по умолчанию считается равной нулю, температура и влажность воздуха задаются постоянными величинами. Такое приближение допустимо для небольших интервалов модельного времени.

337

Если же моделирование проводится для интервалов порядка суток, месяцев и более, то игнорировать суточные и годовые колебания параметров атмосферы становится невозможным. В первом случае задаются суточные колебания температуры и влажности воздухадлямоментоввремени синтерваломвдвачаса.

Во втором случае указываются усредненные параметры атмосферы для каждого из двенадцати месяцев в году.

338

8.3.4. Источники выделения тепла, влаги и газа на схеме

Включение опции «Источники тепловыделения» в пара-

метрах теплофизической модели позволяет моделировать изменение температуры воздуха в выработках вследствие теплообмена с источниками тепловыделения на схеме, а включение опции «Источники газовыделения» – изменение концентрации газа. При этом все такие источники должны быть предварительно размещены на схеме (категория «Тепло- и газовыделение» в галерее объектов на выработках).

Все источники тепловыделения делятся на две категории: обобщенные источники, в которых параметры задаются на основе простых формул, и специальные, в которых итоговые свойства выводятся на основе проектных параметров источников. В категорию простых источников тепловыделения попадают калориферная и холодильная установки, у которых мощность нагревания и остывания проходящего через них воздуха указывается непосредственно.

Для пересчета температуры воздуха на участках, граничащих с калориферной или холодильной установкой, используется следующая зависимость

T n 1

T n L

 

 

Wpart t

,

i

i

portion

 

c S Lpart

 

где Ti n 1 – температура на i-м участке после теплообмена с источником тепловыделения,Ti n – температура на i-м участке до

339

теплообмена с источником тепловыделения, Lportion – доля (от 0

до 1), которую занимает область теплообмена с источником на текущем участке, Wpart – мощность нагревания (остывания) воз-

духа у источника, приходящаяся на длину одного участка, t – время теплообмена, c – удельная теплоемкость воздуха, –

плотность воздуха, S – площадь поперечного сечения выработки, Lpart – длина одного участка.

Для вычисления Lportion сначала определяются границы расположения источника тепловыделения в выработке Lstart и Lend :

1)по умолчанию границы источника указываются совпадающими с границами выработки: Lstart 0, Lend 1;

2)если протяженность источника меньше длины выработ-

ки, то вычисляется местоположение центра источника

Lcenter :

Lcenter Lrelative

Lrib ,

Lrelative

относительное местоположение ис-

точника вдоль длины выработки (от 0 до 1),

Lrib

– длина выра-

ботки;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3) если

 

 

источник

не

помещается

в

выработку

слева

( L

 

 

 

 

Lsource

 

0,

 

 

L

 

 

протяженность

источника):

 

 

 

 

 

 

 

 

 

center

2

 

 

 

 

 

 

 

source

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

Lsource

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

end

 

 

Lrib

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4) если

 

источник

не помещается в выработку справа

( L

 

 

 

Lsource

 

L

):

L

 

1

Lsource

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

center

2

 

 

 

 

rib

 

 

start

 

 

Lrib

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5) если

 

 

источник

помещается в

выработку целиком

L

L

 

Lsource

,

L

 

L

 

Lsource

;

 

 

 

 

 

2 L

 

 

 

start

 

 

relative

 

 

 

2 L

 

 

end

 

relative

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rib

 

 

 

 

 

 

 

rib

 

 

 

 

6) если протяженность источника больше длины выработки, то границы источника обрезаются по размерам выработки.

Протяженность источника Lsource по умолчанию считается равной 2 м. После того как установлены границы Lstart и Lend ,

340