Wolfram mathematica Mathematica中的流线
有人能提出一个更好的方法来解决Mathematica中的流线问题吗? P.S.StreamPlot未提供准确的结果 这里,速度是Wolfram mathematica Mathematica中的流线,wolfram-mathematica,mathematica-8,fluid-dynamics,streamline,Wolfram Mathematica,Mathematica 8,Fluid Dynamics,Streamline,有人能提出一个更好的方法来解决Mathematica中的流线问题吗? P.S.StreamPlot未提供准确的结果 这里,速度是u[x,y]i+v0[x,y]jp0x[x]是压力梯度。代码中的最后一行太长,无法执行更简单的解决方案 代码如下: B1[x_] = 1 + a*Cos[2*Pi*x]; p0x[x_] = -1/B1[x]; u0[x_, y_] = -p0x[x]; v0[x_, y_] = D[p0x[x], x]*y; u1[x_, y_] = \[Lambda]^2*
u[x,y]i+v0[x,y]j
p0x[x]
是压力梯度。代码中的最后一行太长,无法执行更简单的解决方案
代码如下:
B1[x_] = 1 + a*Cos[2*Pi*x];
p0x[x_] = -1/B1[x];
u0[x_, y_] = -p0x[x];
v0[x_, y_] = D[p0x[x], x]*y;
u1[x_, y_] = \[Lambda]^2*D[p0x[x], {x, 2}]/2*y^2 - \[Lambda]^2* D[p0x[x], {x, 2}]/6*(B1[x])^2;
u[x_, y_] = -p0x[x] + \[Delta]^2*\[Lambda]^2*D[p0x[x], {x, 2}]/2*(y^2 - (B1[x])^2/3);
eq = (D[y[x], x] == v0[x, y[x]]/u[x, y[x]]) /. {a -> 0.3, \[Lambda] -> 1, \[Delta] -> 1} ;
sl = DSolve[eq, y[x], x]
可能会在mathematica网站上得到更好的回复哦。谢谢这绝对是一个更好的主意。可能会在mathematica网站上得到更好的回复哦。谢谢这绝对是个更好的主意。