Let $B_R$ be a ball of radius $R>0$ centered at the origin in $\mathbb{R}^2$.
Consider the PDE $u_t + \nabla \cdot (u^2 \nabla \Delta u) = 0$, in $ (0,T)\times B_R,$ with boundary conditions $u =\frac{\partial}{\partial\nu} \Delta u = 0$, in $(0,T)\times B_R,$ and initial data $u(0,\cdot) = g,$ in $ B_R$, where for instance
g[x] = 1
in $B_r$ and $0$ outside $B_r$, with $r < R$. A reference on numerical analysis for this problem is https://dml.cz/handle/10338.dmlcz/134511.
Is it possible to use Mathematica (for example NDSolve
) to obtain and plot solutions to such equation?
NDSolve
and friends? $\endgroup$NDSolve
or if usingNDSolve
is actually the correct approach. $\endgroup$NDSolve
for solving linearized PDEs in order to implement a semi-implicit time integrator... $\endgroup$