Разделы сайта
Читаемое
Обновления Mar-2024
|
Промышленность Ижоры --> Динамика жидкости: уравнения Глава 12 Обобщенные криволинейные координаты при расчете течений в сложных областях, в том числе около тел реальной формы (канал, воздухозаборник, самолет, автомобиль и т. д.), приходится рассматривать расчетные границы, не совпадающие в физическом пространстве с координатными линиями. Для конечно-разностных методов это приводит к тому Область в обобщенных коорлинатах Рис. 12.1. Соответствие между физической областью и областью в обобщенных координатах. ЧТО при постановке граничных условий требуется применять сложную интерполяцию на линии локальной сетки, из-за чего как правило, происходит локальная потеря точности численного решения. Подобные трудности служат поводом для введения отображения или преобразования физического пространства (х, у, z) к пространству (, t,) обобщенных криволинейных координат. Область в обобщенных кородинатах строится таким образом, чтобы границы в физическом пространстве совпадали с координатными линиями в пространстве обобщенных координат. Использование обобщенных координат предполагает, что искривленная в физическом пространстве область преобразуется в обобщенных координатах в прямоугольную (рис. 12.1). Уравнения должны быть записаны с использованием в качестве независимых переменных обобщенных координат, и дискретизация проводится в пространстве этих координат. Таким образом, расчет, проводимый в обобщенных координатах, становится весьма эффективным. Например, при расчете течения в двумерном искривленном канале имеет смысл сделать так, чтобы стенки канала совпадали с линиями постоянного значения г) (рис. 12.2). Перемещение вдоль стенки канала, например от А к В или от D к С, соответствует в этом случае в расчетной области изменению координаты Для точек, лежащих на некоторой ли- КМАХ j=JMAX Рис. 12.2. Двумерный искривленный канал. НИИ г, соединяющей соответствующие точки на АВ и CD, значения li постоянны, а ц меняется (от цх на АВ до т]кмах на CD), Для точки (jyk)y лежащей на этой линии г, = и ц=Цк В физической области ей соответствует точка с координатами л: =x(g/, т];) и у = у{11,щ). Использование криволинейных координат дает дополнительные возможности. Прежде всего расчетная сетка в пространстве обобщенных координат в физическом пространстве может соответствовать движущейся сетке, что имеет место в нестационарных течениях с подвижными границами. Отображение физического пространства на пространство обобщенных координат позволяет произвести сгущение координатных линий в физическом пространстве в областях, где можно ожидать появления больших градиентов. Если области больших градиентов изменяются со временем, например при перемещении ударных волн, физическую сетку можно перестраивать так, что локальная сетка будет достаточно мелкой для получения решений требуемой точности. § 12.1. Преобразование координат В данном параграфе будет установлена связь между физическими (x,y,z) и расчетными (,л,?) координатами. Обобщение на преобразования, зависящие от времени {х, у, г, t-l.r], т), приведено в работах [Sieger, 1978; Thompson et al., 1985]. 12.1,1. Обобщенные координаты Предполагается, что между физическими и обобщенными координатами существует взаимно однозначное соответствие, которое может быть записано в виде 1 = 1{х, у, Z), ц = У]{х, у, 2), Z = Ux, у, z) (12.1) и соответственно х = х{1,цХ) и т. д. Конкретные соотношения определяются при построении сетки в физической области (гл. 13). При заданных функциональных зависимостях l> = l(x,y,z) уравнения могут быть преобразованы к виду, содержащему частные производные по г], . Для стационарных течений возможность перемещения сетки в процессе проведения итераций весьма полезна. Если положение областей больших градиентов заранее неизвестно, использование обобщенных координат позволяет реализовать итерационный процесс с адаптивным построением сетки так, что поле течения во всех частях расчетной области будет получено с одинаковой точностью. В результате получится более эффективное расположение узлов сетки. Использование обобщенных координат приводит к определенным трудностям. Во-первых, необходимо определить вид уравнений в обобщенных координатах. Данный аспект будет рассмотрен в § 12.3. Численное решение типичного уравнения в частных производных будет рассмотрено в § 12.4. В уравнениях, записанных в криволинейных координатах, появляются дополнительные члены, определяющие отображение физической области на пространство обобщенных координат. Эти дополнительные члены (параметры преобразования) имеют форму производных, например дх/д1, и обычно для них также приходится проводить дискретизацию (§ 12.2). Это является до-лолнительным источником ошибок в решении (п. 12.2.3). Вывод параметров преобразования приведен в § 12.1.
|
© 2003 - 2024 Prom Izhora
При копировании текстов приветствуется обратная ссылка |