文件名称:非线性扩散方程的有限差分解:非线性扩散方程的数值解并创建结果影片-matlab开发
文件大小:394KB
文件格式:ZIP
更新时间:2024-06-19 00:08:25
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
【文件预览】:
IsothermalGravityCurrentSolver.zip
IsothermalGravityCurrent.zip