Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
8
Добавлен:
20.02.2016
Размер:
89.6 Кб
Скачать

Програма pre1f77

Mоделювання роботи рідинного пористого невитратного електрода

  1. Обєкт і мета розрахунків

1.1 Фізичний обєкт, який моделюється – пористий електрод, товщиною L і габаритною площею 1 см2, виготовлений із електропровідного інертного матеріалу, який не приймає участь в електрохімічних реакціях. На внутрішній поверхні електрода, який може мати деякий ненульовий питомий опір, здійснюється стаціонарний (безперервний) електрохімічний процес з заданою поляризаційною характеристикою. Враховується лише електрохімічна поляризація, задана апроксимаційною формулою.

Мета роботи-

1-вивчення закономірностей роботи невитратних рідинних пористих електродів, впливу їх структури на характер розподілу інтенсивності електрохімічного процесу в глибині пор і загальну поляризаційну характеристику пористого електрода.

2-ознайомлення з алгоритмом вирішення задачі (системи рівнянь Пуассона)

1.3 Програма написана алгоритмічною мовою ФОРТРАН, має дуже просту структуру, з введенням вхідних даних безпосередньо в текст. Можливість змінювати текст програми дозволяє вирішувати задачі, в яких сам алгоритм змінюється (наприклад, можна змінювати форму локальної поляризаційної характеристики), або доповнюється додатковими елементами.

Програма розраховує стаціонарний режим електрохімічного процесу на внутршній поверхні електрода при постійній заданій поляризації. Розраховуються та виводяться в файл результатів такі (головні) дані:

  • таблиця розподілу густини струму в порах

  • загальні характеристики процесу – габаритний струм, ефективність використаняя внутрішньої поверхні, потенціали фаз фронтальної та тильної поверхні.

C ЭЛЕКТРИЧЕСКОЕ ПОЛЕ В ПОРИСТОМ ЭЛЕКТРОДЕ ПРОГРАММА PRE1f77.for

C...._______________________________________________________________

REAL L

DIMENSION T(40),RP(40),U(40),V(40)

NN=40

DATA RM, RE,R, B1,B2, L, EPS, ALF, T,RP, U, V

* / 0, 10,0.005, 10, 5, 0.2, 0.00001, 1.5,80*1.,40*0.2, 40*0./

C----------------------------------------------------------

OPEN (UNIT=3,FILE='REZ.REZ', STATUS='UNKNOWN')

WRITE (3,7)

7 FORMAT(' UDXTU____PRE1F77___TE_________________________________')

WRITE (3,8)

8 FORMAT

* (' __RM,OM*CM______RE_____U1____B1____B2_____L___ALF_______EPS')

WRITE (3,9) RM, RE, U(1), B1, B2, L, ALF, EPS

9 FORMAT ( E11.2,F8.2,F7.3,4F6.2,F10.6)

SU = 0.5/R

TS = 10000.

DX = L/(NN-1)

CL----------------------------------------------------------------

DO 1 ITER = 1,2000

DO 2 N = 2, NN-1

UN = (U(N+1) + U(N-1) - T(N)*SU*RE*DX*DX)/2.

VN = (V(N+1) + V(N-1) + T(N)*SU*RM*DX*DX)/2.

U(N) = UN*ALF + U(N)*(1.-ALF)

2 V(N) = VN*ALF + V(N)*(1.-ALF)

U(NN)= U(NN-1)

V( 1)= V( 2)

TG = 0.

DO 22 N = 1,NN

T(N) = (U(N) - V(N))/RP(N)

PP = B1*T(N) + B2*(T(N)**2)

PP = B1*T(N) + B2*(T(N)**2)

RP(N)= PP/T(N)

22 TG = TG +T(N)*SU*DX

EFF = TG/(T(1)*SU*DX*(NN-1))

EP=ABS((TG-TS)*2/(TS+TG))

PRINT 52, ITER,EP,TG

52 FORMAT (I5,F10.6,F10.3)

IF(EP.LT.EPS) GOTO 3

1 TS = TG

CL--------------------------------------------------------------------

3 CONTINUE

WRITE (3,6)

6 FORMAT('__ITER_____Rcm__SU.cm2/cm3__GAB,A/CM2___EFF_____EPS')

WRITE (3,4) ITER, R, SU, TG, EFF, EP

4 FORMAT( I6, F8.5, F12.1, F11.6, F6.3, F8.5)

WRITE (3,41)

41 FORMAT (1X,'___________________________________________________')

WRITE (3,50)

50 FORMAT ( '____N_i(N),A/sm2___U(N),B___V(N),B')

DO 51 I=1,NN,3

WRITE (3,5) I, T(I), U(I),V(I)

5 FORMAT( I5,F10.3,2F9.2)

51 CONTINUE

STOP

END

Соседние файлы в папке Материалы по электра