Steady Reaction-Diffusion with Dirichlet-Neumann boundary data. P1 or P2 and with remeshing
\begin{array}{rcl} -\nu \,u_{xx}+a_0\,u & =&0 \quad\mbox{on } (0,1)\quad \nu <<a_0\\ u(0)=1,\;u_x(1) & =&0. \end{array}
These results use P2 finite elements and a row of triangles sitting on the 1d interval to use the freefem tools
![]() |
Results viewed with MATLAB |
![]() | ![]() |
Uniform mesh | computed solution |
![]() | ![]() |
Adapted mesh | computed solution |
download example : SteadyReaction-Diff-1d.edp and its MATLAB companion code plotreactiondiffusion.m