求助：关于满足三个边值条件的二维热传导模型求解（PDE）？
第一张图是求解模型所用的方程组，我想得到在那个矩形区域内的温度分布图（ContourPlot），用动画显示。
第二张图是我在Mathematica上的代码，但是它提示我说方程个数大于了未知量数，删掉1个条件又出现错误，求同志们帮忙。
代码在最下方（可复制）。
谢谢。
![a = 0.5*10^-3; b = 20*10^-3; \[Rho] = 11.3*10^3; c = 130; \[Lambda] = 300; i = 10; R = 0.01; H = 15; fun = \[Rho]*c*D[T[r, z, t], t] - \[Lambda]* Laplacian[T[r, z, t], {r, z}] - (i^2 R)/(\[Pi]*a^2*b); ic = T[r, z, 0] == 300.15; bc1 = fun == NeumannValue[H*(T[r, z, t] - T[r, z, 0]), z == 0]; bc2 = fun == NeumannValue[-H*(T[r, z, t] - T[r, z, 0]), r == a]; bc3 = fun == NeumannValue[-H*(T[r, z, t] - T[r, z, 0]), z == b]; sol = NDSolveValue[{bc2, bc1, bc3, ic}, T[r, z, t], {r, 0, a}, {z, 0, b}] Clear["Global`*"]](https://developer.qcloudimg.com/http-save/yehe-259069/46cca24b2c0c97a73d81897de863747b.png)
a = 0.5*10^-3; b = 20*10^-3; \[Rho] = 11.3*10^3; c = 130; \[Lambda] = 300; i = 10; R = 0.01; H = 15; fun = \[Rho]*c*D[T[r, z, t], t] - \[Lambda]* Laplacian[T[r, z, t], {r, z}] - (i^2 R)/(\[Pi]*a^2*b); ic = T[r, z, 0] == 300.15; bc1 = fun == NeumannValue[H*(T[r, z, t] - T[r, z, 0]), z == 0]; bc2 = fun == NeumannValue[-H*(T[r, z, t] - T[r, z, 0]), r == a]; bc3 = fun == NeumannValue[-H*(T[r, z, t] - T[r, z, 0]), z == b]; sol = NDSolveValue[{bc2, bc1, bc3, ic}, T[r, z, t], {r, 0, a}, {z, 0, b}] Clear["Global`*"]
回答 2
sed是流编辑器
s,%,$,g命令行的一部分是sed要执行的命令。
和开发者交流更多问题细节吧，去 写回答
相关文章
腾讯云开发者
