В.А. Хямяляйнен Установившаяся двумерная фильтрация жидкости вокруг перемычки

20

ПРИЛОЖЕНИЕ

Листинг 1. Программа расчета распределения давления на области фильтрации.

function pdemodel [pde_fig,ax]=pdeinit; pdetool('appl_cb',1); set(ax,'DataAspectRatio',[35 7 1]); set(ax,'PlotBoxAspectRatio',[2 2 1]); set(ax,'XLim',[-30 40]); set(ax,'YLim',[0 14]); set(ax,'XTickMode','auto'); set(ax,'YTickMode','auto');

% Geometry description:

pderect([-20 30 2.2000000000000002 12],'R1'); pderect([-20 30 2.2000000000000002 3],'R2'); pdepoly([ -5.5 -5.5 -2.5 -2.5 2.5 2.5 5.5 5.5],… [2.2 5 5 7 7 5 5 2.2],'P1');

pdepoly([ -2.5 -2.5 0 0 2.5 2.5],...

[ 2.2 3 3.5 3.2 3.5 2.2], 'P2'); set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String',… 'R1+R2+P1+P2')

% Boundary conditions: pdetool('changemode',0) pdesetbd(19,'dir',1,'1','1500000') pdesetbd(18,'dir',1,'1','1500000') pdesetbd(17,'neu',1,'0','0') pdesetbd(16,'dir',1,'1','0') pdesetbd(15,'dir',1,'1','0') pdesetbd(9,'dir',1,'1','0') pdesetbd(8,'dir',1,'1','0') pdesetbd(7,'dir',1,'1','1500000') pdesetbd(6,'dir',1,'1','1500000') pdesetbd(1,'neu',1,'0','0')

21

%Mesh generation: setuprop(pde_fig,'Hgrad',1.3); setuprop(pde_fig,'refinemethod','regular'); pdetool('initmesh')

pdetool('refine')

%PDE coefficients:

pdeseteq(1,'y*0.1*10^-12!y*0.01*10^-12!,…

y*0.01*10^-12!y.*100.*exp(-y.*0.5)*10^-12!,…

y*0.01*10^-12!y*0.01*10^-12!y*0.01*10^-12!y*0.01*10^-12',...

'0.0!0.0!0.0!0.0!0.0!0.0!0.0!0.0','0!0!0!0!0!0!0!0',…

'1.0!1.0!1.0!1.0!1.0!1.0!1.0!1.0','0:10','0.0','0.0','[0 100]')

setuprop(pde_fig,'currparam',...

['y*0.1*10^-12!y*0.01*10^-12!y*0.01*10^-12!,…

y.*100.*exp(-y.*0.5)*10^-12!y*0.01*10^-12!,…

y*0.01*10^-12!y*0.01*10^-12!y*0.01*10^-12';...

'0.0!0.0!0.0!0.0!0.0!0.0!0.0!0.0';...

'0!0!0!0!0!0!0!0';...

'1.0!1.0!1.0!1.0!1.0!1.0!1.0!1.0 '])

%Solve parameters: setuprop(pde_fig,'solveparam',...

str2mat('0','3564','10','pdeadworst',...

'0.5','longest','0','1E-4','','fixed','Inf'))

%Plotflags and user data strings:

setuprop(pde_fig,'plotflags',[1 1 1 1 2 1 1 1 0 0 0 1 1 0 0 1 0 1]); setuprop(pde_fig,'colstring',''); setuprop(pde_fig,'arrowstring',''); setuprop(pde_fig,'deformstring',''); setuprop(pde_fig,'heightstring','');

% Solve PDE: pdetool('solve')

%//// Конец листинга 1.

Составители:

Вениамин Анатольевич Хямяляйнен и др.

УСТАНОВИВШАЯСЯ ДВУМЕРНАЯ ФИЛЬТРАЦИЯ ЖИДКОСТИ ВОКРУГ ПЕРЕМЫЧКИ

Методические указания по выполнению лабораторной работы по теме "Численное моделирование перетока воды и газа через тампонажную завесу" для студентов специальности 070600 "Физические процессы горного производства"

Редактор Е.Л. Наркевич

ИД №06536 от 16.01.02

Подписано в печать 26.09.02. Формат 60×84/16. Бумага офсетная. Печать офсетная. Уч.-изд. л.1,5. Тираж 50 экз. Заказ

ГУ Кузбасский государственный технический университет. 650026, Кемерово, ул. Весенняя, 28.

Типография ГУ Кузбасский государственный технический университет.

650099, Кемерово, ул. Д. Бедного, 4А.