封盖温度 - 使用 ode15 的二维热方程

问题描述

我正在尝试求解二维热方程,我正在使用 ode15 求解,有人指示我必须调整 dT/dt 方程。是否可以获得一些有关我应该如何处理的指导?

我遵循了 Matlab 中的这个示例:https://uk.mathworks.com/matlabcentral/fileexchange/59916-simple-heat-equation-solver?s_tid=srchtitle

怎么可能

2

示例中的 getRHS 函数的一部分要更改以限制其升高的温度?

解决方法

只需将其添加在 ode 之前作为 thr 的阈值 u0,例如:

thr=0.5; % just some # I picked,u0(u0>thr)=thr;