请问这样一个半无限域扩散方程用Dsolve怎么解
mathematica吧
全部回复
仅看楼主
level 4
这是我程序中给出的方程和边界:
eqn=(u^(1,0))[t,y]==\[Nu] (u^(0,2))[t,y];
bc1=u[0, y]==0;
bc2=u[t, 0]==Piecewise[{{U,t>0}},0];
bc3=u[t, y]==0/.y->Infinity;
有个类似的帖子我看了,不过他的是第二类边界,而且没有涉及到无穷处的取值。
我这个问题是不是用Fourier变换更简单一点?用laplace变换怎么做,我边界条件还是处理不好。
2018年11月10日 07点11分 1
吧务
level 15
……你没看懂这帖:https://tieba.baidu.com/p/4126799131
eqn = Derivative[1, 0][u][t, y] == \[Nu]*Derivative[0, 2][u][t, y];
bc1 = u[0, y] == 0;
bc2 = u[t, 0] == Piecewise[{{U, t > 0}}, 0];
eqnlp = LaplaceTransform[{eqn, bc2}, t, s] /. Rule @@ bc1 /.
HoldPattern[LaplaceTransform[a_, t, s]] :> a
solla = FullSimplify[DSolve[eqnlp, u[t, y], y]][[1]]
mid = Collect[u[t, y] /. solla, E^(-((Sqrt[s] y)/Sqrt[\[Nu]]))]
InverseLaplaceTransform[mid /. C[1] -> 0, s, t]
我自以为那帖是的解说已经很详尽了,如果还是有疑问请具体说说哪里不懂。
2018年11月30日 18点11分 8
1