Понятие моделирования Моделирование – построение образа реальной системы, учитывающего только интересующие нас стороны ее поведения (упрощённое представление, принимающее во внимание только главные её особенности) Цель: исследовать объект/явление, понять его и научиться предсказывать результаты будущих наблюдений, а также помочь в управлении (либо показать пути совершенствования структуры и свойств объекта) При построении модели объекта необходимо понять: • Какие стороны объекта нас интересуют (или нам достаточно ими ограничиться), поскольку реальный объект многогранен. • Какие процессы/явления относятся к этой стороне объекта. • Какие соотношения (либо физические/аналоговые модели) описывают данные процессы/явления
Физическое моделирование Проведение исследований над объектом, аналогичным реальному, но уменьшенным или увеличенным и допускающим исследование с помощью последующего перенесения свойств изучаемого процесса или явления с модели на объект на основе теории подобия Это экспериментальный метод изучения Применение физического моделирования: • Исчерпывающе точного математического описания (модели) явления/объекта в данный момент не существует, либо оно слишком громоздко • Воспроизведение физического явления в целях изучения в реальных масштабах невозможно, нежелательно или слишком дорогостояще Метод может дать достоверные результаты лишь в случае соблюдения физического подобия реального явления и модели. Подобие достигается за счёт равенства для модели и реального явления значений критериев подобия – безразмерных чисел, зависящих от физических (в том числе геометрических) параметров, характеризующих явление
Теория подобия – учение об условиях подобия физических явлений Физические явления, процессы или системы подобны, если в сходственные моменты времени в сходственных точках пространства значения переменных величин, характеризующих состояние одной системы, пропорциональны соответствующим величинам другой системы. Коэффициенты пропорциональности для каждой из величин называются критериями подобия. Виды подобия: • Геометрическое (пропорциональность элементов фигур или тел) сходственных геометрических • Кинематическое подобие (пропорциональность полей скоростей) • Динамическое подобие (пропорциональность систем действующих сил или силовых полей различной физической природы) • Механическое подобие (потоков жидкости/газа, упругих систем) предполагает подобие геометрическое+кинематическое+динамическое • Тепловое (пропорциональность полей температур и тепловых полей) • Электродинамическое (пропорциональность мощностей, электромагнитных сил) полей токов, нагрузок, • Другие (процессов трения и износа деталей машин и мех-мов, кинетики физико-химических превращений, …)
Теория подобия – учение об условиях подобия физических явлений Пропорциональность для подобных явлений характеризующих их параметров приводит к тому, что все безразмерные комбинации этих параметров имеют для подобных явлений одинаковые численные значения Безразмерные комбинации, составленные из определяющих параметров рассматриваемых явлений, называются критериями подобия Чаще всего встречается неполное (частичное подобие) – равенство не всех, а лишь некоторых независимых критериев подобия. Важно, чтобы влияние на процесс/явление критериев, равенство которых не соблюдается, было малозначительным. Размерные физические параметры, входящие в критерии подобия, для подобных систем могут различаться на порядки величины, но одинаковыми должны оставаться безразмерные критерии подобия. Это свойство подобных систем составляет основу моделирования.
Теория подобия – учение об условиях подобия физических явлений Определение критериев подобия: • Приведение определяющих уравнений для системы/явления, к безразмерному виду (критерии подобия – безразмерные коэффициенты при некоторых членах безразмерных уравнений) • Из анализа размерностей, определяющих физические параметры Примеры: Число Рейнольдса Отношение инерционных сил при движении ж/г к силам вязкости Число Маха До/сверхзвуковой режим течения Число Фруда Отношение инерционных сил к силам тяжести Число Прандтля Влияние физических свойств теплоносителя на теплопередачу Число Нуссельта Отношение тепловых потоков конвекции и теплопередачи
Теория подобия – учение об условиях подобия физических явлений Примеры: Число Грасгофа Мера соотношения между архимедовой силой, вызванной неравномерным распределением плотности в неоднородном поле T, и силами межмолекулярного трения Число Пекле Хар-ка соотношения между конвекцией и диффузией: между конвективным и молекулярным процессами переноса тепла (примесей, импульса, характеристик турбулентности) в потоке жидкости Число Стэнтона Характеризует интенсивность диссипации энергии в потоке жидкости или газа
Важнейшие аспекты возникновения и функционирования зон поддвига, актуальные для зон субдукции в литосфере 1. Возникновение на существующих ослабленных границах раздела в результате изменения регионального или глобального напряженного состояния среды Давление толщи осадков Gerya et al. , 2008 Regenauer-Lieb et al. , 2001 «Подготовленные» границы раздела могут формироваться в молодом ледовом покрове и функционировать как дивергентные или трансформные в течение длительного времени
Важнейшие аспекты возникновения и функционирования зон поддвига, актуальные для зон субдукции в литосфере 2. Необходимым условием инициирования субдукции является достаточно мощный «импульс сжатия» плит на «подготовленной» границе раздела, характеризующейся низкой прочностью Моделирование методом МСА Hall et al. , 2003 «Импульс сжатия» может быть инициирован в результате исчерпания релаксационной способности деформационных механизмов низких рангов и быстрого роста напряжений в среде
Образование зоны поддвига как деформационной структуры субдукционного типа Эволюция границы раздела Образование пододвигающегося блока Vload
Важнейшие аспекты возникновения и функционирования зон поддвига, актуальные для зон субдукции в литосфере 3. Субдукция «зрелых» ледовых плит, как и океанических литосферных плит, является односторонней Gerya et al. , 2008 Gerya, Yuen, 2007 Причины: относительно высокая прочность пластин и относительно низкая прочность межплитного интерфейса, обводнение поверхности проскальзывания Надводная часть Подводная часть
Важнейшие аспекты возникновения и функционирования зон поддвига, актуальные для зон субдукции в литосфере 4. Крайне низкий «коэффициент и надвигающейся плит: трения» пододвигающейся ~0, 01 в случае литосферных плит (Sobolev S. V. , 2005, 2006) ~0, 01 – 0, 05 в случае ледовых пластин (2 D моделирование методом МСА) Фотографии подводной части зоны поддвига (вид вкрест линии границы раздела) Надвигающаяся плита Пододвигающийся блок Надвигающаяся плита Пятна контакта Пододвигающийся блок
Важнейшие аспекты возникновения и функционирования зон поддвига, актуальные для зон субдукции в литосфере 5. Развивающаяся фрагментация пододвигающихся плит, способствующая уменьшению сопротивления изгибу Ranero et al. , 2003, 2005 В процессе пододвигания ледовая плита испытывает значительный изгиб, что приводит к образованию в них серий вертикальных или диагональных трещин
Анализ подобия деформационных процессов в литосфере и ледовом покрове с использованием критериев подобия В приближении литосферы и астеносферы как Ньютоновских несжимаемых и медленно текущих жидкостей: Подстановка: - опорные значения - безразмерные параметры Системы с одинаковым значением безразмерного параметра описываются одним и тем же уравнением и являются подобными с точки зрения деформационных процессов
Анализ подобия деформационных процессов в литосфере и ледовом покрове с использованием критериев подобия Критерии подобия: Значения для литосферы: Значения для ледового покрова:
Аналоговое моделирование Основано на аналогии процессов и явлений, имеющих различную физическую природу, но одинаково описываемых формально (одинаковыми математическими соотношениями, логическими и структурными схемами) Пример: • Одна система механическая (ось передаёт вращения через пружину и маховик, частично погруженный в вязкую жидкость, валу, жёстко связанному с маховиком) • Вторая система электрическая (источник ЭДС соединён через катушку индуктивности, конденсатор и активное сопротивление со счётчиком электроэнергии) ОДУ 2 -го порядка: Другие примеры: электротепловое и гидротепловое моделирование (для изучения процессов теплопроводности), световое моделирование (для изучения радиационного переноса тепла) До возникновения цифровых компьютеров аналоговое моделирование было основным способом «предметно-математического моделирования» процессов распространения звуковых и электромагнитных волн, движения газов и жидкостей, фильтрации, кручения стержней и т. д.
Интуитивное моделирование Основано на интуитивном (не обоснованном формально) представлении об объекте исследования, не поддающимся формализации или не нуждающемся в ней Проводится на вербальном уровне. Не устанавливает строгие количественные соотношения между моделируемыми явлениями, ограничиваясь лишь анализом качественных обобщённых понятий, отражающих общие тенденции развития явлений, направления изменения свойств изучаемых объектов и т. д. Широко применяется учёными-экспериментаторами с целью выдвижения различного рода гипотез поведения субъектов сложных систем, формирования эвристик относительно взаимоотношений между активными элементами системы и их развития
Научное моделирование Логически обоснованное моделирование, использующее минимальное число предположений, принятых в качестве гипотез на основании наблюдений за объектом моделирования Математическое моделирование Моделирование на основе математического представления выделенных процессов и явлений в объекте в виде совокупности алгебраических, интегральных или дифференциальных уравнений и их совместного решения методами математики Виды: аналитическое и численное (компьютерное) Математические модели: размерные и безразмерные функциональные и структурные детерминистические и вероятностно-статистические
Математическое моделирование Основные этапы: 1. Построение модели. Определяется «нематический» объект исследования (явление или процесс, материал или конструкция). Далее выявляются основные особенности явления и связи между ними, отбрасываются несущественные для поставленной цели. Найденные качественные зависимости формулируются на языке математики – строится математическая модель. 2. Решение математической задачи, к которой приводит построенная модель. На этом этапе разрабатываются алгоритмы и численные методы решения задачи на компьютере, обеспечивающие получение результата с приемлемой точностью и за приемлемое время. 3. Интерпретация полученных решений (цифр, зависимостей, …) на языке, принятом в данной предметной области. 4. Проверка адекватности модели. Проверяется, согласуются ли теоретические следствия из модели с данными эксперимента, теоретическими работами других авторов и т. д. 5. Модификация модели. Производится коррекция модели (усложнение модели с целью получения более адекватного результата, либо упрощение для достижения практически приемлемого решения).
Проблема B. D. Wirth et al. Journal of Nuclear Materials V. 329– 333 (2004) P. 103
Проблема
Проблема Можно ли в рамках единого формализма моделировать материалы/среды во всем спектре масштабов от нано до гео?
Два подхода к описанию реальных сред В 1822 и 1823 году Коши и Навье представили в Парижскую академию наук мемуары, положившие начало двум подходам к описанию материалов и сред различной природы: Cauchy (1822) континуальный подход – механика сплошных сред «дискретный подход» – система взаимодействующих элементов Вплоть до середины ХХ века развивался преимущественно континуальный подход. Он открыл эпоху эффективных аналитических исследований.
Некоторые уравнения механики сплошных сред Пространственное описание (формализм Эйлера) Субстанциальное описание (формализм Лагранжа)
Континуальный и дискретный подходы С возникновением и развитием вычислительных устройств самое широкое применение получили вычислительные методы. Новое научное направление: вычислительная механика Основные численные методы континуального подхода: • Метод конечных разностей • Метод конечных элементов С развитием вычислительной техники был дан толчок к развитию вычислительных методов дискретного подхода. Дискретный подход начал развиваться как отдельное научное направление.
Основные численные методы континуального подхода Метод конечных элементов: численный метод решения дифференциальных уравнений в частных производных или интегральных уравнений, возникающих в прикладной физике Метод конечных разностей: метод интерполяции, заключающийся в замене дифференциальных коэффициентов уравнения их разностными аналогами
РАЗВИТИЕ ЧИСЛЕННЫХ МЕТОДОВ ДИСКРЕТНОГО ПОДХОДА Метод молекулярной динамики пионерный численный метод дискретного подхода в механике Применение метода молекулярной динамики для описания систем на атомном масштабном уровне: 1. Изучение механизмов зарождения пластической деформации в металлах и сплавах при наноиндентировании. 2. Изучение генерации и эволюции дефектов структуры в каскадах атомных соударений в материалах при радиационном воздействии. 3. Моделирование и изучение свойств наноразмерных структур, сформированных на основе тонких металлических пленок. Применение метода молекулярной динамики ограничивается объектами/фрагментами пространственного масштаба ~1 мкм.
Метод частиц: нано- до макро- уровня Редукция размерности фазового пространства На мезо уровне размерность фазового пространства - 6 Nm (Nm=N 0/m; m > 1010 ) мезо-элементы (методы частиц, КА) нано-уровень (методы МД, Монте-Карло) На атомном уровне размерность фазового пространства - 6 N 0 (N 0 – число атомов)
Основы формализма методов дискретного подхода Концепция: моделируемый объект рассматривается как ансамбль взаимодействующих фрагментов (элементов, частиц, автоматов) конечного размера и заданной исходной формы Structure Velocity field DEM model of Lac du Bonnet granite (Potyondy, Cundall, 2004) Базовым элементом математических моделей в методах дискретного подхода является взаимодействие пары фрагментов/частиц/автоматов. Принцип локальности (масштабный фактор)!!!
ОСНОВНЫЕ КЛАССЫ ЧИСЛЕННЫХ МЕТОДОВ ДИСКРЕТНОГО ПОДХОДА В настоящее время численные методы дискретного подхода являются общепризнанными и используются в самых различных отраслях фундаментальной и инженерной науки. Два основных класса методов дискретного описания среды: 1. Методы клеточных автоматов применяются для изучения термодинамических процессов, связанных с изменением фазового состава материала на фиксированной в пространстве сетке автоматов. 2. Методы частиц (DEM, particle-in-cell, SPH etc. ) используются преимущественно для моделирования механического отклика среды. 3. «Гибридные» методы (MCA) применяются для решения связанных задач, включающих в рассмотрение как процессы деформирования среды, так и химические реакции и фазовые превращения.
ОСНОВНЫЕ КЛАССЫ ЧИСЛЕННЫХ МЕТОДОВ ДИСКРЕТНОГО ПОДХОДА Классическая концепция клеточных автоматов - внутренние источники - взаимодействие автоматов Классическая концепция частиц (концепция дискретных элементов): Явные методы (discrete element method) Неявные методы (discontinuous deformation analysis) MD + PM PM Cundall, Strack, 1979 Shi, Goodman, 1984 -1989
УРАВНЕНИЯ ДВИЖЕНИЯ В МДЭ Пространственная эволюция автоматов (частиц) определяется уравнениями движения Ньютона-Эйлера: + уравнения: а) теплопроводности; б) дифузии; в) химические реакции… Силы взаимодействия включают потенциальную и диссипативную составляющие: Диссипативные силы компонентами вектора скорости: определяются относительной Модель гранита Lac du Bonnet
Скачать презентацию Понятие моделирования Моделирование построение образа реальной системы
presentation_introduction_1.ppt