HT HW7

School: Arizona State University, Tempe - Course: MAE 589 - Subject: Accounting

!pip install pyfluids import numpy as np import matplotlib.pyplot as plt from CoolProp import CoolProp as CP from scipy.optimize import fsolve import math Looking in indexes:https://pypi.org/simple,https://us-python.pkg.dev/colab-wheels/p Requirement already satisfied: pyfluids in /usr/local/lib/python3.8/dist-packages (2. Requirement already satisfied: CoolProp==6.4.3.post1 in /usr/local/lib/python3.8/dist Fluid='Air' L=12 W=0.2 Tinf=10 Ti=60 E=0.3 As=4*W*L D=float(4*(W**2))/(4*W) #Dh=4Ac/P x=[] y=[] z=[] sig=5.67*(10**-8) V=np.linspace(4,10,50) Te=45.1 #Assume ho=10 for i in V: Tavg=(Ti+Te)/2 C=CP.PropsSI('C','T',Tavg+273,'P',101325,Fluid) Pr=CP.PropsSI('Prandtl','T',Tavg+273,'P',101325,Fluid) rho=CP.PropsSI('D','T',Tavg+273,'P',101325,Fluid) k=CP.PropsSI('L','T',Tavg+273,'P',101325,Fluid) mu=CP.PropsSI('V','T',Tavg+273,'P',101325,Fluid) nu=float(mu)/(rho) alpha=float(k)/(rho*C) Re=float(i*D)/(nu) m=rho*i*(W**2) Nu=0.023*(Re**0.8)*(Pr**0.3) hi=float (k*Nu)/(D) def equations(p): Ts,Tex,q = p q1=-Tex+Ts-((Ts-Ti)*math.exp(float(-hi*9.6)/(m*C))) q2=-q+(ho*9.6*(Ts-10))+(E*9.6*sig*((Ts+273)**4-(283)**4)) q3=-q+(m*C*(Ti-Tex)) return (q1,q2,q3) Ts,Tx,q=fsolve(equations,(1,1,1)) Te=Tx Untitled4.ipynb - Colaboratoryhttps://colab.research.google.com/drive/1XisYxgkrPH-_Fa4IUsyJoMJ... 1 of 303/05/23, 6:51 PM

x.append(q) y.append(Te) [2631.718421437427, 2669.186428569941, 2704.9824630658795, 2739.8095285660374, 2773.7 [44.94601854102865, 45.188653986948964, 45.417657081136824, 45.63905628194276, 45.853