ПОИСК Статьи Чертежи Таблицы Особенности методики расчета из "Нелинейная теория крыла и ее приложения " Выше было показано, что расчет нелинейных нестационарных характеристик профиля сводится к решению системы линейных алгебраических уравнений. Число этих уравнений равно числу неизвестных циркуляций, а коэффициенты уравнений представляют собой совокупность безразмерных скоростей, вычисленных от дискретных вихрей, моделирующих профиль и его след. [c.77] Указанная система линейных уравнений устойчива. Малые изменения коэффициентов уравнений, вызванные малыми (за расчетный промежуток времени) деформациями вихревого следа, приводят к малым изменениям решений. Это объясняется тем, что элементы главной диагонали матрицы коэ 1к1)ициентов по модулю превьпиаю остальные элементы, так как наибольшую безразмерную скорость вызьшает тот вихрь, ближе к которому находится контрольная точка. [c.77] На рис. 3.3 проводится сопоставление зависимостей (т), полученных авторами при р = 17 и взятых из работы [3.32]. Совпадение этих зависимостей во всем рассмотренном диапазоне изменения т хорошее. [c.78] Исследовалось влияние числа дискретных вихрей п по хорде пластины на практическую сходимость результатов расчета. Чтобы полнее проследить начальный этап обтекания, бралось большое число вихрей ( = 100). При большом т расчеты проводились при меньшем числе вихрей (л = 10,20), поскольку увеличение числа вихрей уже мало влияло на точность результатов. Как показали методические исследования, расчетный шаг безразмерного времени Дт целесообразно задавать обратно пропорциональным числу вихрей на яэрде (Ат /п). [c.78] Этим объясняется и то, что зависимости (т) и (т) в рассмотренном диапазоне также совпадают (см. рис. 3.8.). [c.80] На рис. 3.9 расчетные вихревые структуры пластины сравниваются с полученными экспериментально в гидроканале. ГЪзультаты расчета удовлетворительно согласуются с экспериментальными данными. [c.80] Одним из наиболее простых и удобных способов регуляции расчетов заключается в следующем. Пусть е — расстояние от вихря до бли жайшей контрольной точки. Для выбранной вихревой системы нельзя претендовать на правильное определение поля скоростей внутри малого интервала е. Поэтому всегда свободные вихри в процессе численного расчета подходят друг к другу ближе, чем на интервале е на возмущенные скорости следует ввести ограничения. На оси вихря возмущенные им скорости равны нулю. Внутри окружности с центром на оси вихря и радиусом е поля скоростей следует определять линейной интерполяцией между значениями скорости на границах или просто полагать равным нулю. [c.82] Расчеты показали, что в ряде случаев свободные вихри проскакивают сквозь поверхность профиля, что является следствием дискретности схемы по координатгш и времени. Ликвидировать этот недостаток схемы удалось введением условия, согласно которому проскочивший свободный вихрь на следующем шаге возвращается в первоначальное или близкое к нему положение. [c.82] Вернуться к основной статье