简化Mathematica

问题描述

有人可以提出一种更好的方法解决Mathematica中的流线型问题吗? 附言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*D[p0x[x],{x,2}]/2*y^2 - \[Lambda]^2* D[p0x[x],2}]/6*(B1[x])^2;

u[x_,y_] = -p0x[x] + \[Delta]^2*\[Lambda]^2*D[p0x[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]

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)