NDSolve求解偏微分方程组求助:
mathematica吧
全部回复
仅看楼主
level 6
纯小树 楼主
如下形式的偏微分方程组:
u[x,t]与s[t]的耦合,由于是偏微分方程和常微分方程的耦合,求解会报错。
具体情况如下截图,后附截图对应的代码:
NDSolve[{
\!\(
\*SubscriptBox[\(\[PartialD]\), \(t\)]\(v[\[Xi], t]\)\) -
s'[t]/s[t] \[Xi] \!\(
\*SubscriptBox[\(\[PartialD]\), \(\[Xi]\)]\(v[\[Xi], t]\)\) -
1/s[t]^2 \!\(
\*SubscriptBox[\(\[PartialD]\), \(\[Xi], \[Xi]\)]\(v[\[Xi],
t]\)\) == \[Xi] s[t] E^t + 2,
s'[t] == -(1/s[t])
\!\(\*SuperscriptBox[\(v\),
TagBox[
RowBox[{"(",
RowBox[{"0", ",", "1"}], ")"}],
Derivative],
MultilineFunction->None]\)[1, t],
v[0, t] == 0,
v[1, t] == 0,
v[\[Xi], 0] == \[Xi] (1 - \[Xi]),
s[0] == 1
}, {s, v}, {\[Xi], 0, 1}, {t, 0, 0.1}]
2020年04月13日 13点04分 1
吧务
level 10
要不先令s和v具有相同的参数,作PDE求解?PDE边界条件用NeumannValue和DirichletCondition表示。
2020年04月14日 04点04分 3
感谢回答!我确实先令s和v具有同样的参数,作PDE求解了。边界条件也和你建议的方法一样。但是结果还是报错,报错内容是mathematica不能求解当前这种时滞的偏微分方程。
2020年04月14日 13点04分
@纯小树 哦,你第三个方程和Derivative[0,1][v][1,t]有关,这个量会随计算过程而变化,传统积分方法不适用。软件提示不支持,学界有没有解决方法我也不知道。。。你再找找吧
2020年04月15日 15点04分
@asdasd1dsadsa [大拇指][大拇指]对的,就是这个地方,不是v[x,t],所以是时滞的。这个问题数学上好像还是有解析解的。。[狂汗]
2020年04月16日 00点04分
吧务
level 15
2020年05月02日 04点05分 4
1