a0803e6c-5acc-4a89-86bc-f6e5d8493e06 23 4 5 false false 0 7 decimal

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

m 0.5886 :

t.R 2824 * 60 * 60 * : 2.4192 106^*

ϕ.k 0.86 :

γ.c 1.0 :

γ.p 1.11 :

t.c 30365 * 24 * 60 * 60 * : 9.4608 108^*

C.lim 0.004 426 * : 1.704

β.cl 1.5 :

Cs 6.0 :

D.R 7.516 1012 -^* :

c 85 :

Ci 0.0 :

t D.p t 30 ≥ D.R 1 m - (/ 1 m - (m t.c t 365 * 24 * 60 * 60 * / * + ((* t.R t.c / (m^* D.R 1 m - (/ t.R t 365 * 24 * 60 * 60 * (/ (m^* if 11 line :

t D.k t D.p γ.p ϕ.k / * :

t D.d t D.k γ.c * β.cl * :

t C.d Cs Ci - (1 c 0.001 * (2 t D.d t * 365 * 24 * 60 * 60 * sqrt * / erf - (* Ci + :

t γC.d γ.p t C.d * :

kg/m^3

x f γ.p x C.d * :

t ϕC.lim ϕ.k C.lim * :

x f ϕ.k C.lim * - x 01000 solve #

Optimization Symbolic x f Unknowns #

x f ϕ.k C.lim * - x 1 roots #

x f x ϕC.lim - 100 FindRoot #

x f x ϕC.lim - 100200105 -^BDQRF #

x f x ϕC.lim - 100200105 -^Bisection #

x f x ϕC.lim ≡ 100105 -^NewtonRaphson #

x eq x f x ϕC.lim - :

v g v 1 el eq :

5011 mat 0105 -^g al_nleqsolve 132.714 11 mat

Custom bisection algortihm for f function a# b# goal ε# fsolver Δ 1 : x.a a# : x.b b# : Δ ε# > β x.a x.b + 2 / : S β f eval : goal S - 0 > x.a β : goal S - 0 < x.b β : 0 if Δ x.b f x.a f - abs : 41 line while β 51 line :

x ϕC.lim 1.4654

xf 100200 x ϕC.lim 105 -^fsolver eval :

xf 132.7141

x γC.d xf 0 xf x ϕC.lim 0 x ϕC.lim 3 2 mat 2 1 sys