38f0dda8-9a1b-435f-83b5-d023fd863de6 49 4 5 decimal

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

Colebrook memo m ε 0.00003<ε>0.0001 ≡ m D Pipe ID(m) ≡ f friction factor [dimensionless] ≡ R Reynolds [dimensionless] ≡ ε D 12 mat in the same units system ≡ f 0.02 : initial guess all cases ≡ 61 line

Moody Diagram from Colebrook

Colebrook formula 1 f sqrt / 2.0 ε 3.71 D * / 2.51 R f sqrt * / + log10 * - ≡

t0 1 time :

time 1.2 min * ≡

f R D ε Find 1 f sqrt / 2.0 ε 3.71 D * / 2.51 R f sqrt * / + log10 * - ≡ f 01 solve :

Moody ε 0.00005 : R 10000100000020000 range : D 0.04 0.05 0.06 0.07 14 mat transpose : i 1 R rows range j 1 D rows range f i j el f R i el D j el ε Find : for for R f 1 col augment R f 2 col augment R f 3 col augment R f 4 col augment 41 mat 61 line eval :

Read this: X_Y plot is for publishing, project collection. It is terribly slow compared to Smath QuickPlot. 1. Smath => 5 min 2. X_Y => 27 min. The other interpretation is that X_Y plot is NOT directly compatible to Smath construction. Alternately: spool to file => 1.2 min.

No need for "spool to file", eval() the Moody function instead.

CWD C:\SmathFile CurrentDirectory :

Moody CWD Moody wfile 1

Moody CWD Moody rfile :

Moody 1 el Moody 2 el Moody 3 el Moody 4 el 4 1 sys

1 f sqrt / 2.0 ε 3.71 D * / 2.51 R f sqrt * / + log10 * - ≡

1 time t0 - 60 sec * / 0.3

Next associated Engineering piping stuff: "Inst_Colebrook Pipe Data"

yes moody_plot.xlsx Sheet1 A1 D1 [m] D2 [m] D3 [m] D4 [m] 28 mat Moody 1 el Moody 2 el Moody 3 el Moody 4 el augment stack eval excel_IN 1

Y29weSBzdWNjZXNmdWxs moody_plot.xlsx Sheet1 K6:R22 1 excel_EMF