Variables. psi { define PSI as the system variable }

psi { define PSI as the system variable }

Definitions

far = 5 { size of solution domain }

! psi_far =- 5*arctan(y/x)

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

psi1=-5*arctan(y/x)

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 }

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.0001, 0) { start at lower left corner of airfoil }

value(psi)= -5*arctan((-abs(y))/(-abs(x))) { 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 (-4, -4, 8, 8) as " stream lines analitik "

vector(-dy(psi), dx(psi)) zoom (-4, -4, 8, 8)as " pole skorostey chisleno" norm

vector(-dy(psi1), dx(psi1)) zoom (-4, -4, 8, 8)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