解常微分方程组无法得到正确解，求帮助！！！

sina2532

新手用COMSOL，想尝试解一个常微分方程组，但是解的结果和用Matlab解出来的解不一致，也没有找到原因

附上源文件，word中是常微分方程和在Matlab中模拟的结果

cmut

tao_c 的单位是 秒吗？
我看到 COMSOL 里面显示  m.s ?

sina2532

cmut 发表于 2013-6-26 09:58
tao_c 的单位是 秒吗？
我看到 COMSOL 里面显示  m.s ?

确实错了，单位应该是 (秒)，我更正了一下，重新传上来了，多谢哈！！

不过这个错误不怎么影响结果，我觉得是不是solver需要一些设置？

cmut

个人建议，global parameters 先直接给数值，看ODE公式，似乎没有问题。
solver设置，不太清楚。
祝好运！

TBE_Legend

用 generalized alpha 积分方法，设置 高频放大系数为0.9， 设置不长最大为 1e-9 试试

TBE_Legend

TBE_Legend 发表于 2013-6-27 16:18
用 generalized alpha 积分方法，设置 高频放大系数为0.9， 设置不长最大为 1e-9 试试
...

油条了下

sina2532

TBE_Legend 发表于 2013-6-27 16:49
油条了下

按照版主的提示我改了下，怎么老是给出解不收敛的错误

Nonlinear solver did not converge.
Time : 4.317170537135947e-006
Last time step is not converged.

求指导

TBE_Legend

sina2532 发表于 2013-6-27 21:55
按照版主的提示我改了下，怎么老是给出解不收敛的错误

Nonlinear solver did not converge.

提高迭代次数到1e9

TBE_Legend

TBE_Legend 发表于 2013-6-27 22:42
提高迭代次数到1e9

抽空又调了下，这个看起来好多了。


以下是mathematica代码供参考：

1. Rp = 479 10^25;
   
2. c = 2998 10^5;
   
3. \[Sigma] = 65 10^-24;
   
4. \[Tau]f = 230 10^-6;
   
5. \[Tau]c = 227 10^-11;
   
6. sol = NDSolve[ {n'[t] ==
   
7.      Rp - c \[Sigma] \[Phi][t] n[t] - n[t]/\[Tau]f, \[Phi]'[t] ==
   
8.      c \[Sigma] \[Phi][t] n[t] - \[Phi][t]/\[Tau]c ,
   
9.     n[0] == 0, \[Phi][0] == 1  }, {n, \[Phi]}, {t, 0, 70 10^-6},
   
10.    MaxStepSize -> 10^-8, Method -> "StiffnessSwitching"];
    
11. Plot[n[t] /. sol // Evaluate, {t, 0, 70 10^-6}, PlotRange -> All,
    
12. Frame -> True, PlotStyle -> {Red}]
    
13. Plot[\[Phi][t] /. sol // Evaluate, {t, 0, 70 10^-6}, PlotRange -> All,
    
14.   Frame -> True, PlotStyle -> {Black}]

复制代码

sina2532

前段时间一直忙着其他事情，没有再处理这件事情了。首先感谢楼主的热情帮助，您的方法给了我很多启发！

下面是我用Time Explicit Solver做的解，比Time Dependent Solver的解药更好，速度也更快。

附上求解结果的图，计算的文件用的是Comsol 4.3b，要打开附件中的文件可能需要相同的版本，不好意思哈！

flatime

看看你的结果
对了

zraforgood

学习了，刚学习comsol。

yaogang0209

学习一下

qitiandasheng

瞅瞅看！！！！111

自强厚德

学习一下

g574100239

正在学用comsol求解PDE和ODE，很有帮助

pp792696970

学校学习中

s1302563090

检查微分方程有没有错误