Моделирование структуры ЭГД-течений в несимметричной системе электродов, страница 6

*cfopen,foY              *vwrite,fieldy(1) (F12.6) *cfclos

Аналогично данные из массива fieldy записываются в файл foY.

4. Решение гидродинамической задачи с электростатическими силами в качестве нагрузки осуществляется при помощи программы Flotran. Перед запуском этой программы необходимо загрузить геометрию задачи вместе с построенной сеткой. Ниже приводится описание программы.

et,3,fluid141                             Выбирается новый элемент – fluid141.

esel,all

EMODIF,all,TYPE,3                Замена всех элементов на гидродинамический.

*get,nodemx,node,0,count       Определяется общее количество узлов.

*dim,fox,array,nodemx            Создаются массивы в которые будут записываться значения

*dim,foy,array,nodemx            компонент плотности сил.

*vread,fox(1),fox	Числа из файла fox (X компонента плотности электростатических сил) записываются в массиве fox начиная с первого элемента.
(F12.6)	Формат записи числа
*vread,foy(1),foy (F12.6)	Аналогично Y компонента плотности сил для каждого узла записывается в массив
Lsel,s,,,9,12 dl,all,,VX,0,1 dl,all,,Vy,0,1 Lsel,s,,,2,6,2 lsel,a,,,1 dl,all,,VX,0,1 dl,all,,Vy,0,1	На стенках кюветы и на поверхности электродов задаются условия «прилипания» (нулевые значения компонент скорости).

nsel,all *do,n,1,nodemx,1	Создается цикл по всем узлам
bf,n,forc,fox(n),foy(n)	В каждом узле задаются компоненты сил, которые берутся из массивов fox и foy. Формат команды: в узле с номером n задается плотность сил (ключ forc), fox(n) – значение X компоненты плотности сил, foy(n) - значение Y компоненты плотности сил. Они равны n – м элементам массивов fox и foy.
*enddo	Конец цикла
allsel,all FINISH
/solu	Переход в решение.

Далее следует настройка параметров решения и задание свойств жидкости.

FLDATA1,SOLU,flow,t                        Решение уравнений для скоростей и давления.

FLDATA2,ITER,EXEC,200                 Назначение числа итераций 200

FLDATA2,ITER,OVER,50                   Количество итераций после которых выводятся проме-

                                                               жуточные результаты.

fldata3,term,vx,0.001                            Устанавливается уровень сходимости для X компонен-

                                                              ты скорости.

fldata3,term,vy,0.001                            Аналогично для Y компоненты.

FLDATA12,PROP,DENS,2                

FLDATA13,VARY,DENS,1                Разрешается вариация плотности

FLDATA12,PROP,VISC,0

FLDATA13,VARY,VISC,0                  Запрет вариации вязкости

FLDATA7,PROT,DENS,LIQUID        Определяет тип свойств жидкости как Liquid

FLDATA8,NOMI,DENS,800                Начальное значение плотности 800 кг/м3

FLDATA9,COF1,DENS,273                 Определяет значение температуры при которой

 FLDATA10,COF2,DENS,-0.001  

FLDATA11,COF3,DENS,0  

FLDATA7,PROT,VISC,CONSTANT       Вязкость постоянная и равна 0.01 

FLDATA8,NOMI,VISC,0.01

FLDATA18,METH,PRES,1                      Определяется тип используемого решателя

FLDATA19,TDMA,PRES,100 

SOLVE                                                       Запуск решения гидродинамической задачи

FINISH                                                       

/POST1                                                       Переход в постпроцессор

SET,LAST                                                  Считывание результатов из последнего набора.