ÐÐТÐÐÐТÐЧÐÐ ÐÐÐÐÐЮÐÐÐÐЯ
ÐÐТÐÐÐТÐЧÐÐ ÐÐÐÐÐЮÐÐÐÐЯ
ÐÐТÐÐÐТÐЧÐÐ ÐÐÐÐÐЮÐÐÐÐЯ
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