книги / Методы математического моделирования рудничных аэрологических процессов и их численная реализация в аналитическом комплексе Аэросеть
..pdfРасстояния до центров ближайших участков:
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