LINE TO(fa,-fa) TO (fa,fa) TO (-fa,fa) TO CLOSE. start(a,0)value(psi1)=b*ln(r)

start(a,0)value(psi1)=b*ln(r)

arc(center=0,0)angle=360

MONITORS { show progress }

PLOTS { save result displays }

contour(psi) as "linii toka analit "

vector(Vx,Vy)as "pole skorostei analit" norm

contour(P) as "pole davleniya analit"painted

contour(psi1) as "linii toka "

vector(Vxx,Vyy)as "pole skorostei" norm

contour(P1) as "pole davleniya"painted

! CONTOUR(u)

END

Поток

Пусть в плоскости Z задано течение, определяемое комплексным потенциалом

W= ɑz = ɑx+iɑy, (4.3)

тогда потенциал скорости φ =ɑx, функция тока Ψ=ɑy. Соответственно, эквипонтенциали ɑx=const, прямые параллельные оси oy, а линии тока ɑy= const, прямые параллельные оси ox (рис.4.1), при этом

. (4.4)

Следовательно, комплексный потенциал, задаваемый в виде линейной функции, порождает течение однородного поступательного потока.

{ Fill in the following sections (removing comment marks! if necessary),

and delete those that are unused.}

TITLE" postupatelnoe techenie" { the problem identification }

COORDINATES cartesian2 { coordinate system, 1D,2D,3D, etc }

VARIABLES { system variables }

psi { choose your own names }

DEFINITIONS { parameter definitions }

fa=7

V0=50

ro=1

P0=1

Vx=-dy(psi)

Vy=dx(psi)

V=sqrt(Vx^2+Vy^2)

P=(V0^2/2-V^2/2+P0/ro)*ro

EQUATIONS { PDE's, one for each variable }

div(grad(psi))=0 { one possibility }

BOUNDARIES { The domain definition }

REGION 1 { For each material region }

START(-fa,-fa) VALUE(psi)=V0*y { Walk the domain boundary }

Не нашли, что искали? Воспользуйтесь поиском: