2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 Пересчет функций на границе по Эйлеру
Сообщение24.08.2007, 08:32 


09/07/07
26
Кемерово
В некоторой области задано течение идеальной несжимаемой жидкости. Расчет проводится методом, подобным МКЭ. Для движения используется уравнение Эйлера, которое выписано в переменных Лагранжа (нет конвективного слагаемого, производная по времени свернута в полную), расчетные узлы движутся вместе с сеткой. При таком движениии возникает разреженность узлов на границе. Можно ли граничные узлы оставлять на месте, но пересчитывать в них значения скорости через полную производную (из найденной скорости в узле в Лагранжевых переменных отнимать конвективный член). Не сталкивался ли кто-нибудь с подобной задачей? Как можно обойти подобный момент?

 Профиль  
                  
 
 
Сообщение24.08.2007, 12:51 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
Проще сделать перинтерполяцию значений в узлах и ячейках сетки с деформированной на исходную.

 Профиль  
                  
 
 
Сообщение08.09.2007, 07:22 


09/07/07
26
Кемерово
Zai писал(а):
Проще сделать перинтерполяцию значений в узлах и ячейках сетки с деформированной на исходную.


Двигать узлы по найденным скоростям, строить по ним, например, сплайн и переразбивать границу? Или есть еще какой-нибудь вариант? Сплайн нужно выбрать подходящий. При сближении 2 точек на малое расстояние и сильно отличающихся значениях в них могут быть большие осцилляции после интерполяции.

 Профиль  
                  
 
 
Сообщение08.09.2007, 08:13 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
Сплайн может исказить консервативность алгоритма. Старая и новая сетки можно переразбить на промежуточные элементарные фигуры. Для каждой элементарной фигуры со старой сетки определяются параметры ячеек и узлов на переразбитой. Затем исходя из законов сохранения энергии и импульса определяются параметры в новых ячейках и новых узлах с соответствующим осреднением. Данная процедура называется rezoning. Она вносит дополнительную искусственную вязкость в алгоритмы, но более достаточно устойчива от выворачивания ячеек по сравнению с прямым учетом конвективных членов.

 Профиль  
                  
 
 
Сообщение08.09.2007, 09:00 


09/07/07
26
Кемерово
Zai писал(а):
Сплайн может исказить консервативность алгоритма. Старая и новая сетки можно переразбить на промежуточные элементарные фигуры. Для каждой элементарной фигуры со старой сетки определяются параметры ячеек и узлов на переразбитой. Затем исходя из законов сохранения энергии и импульса определяются параметры в новых ячейках и новых узлах с соответствующим осреднением. Данная процедура называется rezoning. Она вносит дополнительную искусственную вязкость в алгоритмы, но более достаточно устойчива от выворачивания ячеек по сравнению с прямым учетом конвективных членов.

Не подскажите, где более подробно можно об этом почитать? Уже пробовал решать задачи при помощи сплайнов. Возникают трудности. В принципе, вообще не хотелось бы перемещать граничные узлы, а только пересчитывать значения функции в них. Сплайн показался нормальным выходом - сначала узел убежал по Лагранжеву подходу, а потом мы его вернули назад и выставили функцию через сплайн.

 Профиль  
                  
 
 
Сообщение09.09.2007, 16:26 
Заслуженный участник
Аватара пользователя


11/04/07
1352
Москва
Подробности rezoning мне не совсем известны. Я знаю что они исползуются а пакете ABAQUS/Explicit. Может быть на сайте этого пакета Вы найдете Вам необходимые ссылки.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 6 ] 

Модераторы: photon, profrotter, Парджеттер, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group