Математика от А до Я: Справочное пособие (издание третье с дополнениями)
Шрифт:
Как известно [121], любая консервативная субстанция, смешивающаяся с движущейся жидкостью, переносится относительно системы координат, связанной со средним ее движением, путем турбулентного и молекулярного обмена. Общий поток массы выражается в виде
В соответствии с законом Фика
где kv —
Уравнение сохранения консервативной пассивной примеси в предположении постоянства pv и kv по пространству имеет вид [121]
где q — концентрация примеси.
Отметим, что это уравнение имеет весьма общий вид, и им можно пользоваться для определения изменения любой консервативной и пассивной примеси или любого свойства воздуха, заменив q на концентрацию, выраженную отношением массы примеси к единичному объему общей массы воздуха и понимая под ки коэффициент молекулярной диффузии этой примеси.
Записанное выше уравнение диффузии можно решить, выбрав подходящие граничные условия и зная распределение поля скорости. Граничные условия задаются трех типов, а именно на поверхности z = 0 задается либо значение q, либо поток рассматриваемой примеси, либо поток примеси выражается через другие компоненты теплового баланса.
К сожалению, уравнение (4.26) не находит непосредственного применения в практических задачах, так как реальные потоки имеют турбулентный характер. Это означает, что в действительности невозможно определить скорость переноса и концентрации примесей в любой заданной точке пространства и времени, а можно найти только их статистические характеристики.
Для этого рассматривают осредненные величины, и в соответствии с общепринятым подходом, предложенным Рейнольдсом, зависимые переменные представляют в виде сумм не возмущенных величин и возмущений:
Применяя затем обычный метод осреднения по времени с соответствующим периодом осреднения и используя уравнение неразрывности, из уравнения (4.26) получается соотношение для нахождения
Члены в левой части этого уравнения представляют скорость изменения средней массовой доли вещества примеси, перемещающейся с осредненной скоростью движения воздуха. Ковариации пульсаций в правой части уравнения можно назвать турбулентными потоками по аналогии с напряжениями Рейнольдса. Они являются компонентами диффузионного потока, обусловленного турбулентным движением. Последний член представляет перенос средней субстанции за счет молекулярной диффузии.
Это уравнение должно быть дополнено уравнениями неразрывности, количества движения и энергии в терминах средней скорости движения несжимаемой жидкости. Вид этих уравнений для пассивных и консервативных примесей общепринятый и поэтому здесь не приводится.
Записанные таким образом уравнения сохранения имеют незамкнутый вид, и поэтому их решение представляет большую проблему. Уравнения для моментов низших порядков (для осредненных величин) содержат потоки, обусловленные пульсациями метеорологических элементов и содержат моменты более высокого порядка. Таким образом, любой конечный набор уравнений для моментов турбулентных флуктуаций всегда включает
Упрощение этой проблемы достигается несколькими подходами. Во-первых, путем выделения в атмосфере вблизи подстилающей поверхности особого пограничного слоя, в котором вертикальные градиенты значительно больше горизонтальных. Во-вторых, путем использования принципов подобия и полуэмпирической теории турбулентности, выражая моменты второго и более высоких порядков через осредненные переменные и моменты более низких порядков.
В качестве примера использования инженерного подхода для решения задачи распространения консервативных пассивных примесей в атмосфере приведем математическую модель атмосферной диффузии примеси при тумане. Расчет распространения примеси от источников при тумане основывается на решении уравнения турбулентной диффузии, записанном в виде [129]:
Здесь u — скорость ветра; q — концентрация примеси; Ку — горизонтальная составляющая коэффициента обмена; а' — показатель степени поглощения примеси водяными каплями (вне тумана а' — 0).
Начальным условием при х=0 принимается наличие источника на некотором уровне z = Нmр при у = 0 и в качестве граничных условий, как обычно, убывание q до нуля при неограниченном удалении от источника и отсутствие потока примеси на подстилающей поверхности, т. е. при z = 0:
Способ определения а' рассмотрен в работе [122]. Для расчета концентрации примеси q необходимо знать распределение водности в тумане и высоту тумана. Решение приводится в работе [122]. Из анализа решения следует, что поглощение примеси, содержащейся в газообразном виде в атмосфере, происходит в основном в верхнем слое тумана; вблизи земли ее концентрация близка к нулю. Причем на расстояниях х > 0,5 км от источника практически вся газообразная примесь в тумане растворена в каплях.
4.6. Взрывной разлет твердых и жидких частиц
Жидкие и твердофазные выбросы являются важной загрязняющей компонентой при авариях на промышленных объектах. Сносящий ветровой поток приводит к переносу частиц на большие удаления от места аварии и загрязнению обширных ареалов.
Исследованию процессов разлета частиц и фрагментов взрываемых объектов разного размера, а также изучению загрязнения атмосферного воздуха и поверхности земли твердофазными и жидкими продуктами взрыва посвящено большое количество работ [71–72,74-85], основная часть которых описывает возникновение и разлет частиц и осколков при взрывах емкостей, снаряженных газами и конденсированными твердыми топливами.
Авторы большинства работ ограничиваются рассмотрением движения массивных тел по баллистическим траекториям в пренебрежении воздействия ветра. Согласно упрощенному анализу [74] движение тела предполагается в одной плоскости, причем допускается, что оно может вращаться вокруг продольной оси, что придает осколку необходимую устойчивость и позволяет считать, что тело не сносится ветром.
В действительности фрагменты разрушенного объекта, жидкие и твердые частицы при их взрывном разлете в ветровом потоке заметно отклоняются от первоначальной плоскости. Причем эти отклонения тем больше, чем мельче частицы. Пространственный характер движения частиц при наличии возмущающего воздействия внешних сил может быть учтен в предположении независимости их движения в горизонтальной плоскости и в плоскости разлета [87,76].