Лекция_2 Магистры

Содержание

Слайд 2

Структурированная и неструктурированные сетки Структурированная сетка (упорядоченная): накладывает ограничения на форму

Структурированная и неструктурированные сетки

Структурированная сетка (упорядоченная):
накладывает ограничения на форму элементов, структуру

их размещения и на геометрическую модель;
требует затрат большего количества «человеческого» времени;
требует меньших затрат «машинного времени» счета;
требует меньшего количества элементов для получения независящего от числа элементов решения.
Слайд 3

Получение независящего от сетки решения для структурированной и неструктурированной сетки Ф

Получение независящего от сетки решения для структурированной и неструктурированной сетки

Ф

– функция, характеризующая зависимость сходимости решения дифференциальных уравнений от количества элементов сетки N

1 – структурированная гексаэдрическая сетка;
2 – неструктурированная сетка

Следует понимать, что для получения независящего от сетки решения (расчетные значения не изменяются при дальнейшем измельчении сетки), количество элементов и, следовательно, точек интегрирования решаемых уравнений, для структурированной гексаэдрической сетки будет меньше, чем для неструктурированной. То есть на гексаэдрической сетке решение, не зависящее от сетки, будет получено за меньшее время счета.

Ф

N

Независящее от количества элементов решение

2

1

Слайд 4

Структурированную сетку следует строить в случаях, когда лимитирующим фактором для получения

Структурированную сетку следует строить в случаях, когда лимитирующим фактором для получения

решения является продолжительность вычислений, зависящая от возможностей процессора ПК, размера оперативной памяти, объема свободного пространства на жестком диске и пр.
Если необходимо провести разовое моделирование геометрически сложного объекта, и при этом нет ограничений по ресурсам, то целесообразно применение неструктурированной сетки для экономии затрат времени на её создание.
Слайд 5

Примеры структурированной и неструктурированной сетки Структурированная гексагональная сетка Неструктурированная сетка

Примеры структурированной и неструктурированной сетки

Структурированная гексагональная сетка

Неструктурированная сетка

Слайд 6

Пример смешения нескольких типов сеток в одной модели Призматический подслой в

Пример смешения нескольких типов сеток в одной модели

Призматический подслой в тетраэдрической

сетке – компромисс между структурированной и неструктурированной сеткой
Слайд 7

Примеры использования призматических слоев Призматические структурированные подслои чаще всего используют в

Примеры использования призматических слоев

Призматические структурированные подслои чаще всего используют в пограничной

области; у поверхности контакта разнородных материалов; у поверхности раздела фаз.
Слайд 8

II. Методы построения сеток 1. Метод граничной коррекции Разработан Марком Шепардом

II. Методы построения сеток

1. Метод граничной коррекции
Разработан Марком Шепардом и Марком

Йери в политехническом университете Ренсселера (США) в 1984 году и получил название метод «Octree / Quadtree»

Основные этапы построения расчетной сетки с использованием
метода «Octree / Quadtree»

Идея метода: на заданную область «накладывается» предварительно сгенерированная сетка. Далее от этой сетки отсекаются все лишние фрагменты, которые выходят за пределы расчетной области; затем корректируется положение узлов, лежащих вблизи границы, так, чтобы они попали в «углы», на «ребра» и на «грани» области.

1) построение "первичной" сетки; 2) коррекция сетки.

Слайд 9

Пошаговая визуализация метода граничной коррекции: 1. Построение "первичной" сетки 1.1. Разбиение на квадраты / кубы

Пошаговая визуализация метода граничной коррекции:
1. Построение "первичной" сетки
1.1. Разбиение на

квадраты / кубы
Слайд 10

1.2. Разбиение на треугольники / тетраэдры "Классический" набор шаблонов для разбиения квадратов

1.2. Разбиение на треугольники / тетраэдры

"Классический" набор шаблонов для разбиения квадратов

Слайд 11

2. Корректировка "первичной" сетки 2.1. Удаление элементов, полностью выходящих за границу

2. Корректировка "первичной" сетки
2.1. Удаление элементов, полностью выходящих за границу

геометрии
2.2. Смещение выступающих
элементов сетки
Достоинства и недостатки метода:

+ Обладает достаточно высокой скоростью работы
+ Сравнительная простота реализации
– построенные сетки необходимо проверять на правильность структуры для исключения вырожденных тетраэдров
– низкое качество элементов вблизи границы, поэтому необходим этап оптимизации
– низкая "чувствительностью", при недостаточно малом шаге первичной сетки некоторые особенности области могут быть потеряны.

Типичные ошибки – это "слипание" узлов (когда два узла первичной сетки оказываются перемещены в одну точку), что приводит к появлению вырожденных тетраэдров, а также появлению тетраэдров с нулевым объемом (все четыре узла оказываются в одной плоскости).

Слайд 12

Метод граничной коррекции используется приложением Meshing в методе построения сетки Mesh Control

Метод граничной коррекции используется приложением Meshing в методе построения сетки Mesh

Control < Method < Tetrahedrons (опция Patch Independent [независимое разбиение]).
Слайд 13

2. Метод «продвижения фронта» (Advancing Front) Впервые идея метода предложена Рейнальдом

2. Метод «продвижения фронта» (Advancing Front)

Впервые идея метода предложена Рейнальдом Лонером

(R. Lohner), а его трехмерный вариант разработал профессор Гонконгского университета С.Х. Ло.
Общая идея этого класса методов заключается в последовательном изымании из заданной области фрагментов тетраэдрической или гексаэдрической формы до тех пор, пока вся область не окажется "исчерпана“ (заполнена конечными объемами).
Слайд 14

Пошаговая визуализация метода «продвижения фронта» 1. Отправной точкой алгоритма является начальная

Пошаговая визуализация метода «продвижения фронта»

1. Отправной точкой алгоритма является начальная дискретизация

границы заданной области, наиболее точно описывающая исходную геометрию.
2 – 3. Формирование фронта дискретизации – набора отрезков/ребер или треугольных граней.

1

2

3

Слайд 15

4 – 5. Поиск оптимального узла С для каждой грани/ребра АВ.

4 – 5. Поиск оптимального узла С для каждой грани/ребра АВ.
6.

Попытка образовать треугольник АВС, соединив концы ребер.

4

5

6

Слайд 16

7- 8. Задание радиуса поиска узлов r. Поиск других узлов в

7- 8. Задание радиуса поиска узлов r. Поиск других узлов в

пределах радиуса r для формирования второго элемента на базе отрезка AB. Вместо узла С – выбор узла D.
9-10. Добавление в сетку нового узла и образование с ним элемента. Добавление нового элемента в сетку и обновление фронта.

7

8

9

10

Слайд 17

11. Переход к следующей итерации. 11 Методы «продвижения фронта» универсальны и

11. Переход к следующей итерации.

11

Методы «продвижения фронта» универсальны и могут быть

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

Триангуляция пространства вокруг модели самолета для решения задачи аэродинамики

Слайд 18

Метод «продвижения фронта» используется приложением Meshing в методе построения сетки Mesh Control

Метод «продвижения фронта» используется приложением Meshing в методе построения сетки Mesh

Control < Method < Sweep.
Слайд 19

+ Сетки обладают неплохим качеством. + Наиболее эффективны, если изначально задана

+ Сетки обладают неплохим качеством.
+ Наиболее эффективны, если изначально задана дискретизация

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

Достоинства и недостатки метода

Слайд 20

3. Сетки на основе критерия Делоне (Delaunay Mesh) Впервые предложены Чарльзом

3. Сетки на основе критерия Делоне (Delaunay Mesh)
Впервые предложены Чарльзом Лоусоном

(Charles Lawson) и Дэвидом Уотсоном (Dave Watson), 1977-1981 гг.

Метод основе критерия Делоне используется приложением Meshing в методе построения сетки Mesh Control < Method < Tetrahedrons (опция Patch-Conforming [согласованное разбиение]).

Слайд 21

Критерий Делоне: треугольная сетка на плоскости удовлетворяет критерию Делоне, если внутрь

Критерий Делоне: треугольная сетка на плоскости удовлетворяет критерию Делоне, если внутрь

окружности, описанной вокруг любого треугольника, не попадают никакие другие узлы этой сетки.
Этот термин также употребляется и по отношению к треугольнику сетки: треугольник удовлетворяет критерию Делоне, если критерию Делоне удовлетворяет сетка, составленная только из самого треугольника и соседних с ним треугольников.

Сетка, удовлетворяющая критерию Делоне (слева), и не удовлетворяющая ему (справа).

Слайд 22

Пошаговая визуализация метода на основе критерия Делоне 1. Формирование множества U

Пошаговая визуализация метода на основе критерия Делоне
1. Формирование множества U -

набора заданных узлов
2. Создание "суперструктуры", представляющей собой произвольный выпуклый многогранник с треугольными гранями, такой, что все заданные узлы лежат внутри него.

1

2

Слайд 23

3. Разбиение на треугольники с учетом критерия Делоне с вершинами на заданном множестве точек 3.1 3.2

3. Разбиение на треугольники с учетом критерия Делоне с вершинами на

заданном множестве точек

3.1

3.2

Слайд 24

3.3 3.4

3.3

3.4

Слайд 25

4. Удаление из сетки всех треугольников, среди вершин которых были вспомогательные узлы суперструктуры

4. Удаление из сетки всех треугольников, среди вершин которых были вспомогательные

узлы суперструктуры
Слайд 26

Достоинства и недостатки метода на основе критерия Делоне: + оптимальное разбиение

Достоинства и недостатки метода на основе критерия Делоне:
+ оптимальное разбиение расчетной

об­ласти по заданному множеству узлов с точки зрения получения максимально возможных значений малых углов треугольников (стремление полученных треугольников к равноугольным).
- крайне высокая чувствительность к точности машинных вычислений.
- снижение скорости вычислений.
Слайд 27

Архитектура комплекса ANSYS CFX CAD DesignModeler (Создание геометрии) ANSYS Meshing ANSYS

Архитектура комплекса ANSYS CFX

CAD
DesignModeler
(Создание геометрии)

ANSYS Meshing
ANSYS ICEM CFD
(Генерация сетки)

CFX - Pre
(Препроцессор
физики)

CFX

- Solver
(Решатель)

CFX - Post
(Постпроцессор)

CFX – Solver Manager
(Менеджер задач)

ANSYS CFX

Geometry (Геометрия)

Mesh (Сетка)

Solver (Решение)

Setup (Ввод исх. данных)

Results (Результаты)

ANSYS Workbench

Слайд 28

Структура взаимодействия информационных потоков в ANSYS CFX ANSYS CFX Начальные условия

Структура взаимодействия информационных потоков в ANSYS CFX

ANSYS CFX

Начальные
условия

Граничные
условия

Выходные данные:
поля скоростей,
давлений и

пр.

Сеточная
модель

Геометрическая
модель

Дополнительная
информация

Слайд 29

Общий подход к моделированию гидро-, газодинамики Гидро-, газодинамика – это раздел

Общий подход к моделированию гидро-, газодинамики

Гидро-, газодинамика – это раздел механики,

описывающий движение жидкостей и газов в рамках модели сплошной среды.
Сплошная среда (СС) – механическая система, обладающая бесконечным числом внутренних степеней свободы, и которую можно рассматривать как непрерывную, пренебрегая её дискретным атомно-молекулярным строением.
Различают:
1) однородную СС – в разных точках среды физические свойства одинаковы при одинаковых температуре и давлении;
2) неоднородную СС – в разных точках среды физические свойства неодинаковы при одинаковых температуре и давлении;
3) изотропная и анизотропная СС (одинаковые/разные физические свойств среды по всем направлениям);
4) сжимаемая и несжимаемая СС (плотность СС не равна/равна константе).
Слайд 30

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

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

которые описывают процессы переноса импульса, массы и энергии и известны как уравнения Навье-Стокса. Эти уравнения были получены в начале XIX века и не имеют общего аналитического решения, но могут быть дискретизированы и решены численно.
Дискретизация – преобразование непрерывной функции в дискретную.
Система уравнений гидро-, газодинамики в общем виде состоит из следующих зависимостей:

МЕТОДОЛОГИЯ ЧИСЛЕННЫХ РАСЧЕТОВ

Слайд 31

1. Уравнение неразрывности (уравнение сохранения масс): где – плотность среды; – время; – вектор скорости.

1. Уравнение неразрывности (уравнение сохранения масс):

где – плотность среды;
– время;

– вектор скорости.

 

Слайд 32

 

Слайд 33

2. Уравнение движения (сохранения количества движения): Это три уравнения для компонент

2. Уравнение движения (сохранения количества движения):

Это три уравнения для компонент вектора

скорости, которые являются конкретизацией второго закона Ньютона (a = F/m) для контрольного объема V = dx∙dy∙dz, в котором происходит изменение количества движения ρ∙v = G (импульса).
Векторная сумма импульсов всех тел системы есть величина постоянная, если векторная сумма внешних сил, действующих на систему тел, равна нулю.

– давление в среде;
– тензор напряжений;
– вектор массовых сил.

– динамическая вязкость среды;
– индекс, обозначающий оператор транспонирования матрицы;
– дельта-функция Кронекера.

Слайд 34

Проекция на ось х:

 

Проекция на ось х:

 

Слайд 35

3. Полное уравнение сохранения энергии (учитывает теплообмен и сжимаемость среды): –

3. Полное уравнение сохранения энергии
(учитывает теплообмен и сжимаемость среды):

– полная энтальпия,

равная сумме статической энтальпии и скоростного напора;
– температура;
– теплопроводность;
– источниковый член.

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

Слайд 36

 

 

 

Слайд 37

4. Уравнение состояния: – удельная теплоёмкость

4. Уравнение состояния:

 

– удельная теплоёмкость

Слайд 38

Методы решения систем уравнений Метод конечных разностей; Метод контрольных объемов; Метод

Методы решения систем уравнений

Метод конечных разностей;
Метод контрольных объемов;
Метод конечных элементов;
Метод сглаженных

частиц;
Метод с использованием функции распределения вероятностей.
Слайд 39

Метод контрольных объемов Дискретизация – преобразование непрерывной функции в дискретную. ANSYS

Метод контрольных объемов

Дискретизация – преобразование непрерывной функции в дискретную.
ANSYS CFX использует

метод конечных объемов на основе элементов дискретизации пространственной области с использованием сетки. Сетка нужна для построения конечных объемов, которые используются для применения законов сохранения соответствующих величин, таких как масса, импульс и энергия. Сетка трехмерна, но для простоты рассмотрим двухмерную.
Построение сеточной модели – дискретизация пространства.
Слайд 40

Все переменные решения и свойства текучей среды хранятся в узлах Node

Все переменные решения и свойства текучей среды хранятся в узлах Node

(вершины сетки). Контрольный объем Control Volume (заштрихованная область) строится вокруг каждого узла сетки следующим образом: контрольный объем ограничивается линиями, соединяющими центры ребер (т. 1) и центры граней Element Center (т. 2) сеточных элементов Element, окружающих узел Node (т. 0).

Типичная двумерная сетка

1

2

0

Слайд 41

Методология метода конечного объёма Для иллюстрации методологии метода конечного объема рассмотрим

Методология метода конечного объёма

Для иллюстрации методологии метода конечного объема рассмотрим уравнения

сохранения массы, импульса, выраженные в декартовых координатах:

 

 

Эти уравнения интегрируются по каждому контрольному объему с использованием теоремы Гаусса о преобразовании объемных интегралов объема в поверхностные интегралы.

Слайд 42

Методология метода конечного объёма где V и S соответственно, объемные и

Методология метода конечного объёма

 

 

где V и S соответственно, объемные и поверхностные

области интегрирования, а dni - дифференциальные декартовы компоненты внешнего нормального поверхностного вектора.
Следующим шагом в численном алгоритме является дискретизация объемных и поверхностных интегралов.
Слайд 43

Объемные интегралы дискретизируются в каждом секторе Sector сеточного элемента Element и

Объемные интегралы дискретизируются в каждом секторе Sector сеточного элемента Element и

накапливаются в контрольном объеме Control Volume, к которому принадлежит сектор.
Поверхностные интегралы дискретизируются в точках интегрирования (ipn), расположенных в центре грани каждого сегмента сеточного элемента.

Методология метода конечного объёма