- •Моделирование физических процессов Учебное пособие
- •Содержание
- •3. 8. Математическая модель массообмена ……………………………. 54
- •5.1.1. Формулы конечного дифференцирования …………………… 74
- •5.6.2. Метод Ньютона ………………………………………………. 95
- •Введение
- •VII этап
- •2. Структурные математические модели
- •Основы структурного моделирования
- •2.2. Типовые элементы
- •2.2.1. Способы математического описания звеньев и систем.
- •2.2.2. Типовые динамические звенья.
- •2.3. Схемы соединения типовых элементов
- •2.4. Структурные модели теплообменного оборудования
- •2.4.1. Составление математических моделей тепловых объектов.
- •3. Физические математические модели
- •3.1 Использование математических моделей
- •Под углом 57 к горизонту со скоростью 20 м/с
- •Изменение температуры при остывании тела
- •3.2 Математическая модель теплопроводности
- •. Математическая модель гидродинамики
- •3.4 Математическая модель теплообменника
- •3.5 Математическая модель парогенератора
- •С естественной циркуляцией
- •Математическая модель горения
- •Математическая модель образования токсичных
- •3. 8. Математическая модель массообмена
- •Дифференциальные уравнения диффузионного пограничного слоя
- •Моделирование теплоотдачи диффузией.
- •3.9. Математическая модель парового котла
- •3.9.1. Математическая модель динамических процессов.
- •3.9.2. Регулирование давления пара и тепловой нагрузки
- •3.10. Пример построения математической модели объекта
- •4. Планирование эксперимента
- •Полный трехфакторный эксперимент
- •Численные методы
- •5.1. Интерполяция
- •5.1.1. Формулы конечного дифференцирования
- •5.1.2. Сплайны
- •Способы получения краевых условий:
- •5.2. Решение систем алгебраических линейных уравнений
- •Метод прогонки. Метод прогонки используется для решения систем специального вида
- •5.3. Решение систем алгебраических нелинейных уравнений
- •5.3.1. Метод итераций
- •5.3.2. Метод Ньютона
- •5.4. Решение дифференциальных уравнений
- •5.4.1. Постановка задачи. Разностные схемы.
- •5.4.2. Разностные схемы. Метод Рунге – Кутта.
- •Пример решения дифференциального уравнения методом Рунге-Кутта.
- •Результаты решения уравнения методом Рунге-Кутта
- •5.5. Решение краевых задач методом прогонки
- •5.5.2. Распространение тепла в пластине радиатора.
- •Сравнение метода прогонки с точным решением
- •5.6. Решение краевых задач методом пристрелки
- •5.6.1. Метод пристрелки.
- •5.6.2. Метод Ньютона. Рассмотрим граничную задачу, определяемую дифференциальным уравнением второго порядка
- •5.7. Решение краевых задач методом конечных разностей
- •5.7.1. Метод конечных разностей.
- •5.7.2. Линейные дифференциальные уравнения.
- •Моделирующие программы
- •Смеси газов в реакторе
- •Перемещение свободной поверхности жидкости
- •Скоростей течения испаряющегося метана
- •Рекомендованный библиографический список
- •Пасько Петр Иванович «моделирование физических процессов»
- •346360, Г. Волгодонск, ул. Ленина 73/94
VII этап
Рис. 1.1. Основные этапы математического моделирования
Второй этап – это формальное, математическое описание РС. Это описание в виде математических соотношений, устанавливающих связь между параметрами, характеризующими РС ТО, и называют математической моделью (ММ).
На третьем этапе проводят качественный и оценочный количественный анализ построенной ММ. При этом могут быть выявлены противоречия, ликвидация которых потребует уточнения или пересмотра РС (штриховая линия на рис. 1.1). Количественные оценки могут дать основания упростить модель, исключив из рассмотрения некоторые параметры, соотношения или их отдельные составляющие, несмотря на то, что влияние описываемых ими факторов учтено в РС. В большинстве случаев, принимая дополнительные по отношению к РС допущения, полезно построить такой упрощенный вариант ММ, который позволял бы получить или привлечь известное точное решение. Это решение затем можно использовать для сравнения при тестировании результатов на последующих этапах. В некоторых случаях удается построить несколько ММ для одного и того же ТО, отличающихся различным уровнем упрощения. В этом случае говорят об иерархии ММ, т.е. об упорядочении ММ по признаку их сложности и полноты. Построение иерархии ММ связано с различной детализацией свойств изучаемого ТО. Сравнение результатов исследования различных ММ может существенно расширить и обогатить знания об этом ТО. Кроме того, такое сравнение позволяет оценить достоверность результатов последующего вычислительного эксперимента: если более простая ММ правильно отражает некоторые свойства ТО, то результаты исследования этих свойств должны быть близки к результатам, полученным при использовании более полной и сложной ММ. Итог анализа на рассматриваемом этапе – это обоснованный выбор рабочей ММ ТО, которая подлежит в дальнейшем детальному количественному анализу. Успех в проведении третьего этапа зависит, как правило, от глубины понимания связи отдельных составляющих ММ со свойствами ТО, нашедшими отражение в его РС, что предполагает органическое сочетание владения математикой и инженерными знаниями в конкретной предметной области.
Четвертый этап состоит в обоснованном выборе метода количественного анализа ММ, в разработке эффективного алгоритма вычислительного эксперимента, а пятый этап – в создании работоспособной программы, реализующей этот алгоритм средствами вычислительной техники.
Получаемые на шестом этапе (в итоге работы программы) результаты вычислений должны пройти тестирование путем сопоставления с данными количественного анализа выбранного упрощенного варианта ММ рассматриваемого ТО. Тестирование может выявить недочеты как в программе, так и в алгоритме, и потребовать доработки программы или же модификации и алгоритма, и программы. Анализ результатов вычислений и их инженерная интерпретация могут вызвать необходимость в корректировке РС и соответствующей ММ. После устранения всех выявленных недочетов триаду „модель – алгоритм – программа" можно использовать в качестве рабочего инструмента для проведения вычислительного эксперимента и выработки на основе получаемой количественной информации практических рекомендаций, направленных на совершенствование ТО, что составляет содержание седьмого, завершающего „технологический цикл" этапа математического моделирования.
Свойства математической модели:
Полнота ММ – это подробность описания необходимых с точки зрения модели параметров ТО.
Точность ММ – это допустимая величина несовпадения реальных и найденных при помощи ММ значений выходных параметров ТО.
Адекватность ММ – это способность ММ отражать характеристики ТО с относительной погрешностью не более некоторого заданного значения . В более общем смысле под адекватностью ММ понимают правильное качественное и достаточно точное количественное описание именно тех характеристик ТО, которые важны в данном конкретном случае. Модель, адекватная при выборе одних характеристик, может быть неадекватной при выборе других характеристик того же ТО.
Экономичность ММ оценивают затратами на вычислительные ресурсы, необходимые для реализации ММ на ЭВМ.
Робастность ММ характеризует способность ММ не допускать чрезмерного влияния погрешностям исходных данных на результат вычислительного эксперимента. Причинами низкой робастности ММ могут быть необходимость при ее количественном анализе вычитания близких друг к другу приближенных значений величин или деления на малую по модулю величину, а также использование в ММ функций, быстро изменяющихся в промежутке, где значение аргумента известно с невысокой точностью.
Продуктивность ММ связана с возможностью располагать достаточно достоверными исходными данными, точность определения которых должна быть выше, чем для тех параметров, которые получаются при использовании ММ. В противном случае ММ будет непродуктивной и ее применение для анализа конкретного ТО теряет смысл. Ее можно будет использовать лишь для оценки характеристик некоторого класса ТО с гипотетическими исходными данными.
Наглядность ММ является ее желательным, но необязательным свойством. Тем не менее использование ММ упрощается, если ее составляющие имеют ясный содержательный смысл.
Структурные и функциональные математические модели. Различные особенности и признаки математических моделей лежат в основе их классификации. Среди таких признаков выделяют характер отображаемых свойств технического объекта, степень их детализации, способы получения и представления ММ.
Один из существенных признаков классификации связан с отражением в ММ тех или иных особенностей ТО. Если ММ отображает устройство ТО и связи между составляющими его элементами, то ее называют структурной математической моделью. Если же ММ отражает происходящие в ТО физические, механические, химические или информационные процессы, то ее относят к функциональным математическим моделям. Существуют и структурно-функциональные математические модели.
Структурные ММ делят на топологические и геометрические, составляющие два уровня иерархии ММ этого типа. Первые отображают состав ТО и связи между его элементами. Топологическую ММ целесообразно применять на начальной стадии исследования сложного по структуре ТО, состоящего из большого числа элементов, прежде всего для уяснения и уточнения их взаимосвязи. Геометрическая ММ дополнительно к информации, представленной в топологической ММ, содержит сведения о форме и размерах ТО и его элементов, об их взаимном расположении. В геометрическую ММ обычно входят совокупность уравнений линий и поверхностей и алгебраические соотношения, определяющие принадлежность областей пространства телу ТО или его элементам. Такую ММ иногда задают координатами некоторого множества точек, по которым интерполированием можно построить ограничивающие область линии или поверхности. Границы области задают и кинематическим способом: линию – как траекторию движения точки, а поверхность – как результат перемещения линии. Представление формы и размеров области совокупностью типовых фрагментов достаточно простой конфигурации. Такой способ характерен, например, для метода конечных элементов, широко используемого в математическом моделировании.
Функциональные ММ состоят из соотношений, связывающих между собой фазовые переменные, т.е. внутренние, входные и выходные параметры ТО. Функционирование сложных ТО нередко удается описать лишь при помощи совокупности его реакций на некоторые известные (или заданные) входные воздействия (сигналы). Такую разновидность функциональной ММ относят к типу черного ящика и обычно называют имитационной математической моделью, имея в виду, что она лишь имитирует внешние проявления функционирования ТО, не раскрывая и не описывая существа протекающих в нем процессов. По форме представления имитационная ММ является примером алгоритмической математической модели, поскольку связь в ней между входными и выходными параметрами ТО удается описать лишь в форме алгоритма, пригодного для реализации в виде ЭВМ-программы. По этому признаку к типу алгоритмических относят более широкий класс как функциональных, так и структурных ММ. Если связи между параметрами ТО можно выразить в аналитической форме, то говорят об аналитических математических моделях. При построении иерархии ММ одного и того же ТО обычно стремятся к тому, чтобы упрощенный вариант ММ был представлен в аналитической форме, допускающей точное решение, которое можно было бы использовать для сравнения при тестировании результатов, полученных при помощи более полных и поэтому более сложных вариантов ММ.
Теоретические и эмпирические модели. По способу получения ММ делят на теоретические и эмпирические. Первые получают в результате изучения свойств ТО и протекающих в нем процессов, а вторые являются итогом обработки результатов наблюдений внешних проявлений этих свойств и процессов. Один из способов построения эмпирических ММ заключается в проведении экспериментальных исследований, связанных с измерением фазовых переменных ТО, и в последующем обобщении результатов этих измерений в алгоритмической форме или в виде аналитических зависимостей. Поэтому эмпирическая ММ по форме представления может содержать признаки как алгоритмической, так и аналитической математической модели.
При построении теоретических ММ прежде всего стремятся использовать известные фундаментальные законы сохранения таких субстанций, как масса, электрический заряд, энергия, количество движения и момент количества движения. Кроме того, привлекают определяющие соотношения, (называемые также уравнениями состояния), в роли которых могут выступать так называемые феноменологические законы (например, уравнение Клапейрона – Менделеева, закон Ома, закон Гука, закон Фурье и т.п.).
Сочетание теоретических соображений качественного характера с обработкой результатов наблюдения внешних проявлений свойств изучаемого ТО приводит к смешанному типу ММ, называемых полуэмпирическими. При построении таких ММ используют основные положения теории размерностей, в том числе так называемую Пи-теорему: если между параметрами, характеризующими изучаемый объект, существует зависимость, имеющая физический смысл, то эту зависимость можно представить в виде – k зависимостей их безразмерных комбинаций, где k – число независимых единиц измерения, через которые можно выразить размерности этих параметров. При этом величина – k определяет число независимых (не выражаемых друг через друга) безразмерных комбинаций, обычно называемых критериями подобия.