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

Ошибка Floating point exception

combo.o
Аватар пользователя combo.o

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

ВложениеРазмер
Иконка изображения error.jpg81.92 КБ
Кирилл
Аватар пользователя Кирилл

Добрый день! Какая в данном случае сетка, граничные условия, методы расчета?

combo.o
Аватар пользователя combo.o

Качество сетки хорошее, искривленных элементов нет, вырожденных тоже.

Задачу решаю в 2D, осесимметричный случай. Решатель по плотности, расчитываю нестационар. В качестве модели турбулентности SST модель.

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

По граничным условиям - на входе pressure inlet, на выходе - pressure outlet, wall - поверхность тела, стенки форкамеры, сопло и стенки рабочей части.

Кирилл
Аватар пользователя Кирилл

Это хорошо. Но сетку все равно бы посмотреть. Погранслой разрешается? Сколько элементов получается поперек погранслоя? Должно быть как минимум десять.

 

А если попробовать ламинарный или даже невязкий расчет?

 

Попробуйте неявную схему - она поустойчивее. 

combo.o
Аватар пользователя combo.o

Так, сетку все же проверил - пограничный слой состоит из двенадцати элементов, устранил искривленные элементы, повысил качество.

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

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

Проблемы начинаются, когда поток воздуха переходит за критическое сечение сопла - скорость потока возрастает до 70 Махов и выше. Никакого физического объяснения для ТАКОГО роста не могу найти.

Делюсь графиками - невязки и Сх.

Кирилл
Аватар пользователя Кирилл

Покажите картинку сетки. Кто знает - может там проблема. Кстати, на стенке задана фиксированная температура или что-то ещё?

Физического объяснения тут и нет. Косяк явно численный

Добавить комментарий

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

Максимальный размер файла: 128 МБ.
Допустимые типы файлов: txt doc docx xls xlsx pdf rar zip 7zip tar.