Question

write a program in Fortran 77 for moving magnet linear compressor operated under dynamic pressure and Tesla motor....

Answer 1

Stuff$='mercury-liquid'
T=350   [K]
k=k_(Stuff$,T)          
rho=rho_(Stuff$,T)      
c=c_(Stuff$,T)*Convert(kJ, J)  
mu=mu_(stuff$,T)      
alpha=k/(rho*c)      
nu=mu/rho          
Pr=nu/alpha
beta=beta_(stuff$, T)      

"Note that the substance can be a string (e.g., substance$) or a case insensitive string constant (e.g., 'mercury') and, for ease of use, can be copied and pasted into EES from the list of materials in the help file. The temperature can be a variable (e.g., T) or a number (e.g., 350). Unlike the EES built-in functions, the format T=T or T=350 is not used.

A complete list of substances is available in the Help provided for any of the functions."

"! Error and Warning messages. Error messages will stop EES but Warning messages are displayed after the calculations are complete even if the 'display warning messages' option is not checked in the menu Options\Preferences\Options."

"If the Solid-Liquid_Props library krhoc is not loaded and the functions, k_ (or any of the other functions: rho_, c_, mu_ or beta_) is used, the error message 'The function k_ is not defined' will be displayed. The library can be automatically loaded if the folder 'Solid-Liquid_Props' is placed in the folder Userlib.
"
"If a property does not exist, e.g., beta_('soil', 300), an error message 'A column with the name beta was not found' will be displayed."
{beta_soil=beta_('soil', 300)}

"If a material does not exist, e.g., 'mud', an error message 'A lookup table or file with the name 'Solid-Liquid_Props/mud' was not found' will be displayed."
{rho_mud=rho_('mud', 300)}

"If the temperature specified for a material is out of the allowable range, a warning message such as 'The temperature supplied to the solid density lookup table for mercury is 100 K. The allowable range is 273 to 600 K. ' will be displayed."
{rho_coldMercury=rho_('mercury', 100)}

"!Supermarket cooling system with bypass"

Select the Solve Table command in the Calculate menu (or press F3) to initiate the calculations. The Parametric table will appear showing the results. Select Plot from the Windows menu to view a plot of the results in the table."

"!Supply Air"
Vol[5]=4000 [cfm]
T[5]=62 [F]
rh[5]=0.55
P_atm=14.7 [psia]
v[5]=volume(AIRH2O,P=P_atm,T=T[5],r=rh[5])
m_dot[5]=Vol[5]/v[5]        "mass flowrate of dry air"
h[5]=enthalpy(AIRH2O,P=P_atm,T=T[5],r=rh[5])
w[5]=humrat(AIRH2O,P=P_atm,T=T[5],r=rh[5])

"!Return air"
T[6]=74 [F]
rh[6]=0.54       
h[6]=enthalpy(AIRH2O,P=P_atm,T=T[6],r=rh[6])
w[6]=humrat(AIRH2O,P=P_atm,T=T[6],r=rh[6])
m_dot[1]=0.15*m_dot[5]

"!Outside air"
T[1]=82 [F]
rh[1]=0.48;
h[1]=enthalpy(AIRH2O,P=P_atm,T=T[1],r=rh[1])
w[1]=humrat(AIRH2O,P=P_atm,T=T[1],r=rh[1])

"!Mix return and outdoor"
m_dot[7]=0.85*m_dot[5]*(1-ByPass)
m_dot[2]=m_dot[1]+m_dot[7]
w[2]*m_dot[2]=m_dot[1]*w[1]+m_dot[7]*w[6]   {water balance}
m_dot[1]*h[1]+m_dot[7]*h[6]=m_dot[2]*h[2]   {energy balance}

"!Cooling Coil"
rh[3]=1
m_dot[3]=m_dot[2]
Q_C=m_dot[3]*(h[3]-h[2])
h[3]=enthalpy(AIRH2O,P=P_atm,w=w[3],r=rh[3])
T[3]=temperature(AIRH2O,P=P_atm,w=w[3],r=rh[3])

"!Mix coil outlet with bypass"
m_dot[8]=0.85*m_dot[5]*(ByPass)
m_dot[4]=m_dot[3]+m_dot[8]
w[6]*m_dot[8]+w[3]*m_dot[3]=m_dot[4]*w[4]
m_dot[8]*h[6]+m_dot[3]*h[3]=m_dot[4]*h[4]

"!Reheat coil"
w[5]=w[4]
Q_H=(h[5]-h[4])*m_dot[5]

"!Get missing Temps for plotting states on the Psych chart"
T[2]=Temperature(AIRH2O,P=P_atm,w=w[2],h=h[2])
T[4]=Temperature(AIRH2O,P=P_atm,w=w[4],h=h[4])

End

Stop

End