Al discretizar con FEM y resolver un problema de reacción-difusión, por ejemplo, con 0 < ε ≪ 1 (perturbación singular), la solución del problema discreto exhibirá típicamente capas oscilatorias cerca del límite. Con Ω = ( 0 , 1 ) , ε = 10 - 5 y elementos finitos lineales, la solución u h parece
Veo que hay mucha literatura sobre tales efectos no deseados cuando son causados por convección (p. Ej., Discretizaciones de viento), pero cuando se trata de reacción, las personas parecen centrarse en mallas refinadas (Shishkin, Bakhvalov).
¿Existen discretizaciones que eviten tales oscilaciones, es decir, que preserven la monotonicidad? ¿Qué más puede ser útil en este contexto?