Моделирование течения вихревая точка в среде FLEX.PDE

title " Vichrevaya tochka"

Variables

psi { define PSI as the system variable }

Definitions

far =10 { size of solution domain }

psi_far = 5*ln(r) { solution at large x, y }

V1=1 P1=1 a=0.5 ro=1

psi1=50*ln(r)

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

div(grad(psi)) = 0

Boundaries

region 1 { define the domain boundary }

start(-far, -far) { start at the lower left }

value(psi) =0

line to (far, -far) { walk the boundary Counter-Clockwise }

Line to (far, far)

Line to (-far, far)

line to close { return to close }

start(0.001, 0) { start at lower left corner of airfoil }

value(psi)= 50*ln(r) { 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 }

contour(psi) zoom (-0.001, -0.001, 0.002, 0.002) as " stream lines shisleno" { show the flow vectors: }

contour(psi1) zoom (-0.001, -0.001, 0.002, 0.002) as " stream lines analitik "

vector(-dy(psi), dx(psi)) zoom (-0.001, -0.001, 0.002, 0.002)as " pole skorostey chisleno" norm

vector(-dy(psi1), dx(psi1)) zoom(-0.001, -0.001, 0.002, 0.002)as " pole skorostey analitik" norm

contour(Pt) zoom(-0.001, -0.001, 0.002, 0.002) as " Davlenie analitik" painted

contour(P) zoom(-0.001, -0.001, 0.002, 0.002) as " Davlenie shisleno" painted

End