Variables

temp { bezrazmerniy kvadrat davleniy v zazore }

Select

Definitions

Par=0.98e5 {davlenie na kromke podshipnika, n/m2 }

Psr=5*Par {davlenie nadduva, n/m2 }

Rnr=0.112 {narushny radius podshipnika, m}

Rkr= 0.03 {vnutrenniy radius podshipnika, m}

Rdr=0.02 {radius okrushnosti raspoloshniy pitateley, m}

dpr=0.8e-3! 0.8e-3 { diametr pitateley, m}

H=30e-6! staged(20e-6, 25e-6, 30e-6, 35e-6, 40e-6, 45e-6, 50e-6, 60e-6) { velichina zazora podchipnika, m }

k1=1.4 { pokazatel adiabaty gaza }

mu=179.2e-7! 179.2e-7{nc/m2}

ass=341.7 {skorost zvuka, m/c }

AL=0.7

N=3 {chislo pitateley}

A= pi*dpr*dpr/4 {plochad pitately }

gamma=((12*AL*ass*mu*N*A*Rnr/Rkr)/(3.14*(H^3)*Psr))*(2/(k1+1))^((k1+1)/(2*(k1-1))) {konstruktivniy parametr }

Pa=Par/Psr

Rd=Rdr/Rnr

Rn=1

Rk=Rkr/ Rnr

dp=dpr/(Rnr*2)

p=(abs(temp))^(1/2)

omega=10000! staged(100, 500, 1000, 3000, 4000, 5000, 5500, 6000, 6500, 8000, 10000, 15000, 20000, 35000, 50000)

sigma=(12*mu*omega*Rnr*Rnr)/(H*H*Psr)

nOM=omega*30/pi

Qfi=-dx(p)*y+dy(p)*x

F1=(2/(k1+1))^(2*k1/(k1-1))

F2=sqrt(2/(k1-1)*((k1+1)/2)^((k1+1)/(k1-1)))

Pkn2=val(temp, 0, Rk)

Pk=sqrt(Pkn2)

Initial values

temp =0.999

equations { define the heatflow equation }

div(grad(temp)) =sigma* Qfi

Boundaries

Region 1 { the outer boundary defines the copper region }

Start (0, -Rn)

value (temp) =Pa ^2 {convection boundary}

arc(center=0, 0) angle 360

Region 2

Start (0, -Rk)

natural(temp)=gamma*(((abs(temp))^(1/k1)-(abs(temp))^((k1+1)/(2*k1)))^(1/2))*F2

arc(center=0, 0) angle 360

repeat i=0 to N-1{ an indexed loop on X-position } { an array of circular dots at the tabulated coordinates: }

start (Rd*sin(2*pi/N*i)+ dp, Rd*cos(2*pi/N*i))

natural(temp)=gamma*(((abs(temp))^(1/k1)-(abs(temp))^((k1+1)/(2*k1)))^(1/2))*F2

arc(center=Rd*sin(2*pi/N*i), Rd*cos(2*pi/N*i)) angle=360 endrepeat

Plots grid(x, y)

contour(p) zoom(-1, -1, 2, 2)as " pole davlenii" painted

report (Rnr^2*Psr*INTEGRAL(p-Pa))report(p)

surface(p)zoom(-1, -1, 2, 2) as " pole davlenii" painted report(sigma)

Elevation(p) from (Rk, 0) to (Rn, 0) report(sigma)

Summary

report (Psr*(Rnr^2)*(INTEGRAL(p-Pa)))as " W nesuchay sposobnost"

report (Psr*(Rnr^2)*(INTEGRAL(p-Pa))/(pi*(Psr-Par)*(Rnr^2))) as " Cw koefichent nesuchay sposobnost"

{histories

! history (Psr*(Rnr^2)*(INTEGRAL(p-Pa))) versus nOM as " W nesuchay sposobnost"

! history (Psr*(Rnr^2)*(INTEGRAL(p-Pa))) versus omega as " W nesuchay sposobnost"

! history (Psr*(Rnr^2)*(INTEGRAL(p-Pa))/(pi*(Psr-Par)*(Rnr^2))) versus omega as " Cw koefichent nesuchay sposobnost"

history (Psr*(Rnr^2)*(INTEGRAL(p-Pa))/(pi*(Psr-Par)*(Rnr^2)))versus H as " Cw koefichent nesuchay sposobnost"

history (Psr*(Rnr^2)*(INTEGRAL(p-Pa))) versus H as " W nesuchay sposobnost"

history(-((pi*k1*(H^3)*(Psr^2))/(12*(ass^2)*mu))*(Pk^2-Pa^2)/ln(Rkr/Rnr)) versus H as " Q rasxod smazkit" }

End