求解非线性扩散方程,该方程可以线性化,如 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
2023-05-05 10:22:48
394KB
matlab
1