8 votos

minimizar la norma de un enrollamiento sobre un dominio

Según mis cálculos:

La función que minimiza $$\int_\Omega \|\operatorname{curl} f\|^2\,dx$$ should satisfy $$\operatorname{curl}(\operatorname{curl}f) = 0$$ everywhere on $\Omega$, provided $\operatorname{curl} f = 0$ on $\partial \Omega$.

He seguido el mismo tipo de cálculo que el uno que demuestra la argmin de $\int_\Omega \|\nabla f\|^2~dx$ debe satisfacer $\Delta f = 0$. Sin embargo, no estoy seguro de si mis cálculos son de derecha... nadie Podrá comprobar que, por favor ?

1) primero vamos a empezar con un funcional $$G(f) = \int_\Omega \|\operatorname{curl} f\|^2\,dx.$$

2) calculamos el $$V(f,h) = \lim_{\epsilon\rightarrow 0} \frac{G(f+\epsilon h)-G(f)}{\epsilon} = 2\int_\Omega \operatorname{curl} f\cdot\operatorname{curl} h\,dx.$$

3) utilizamos la identidad : $$\operatorname{div}(A\times B) = -A\cdot\operatorname{curl} B + B\cdot\operatorname{curl} A,$$ with $A=\operatorname{curl} f$ and $B=h$.

4) obtenemos $$\int_\Omega \operatorname{curl} f\cdot \operatorname{curl} h\:dx = -\int_\Omega \operatorname{div}(\operatorname{curl} f\times h)\,dx + \int_\Omega h\cdot \operatorname{curl}(\operatorname{curl}f)\,dx.$$

5) utilizamos el teorema de la divergencia para obtener : $$\int_\Omega \operatorname{curl} f\cdot \operatorname{curl} h\,dx = -\int_{\partial\Omega} \operatorname{curl} f\times h\,ds + \int_\Omega h\cdot\operatorname{curl}(\operatorname{curl}f)\,dx$$

6) Hemos asumido $\operatorname{curl} f = 0$$\partial\Omega$, por lo que el primer término es cero.

7) $V(f,h)$ debe ser igual a cero para todos los $h$ para la función que se reduce al mínimo, por lo que $$\int_\Omega h\cdot\operatorname{curl}(\operatorname{curl}f) dx = 0\quad\forall h,$$ which implies $\operatorname{curl}(\operatorname{curl}f) = 0$ localmente.

Supongo que esto reasonning puede estar equivocado en varios lugares..... o ¿es lo correcto?? Gracias!

4voto

Evan Anderson Puntos 118832

OP es principalmente perfecto hasta después de algunos partes.

Para el curl operador tenemos una similar de integración por partes de la fórmula del teorema de la divergencia: $$ \int_{\Omega}\nabla \cdot (u\nabla v) = \int_{\Omega}(u\Delta v + \nabla u\cdot \nabla v) = \int_{\partial \Omega} u\nabla v\cdot \nu \,dS, $$ donde $\nu$ es el exterior de la unidad normal a la frontera. El curl de integración por partes fórmula es: $$ \int_{\Omega}\nabla \cdot (\phi\times \psi) = \int_{\Omega}(\psi\cdot\nabla \times \phi \phi\cdot \nabla\times \psi) = \int_{\partial \Omega} \phi\times \psi\cdot \nu \,dS = \int_{\partial \Omega} \nu\times \phi\cdot \psi \,dS.\la etiqueta{1} $$

Dicen que ahora quieres minimizar (aquí utilizo $u$ $v$ en lugar de $f$$h$): $$ \mathcal{G}(v) = \int_{\Omega}|\nabla\times v|^2\,dx, $$ si el minimizer es $u$ y el variacional límite es también no hay ningún problema: $$ 0= \frac{d}{d\epsilon}\mathcal{G}(u+\epsilon v)\Bigg|_{\epsilon= 0} = 2\int_{\Omega} \nabla \times u\cdot \nabla \times v. $$

Ahora vamos a utilizar la fórmula (1), dejando $\phi = \nabla \times u$, e $\psi = v$: $$ \int_{\Omega}(v\cdot\nabla \times (\nabla \times u)- \nabla \times u\cdot \nabla\times v) =\left\{\begin{aligned} \int_{\partial \Omega} \big(\nu\times (\nabla \times u)\big)\cdot v\,dS \\ \\ \int_{\partial \Omega} (v\times \nu) \cdot (\nabla \times u)\,dS \end{aligned}\right. $$ Hay dos opciones de modo de escribir límite de los términos a ser cero:

  1. Se pueden establecer para el minimizer $\color{blue}{\nu\times (\nabla \times u)|_{\partial \Omega} = 0}$, y sólo tangencial parte es suficiente. $\nabla \times u|_{\partial \Omega = 0}$ es demasiado fuerte.

  2. Puede establecer para la función de prueba de espacios de $\color{red}{\nu\times v = 0}$, $\nu \times u$ puede ser prescrito como cualquier permitido datos de límite de $g$ (aquellos que tienen 0 de la superficie de la divergencia).

La elección de azul, usted tiene que modular de un gradiente de una función adecuada de espacio $H$ a configurar una bien definida problema: $$ \left\{\begin{aligned} \nabla\times (\nabla \times u) &=0 \quad \text{ in } \Omega, \\ \nu \times (\nabla \times u) &= 0\quad \text{ on } \partial \Omega. \end{aligned}\right.\la etiqueta{2} $$ Para $u \in H/\nabla p$ el cociente del espacio, lo que significa que no necesitamos decirle a la diferencia entre el$u $$u+\nabla p$, por encima de problema tiene una solución única. La única solución de (2) es cero, un.k.una. $u = 0$ o decir $u = \nabla p$. Una vez que añadimos un indicador de condición de $\nabla \cdot u = f$, (2) es equivalente al problema de Neumann para $p$: $$ \left\{\begin{aligned} \Delta p &= f \quad \text{ in } \Omega, \\ \nabla p\cdot \nu &= 0\quad \text{ on } \partial \Omega. \end{aligned}\right. $$

La elección de rojo, el problema es: $$ \left\{\begin{aligned} \nabla\times (\nabla \times u) &=0 \quad \text{ in } \Omega, \\ \nu \times u &= g\quad \text{ on } \partial \Omega. \end{aligned}\right.\la etiqueta{3} $$ De nuevo, este problema tiene una única solución en $H/\nabla p$, excepto que esta vez $p$ es constante en el límite que implica la $\nabla p\times \nu =0$. Si establecemos $g$ 0, luego por encima de problema tiene una solución cero, lo que significa que $u = \nabla p$. Ahora un indicador de condición de $\nabla \cdot u = f$, a continuación, (3) puede ser transformado a un Dirichlet fronteras problema para $p$: $$ \left\{\begin{aligned} \Delta p &= f \quad \text{ in } \Omega, \\ p &= \mathrm{Constant}\quad \text{ on } \partial \Omega. \end{aligned}\right. $$


Resumen: el minimizer $u = \mathrm{arg}\min\mathcal{G}(v)$ es un gradiente de campo $\nabla p$, la condición de contorno de $p$ se basa en la condición límite de $u$. Tanto problema (2) y (3) no requieren de la componente normal de la $\nabla \times u$ o $u$ a imponerse. Si prescribe $\nabla\times u = 0$, en el límite, lo que implica $\nabla\times u\cdot \nu = 0$, el problema puede ser sobre-determinado.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X