wolfram-mathematica - Mathematica 中的流线型

Question

有人可以建议一种更好的方法来解决 Mathematica 中的流线问题吗？PS StreamPlot 没有提供准确的结果。

在这里，速度是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]