Matlab practice

The flow rate ranges from $20$ at $200\frac{L}{m}in$ for the pipe system below containing $50mm$ ID

PVC pipes expands up to $90mm,$ then reduces down to a $35mm$ ID pipe. The figure below is

a side view. The change in height from the entrance to the exit is marked on the diagram. The

pressure at $P_t$ is atmospheric pressure and you may use the large tank approximation to

estimate $v_t.$ You are asked to calculate and then plot the pressure $($m$)$ at the exit as a function

of flow rate $(\frac{m^3}{s}).$ Diagram not to scale

$lo{n}_{PVC}=2\mathrm{.}5*{10}^{-6}m,g=9\mathrm{.}81\frac{m}{s^2},P_{\text{a}\text{t}\text{n}}=101,325Pa$

${}_{\text{w}\text{a}\text{t}\text{e}\text{r}}=998k\frac{g}{m^3},{}_{\text{w}\text{a}\text{t}\text{e}\text{r}}=1\mathrm{.}0*{10}^{-3}Pa*s$ Calculate the following:

the flow rate at the entrance as a vector from $20$ to $200\frac{L}{m}in$ in increments of $20$

L$/$min$,$ Q$\_$L$\_$min

convert the above flow rate to SI units, $\frac{m^3}{s},Q$

the fluid velocity as vectors in $\frac{m}{s}$ for the first, second and third section of pipe, v$1,$ v$2$

and $v3$

Reynolds number for the first, second and third section of pipe as vectors, $Re1,Re2$

and $Re3$

friction factor for the first, second and third section of pipe as vectors, $f1,f2$ and $f3$

the total pressure drop for the pipes in $m$ for the first, second and third section of pipe

as vectors, h$\_$pipe$1,$ h$\_$pipe$2$ and h$\_$pipe$3$

the total pressure drop for the pipes as a vector in $m,$ h$\_$pipe

the K$-$value of reducer as a vector, K$\_$red

the K$-$value of expander as a vector, K$\_$exp

the total pressure drop for the fittings as a vector in $m,h_{-}$fittings

the absolute pressure at the outlet in $m,P_{{?}_{}}$ out $($note: in METERS$)$

plot the absolute pressure drop, P$\_$out, in as $m$ a function of flow rate, $Q,$ in $\frac{m^3}{s}.$On the plot,

label the x$-$axis Flow Rate

label the y$-$axis Pressure $($m$)$

add a title Graph of pressure versus flow rate

You need to use if$/$else statement to check the flow regime and you must us a for loops

to calcualte the vectors for f$1,$ f$2,$ f$3,$ K$\_$red and K$\_$exp with the if$/$else statements

Caution when performing calculations:

All calculations to be performed in Matlab. Manually calculated values that are entered

into Matlab are not accepted.

DO NOT use any 'clear' or $\text{'}$clc$\text{'}$ commands in your script.

For constants, use applied Matlab functions e$.$g$.$ pi instead of $3\mathrm{.}142.$

You must use the variable names specified in the question.

Parameters:

X$=12$

Y$=34$

Z$=3$

Laminar flow, f$=64/$Re

Turbulent flow, f$=1/(-1\mathrm{.}8$log$10(6\mathrm{.}9/$Re$+(($epsilon$/$D$)/3\mathrm{.}7)^1\mathrm{.}11)^2$