Вадим Романов - Прикладные аспекты аварийных выбросов в атмосферу

Поскольку истечение из ракетных сопел происходит с большими скоростями, то в первом приближении может быть оправданным подход при котором считается формирование полусферического выброса происходящим за первый шаг интегрирования задачи. Далее выброс представляется суммой полусферы и удлиняющегося усеченного конуса (Рис. 3.4).

Для определения координаты центра масс полусферического выброса х *радиуса R = d 0(Рис. 3.4а) приравняем массы газа в части выброса при х ≤ х *массе газа в части выброса при х > х *.

Получаем:

В этом выражении:

— радиус сопла;

Рис. 3.4. Схема формирования кратковременного выброса при истечении газа из сопла: а) переходный процесс возникновения выброса в окрестности сопла; б) развитый самоподобный выброс.

— уравнение образующей полусферической поверхности выброса;

ρ 1и ρ 2— плотности газа в левой (х ≤ х *) и правой (х > х *) части выброса, соответственно.

Если предположить, что вещество выброса имеет одинаковую плотность в разных его частях, т. е. ρ 1= ρ 2, то приходим к уравнению относительно искомой координаты центра массы х *. Получаем:

В уравнении (3.21): — безразмерная продольная координата центра масс.

Решением уравнения (3.21) является

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

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

Как следует из него в предложении однородности вещества выброса объем усеченной части выброса до координаты х„должен быть равен сумме объемов остальной части выброса.

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

Координаты сопряжения х сконической части выброса со сферической определяется приравнивания объемов этих частей выброса.

Получаем:

где

у 1= кх — уравнение образующей конической поверхности выброса;

— управление поверхности сферической его части.

Подставив значения у 1и у 2в это соотношение, получаем:

Вещественный корень этого уравнения может быть определен по формуле Кардана [172]:

Окончательное выражение для безразмерной продольной координаты сопряжения конической и сферической частей выброса может быть получено при подстановке в соотношение (3.24) вместо р и q их значений. Из-за громоздкости мы его не приводим.

Если известен радиус полусферической «шапки» выброса R, то выражение для продольной координаты сопряжения может быть записано в виде компактного соотношения. Приравниваем объем цилиндрической части выброса