Oscilaciones en problemas de reacción-difusión singularmente perturbados con elementos finitos

12

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

εΔu+u=1 on Ωu=0 on Ω
0<ε1Ω=(0,1)ε=105uh

solución de un problema singularmente perturbado

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?

Nico Schlömer
fuente
1
¿No es el esquema de diferencia central preservar la monotonicidad porque conduce a una matriz M ?
Hui Zhang
1ϕi,ϕj>0
@HuiZhang Tienes razón, por supuesto, en el caso de diferencias finitas (y también en volúmenes finitos). Adaptaré la respuesta para indicar más claramente que estoy interesado en elementos finitos.
Nico Schlömer
1
Los métodos discontinuos de Galerkin se han vuelto bastante populares para tales problemas. ¿Has mirado el libro de Di Pietro y Ern?
Christian Clason

Respuestas:

6

En el caso que muestra, la solución tiene una capa límite. Si no puede resolverlo porque su malla es demasiado gruesa, entonces, para todos los asuntos prácticos, la solución es discontinua al esquema numérico.

N

εh0

Wolfgang Bangerth
fuente
4

TL; DR: sus opciones son limitadas 1) ir adaptable a la fuerza bruta para una solución precisa y costosa 2) usar la difusión numérica para una solución menos precisa pero estable o (mi favorita) 3) aprovechar el hecho de que este es un problema de perturbación singular y resolver ¡Dos problemas internos / externos económicos y dejar que las asintóticas coincidentes hagan su magia!


δ=O(ϵ)

x=O(δ)η=x/δ

Δui+ui=1

u(0)=0ui(η)=uo(x0)uox=O(1)u1u0=1 solución interna con facilidad, en este caso incluso analíticamente.

De hecho, esta es la técnica que fue (y sigue siendo) muy popular para resolver problemas de capa límite laminar en mecánica de fluidos en el pasado. De hecho, si observa las ecuaciones de Navier-Stokes, en números altos de Reynolds, se enfrenta efectivamente a un problema de perturbación singular que, como el que mencionó aquí, desarrolla una capa límite (hecho divertido: los términos "capa límite" en perturbación el análisis en realidad proviene del problema de la capa límite fluida que acabo de describir).

u0=1

GradGuy
fuente