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

title " Chrkulychonnoe obtekanie chlindra"

Variables

psi { define PSI as the system variable }

Definitions

far = 5 { size of solution domain }

{ solution at large x, y }

V0=15 P0=1 a=0.5 ro=1 D=-15

GA=-4*pi*a*V0+D

psi1=V0*y*(1-a^2/r)+GA/(2*pi)*LN(r)

Pt=P0+ro/2*(V0^2-(2*V0*(y/sqrt(x^2+y^2)+GA/(2*pi*a)))^2)

Vx=-dy(psi1) Vy=dy(psi1)

psi_far =V0*y

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 }

to (far, far) line to (-far, far) line to close { return to close }

start " kontur" (a, 0)

value(psi)=1.*(GA/(2*pi))*ln(a)

arc(center=0, 0)angle=360 plots { write hardcopy files at termination }

contour(psi) zoom(-1.5, -1.5, 3, 3)as " linii toka shisleno" { show the flow vectors }

contour(psi1) zoom (-1.5, -1.5, 3, 3) as " linii toka analitishski"

vector(-dy(psi), dx(psi)) zoom (-1.5, -1.5, 3, 3) as " pole skorosty shislo " norm

vector(-dy(psi1), dx(psi1)) zoom (-1.5, -1.5, 3, 3)as" pole skorosty analitik " norm

contour(Pt) zoom (-2, -2, 4, 4) as " Pole davleniy analitishski" norm! painted

End

Движение вязкой жидкости