Историческое развитие методов моделирования полей
Метод электрогидродинамической аналогии (метод сплошных сред)
Рассмотрим вопрос развития методов моделирования полей в историческом плане. Это поможет более глубоко понять сущность проблемы и области применимости тех или иных методов.
В качестве первого, наиболее выдающегося примера мышления по аналогии, следует отметить предложенный в 1911 году Павловским метод электрогидродинамической аналогии (ЭГДА), который известен также как метод сплошных сред (рис.7.7 ).
Рис. 7.7 Метод электрогидродинамической аналогии
Идея метода, достаточно наглядно иллюстрируемая рисунком на примере исследования фильтрации воды под плотиной, заключается в следующем. Интенсивность фильтрации воды под плотиной (что может привести к ее осадке и разрушению) зависит прежде всего от разности уровней верхнего и нижнего бьефов, свойств грунта и глубины заложения шпунта (рис.7.7, а). Авторы проектов при расчете плотин, стремясь обезопаситься, как правило, завышают глубину заложения шпунта, что выражается в перерасходе десятков и даже сотен тысяч тонн бетона.
Павловский предложил (рис.7.7,б) в качестве аналога грунта использовать лист станиоля (электропроводной бумаги), форма выреза в котором воспроизводит нижний профиль плотины, а разность потенциалов напряжений – аналог перепада высот и . При этом линии равного потенциала становятся аналогами линий равного давления, от расположения которых зависит уровень фильтрации. Естественно, что в первоначальном виде метод имел некоторые недостатки. Например, для каждого варианта нижнего профиля плотины нужно было иметь отдельный лист электропроводящего материала. Возникали также трудности при моделировании грунтов с неоднородными по плотности слоями.
Метод электрических сеток
Основа метода электрических сеток
Качественно новый уровень этот метод приобрел в 60 – 70 годы, когда на основе этой идеи и нового технического уровня было создано целое семейство, так называемых, сеточных интеграторов. Метод электрических сеток, основанный на математическом описании полей с помощью уравнений в конечных разностях, позволил значительно расширить круг решаемых задач. С известным приближением модели в виде сплошной среды можно заменить электрической сеткой, состоящей из схем замещения элементарных объемов или площадей (рис.7.8), на которые может быть разделена сплошная среда, например, ранее представленная на рис.7.7.
Рис. 7.8 Схема замещения элементов и плоскопараллельного слоя
Преимущество метода электрических сеток
Важным преимуществом электрических сеток является возможность моделирования трехмерных полей, в том числе описываемых уравнениями с правой частью. Решение такого рода задач даже на достаточно быстродействующих современных цифровых ЭВМ представляет значительные трудности.
Решение задач, описываемых уравнениями с правой частью
Моделирование уравнения Лапласа
На рис.7.9 представлена схема замещения элементарного объема потенциального поля, описываемого уравнением Лапласа
Рис. 7.9 Схема замещения элементарного объема
- – потенциал поля,
- – удельная проводимость.
В декартовых координатах при .
Это уравнение принимает следующий вид
Моделирование уравнения Пуассона
Если в узлы электрической сетки включить источники тока в соответствии с некоторой функцией распределения , то такая сетка (рис.7.10) будет моделировать уравнение Пуассона.
Рис. 7.10 Схема замещения элементарного объема со стоком и источником
Моделирование уравнений при включении в узлы сетки стыков
Стоки – резистивные сопротивления
Включение в узлы сетки нагрузочных проводимостей на землю (стыков) открывает возможность моделирования уравнений с другими формами правой части. В зависимости от вида проводимости стока ток через него оказывается пропорциональным напряжению узла сетки относительно земли или его производной по времени, что позволяет задавать граничные условия первого или второго рода.
Если стоки представляют собой резистивные сопротивления (рис.7.10), получаем модель для уравнения вида
Стоки – конденсаторы
Если для той же сетки сделать стоки в виде конденсаторов (рис.7.11), то будет моделироваться уравнение диффузии (или уравнение теплопроводности Фурье)
Рис. 7.11 Схема замещения элементарного объема с емкостным стоком
Стоки – катушка индуктивности
Если заменить резисторную сетку на реактивную катушку индуктивности вместо резисторов, получаем модель для волнового уравнения
Таким образом возможность введения стоков и источников значительно расширяет круг задач, решаемых методом электрических сеток, по сравнении с методом сплошных сред.
Решение задач подземной гидравлики нефтяных месторождений
Очень показательным примером широкого практического применения сеточного моделирования явилось решение задач подземной гидравлики нефтяных месторождений [55]. Нефтяной пласт представляет собой пористую среду пропитанную нефтью, которая движется под действием перепада давлений между внешним контуром питания и скважинами. Величина дебитов отдельных скважин зависит от распределения давлений, которое, в свою очередь, зависит от расположения и взаимного влияния скважин. Ставится задача определения дебитов скважин при их оптимальном расположении, которая может быть представлена в виде решения уравнения электрического типа.
Моделируемое поле представляется по форме нефтяной залежи, на контуре питания и на скважинах устанавливаются потенциалы, пропорциональные давлениям. Относительно малая (по сравнению с площадью) толщина пласта позволяет рассматривать задачу как двумерную, выполняя модель в виде слоя переменной толщины. На рис.7.12 представлен фрагмент такой модели в виде электрической сетки. Гидравлические сопротивления скважин моделируются с помощью добавочных сопротивлений в цепи тока. Токи в этих сопротивлениях являются аналогами дебита скважин.
Рис. 7.12 Сеточная модель нефтяного месторождения
Моделирование одного из действующих нефтяных месторождений, на котором было набурено около ста скважин, показало, что половину из этих скважин не нужно было бурить, то есть такой же суммарный дебит могли обеспечить 50 скважин. Многомиллионные затраты на бесполезные и даже “вредные” скважины, просто снижающие давление пласта, очень наглядно подчеркнули эффективность использования методов математического моделирования и невозможность решения такого рода задач по “наитию”. Следует подчеркнуть, что упомянутая выше модель представляла собой весьма сложную систему из восьмисот уравнений, которые решились на гибридном (аналого-цифровом) комплексе. Этот подход в дальнейшем получил достаточно широкое развитие на новом научном и техническом уровне применительно к современным цифровым ЭВМ, быстродействие которых за последние годы сильно возросло.
Метод конечных элементов
Одной из широко распространенных направлений такого подхода (также основанного на принципе прямой аналогии) является метод конечных элементов. Не останавливаясь здесь на описании этого метода, мы несколько ниже рассмотрим пример его реализации с использованием системы компьютерной математики ”MATLAB”.
Модели полей тепло и массопереноса на табличном процессоре Excel
Использование табличного процессора Excel
Большинство пользователей персональных ЭВМ знакомо с замечательными свойствами электронных таблиц, в которых с помощью определенных математических соотношений можно связать ячейки таким образом, чтобы при изменении числовых данных в одних ячейках автоматически менялись результаты также в других связанных с ними ячейках. Во многих случаях применение табличного процессора Excel, входящего в состав основного программного обеспечения Microsoft Office, ограничивается сравнительно простыми вычислениями типа калькуляционных. Однако на самом деле в этом процессоре кроются огромные возможности, которые представляются при квалифицированном использовании встроенного в него внутреннего оптимизатора. При этом становится возможным превратить ячейки электронных таблиц в цифровой аналог рассмотренных выше элементов замещения поля и таким образом решать задачи, описываемые уравнениями в частных производных. Пример такой задачи приведен ниже.
Моделирование процесса теплообмена в трехмерном приближении
Для математического описания теплообмена трехмерного тела с геометрической областью и кусочно-гладкой границей использованы модели в виде краевых задач математической физики III рода для дифференциального уравнения с частными производными параболического типа в случае нестационарного теплообмена:
Особенности методики расчета и реализации на ЭВМ трехмерных задач теплообмена
Особенности
Особенность методики решения задачи заключается в разработке соответствующих Excel-приложений на языке визуального и событийно-управляемого программирования VBA (Visual Basic for Application) в составе Microsoft Office 97 операционной системы Windows, реализующих алгоритмы расчета, соответствующие математическим моделям (7.22) и (7.23).
Основные этапы методики
Основные этапы методики сводятся к следующему:
- Анализ геометрии исследуемого трехмерного тела и выделение на основе принципов симметрии расчетной области для постановки краевых задач.
- Разбиение расчетной области на элементарные непересекающиеся однотипные подобласти с учетом особенностей исходной геометрии .
- Аппроксимация задач (7.22) и (7.23) соответствующей конечно-разностной задачей с помощью метода конечных разностей.
- Разработка и реализация в модулях Excel-приложений на VBA итерационных типовых алгоритмов решения конечно-разностных уравнений для областей . Расчет распределений температуры в трехмерной области производится набором процедур в модулях VBA с использованием трехмерных структур данных (например, массивов).
- Проектирование и реализация интерфейса пользователя путем создания соответствующих диалоговых окон, элементов управления и диаграмм для оперативного управления процессом вычислительного эксперимента и анализа его результатов в виде как числовой, так и графической информации. Параметрами оперативного управления являются как теплофизические характеристики горячих сред и охладителей, так и конструктивные параметры охлаждаемых элементов исследуемых объектов. Данные для визуализации распределения температуры формируются на листах табличного процессора Excel с помощью процедур извлечения двумерных структур данных в соответствии с интересующими сечениями трехмерной расчетной области. Эти процедуры так же находятся в модулях VBA. На основе двумерных структур данных формируются диаграммы распределения температурных полей в сечениях, которые компонуются в виде документа, характеризующего в целом распределение температуры по объему тела.
- Формирование критерия оптимизации в соответствии с требуемой естественно-физической постановкой задачи и реализация поставленной оптимизационной задачи с помощью оптимизатора (Solver) табличного процессора Excel.
Исследование теплообменных процессов в плоской стенке при охлаждении с помощью специальных охлаждаемых элементов
Исследование температурных полей в медной плоской стенке
Применение описанной выше методики рассмотрено [57] на примере исследования температурных полей в медной плоской стенке толщиной охлаждаемой с помощью специальных элементов охлаждения, в которых циркулирует в качестве охладителя вода с температурой и которые конструктивно вмонтированы в специальные посадочные гнезда глубиной , расположенные на стенке с охлаждаемой стороны равномерно, т.е. в вершинах квадратов, на расстоянии между центрами элементов. Нагрев стенки осуществляется с противоположной стороны за счет контактирования с горячим газом, температура которого равна . Задача решается в декартовом приближении. Сечение элемента водоохлаждения представляет собой квадрат, сторона которого равна . Фрагмент стенки с выделенной расчетной областью представлен на рис.7.13.
Рис. 7.13 Фрагмент плоской стенки с посадочными гнездами для элементов охлаждения
Для конечно-разностной аппроксимации краевых задач область покрыта равномерной по каждой кординате сеткой с шагами,и индексацией узлов с помощью индексов ,по соответствующей координате. По времени выбран шаг дискретизации.
Рис. 7.14 Температурные поля в сечении , проходящем через центры элементов водоохлаждения при значениях глубины посадки: а) 0,02 м; б) 0,03 м; в) 0,04 м
На рис.7.14 представлены распределения температуры в сечении (рис.7.13), проходящем через центры элементов водоохлаждения, для случаев трех значений глубины посадки элементов водоохлаждения: 0,02 м, 0,03 м, 0,04 м. При этом основные конструктивные размеры имели следующие значения:,,.Значения температур, приведенные на диаграммах, соответствуют значениям температур на изолиниях. Температура газа равна , воды – . Коэффициент теплоотдачи со стороны газа определялся на основании теории лучистого теплообмена и принимал значения в рассматриваемых трех случаях от 70 до 75, со стороны воды коэффициент теплоотдачи задавался равным 2000 . Распределение температуры в сечении Б (рис.7.13), проходящем через середину толщины стенки, для тех же случаев представлено на рис. 7.15.
Рис. 7.15 Температурные поля в сечении , проходящем через середину толщины плоской стенки при значениях глубины посадки элементов водоохлаждения: а) 0,02 м; б) 0,03 м; в) 0,04 м
Анализ приведенной графической информации в серии расчетов в совокупности с решением оптимизационных задач позволяет правильно выбрать как режимы охлаждения, так и конструктивные решения при проектировании систем охлаждения технологическими агрегатами.