Question

write a FORTRAN 77 software programme for moving coil type linear compressor using R407c refrigerant in pulse tube refrigerator

Answer 1

"Test of the genetic optimization method

This function has two peaks. A traditional optimization algorithm would likely be attracted to the broader peak, although it represents a local optimum. The genetic algorithm correctly identifies the global optimum.

Press F4 to start the optimization"

D=2
f=A_1*exp(-r_1^2/sigma_1^2)+A_2*exp(-r_2^2/sigma_2^2)
r_1^2=sum((x[j]-0.5)^2,j=1,D)
r_2^2=sum((x[j]-0.2)^2,j=1,D)
A_1=0.7
A_2=1-0.7*exp(-r_1^2/sigma_1^2)
sigma_1^2=0.15
sigma_2^2=0.005

"!This program demonstrates the use of the Integral functions to solve second order equations. "

"Here EES is used to calculate the velocity and position of a freely falling sphere, subject to aerodynamic drag. The unit system is set to English. The graph is set to automatic update - change v_o to -50 to see the impact of an initial upward velocity.

Note how the Integral function displays on the Formatted Equations Window."

D=0.25 [ft]
m=1.0 [lb_m]    "mass of sphere"
v_o=0 [ft/s]    "initial velocity."
z_o=0 [ft]     "initial position"
time=5 [s]   "time period for analysis"
g=32.17 [ft/s^2]    "gravitational acceleration"

F=m*g*Convert(lbm-ft/s^2,lbf)         "Newton's Law"
m*a*Convert(lbm-ft/s^2,lbf)=F-F_d        "force balance"
Area=pi*D^2/4           "frontal area of sphere"
F_d=Area*C_d*(1/2*rho*v^2)*Convert(lbm-ft/s^2,lbf)    "definition of drag coefficient"

"Find Reynolds number"
mu=viscosity(air, T=70)*Convert(1/hr,1/s)
rho=density(Air,T=70,P=14.7)
Re=rho*abs(v)*D/mu

"Find drag coefficient from the Reynolds number. The Lookup table contains ln(Re) and ln(C_d). The max function is used to prevent attempting to find the log of zero (i.e., when the velocity is zero use a small value of Re)"
C_d=exp(interpolate1( 'LnRe', 'LnCd', LnRe=Ln(max(.01, Re))))

"As a test of the need for tthe variable drag coefficient, set C_d to a constant value, say C_d=0.4. Turn off automatic update on the plots (click on the plot window) and overlay the new plots, using the left scale."
{C_d=0.4}

"Use EES integral function to determine velocity and position given the acceleration."
v=v_o+integral(a,t,0,time)   "velocity after 5 seconds"
z=z_o+integral(v,t,0,time)   "vertical position after 5 seconds"

"The following directive instructs EES to store values of v (velocity), z (elevation) and C_d (drag coefficient) as a function of t (time) at increments of 0.2 sec.
"
$integraltable t:0.2, v,z, C_d
$tabstops 1 in