Программа для моделирования энергомассопереноса
Рассмотрим подробно блок-схему программы (см. рис. 13.10) для моделирования процессов энергомассопереноса в связной постановке.
1. Первым этапом работы пользователя с программой является ввод размеров для построения геометрической модели. В современных программных комплексах МКЭ этот этап проходит в интерактивном режиме с использованием ЗБ-графики, позволяющей рассмотреть построенные части сложной модели и убедиться в правильности введенных данных. Из построенной геометрической модели (чертежа) программа извлекает координаты узлов и размеры элементов, из которых состоит модель.
2. Программа обеспечивает автоматическую генерацию сетки конечных элементов. Геометрические параметры модели включают две группы данных:
- размеры Ins каждого конечного элемента, позволяющие рассчитать его электрическое сопротивление R по формуле (13.4), а также аналогичные параметры для тепловой и диффузионной моделей;
- данные для включения модели элемента в общую модель (для каждого элемента R это пара номеров узлов (см. рис. 13.5), к которым он подключен).
3. Для моделирования нестационарных процессов пользователь должен задать начальное состояние, т. е. температуру, фазовый состав материала и содержание водорода (они могут быть различными в разных частях модели). Для моделирования стационарных процессов этот этап не нужен.
4. После завершения подготовительных этапов 1-3 программа переходит к моделированию комплекса физических процессов. Весь период времени разбивается на шаги, обеспечивающие необходимую точность моделирования, и начинается цикл по шагам.
5. В начале каждого шага пользователь вводит граничные условия (если они изменялись) для каждого процесса.
6. По известным значениям температуры и остальных параметров в начале шага программа обеспечивает расчет свойств материала, необходимых для моделирования всех процессов. Затем на -
уравнению (2.50) для Na дает % ~ 0,21; для Аг получаем % ~
_4
- 0,23 10 . Следовательно, степень ионизации Аг по сравнению с
4
Na меньше в 10 раз.
