Windows – 运行NDSolve时内存耗尽

发布时间：2020-12-14 04:31:58 所属栏目：Windows 来源：网络整理

导读：我碰到数学中的“无可用内存”错误消息.我明白“并行[]”不是(显然)会帮助我.也没有“ClearSystemCache []”. 是什么赋予了？我只需要更多的RAM？我的代码 Needs["VectorAnalysis`"]Needs["DifferentialEquations`InterpolatingFunctionAnatomy`"];Clear[Eq

我碰到数学中的“无可用内存”错误消息.我明白“并行[]”不是(显然)会帮助我.也没有“ClearSystemCache []”.

是什么赋予了？我只需要更多的RAM？

我的代码

Needs["VectorAnalysis`"]
Needs["DifferentialEquations`InterpolatingFunctionAnatomy`"];
Clear[Eq4,EvapThickFilm,h,S,G,E1,K1,D1,VR,M,R]
Eq4[h_,{S_,G_,E1_,K1_,D1_,VR_,M_,R_}] := !(
*SubscriptBox[([PartialD]),(t)]h) + 
    Div[-h^3 G Grad[h] + 
      h^3 S Grad[Laplacian[h]] + (VR E1^2 h^3)/(D1 (h + K1)^3)
        Grad[h] + M (h/(1 + h))^2 Grad[h]] + E1/(
    h + K1) + (R/6) D[D[(h^2/(1 + h)),x] h^3,x] == 0;
SetCoordinates[Cartesian[x,y,z]];
EvapThickFilm[S_,R_] := 
  Eq4[h[x,t],{S,R}];
TraditionalForm[EvapThickFilm[S,R]];

L = 318; TMax = 10;
Off[NDSolve::mxsst];
Clear[Kvar];
Kvar[t_] :=  Piecewise[{{1,t <= 1},{2,t > 1}}]
(*Ktemp = Array[0.001+0.001#^2&,13]*)
hSol = h /. NDSolve[{
     (*S,E,K,D,M*)

     EvapThickFilm[1,3,0.1,7,0.01,160],h[0,t] == h[L,h[x,t] == h[x,L,(*h[x,0] == 1.1+Cos[x] Sin[2y] *)
     h[x,0] == 
      1 + (-0.25 Cos[2 [Pi] x/L] - 0.25 Sin[2 [Pi] x/L]) Cos[
         2 [Pi] y/L]
     },{x,L},{y,{t,TMax},MaxStepSize -> 0.1
    ][[1]]

hGrid = InterpolatingFunctionGrid[hSol];

错误信息

No more memory available.
Mathematica kernel has shut down.
Try quitting other applications and then retry.

我的操作系统规格

英特尔酷睿2双核与4.00 GB RAM,64位操作系统(Windows 7)

在这里,您可以了解发生的情况：

更换

MaxStepSize -> 0.1

通过

MaxStepFraction -> 1/30

并运行你的代码.

然后：

p = Join[#,Reverse@#]&@ 
     Table[Plot3D[hSol[x,i],PlotRange -> {All,All,{0,4}}],{i,8,.1}]

Export["c:plot.gif",p]

所以,妈妈正在努力改善这些高峰的解决方案,无济于事.

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!