Моделирование течения бесциркуляционного обтекания цилиндра поступательным потоком в среде FLEX.PDE

title " Bezchirkulechionnoe obtekanie chilindra"

Variables

psi { define PSI as the system variable }

Definitions

far =15 { size of solution domain }

{ solution at large x, y }

V1=10 P1=1 a=1.5 ro=1

psi_far = -y* V1

psi1=-y* V1*(1-(a/r)^2)

Vx=-dy(psi1) Vy=dx(psi1) Vxx=-dy(psi) Vyy=dx(psi)

Pt=P1+ro/2*(V1^2-(Vx^2+Vy^2))

P=P1+ro/2*(V1^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

plots { write hardcopy files at termination }

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 "

vector(-dy(psi), dx(psi)) zoom (-4, -4, 8, 8)as " flow shisleno" norm

vector(-dy(psi1), dx(psi1)) zoom (-4, -4, 8, 8)as " flow analitik" norm

contour(Pt) zoom (-4, -4, 8, 8) as " Davlenie analitik" painted

contour(P) zoom (-4, -4, 8, 8) as " Davlenieshisleno" painted

surface(psi) zoom (-4, -4, 8, 8) as " stream lines"

End