非线性扩散方程的有限差分解：非线性扩散方程的数值解并创建结果影片-matlab开发

求解非线性扩散方程，该方程可以线性化，如 Richtmyer & Morton [1] 中的一般非线性扩散方程所示。 该方法是将 pde 线性化并应用 Crank-Nicolson 隐式有限差分格式以数值方式求解方程。 Matlab 运行命令-------------------------- 类型： IsoFreeSurfaceSolver 解决 pde： ------------------- \frac{\partial h}{\partial t}=\frac{1}{12}\frac{\partial^2 h^4}{\partial x^2} pde 可应用于在重力作用下在水平基底上扩散的等温粘性流体流动 - Huppert [2]。 请注意，PDE 已无量纲化。 初始条件： t=0: h = (1 - x^2)_{+} + 10^-6（有预湿膜） 考虑到 x = 0