Given the following parameters, use the Colebrook equation, the formula for Reynold's number, and v $=$ Q$/$A as well as the following formula for NPSH$\_$avail with the end goal of creating a plot in Python that graphs NPSH$\_$avail $($in m$)$ vs$.$ volumetric flow rate Q $($in Liters per minute$).$ NPSH$\_$avail $=7\mathrm{.}836-(5\mathrm{.}946*$ f $+0\mathrm{.}059)*$v$**2($all in SI units$)$

$**$Note: I have already tried to figure this out but my code keeps not working. I have attached my code below:$**$

import numpy as np

import scipy.optimize as opt

from scipy.optimize import fsolve

import matplotlib.pyplot as plt

from math import pi$,$ sqrt$,$ log$10$

# Parameters

rho $=997\mathrm{.}0$ #kg$/$m$^3$

mu $=8\mathrm{.}91$e$-4$ #kg$/$m$*$s

Pv $=3169$ #Pa

z $=2\mathrm{.}2$ #m

KL$\_$tot $=1\mathrm{.}15$

L $=2\mathrm{.}8$ #m

D $=0\mathrm{.}024$ #m

epsilon $=0$

P$\_$atm $=101325$ #Pa

g $=9\mathrm{.}81$ #m$/$s$^2$

df # Colebrook equation

def colebrook$\_$eq$($f$,$ Re$,$ epsilon, D$)$:

try:

return $1/$ sqrt$($f$)+2*$ log$10(2\mathrm{.}51/($Re $*$ sqrt$($f$)))$

except $($ZeroDivisionError$,$ ValueError$)$:

return float$(\text{'}$inf$\text{'})$

# NPSH Available

def calculate$\_$NPSH$\_$avail$($f$,$ v$)$:

return $7\mathrm{.}836-(5\mathrm{.}946*$ f $+0\mathrm{.}059)*$ v$**2$

# Generate flow rates

Q$\_$m$3$s $=$ np$.$linspace$(0\mathrm{.}1,200,100)$

# Calculate NPSH available for each flow rate

NPSH$\_$avail$\_$values $=[]$

for Q in Q$\_$m$3$s:

v $=$ Q $/($pi $*($D $/2)**2)$

Re $=($rho $*$ v $*$ D$)/$ mu

initial$\_$guess $=$ max$(0\mathrm{.}02,0\mathrm{.}05*$ Re $*$ D $/$ Q$)$

f $=$ fsolve$($colebrook$\_$eq$,0\mathrm{.}02,$ args$=($Re$,$ epsilon, D$))[0]$

NPSH$\_$avail $=$ calculate$\_$NPSH$\_$avail$($f$,$ v$)$

NPSH$\_$avail$\_$values.append$($NPSH$\_$avail$)$

# Convert to Lpm

Q$\_$Lpm $=$ Q$\_$m$3$s $*60*1000$

# Print the first few values of Q$\_$Lpm and NPSH$\_$avail$\_$values

print$("$Flow Rates $($L$/$min$)$:$",$ Q$\_$Lpm$[$:$5])$

print$("$NPSH Available $($m$)$:$",$ NPSH$\_$avail$\_$values$[$:$5])$

# Check for NaN or infinite values in NPSH$\_$avail$\_$values

print$("$NaN or Infinite values:", np$.$isnan$($NPSH$\_$avail$\_$values$).$any$()$ or np$.$isinf$($NPSH$\_$avail$\_$values$).$any$())$

# Plot

plt$.$plot$($Q$\_$Lpm$,$ NPSH$\_$avail$\_$values$)$

plt$.$xlabel$(\text{'}$Flow Rate $($L$/$min$)\text{'})$

plt$.$ylabel$(\text{'}$NPSH Available $($m$)\text{'})$

plt$.$title$(\text{'}$NPSH Available vs Flow Rate'$)$

plt$.$grid$($True$)$

plt$.$xlim$(0,200)$

plt$.$ylim$(0,20)$

plt$.$show$()$