Windows – 运行NDSolve时内存耗尽

前端之家收集整理的这篇文章主要介绍了Windows – 运行NDSolve时内存耗尽前端之家小编觉得挺不错的,现在分享给大家,也给大家做个参考。
我碰到数学中的“无可用内存”错误消息.我明白“并行[]”不是(显然)会帮助我.也没有“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]

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

猜你在找的Windows相关文章