30.10.2014 Views

МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ

МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ

МАТЕМАТИЧНЕ МОДЕЛЮВАННЯ

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

В.Є. Бахрушин. Математичне моделювання<br />

Розглянемо як приклад математичну модель дифузії домішки у твердому<br />

тілі. Вихідну модель у разі вільної дифузії можна записати в такому<br />

вигляді:<br />

<br />

⎧ jn<br />

=−D(r,t) ⋅grad n(r,t);<br />

⎪ ∂ n <br />

⎪ = div( D(r,t)grad n(r,t) ) ;<br />

⎨ ∂ t<br />

<br />

⎪ n(r<br />

0,t)<br />

=ϕ(t);<br />

⎪ <br />

⎪ ⎩n(r,t) 0<br />

=ψ(r) ,<br />

(1.13)<br />

де r – радіус-вектор.<br />

Модель 1.13 не має аналітичного розв’язку. Але за деяких умов вона<br />

може бути перетворена до моделей, що можна досліджувати аналітично.<br />

Наприклад, однією з типових задач дифузії є одновимірна дифузія з необмеженого<br />

однорідного плоского джерела в напівнескінченне однорідне тіло<br />

при сталому коефіцієнті дифузії. Така модель дифузії широко використовується<br />

при аналізі багатьох процесів хіміко-термічної обробки металів,<br />

процесів створення твердотілих електронних приладів тощо. Для цього<br />

випадку, якщо можна нехтувати присутністю домішки в тілі при t = 0, модель<br />

1.13 набирає вигляду:<br />

⎧ dn<br />

⎪<br />

j = −D<br />

;<br />

dx<br />

⎪<br />

2<br />

⎪dn<br />

d n<br />

⎪ = D ;<br />

2<br />

⎨ dt dx<br />

⎪n( 0,t)<br />

= n 0;<br />

⎪<br />

⎪ n<br />

⎪<br />

n(x,0) =<br />

⎩ 0<br />

0 ( x = 0)<br />

( x ≠ 0)<br />

.<br />

(1.14)<br />

Із системи 1.14 можна одержати таке рівняння, що описує концентрацію<br />

в довільний момент часу в довільній точці тіла:<br />

де<br />

⎛ x ⎞<br />

n(x, t) = n 0 erfc⎜<br />

⎟ , (1.15)<br />

⎝ 2 Dt ⎠<br />

2<br />

erfc (z) 1 e dx .<br />

z<br />

2<br />

−x<br />

= −<br />

π<br />

∫ (1.16)<br />

0<br />

12

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!