Vector(-dy(psi),dx(psi)) zoom (-4,-4,8,8)as "flow shisleno" norm
contour(P) zoom (-0.004,-0.004,0.008,0.008) as "Davlenieshisleno" painted
surface(psi) zoom (-0.004,-0.004,0.008,0.008) as "stream lines"
contour(psi) zoom (-4,-4,8,8) as "stream lines" painted
elevation(P) from(0,-1) to (0,1) report(P)
history(P) versus V0 as "P'
end
title "chirkulechionnoe obtekanie chilindra"
!variables
!psi { define PSI as the system variable }
definitions
far = 15 { size of solution domain }
V0=10
M=1
a=1.5
ro=1
b=5
psi=y* V0*(1-(a/r)^2)-b*ln(r)
Vx=-dy(psi)
Vy=dx(psi)
P0=1
V=sqrt(Vx^2+Vy^2)
P=P0+(ro*(V0^2))/2-V
Boundaries
region 1 { define the domain boundary }
start(-far,-far) { start at the lower left }
line to (far,-far) { walk the boundary Counter-Clockwise }
line to (far,far)
line to (-far,far)
line to close { return to close }
start(-a,0) { start at lower left corner of airfoil }
arc(center=0,0)angle=360
monitors
contour(psi) zoom (-4,-4,8,8)as "stream lines"
plots { write hardcopy files at termination }
!grid(x,y) zoom (-4,-4,8,8)
contour(psi) zoom (-4,-4,8,8) as "stream lines shisleno" { show the flow vectors: }
Vector(-dy(psi),dx(psi)) zoom (-4,-4,8,8)as "flow shisleno" norm
contour(P) zoom (-4,-4,8,8) as "Davlenieshisleno" !painted
surface(psi) zoom (-4,-4,8,8) as "stream lines"
end
1-ый случай
2-ой случай
3-ий случай
title "Bezchirkulechionnoe obtekanie chilindra"
variables
psi { define PSI as the system variable }
definitions
far = 15 { size of solution domain }
psi_far = -y { solution at large x,y }
V0=10
P1=1
a=1.5
ro=1
psi1=-y* V0*(1-(a/r)^2)
Vx=-dy(psi1)
Vy=dx(psi1)
Vxx=-dy(psi)
Vyy=dx(psi)
Pt=P1+ro/2*(V0^2-(Vx^2+Vy^2))
P=P1+ro/2*(V0^2-(Vxx^2+Vyy^2))
equations { the equation of continuity: }
div(grad(psi)) = 0
Boundaries
region 1 { define the domain boundary }
start(-far,-far) { start at the lower left }
value(psi) = psi_far
line to (far,-far) { walk the boundary Counter-Clockwise }
line to (far,far)
line to (-far,far)
line to close { return to close }
start(-a,0) { start at lower left corner of airfoil }
value(psi)=0 { specify no flow through the airfoil surface }
arc(center=0,0)angle=360
monitors
contour(psi) zoom (-4,-4,8,8)as "stream lines"
plots { write hardcopy files at termination }
!grid(x,y) zoom (-4,-4,8,8)
contour(psi) zoom (-4,-4,8,8) as "stream lines shisleno" { show the flow vectors: }
contour(psi1) zoom (-4,-4,8,8) as "stream lines analitik "