Problem 5: O restart; b r O T r d Tb K Tb KTa $ b ln a ln ; b r

Tb KTa ln T := r/Tb K

b ln a 1 u r O eq1 d diff u r , r$2 C $diff u r , r K 2 = 1 C $$diff T r , r ; r r d u r 2 1 C Tb KTa u r d dr eq1 := 2 u r C K 2 = r b dr r r ln a O solve1 d dsolve eq1, u r ; Tb _C2 1 r 1 C K CTa 2 ln r K1 solve1 := u r = r _C1 C K r b 4 ln a 1 C1$r C2 $int T r $r, r = a ..r C C ; O u d 1 C $$ r 2 r 1 1 b b 2 2 2 2 u := 1 C K Ta Ca Tb K2 a Ta ln a C2 Tb r ln b 4 a a r ln a b b 1 C2 2 2 2 2 K2 Tb r ln KTb r C2 Ta r ln CTa r C C1 r C r r 2 r O er d diff u, r ; 1 b b 1 er := K 1 C K 2 Ta Ca2 Tb K2 a2 Ta ln a C2 Tb r2 ln a a b 4 2 r ln a b b 2 2 2 2 K2 Tb r ln KTb r C2 Ta r ln CTa r r r b b b 1 C 4 Tb r ln K4 Tb r ln C4 Ta r ln C2 1 1 a r r C C C1 K 2 b 4 2 r r ln a O eh d eh := 1 r 1 4 u r ; 1 r ln b a 1 C K 2 Ta Ca2 Tb K2 a2 Ta ln a b a C2 Tb r2 ln b a

(1)

(2)

(3)

(4)

(5)

(6)

K2 Tb r ln

b r

KTb r C2 Ta r ln

b r

CTa r

1 C2 C1 r C 2 r

This is the equation for the single differential equation for sigmar: E $ erC$eh K 1 C $$T r ; O sigmar d 2 1 K

sigmar :=

1 1 K
2

E K

1 4

1 b r ln a
2

1 C K 2 Ta Ca2 Tb K2 a2 Ta ln a

b a

(7)

C2 Tb r ln

b a

K2 Tb r ln

b r

KTb r C2 Ta r ln

b r

CTa r

1 4

1 C 4 Tb r ln

b a

K4 Tb r ln r ln b a

b r

C4 Ta r ln

b r

1 C2 C1 K 2 2 r

1 r

1 4

1 r ln b a b a
2

1 C K 2 Ta Ca2 Tb K2 a2 Ta ln a

b a

C2 Tb r ln

K2 Tb r ln

b r

KTb r C2 Ta r ln b r

b r

CTa r

1 C1 r 2

C2 C r O sigmah d 1

Tb KTa ln K 1 C Tb K ln E 1 K
2

b a ;

eh C$er K 1 C $$T r

sigmah :=

1 K

1 r

1 4

1 b r ln a

1 C K Ta Ca Tb K2 a Ta ln a

b a

(8)

C2 Tb r ln

b a

K2 Tb r ln

b r

KTb r C2 Ta r ln

b r

CTa r

1 C1 r 2

C2 r

C K

1 4

1 b 2 r ln a
2

1 C K Ta Ca Tb K2 a Ta ln a

b a

C2 Tb r ln

b a

K2 Tb r ln b a

b r

KTb r C2 Ta r ln b r

b r

CTa r b r

1 C 4

1 C 4 Tb r ln

K4 Tb r ln r ln b a

C4 Ta r ln

1 C2 C1 K 2 2 r

Tb KTa ln K 1 C Tb K ln O sigmaz d $E$ b a

b r

$int T r $r, r = a ..b ; a b2K 2 b Tb KTa ln 1 r sigmaz := E K C Tb b 1 K ln a b b 2 2 2 a2 Ta Ka2 Tb C2 a2 Ta ln K2 b Tb ln Cb Tb Kb Ta 1 a a K b 2 b2 Ka2 ln a O C1 d 1 C1 := K 2 2$ 1 K $ b2 Ka2 $int T r $r, r = a ..b ; b a
2

1 $ K r C T 1 K

2$

(9)

1 K a2 Ta Ka2 Tb C2 a2 Ta ln b ln a 1 C $$a2 b2 Ka2 1 C a

K2 b Tb ln
2

b a

Cb Tb Kb Ta (10)

b Ka

O C2 d C2 := K 1 4

$int T r $r, r = a ..b ; (11) b a ln b a

a Ta Ka Tb C2 a Ta ln b Ka
2 2

K2 b Tb ln b a

Cb Tb Kb Ta

O sigmar;

1 1 K
2

E K

1 4

1 b r ln a
2

1 C K 2 Ta Ca2 Tb K2 a2 Ta ln a

b a

(12)

C2 Tb r ln

b a

K2 Tb r ln

b r

KTb r C2 Ta r ln

b r

CTa r

1 4

1 C 4 Tb r ln

b a

K4 Tb r ln r ln b a

b r

C4 Ta r ln

b r

1 K 4

1 K a2 Ta Ka2 Tb C2 a2 Ta ln ln b a

b a

K2 b Tb ln

b a

Cb Tb Kb Ta

b2 Ka2

1 4

1 b Ka
2 2

ln

b a
2

1 C a

a Ta Ka Tb C2 a Ta ln

b a

K2 b Tb ln

b a

Cb Tb Kb Ta

1 r

1 4

1 b r ln a K2 Tb r ln
2

1 C

K Ta Ca Tb K2 a Ta ln a

b a

C2 Tb r ln

b a

b r

KTb r

C2 Ta r ln

b r

CTa r

1 4

1 K a Ta Ka Tb C2 a Ta ln ln 1
2 2 b Ka

b a
2

K2 b Tb ln
2

b a

Cb Tb Kb Ta r

b a
2 2

b Ka

1 4

b ln a
2

1 C a r

a Ta Ka Tb C2 a Ta ln

b a b r

K2 b Tb ln O sigmah; 1 1 K
2

b a

Cb Tb Kb Ta

Tb KTa ln K 1 C Tb K ln b a

1 r

1 4

1 r ln b a

1 C K 2 Ta Ca2 Tb K2 a2 Ta ln a

b a

(13)

C2 Tb r2 ln

b a

K2 Tb r2 ln

b r

KTb r2 C2 Ta r2 ln b a
2

b r

CTa r2 b a

1 K 4

1 K a2 Ta Ka2 Tb C2 a2 Ta ln ln 1 b Ka
2 2

K2 b2 Tb ln
2

Cb2 Tb Kb2 Ta r

b a

b Ka

1 4

ln

b a
2

1 C a2 a2 Ta Ka2 Tb C2 a2 Ta ln r 1 4 1 r ln
2

b a

K2 b Tb ln

b a

Cb Tb Kb Ta

C K

b a b r b r

1 C K Ta a

Ca Tb K2 a Ta ln

b a

C2 Tb r ln

b a b a

K2 Tb r ln

KTb r C2 Ta r ln b r

b r

CTa r

1 C 4

1 C 4 Tb r ln

K4 Tb r ln r ln b a

C4 Ta r ln

1 4

1 K a Ta Ka Tb C2 a Ta ln ln b a

b a

K2 b Tb ln

b a

Cb Tb Kb Ta

b2 Ka2

1 4

1 b Ka
2 2

b ln a
2

1 C a

a Ta Ka Tb C2 a Ta ln

b a b r

K2 b Tb ln O sigmaz;

b a

Cb Tb Kb Ta

Tb KTa ln K 1 C Tb K ln b r b a

1 E K C Tb 1 K

Tb KTa ln ln b a

(14) b a b a

1 2

a2 Ta Ka2 Tb C2 a2 Ta ln

K2 b2 Tb ln b a b a

Cb2 Tb Kb2 Ta

b2 Ka2 ln

O u; 1 1 b 4 r ln a K2 Tb r2 ln 1 K 4 K 1 4
2

1 C K 2 Ta Ca2 Tb K2 a2 Ta ln a b r b r

C2 Tb r2 ln

b a

(15)

KTb r2 C2 Ta r2 ln

CTa r2 b a
2

1 K a2 Ta Ka2 Tb C2 a2 Ta ln ln 1
2 2 b Ka

K2 b Tb ln
2

b a

Cb Tb Kb Ta r

b a
2 2

b Ka
2

b ln a
2

1 C a r
2

a Ta Ka Tb C2 a Ta ln

b a

K2 b Tb ln

b a

Cb Tb Kb Ta

O a d 1; b d 1.05; E d 1e11; d .3; d 10eK5; Ta d 190; Tb d 80; a := 1 b := 1.05 E := 1. 1011 := 0.3 := 0.00010 Ta := 190 Tb := 80 O plot sigmar, r = a ..b ;

(16)

0 1.01 K # 106 1. 1.02 r 1.03 1.04 1.05

K # 106 2.

K # 106 3.

K # 106 4.

K # 106 5.

K # 106 6.

O plot sigmah, r = a ..b ;

4. # 108

2. # 108

0 1.01 1.02 r K # 108 2. 1.03 1.04 1.05

K # 108 4.

O plot sigmaz, r = a ..b ;

r 1 K 108 6.# K 108 8.# K 109 1.# K # 109 1.2 K # 109 1.4 K # 109 1.6 K # 109 1.8 K 109 2.# 1.01 1.02 1.03 1.04 1.05

O plot u, r = a ..b ;

0.0140

0.0139

0.0138

0.0137

0.0136

0.0135

1
O plot er, r = a ..b ;

1.01

1.02 r

1.03

1.04

1.05

0.020

0.018

0.016

0.014

0.012

0.010

0.008

1
O plot eh, r = a ..b ;

1.01

1.02 r

1.03

1.04

1.05

0.01349 0.01348 0.01347 0.01346 0.01345 0.01344 0.01343 0.01342

1
O

1.01

1.02 r

1.03

1.04

1.05