(нижче джерела) розподіл концентрації виглядає наступним чином.
Рисунок А.1? Розподіл концентрації при y =0
При у =0.6, тобто вище джерела розподіл концентрації наступне:
Рисунок А.2? Розподіл концентрації при y =0.6
Другий приклад. В області D={0 lt; x lt; 4, 0 lt; y lt; 2}. Джерело розташовується в точці М (0,1; 0,8). Так само належало: U=V=0.5, k1=0.7, k2=0.4, my=0.6, ny=0.3. Крок по просторовій змінної вибирався рівний 0,2, а з тимчасової 0,01.
При вирішенні цього завдання методом покомпонентного розщеплення були отримані наступні результати при Т=3. При y=0.4 (нижче джерела) розподіл концентрації виглядає наступним чином:
Малюнок А.3? Розподіл концентрації при y=0.4
При у=0.6, тобто вище джерела розподіл концентрації наступне:
Малюнок А.4? Розподіл концентрації при y=0.6
Третій приклад. В області D={0 lt; x lt; 2, 0 lt; y lt; 2}. Джерело розташовується в точці М (1; 0,5). Так само належало: U=- 0.25, V=0.25, k1=1, k2=0.5, my=1, ny=0.5. Крок по просторовій змінної вибирався рівний 0.1, а з тимчасової 0.01.
При вирішенні відмінюванні завдання були отримані наступні результати при Т=2. При y=0,1 і при y=1.
Малюнок А.5? Функція чутливості при y=0.1
Малюнок А.6? Функція чутливості при y=1
Четвертий приклад. В області D={0 lt; x lt; 2, 0 lt; y lt; 2}. Джерело розташовується в точці М (1; 0,5). Так само належало: U=0, V=0 (штиль), k1=10, k2=5, my=2, ny=1. Крок по просторовій змінної вибирався рівний 0.1, а з тимчасової 0.01.
При вирішенні цього завдання були отримані наступні результати при Т=1. При y=0,1 і при y=0,5.
Малюнок А.7? Розподіл концентрації при y=0.1
Малюнок А.8? Розподіл концентрації на рівні джерела
Додаток Б
Лістинг ПРОГРАМИ
const
n1=400;
delta=0.001;=201;=array [0. n1] of Extended;=array [0. n1,0. n1] of Extended;=array [0. nr] of Extended;, M, T, k1, k2, U, V, my, ny, h, tao, cons, tc: Extended;
i, j, k, n, z: integer;, fipr, fi1: matr;
a1, a2, a3, b, pr: vect;, xr, ydr, fr, farr, ftr, ar, br, cr, dr, er, ffr: vectr; func (o, s:integer): real; (o=strtoint (frmmain. LabeledEdit14. Text)) and
(s=strtoint (frmmain. LabeledEdit15. Text)) then:=cons * consfunc:=0 ;;
procedure Progonka (var a1, a2, a3, b, pr: vect) ;, pr2: vect;
i1: integer; [n - 1]:=-a1 [n]/a2 [n]; [n - 1]:=b [n]/a2 [n]; i1:= n - 1 downto 1 do [i1-1]:=-a1 [i1]/(pr1 [i1] * a3 [i1] + a2 [i1]); [i1-1]:=(b [i1] - pr2 [i1] ??* a3 [i1])/(pr1 [i1] * a3 [i1] + a2 [i1]);
end; [0]:=(b [0] - a3 [0] * pr2 [0])/(a2 [0] + a3 [0] * pr1 [0]);
for i1:=1 to n do [i1]:=pr1 [i1-1] * pr [i1-1] + pr2 [i1-1];
end; Prog5 (var nr: integer; ar, br, cr, dr, er, fr, yr: vectr; itask: integer);
var: integer; itask=1 then
begin [1]:=dr [1]/cr [1]; [1]:=er [1]/cr [1]; [2]:=cr [2] - dr [1] * br [2]; [2]:=(dr [2] - er [1] * br [2])/cr [2];
er [2]:=er [2]/cr [2]; i:=3 to nr do [i]:=er [i - 2] * ar [i] + dr [ i - 1] * (dr [i - 2] * ar [i] - br [i]); [i]:=(dr [i] + er [i - 1] * (dr [i - 2] * ar [i] - br [i]))/cr [i]; [i]:=er [i]/cr [i] ;;:=2 ;; [1]:=fr [1]/cr [ 1]; [2]:=(fr [2] + fr [1] * br [2])/cr [2]; i:=3 to n do [i]:=(fr [i] - fr [ i - 2] * ar [i] - fr [i - 1] * (dr [i - 2] * ar [i] - br [i]))/cr [i]; [nr]:=fr [nr ]; [nr - 1]:=dr [nr - 1] * yr [nr] + fr [nr - 1]; i:=nr - 2 downto 1 do [i]:=dr [i] * yr [i +1] - er [i] * yr [i + 2] + fr [i] ;; Prog3 (var nr: integer; ar, cr, br, fr, yr: vectr; itask: integer) ;: integer; itask =1 then [1]:=br [1]/cr [1]; i:=2 to nr do [i]:=cr [i] - br [i - 1] * ar [i]; [i] :=br [i]/cr [i] ;;:=2 ;; [1]:=fr [1]/cr [1]; i:=2 to n do [i]:=(fr [i] + fr [i - 1] * ar [i])/cr [i]; [nr]:=fr [nr]; i:=nr - 1 downto 1 do [i]:=br [i] * yr [ i + 1] + fr [i] ;; TfrmMain. cmdRunClick (Sender: TObject) ;:=StrToFloat (LabeledEdit1. Text) ;:=StrToFloat (LabeledEdit2. Text) ;:=StrToFloat (LabeledEdit3. Text) ;:=StrToFloat (LabeledEdit4. Text) ;:=StrToFloat (LabeledEdit5. Text );
U:=StrToFloat (LabeledEdit6. Text) ;:=StrToFloat (LabeledEdit7. Text);
my:=StrToFloat (LabeledEdit8. Text) ;:=StrToFloat (LabeledEdit9. Text) ;:=StrToInt (LabeledEdit10. Text) ;:=StrToFloat (LabeledEdit11. Text) ;:=StrToFloat (LabeledEdit16. Text) ;:=L/n cmdRun. Enabled:=false; i:=1 to n - 1 for j:=1 to n - 1 do [i, j]:...
