bab113a1-1a40-4c8c-952c-884bf59a3e70 2 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

Given:

k 15 W m K * / * :

q' 20 kW m 3^/ * :

D.1 50 mm * :

T.2 773 K * :

T.2 °C 499.85

Stefan Boltzmann constant:

Stefan Boltzmann constant: σ 5.670400 108 -^* W m 2^K 4^* / * :

A.1 π D.1 * : 0.157 m *

D.2 A.2 π D.2 * :

μ 1 :

T ε D2 eq σ T 4^T.2 4^- (* 1 ε - ε A.1 * / 1 μ A.1 * / + 1 ε - ε π D2 * * / (+ / q' π D.1 2^* 4 / (* - :

Numerical methods can't handle equals without some symbolic help. Use rhs-lhs = 0 instead.

D2 60 mm * 100 mm * 1 m * 13 mat :

initialise wisely T1I 500 °C * : T2I 1100 K * : 21 sys

t0 1 time :

ε .01 :

about insensitive ... so fast MAXITERS 10 :

j 1 D2 length range T1 T1I : T2 T2I : k 1 MAXITERS range T3 j el T1 T2 + 2 / : T3 j el ε D2 j el eq 0 ≡ break T1 ε D2 j el eq T3 j el ε D2 j el eq * 0 > T1 T3 j el : T2 T3 j el : if if 21 line for 31 line for

T3 1038.4 K * 1003.9 K * 951.6 K * 31 mat

Command-line define the plots.
Accelerates the project. Very
convenient export for "Genfit"
Rejects D2 indexed ? algo_1 ε 0.01 1 0.015 range : T3 1 el ε D2 1 el eq vectorize augment : algo_2 ε 0.01 1 0.015 range : T3 2 el ε D2 2 el eq vectorize augment : algo_3 ε 0.01 1 0.015 range : T3 3 el ε D2 3 el eq vectorize augment : 31 sys

algo_1 algo_2 algo_3 3 1 sys

This setup simulates the parametrization in the Mathcad solve blocks [Alvaro]

ε D2 fT # f x a b Bisection k 1 MAXITERS range c a b + 2 / : c f 0 ≡ break a f c f * 0 > a c : b c : if if 21 line for c 21 line : c f c ε D2 eq : c f T1 T2 Bisection 31 line :

D2 60 mm * 100 mm * 1 m * 13 mat :

end low values ε D2 1 el fT ε D2 2 el fT ε D2 3 el fT 31 mat #