Доброго времени суток, товарищи!Записал уравнения Навье-Стокса в сферической системе координат, забил в мэйпл.
Он даже выдаёт мне аналитические решения, однако я не могу ими воспользоваться - команда assign не "выдирает решения" и после неё выдаётся сообщение об ошибке. Вектор решений содержит только две компоненты скорости, и на этом дело глохнет.
Кроме того, не могу выражения для граничных и начальных условий запихнуть в команду pdsolve
Если кто может - помогите.
Вот код, достаточно простой.![$restart;with(PDEtools);with(DEtools);
SetCoordinates('spherical'[r,theta,phi]);declare((Vr,Vtheta,Vphi)(r,theta,phi));
ro:=1040;etha:=0.001;k:=etha/ro;R:=5e-2;omeg:=100;Vtan:=3;#Параметры
#Уравнения, расписанные покомпонентно
eq1:=diff(Vr(r,theta,t),t)+(Vr(r,theta,t)*diff(Vr(r,theta,t),r)+Vtheta(r,theta,t)*diff(Vr(r,theta,t),theta)/R)-k*(diff(Vr(r,theta,t),r,r)+cot(theta)*diff(Vr(r,theta,t),theta)/R^2+diff(Vr(r,theta,t),theta,theta)/R^2)=0;
eq2:=diff(Vtheta(r,theta,t),t)+(Vr(r,theta,t)*diff(Vtheta(r,theta,t),r)+Vtheta(r,theta,t)*diff(Vtheta(r,theta,t),theta)/R)-k*(diff(Vtheta(r,theta,t),r,r)+cot(theta)*diff(Vtheta(r,theta,t),theta)/R^2+diff(Vtheta(r,theta,t),theta,theta)/R^2)=0;
eq3:=diff(Vphi(r,theta,t),t)+(Vr(r,theta,t)*diff(Vphi(r,theta,t),r)+Vtheta(r,theta,t)*diff(Vphi(r,theta,t),theta)/R)-k*(diff(Vphi(r,theta,t),r,r)+cot(theta)*diff(Vphi(r,theta,t),theta)/R^2+diff(Vphi(r,theta,t),theta,theta)/R^2)=0;
sys1:=[eq1,eq2,eq3]; # Система
#Граничные и начальные условия
ICS:={Vr(R,Pi,t)=0,Vtheta(R,Pi,t)=Vtan,Vphi(R,theta,t)=omeg*R*sin(theta),Vr(R,theta,0)=0,Vtheta(R,theta,0)=0,Vphi(R,theta,0)=0};
sol:=pdsolve(sys1);$ $restart;with(PDEtools);with(DEtools);
SetCoordinates('spherical'[r,theta,phi]);declare((Vr,Vtheta,Vphi)(r,theta,phi));
ro:=1040;etha:=0.001;k:=etha/ro;R:=5e-2;omeg:=100;Vtan:=3;#Параметры
#Уравнения, расписанные покомпонентно
eq1:=diff(Vr(r,theta,t),t)+(Vr(r,theta,t)*diff(Vr(r,theta,t),r)+Vtheta(r,theta,t)*diff(Vr(r,theta,t),theta)/R)-k*(diff(Vr(r,theta,t),r,r)+cot(theta)*diff(Vr(r,theta,t),theta)/R^2+diff(Vr(r,theta,t),theta,theta)/R^2)=0;
eq2:=diff(Vtheta(r,theta,t),t)+(Vr(r,theta,t)*diff(Vtheta(r,theta,t),r)+Vtheta(r,theta,t)*diff(Vtheta(r,theta,t),theta)/R)-k*(diff(Vtheta(r,theta,t),r,r)+cot(theta)*diff(Vtheta(r,theta,t),theta)/R^2+diff(Vtheta(r,theta,t),theta,theta)/R^2)=0;
eq3:=diff(Vphi(r,theta,t),t)+(Vr(r,theta,t)*diff(Vphi(r,theta,t),r)+Vtheta(r,theta,t)*diff(Vphi(r,theta,t),theta)/R)-k*(diff(Vphi(r,theta,t),r,r)+cot(theta)*diff(Vphi(r,theta,t),theta)/R^2+diff(Vphi(r,theta,t),theta,theta)/R^2)=0;
sys1:=[eq1,eq2,eq3]; # Система
#Граничные и начальные условия
ICS:={Vr(R,Pi,t)=0,Vtheta(R,Pi,t)=Vtan,Vphi(R,theta,t)=omeg*R*sin(theta),Vr(R,theta,0)=0,Vtheta(R,theta,0)=0,Vphi(R,theta,0)=0};
sol:=pdsolve(sys1);$](https://dxdy-02.korotkov.co.uk/f/d/8/4/d84f2c567ac1c9fc441c76ba8d64fe9282.png)