f1cfddea-23de-470d-ada5-f6e0e45b2b57 30 4 5 decimal

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

BOMBAS CENTRÍFUGAS

Viscocidad Dinámica del agua

Densidad del Agua

Diámetro de la Tubería

μ 1103 -^* Pa * s * :

ρ 1000 kg m 3^/ * :

D 2.5 in * :

Rugosidad de la Tubería

Longitud de Tubería

Constante Gravitacional

g.e 9.8066 m s 2^/ *

ε 0.005 :

L 1500 ft * :

V Q π D 2^* 4 / V * :

V H 20 ft * 5103 -^* ft gal min 1 -^* (/ * V Q * - :

En mathcad los valores semillas llevan el mismo nombre que la variable, pero en SMath tienen que estar "limpios".

V Rn D V * ρ * μ / :

f.tur.guess 0.1 :

V f.tur 1 f.tur sqrt / 2 ε 3.7 / 2.51 V Rn f.tur sqrt * / + 10 log * - ≡ 0.1 FindRoot abs :

Parece que FindRoot devuelve una pequeña parte imaginaria, así que tomo el valor rea (y renombro el numero de Reynolds)

1000000 f.tur 0.0304

V f.lam 64 V Re / :

V f.f V Rn 3000 ≥ V f.tur V f.lam if 11 line :

V Rn 2300 >

1 m s / * f.f 0.0318

Funciona

Ahora probablemente FindRoot o cualquier otro solver del plugin NonlinearSolver no funcuione, pues no puede encontrar las incóngnitas, al no funcionar Unknows( ) en la nueva función, pues usa también un solver. No encuentro ninguno que funciones, ni en maple (que tiene un bug y no puedo usar fsolve) ni en maxima ni en ningún otro.

La raíz ahí está, pero no sé cómo calcularla. Quizás con bisección, pero programada especial para este caso.

V eqV V f.f L D / * V 2^2 / * g.e V H * - :

t0 1 time :

9 min

x m s / * eqV

1 time t0 - min 8.4

The gray (symbol Gy) is the SI unit of absorbed dose and is defined as the absorption of one joule of energy, in the form of ionizing radiation, per kilogram of matter, i.e. one gray = 1 J/kg 2.Feb. 24, 2018

1 m s / * eqV 55.4965 Gy *

3 m s / * eqV 942.6666 Gy *

Deshabilitada la gráfica porque es muy lenta.

u 1 m s / * (3 m s / * (1.25 m s / * range :

t0 1 time :

U u m s / / :

V u eqV vectorize m 2^s 2^/ / eval :

1 time t0 - 53 s *

UV U V augment :

x f U V x cinterp :

X x f x 0.25 1 solve : 0.7133

x U 5120 x * + :

UV 1 55.4965 1.25 118.5615 1.5 195.289 1.75 285.686 2 389.7465 2.25 507.4816 2.5 638.8775 2.75 783.9392 3 942.6666 92 mat

x λ 10 x ≤ cases :

Gy x f x U - 0 ≡ x 13 solve : 2.2591

Gy f 512

x f x λ * 1 - x U Gy x U Gy 100 - 3 2 mat X 0 + 40 red 1 5 mat 3 1 sys