Flujo de acortamiento de curvas - Curve-shortening flow

Convergencia de una curva convexa a un círculo bajo el flujo de acortamiento de la curva. Las curvas interiores (color más claro) son versiones fluidas de las curvas exteriores. Los pasos de tiempo entre curvas no son uniformes.

En matemáticas, el flujo de acortamiento de curvas es un proceso que modifica una curva suave en el plano euclidiano moviendo sus puntos perpendicularmente a la curva a una velocidad proporcional a la curvatura . El flujo de acortamiento de curvas es un ejemplo de flujo geométrico y es el caso unidimensional del flujo de curvatura media . Otros nombres para el mismo proceso incluyen flujo de acortamiento euclidiano , flujo de calor geométrico y evolución de la longitud del arco .

A medida que los puntos de cualquier curva cerrada simple suave se mueven de esta manera, la curva permanece simple y suave. Pierde área a una tasa constante y su perímetro disminuye lo más rápido posible para cualquier evolución continua de la curva. Si la curva no es convexa, su curvatura absoluta total disminuye monótonamente, hasta que se vuelve convexa. Una vez convexa, la relación isoperimétrica de la curva disminuye a medida que la curva converge a una forma circular, antes de colapsar en un solo punto de singularidad. Si dos curvas cerradas lisas simples disjuntas evolucionan, permanecen disjuntas hasta que una de ellas colapsa en un punto. El círculo es la única curva cerrada simple que mantiene su forma bajo el flujo de acortamiento de la curva, pero algunas curvas que se cruzan o tienen una longitud infinita mantienen su forma, incluida la curva de la Parca, una curva infinita que se traslada hacia arriba y espirales que giran. sin dejar de tener el mismo tamaño y forma.

Se puede calcular numéricamente una aproximación al flujo de acortamiento de la curva, aproximando la curva como un polígono y utilizando el método de diferencias finitas para calcular el movimiento de cada vértice del polígono. Los métodos alternativos incluyen calcular una convolución de vértices de polígono y luego volver a muestrear los vértices en la curva resultante, o aplicar repetidamente un filtro mediano a una imagen digital cuyos píxeles en blanco y negro representan el interior y el exterior de la curva.

El flujo de acortamiento de curvas se estudió originalmente como un modelo para el recocido de láminas de metal. Posteriormente, se aplicó en el análisis de imágenes para dar una representación de formas a múltiples escalas. También puede modelar sistemas de reacción-difusión y el comportamiento de autómatas celulares . El flujo de acortamiento de curvas se puede utilizar para encontrar geodésicas cerradas en variedades de Riemann y como modelo para el comportamiento de flujos de dimensiones superiores.

Definiciones

Un flujo es un proceso en el que los puntos de un espacio cambian continuamente sus ubicaciones o propiedades a lo largo del tiempo. Más específicamente, en un flujo geométrico unidimensional como el flujo de acortamiento de curvas, los puntos que experimentan el flujo pertenecen a una curva , y lo que cambia es la forma de la curva, su incrustación en el plano euclidiano determinado por las ubicaciones de cada uno. de sus puntos. En el flujo de acortamiento de curvas, cada punto de una curva se mueve en la dirección de un vector normal a la curva, a una tasa proporcional a la curvatura . Para una curva evolutiva representada por una función de dos parámetros C ( s , t ) donde s parametriza la longitud del arco a lo largo de la curva y t parametriza un tiempo en la evolución de la curva, el flujo de acortamiento de la curva se puede describir mediante el parcial parabólico ecuación diferencial

una forma de la ecuación de calor , donde κ es la curvatura yn es el vector normal unitario.

Debido a que los ingredientes de esta ecuación, la longitud del arco, la curvatura y el tiempo, no se ven afectados por las traslaciones y rotaciones del plano euclidiano, se deduce que el flujo definido por esta ecuación es invariante bajo traslaciones y rotaciones (o más precisamente, equivariante ). . Si el avión se escala por un factor de dilatación constante, el flujo permanece esencialmente sin cambios, pero se ralentiza o acelera por el mismo factor.

Curvas no suaves

Para que el flujo esté bien definido, la curva dada debe ser lo suficientemente suave como para tener una curvatura continua. Sin embargo, una vez que comienza el flujo, la curva se vuelve analítica y permanece así hasta alcanzar una singularidad en la que la curvatura estalla. Para una curva suave sin cruces, la única singularidad posible ocurre cuando la curva colapsa a un punto, pero las curvas sumergidas pueden tener otros tipos de singularidad. En tales casos, con cierto cuidado es posible continuar el flujo más allá de estas singularidades hasta que toda la curva se contraiga a un solo punto.

Para una curva cerrada simple, utilizando una extensión del flujo a curvas no suaves según el método de ajuste de nivel , solo hay dos posibilidades. Las curvas con medida de Lebesgue cero (incluidos todos los polígonos y las curvas suaves por partes) evolucionan instantáneamente a curvas suaves, después de lo cual evolucionan como lo haría cualquier curva suave. Sin embargo, las curvas de Osgood con medidas distintas de cero evolucionan inmediatamente hacia un anillo topológico con un área distinta de cero y límites suaves. La curva sinusoidal del topólogo es un ejemplo que instantáneamente se vuelve suave, a pesar de que ni siquiera está conectado localmente ; ejemplos como este muestran que la evolución inversa del flujo de acortamiento de curvas puede llevar curvas de buen comportamiento a singularidades complicadas en una cantidad de tiempo finita.

Superficies no euclidianas

El flujo de acortamiento de la curva, y muchos de los resultados sobre el flujo de acortamiento de la curva, se pueden generalizar desde el plano euclidiano a cualquier variedad de Riemann de dos dimensiones . Para evitar tipos adicionales de singularidad, es importante que la variedad sea convexa en el infinito ; esto se define en el sentido de que todo conjunto compacto tiene un casco convexo compacto , como se define mediante convexidad geodésica . El flujo de acortamiento de la curva no puede hacer que una curva se aparte de su casco convexo, por lo que esta condición evita que partes de la curva alcancen el límite del colector.

Curvas espaciales

El flujo de acortamiento de curvas también se ha estudiado para curvas en el espacio euclidiano tridimensional . El vector normal en este caso se puede definir (como en el plano) como la derivada del vector tangente con respecto a la longitud del arco, normalizado para ser un vector unitario; es uno de los componentes del marco Frenet-Serret . No está bien definido en los puntos de curvatura cero, pero el producto de la curvatura y el vector normal permanece bien definido en esos puntos, lo que permite definir el flujo de acortamiento de la curva. Las curvas en el espacio pueden cruzarse entre sí o entre sí de acuerdo con este flujo, y el flujo puede conducir a singularidades en las curvas; cada singularidad es asintótica con respecto a un plano. Sin embargo, se sabe que las curvas esféricas y las curvas que pueden proyectarse ortogonalmente en una curva plana convexa regular siguen siendo simples. El flujo de acortamiento de curvas para curvas espaciales se ha utilizado como una forma de definir el flujo pasado singularidades en curvas planas.

Más allá de las curvas

Es posible extender la definición del flujo a entradas más generales que las curvas, por ejemplo, usando varifolds rectificables o el método de ajuste de nivel . Sin embargo, estas definiciones extendidas pueden permitir que partes de curvas desaparezcan instantáneamente o engorden en conjuntos de áreas distintas de cero.

Para redes de curvas, extender el flujo de acortamiento de curvas más allá de una singularidad puede resultar en ambigüedad o engorde.

Una variación del problema comúnmente estudiada involucra redes de curvas suaves interiores disjuntas, con uniones donde se encuentran tres o más de las curvas. Cuando todas las uniones tienen exactamente tres curvas que se encuentran en ángulos de 2 π / 3 (las mismas condiciones que se observan en un árbol de Steiner óptimo o una espuma bidimensional de pompas de jabón ), el flujo está bien definido a corto plazo. Sin embargo, eventualmente puede alcanzar un estado singular con cuatro o más curvas que se encuentran en un cruce, y puede haber más de una forma de continuar el flujo más allá de dicha singularidad.

Comportamiento

Principio de evitación, radio y factor de estiramiento

Si dos curvas cerradas simples lisas disjuntas se someten al flujo de acortamiento de la curva simultáneamente, permanecen disjuntas a medida que avanza el flujo. La razón es que, si dos curvas suaves se mueven de una manera que crea un cruce, entonces, en el momento del primer cruce, las curvas necesariamente serían tangentes entre sí, sin cruzarse. Pero, en tal situación, las curvaturas de las dos curvas en el punto de tangencia necesariamente las separarían en lugar de juntarlas en un cruce. Por la misma razón, una sola curva cerrada simple nunca puede evolucionar para cruzarse a sí misma. Este fenómeno se conoce como principio de evitación.

El principio de evitación implica que cualquier curva cerrada suave debe eventualmente alcanzar una singularidad, como un punto de curvatura infinita. Porque, si una curva suave C dada está rodeada por un círculo, ambos permanecerán separados mientras existan. Pero el círculo circundante se contrae bajo el flujo de curvatura, permaneciendo circular, hasta que colapsa, y por el principio de evitación C debe permanecer contenido dentro de él. Entonces, si C nunca alcanzara una singularidad, quedaría atrapado en un solo punto en el momento en que el círculo colapsa, lo que es imposible para una curva suave. Esto se puede cuantificar observando que el radio del círculo más pequeño que encierra a C debe disminuir a una velocidad que sea al menos tan rápida como la disminución del radio de un círculo que experimenta el mismo flujo.

Huisken (1998) cuantifica el principio de evitación para una sola curva en términos de la relación entre la longitud del arco (del más corto de dos arcos) y la distancia euclidiana entre pares de puntos, a veces llamada factor de estiramiento . Muestra que el factor de estiramiento es estrictamente decreciente en cada uno de sus máximos locales, excepto en el caso de los dos extremos de un diámetro de un círculo, en cuyo caso el factor de estiramiento es constante en π . Esta propiedad de monotonicidad implica el principio de evitación, ya que si la curva llegara a tocarse a sí misma, el factor de estiramiento se volvería infinito en los dos puntos de contacto.

Largo

A medida que una curva experimenta el flujo de acortamiento de la curva, su longitud L disminuye a una tasa dada por la fórmula

donde el intervalo se toma sobre la curva, κ es la curvatura y s es la longitud del arco a lo largo de la curva. El integrando siempre es no negativo, y para cualquier curva cerrada suave existen arcos dentro de los cuales es estrictamente positivo, por lo que la longitud disminuye monótonamente. De manera más general, para cualquier evolución de curvas cuya velocidad normal es f , la tasa de cambio de longitud es

que puede interpretarse como un producto interno negado entre la evolución dada y el flujo de acortamiento de la curva. Por lo tanto, el flujo de acortamiento de la curva se puede describir como el flujo de gradiente por longitud, el flujo que (localmente) disminuye la longitud de la curva lo más rápidamente posible en relación con la norma L 2 del flujo. Esta propiedad es la que da nombre al flujo acortador de curvas.

Zona

Para una curva cerrada simple, el área encerrada por la curva se contrae, a una tasa constante de 2 π unidades de área por unidad de tiempo, independientemente de la curva. Por lo tanto, el tiempo total para que una curva se contraiga a un punto es proporcional a su área, independientemente de su forma inicial. Debido a que el área de una curva se reduce a una tasa constante, y (por la desigualdad isoperimétrica ) un círculo tiene la mayor área posible entre las curvas cerradas simples de una longitud dada, se deduce que los círculos son las curvas más lentas para colapsar a un punto bajo el flujo de acortamiento de la curva. Todas las demás curvas tardan menos en colapsar que un círculo de la misma longitud.

La tasa constante de reducción del área es la única ley de conservación que cumple el flujo que acorta la curva. Esto implica que no es posible expresar el "punto de fuga" donde la curva eventualmente colapsa como una integral sobre la curva de cualquier función de sus puntos y sus derivadas, porque tal expresión conduciría a una segunda ley de conservación prohibida. Sin embargo, al combinar la tasa constante de pérdida de área con el principio de evitación, es posible probar que el punto de fuga siempre se encuentra dentro de un círculo, concéntrico con el círculo circundante mínimo, cuya área es la diferencia de áreas entre el círculo circundante y el círculo circundante. curva dada.

Curvatura absoluta total

La curvatura absoluta total de una curva suave es la integral del valor absoluto de la curvatura a lo largo de la longitud del arco de la curva,

También se puede expresar como una suma de los ángulos entre los vectores normales en pares consecutivos de puntos de inflexión . Es 2 π para curvas convexas y mayor para curvas no convexas, y sirve como medida de la no convexidad de una curva.

El flujo de acortamiento de curvas no puede crear nuevos puntos de inflexión. Cada uno de los ángulos en la representación de la curvatura absoluta total como una suma disminuye monótonamente, excepto en los instantes en que dos puntos de inflexión consecutivos alcanzan el mismo ángulo o posición y ambos son eliminados. Por lo tanto, la curvatura absoluta total nunca puede aumentar a medida que evoluciona la curva. Para curvas convexas es constante en 2 π y para curvas no convexas disminuye monótonamente.

Teorema de Gage-Hamilton-Grayson

Si una curva cerrada simple y suave se somete al flujo de acortamiento de la curva, permanece incrustada sin problemas sin auto-intersecciones. Eventualmente se volverá convexo , y una vez que lo haga, seguirá siendo convexo. Después de este tiempo, todos los puntos de la curva se moverán hacia adentro y la forma de la curva convergerá en un círculo a medida que toda la curva se contraiga a un solo punto. Este comportamiento a veces se resume diciendo que cada curva cerrada simple se reduce a un "punto redondo".

Este resultado se debe a Michael Gage , Richard S. Hamilton y Matthew Grayson. Gage ( 1983 , 1984 ) demostró la convergencia a un círculo para curvas convexas que se contraen en un punto. Más específicamente, Gage mostró que la relación isoperimétrica (la relación entre la longitud de la curva al cuadrado y el área, un número que es 4 π para un círculo y más grande para cualquier otra curva convexa) disminuye monótona y rápidamente. Gage y Hamilton (1986) demostraron que todas las curvas convexas suaves eventualmente se contraen hasta un punto sin formar ninguna otra singularidad, y Grayson (1987) demostró que toda curva no convexa eventualmente se volverá convexa. Andrews y Bryan (2011) proporcionan una prueba más sencilla del resultado de Grayson, basada en la monotonicidad del factor de estiramiento.

La forma límite para todas las redes de dos rayos colineales y dos curvas que conectan los puntos finales de los dos rayos. La lente central tiene la forma de una vesica piscis .

Se pueden extender resultados similares desde curvas cerradas a curvas ilimitadas que satisfagan una condición de Lipschitz local . Para tales curvas, si ambos lados de la curva tienen un área infinita, entonces la curva evolucionada permanece suave y libre de singularidades para siempre. Sin embargo, si un lado de una curva ilimitada tiene un área finita, y la curva tiene una curvatura absoluta total finita, entonces su evolución alcanza una singularidad en el tiempo proporcional al área en el lado de área finita de la curva, con una curvatura ilimitada cerca de la singularidad. . Para curvas que son gráficas de funciones suficientemente bien comportadas, asintóticas a un rayo en cada dirección, la solución converge en forma a una forma única que es asintótica a los mismos rayos. Para las redes formadas por dos rayos disjuntos en la misma línea, junto con dos curvas suaves que conectan los puntos finales de los dos rayos, se cumple un análogo del teorema de Gage-Hamilton-Grayson, según el cual la región entre las dos curvas se vuelve convexa y luego converge a una forma de vesica piscis .

Singularidades de las curvas autocruzadas

Las curvas que tienen autocruces pueden alcanzar singularidades antes de contraerse hasta un punto. Por ejemplo, si una lemniscata (cualquier curva suave sumergida con un solo cruce, que se asemeja a una figura de 8 o símbolo de infinito ) tiene áreas desiguales en sus dos lóbulos, entonces eventualmente el lóbulo más pequeño colapsará en un punto. Sin embargo, si los dos lóbulos tienen áreas iguales, permanecerán iguales a lo largo de la evolución de la curva y la relación isoperimétrica divergerá a medida que la curva colapsa a una singularidad.

Cuando una curva autocruzada localmente convexa se acerca a una singularidad a medida que uno de sus bucles se contrae, se contrae de manera auto-similar o se acerca asintóticamente a la curva de la Parca (descrita a continuación) a medida que se contrae. Cuando un bucle colapsa a una singularidad, la cantidad de curvatura absoluta total que se pierde es al menos 2 π o exactamente π .

Sobre variedades de Riemann

En una variedad de Riemann, cualquier curva cerrada simple y suave seguirá siendo suave y simple a medida que evoluciona, al igual que en el caso euclidiano. Se colapsará hasta un punto en un período de tiempo finito, o permanecerá suave y simple para siempre. En el último caso, la curva necesariamente converge a una geodésica cerrada de la superficie.

Las curvas sumergidas en las variedades riemannianas, con un número finito de autocruces, se vuelven auto-tangentes solo en un conjunto discreto de momentos, en cada uno de los cuales pierden un cruce. Como consecuencia, el número de puntos de autocruzamiento no aumenta.

El acortamiento de la curva en una esfera se puede utilizar como parte de una demostración del teorema de la pelota de tenis . Este teorema establece que cada curva cerrada simple y suave en la esfera que divide la superficie de la esfera en dos áreas iguales (como la costura de una pelota de tenis ) debe tener al menos cuatro puntos de inflexión . La prueba proviene de la observación de que el acortamiento de la curva conserva las propiedades de suavidad y bisección del área de la curva y no aumenta su número de puntos de inflexión. Por lo tanto, permite que el problema se reduzca al problema de las curvas cercanas a la forma límite del acortamiento de la curva, un círculo máximo .

La fórmula de la monotonicidad de Huisken

De acuerdo con la fórmula de monotonicidad de Huisken , la convolución de una curva en evolución con un núcleo de calor invertido en el tiempo no aumenta. Este resultado se puede utilizar para analizar las singularidades de la evolución.

Curvas específicas

Curvas con evolución auto-similar

La curva de la Parca y sus copias traducidas producidas por el flujo de acortamiento de la curva

Debido a que todas las demás curvas cerradas simples convergen en un círculo, el círculo es la única curva cerrada simple que mantiene su forma bajo el flujo de acortamiento de la curva. Sin embargo, hay muchos otros ejemplos de curvas que no son simples (incluyen autocruces) o no cerradas (se extienden hasta el infinito) y mantienen su forma. En particular,

  • Cada línea permanece sin cambios por el flujo de acortamiento de la curva. Las líneas son las únicas curvas que no se ven afectadas por el flujo de acortamiento de curvas, aunque existen redes de curvas estables más complejas, como el mosaico hexagonal del plano.
  • La curva de la Parca y = - log cos x se mueve hacia arriba sin cambiar su forma. De la misma manera, cualquier curva similar a la parca es traducida por el flujo de acortamiento de la curva, desplazado en la dirección del eje de simetría de la curva sin cambiar su forma u orientación. La parca es la única curva con esta propiedad. También se le llama modelo de horquilla en la literatura de física.
  • Una familia de curvas cerradas que se autocruzan, derivadas de las proyecciones de los nudos del toro , se encogen homotéticamente pero siguen siendo auto-similares bajo el flujo de acortamiento de la curva. Éstas han llegado a conocerse como curvas de Abresch-Langer , después del trabajo de Abresch y Langer (1986) , aunque fueron mencionadas anteriormente por Mullins (1956) y redescubiertas de forma independiente por Epstein y Weinstein (1987) . Estas curvas son localmente convexas y, por lo tanto, pueden describirse por sus funciones de soporte . Las versiones adecuadamente escaladas de estas funciones de soporte obedecen a la ecuación diferencial
que tiene soluciones periódicas positivas (correspondientes a curvas con evolución auto-similar) para cualquier período que esté estrictamente entre π y .
  • Otras curvas, incluidas algunas espirales infinitas , siguen siendo auto-similares con movimientos más complicados que incluyen rotación o combinaciones de rotación, contracción o expansión y traslación.
  • Para redes de curvas suaves, que se encuentran de tres en tres en uniones con ángulos de 2 π / 3, las soluciones de contracción auto-similares incluyen una burbuja doble que rodea dos áreas iguales, una forma de lente ( vesica piscis ) limitada por dos arcos de círculos congruentes junto con dos rayos colineales que tienen sus vértices en las esquinas de la lente, y una red "en forma de pez" delimitada por un segmento de línea, dos rayos y una curva convexa. Cualquier otra red de contracción auto-similar implica un mayor número de curvas. Otra familia de redes crece homotéticamente y sigue siendo auto-similar; se trata de redes de curvas en forma de árbol, que se encuentran en ángulos de 2 π / 3 en uniones triples, asintóticas a un abanico de dos o más rayos que se encuentran en un punto final común. El caso de dos rayos de estas formas es una curva suave ilimitada; para tres o más rayos, la evolución de estas formas puede definirse utilizando variantes generalizadas del flujo acortador de curvas como el de los pliegues variables. Un abanico dado de cuatro o más rayos puede ser asintótico a más de una solución diferente de este tipo, por lo que estas soluciones no proporcionan una definición única para el flujo de acortamiento de curvas a partir de un abanico de rayos.

Soluciones antiguas

Una antigua solución a un problema de flujo es una curva cuya evolución se puede extrapolar hacia atrás para todos los tiempos, sin singularidades. Todas las soluciones auto-similares que se encogen o permanecen del mismo tamaño en lugar de crecer son soluciones antiguas en este sentido; se pueden extrapolar hacia atrás invirtiendo la transformación de auto-similitud que sufrirían por el flujo hacia adelante de acortamiento de la curva. Así, por ejemplo, el círculo, la parca y las curvas de Abresch-Langer son soluciones antiguas.

También hay ejemplos que no son auto-similares. Un ejemplo explícito es la solución ovalada de Angenent después del trabajo de Angenent (1992) . Esta familia de curvas se puede parametrizar especificando la curvatura en función del ángulo tangente mediante la fórmula

y tienen como forma limitante en evolución inversa un par de curvas de la Parca que se acercan entre sí desde direcciones opuestas. En el sistema de coordenadas cartesianas , pueden estar dadas por la ecuación de curva implícita

En la literatura de física, las mismas formas se conocen como modelo de clip .

Las soluciones de Angenent ovalada y de círculo que se contrae son las únicas soluciones antiguas cuyas ranuras de tiempo limitaban conjuntos convexos acotados. Las soluciones de Grim Reaper, el medio espacio estacionario y las tiras estacionarias son los únicos ejemplos cuyos cortes de tiempo unen conjuntos convexos ilimitados. Existen muchos ejemplos adicionales (no incrustados ) localmente convexos, así como muchos otros ejemplos incrustados (no convexos).

Aproximaciones numéricas

Para calcular el flujo de acortamiento de la curva de manera eficiente, tanto una curva continua como la evolución continua de la curva deben reemplazarse por una aproximación discreta.

Seguimiento frontal

Los métodos de seguimiento frontal se han utilizado durante mucho tiempo en dinámica de fluidos para modelar y rastrear el movimiento de los límites entre diferentes materiales, de gradientes pronunciados en las propiedades de los materiales, como los frentes climáticos , o de las ondas de choque dentro de un solo material. Estos métodos implican derivar las ecuaciones de movimiento del límite y usarlas para simular directamente el movimiento del límite, en lugar de simular el fluido subyacente y tratar el límite como una propiedad emergente del fluido. Los mismos métodos también se pueden utilizar para simular el flujo de acortamiento de la curva, incluso cuando la curva que experimenta el flujo no es un límite o un choque.

En los métodos de seguimiento frontal para el acortamiento de curvas, la curva que experimenta la evolución se discretiza como un polígono. El método de diferencias finitas se usa para derivar fórmulas para el vector normal aproximado y la curvatura en cada vértice del polígono, y estos valores se usan para determinar cómo mover cada vértice en cada paso de tiempo. Aunque el flujo de acortamiento de la curva se define por el movimiento de una curva perpendicularmente a sí misma, algunas parametrizaciones del flujo de acortamiento de la curva pueden permitir que los vértices que se aproximan a la curva se muevan de forma no perpendicular. En efecto, esto permite que los vértices se muevan a lo largo de la curva, a medida que la curva evoluciona. La elección de una reparametrización cuidadosa puede ayudar a redistribuir los vértices de manera más uniforme a lo largo de la curva en situaciones en las que el movimiento perpendicular haría que se agruparan. Merriman, Bence y Osher (1992) escriben que estos métodos son rápidos y precisos, pero que es mucho más complicado extenderlos a versiones del flujo de acortamiento de curvas que se aplican a entradas más complicadas que las curvas cerradas simples, donde es necesario lidiar con singularidades y cambios de topología.

Para la mayoría de estos métodos, Cao (2003) advierte que "las condiciones de estabilidad no se pueden determinar fácilmente y el intervalo de tiempo debe elegirse ad hoc". Otro método de diferenciación finita de Crandall & Lions (1996) modifica la fórmula de la curvatura en cada vértice añadiéndole un término pequeño basado en el operador de Laplace . Esta modificación se denomina regularización elíptica y puede utilizarse para ayudar a probar la existencia de flujos generalizados, así como en su simulación numérica. Utilizándolo, se puede demostrar que el método de Crandall y Lions converge y es el único método numérico enumerado por Cao que está equipado con límites en su tasa de convergencia. Para una comparación empírica de los métodos de diferencia finita de Euler hacia adelante , Euler hacia atrás y Crank-Nicolson más precisos , consulte Balažovjech & Mikula (2009) .

Convolución remuestreada

Mokhtarian y Mackworth (1992) sugieren un método numérico para calcular una aproximación al flujo de acortamiento de la curva que mantiene una aproximación discreta a la curva y alterna entre dos pasos:

  • Vuelva a muestrear la curva actual colocando nuevos puntos de muestra con un espaciado uniforme, medido por la longitud de arco normalizada.
  • Convierta las ubicaciones de los puntos con una función gaussiana con una pequeña desviación estándar, reemplazando de hecho la ubicación de cada punto con un promedio ponderado de las ubicaciones de los puntos cercanos a lo largo de la curva, con pesos gaussianos. La desviación estándar del gaussiano debe elegirse para que sea lo suficientemente pequeña como para que, después de este paso, los puntos de muestra todavía tengan un espaciado casi uniforme.

Como muestran, este método converge a la distribución de acortamiento de la curva en el límite a medida que aumenta el número de puntos de muestra y se contrae la longitud del arco normalizado del radio de convolución.

Filtrado de la mediana

Merriman, Bence y Osher (1992) describen un esquema que opera en una cuadrícula cuadrada bidimensional, efectivamente una matriz de píxeles . La curva a evolucionar se representa asignando el valor 0 (negro) a los píxeles exteriores a la curva, y 1 (blanco) a los píxeles interiores a la curva, dando la función de indicador para el interior de la curva. Esta representación se actualiza alternando dos pasos:

  • Convierta la imagen pixelada con un núcleo de calor para simular su evolución bajo la ecuación de calor durante un breve paso de tiempo. El resultado es un desenfoque gaussiano de la imagen o, de manera equivalente, la transformada de Weierstrass de la función indicadora, con un radio proporcional a la raíz cuadrada del paso de tiempo.
  • Establezca cada píxel con valor numérico menor que 1/2 a 0, y cada píxel con valor numérico mayor que 1/2 a 1, limitando la imagen a sus valores originales en nuevas posiciones.

Para que este esquema sea preciso, el intervalo de tiempo debe ser lo suficientemente grande como para que la curva se mueva al menos un píxel incluso en puntos de curvatura baja, pero lo suficientemente pequeño como para que el radio de desenfoque sea menor que el radio mínimo. de curvatura. Por lo tanto, el tamaño de un píxel debe ser O (min κ / max κ 2 ) , lo suficientemente pequeño como para permitir que se elija un intervalo de tiempo intermedio adecuado.

El método se puede generalizar a la evolución de redes de curvas, encontrándose en los cruces y dividiendo el plano en más de tres regiones, aplicando el mismo método simultáneamente a cada región. En lugar de difuminar y establecer umbrales, este método se puede describir alternativamente como la aplicación de un filtro mediano con pesos gaussianos a cada píxel. Es posible utilizar núcleos distintos al núcleo de calor, o refinar adaptativamente la cuadrícula para que tenga una alta resolución cerca de la curva pero no pierda tiempo y memoria en píxeles alejados de la curva que no contribuyen al resultado. En lugar de usar solo los dos valores en la imagen pixelada, una versión de este método que usa una imagen cuyos valores de píxeles representan la distancia firmada a la curva puede lograr una precisión de subpíxeles y requerir una resolución más baja.

Aplicaciones

Recocido de láminas de metal

Una referencia temprana al flujo de acortamiento de curvas de William W. Mullins  ( 1956 ) lo motiva como modelo para el proceso físico de recocido , en el que el tratamiento térmico hace que los límites entre los granos de metal cristalizado se desplacen. A diferencia de las películas de jabón , que se ven obligadas por las diferencias en la presión del aire a convertirse en superficies de curvatura media constante , los límites de grano en el recocido están sujetos solo a efectos locales, que hacen que se muevan de acuerdo con el flujo de curvatura media. El caso unidimensional de este flujo, el flujo de acortamiento de curvas, corresponde a hojas de recocido de metal que son lo suficientemente delgadas como para que los granos se vuelvan efectivamente bidimensionales y sus límites unidimensionales.

Análisis de formas

En procesamiento de imágenes y visión por computadora , Mokhtarian y Mackworth (1992) sugieren aplicar el flujo de acortamiento de curvas al contorno de una forma derivada de una imagen digital, para eliminar el ruido de la forma y proporcionar un espacio de escala que proporcione una descripción simplificada. de la forma en diferentes niveles de resolución. El método de Mokhtarian y Mackworth implica calcular el flujo de acortamiento de la curva, rastrear los puntos de inflexión de la curva a medida que avanzan a través del flujo y dibujar un gráfico que traza las posiciones de los puntos de inflexión alrededor de la curva contra el parámetro de tiempo. Los puntos de inflexión normalmente se eliminarán de la curva en pares a medida que la curva se vuelve convexa (según el teorema de Gage-Hamilton-Grayson) y la vida útil de un par de puntos corresponde a la prominencia de una característica de la forma. Debido al método de convolución remuestreada que describen para calcular una aproximación numérica del flujo de acortamiento de la curva, llaman a su método el espacio de escala de curvatura remuestreada . Observan que este espacio de escala es invariante bajo las transformaciones euclidianas de la forma dada, y afirman que determina de forma única la forma y es robusto contra pequeñas variaciones en la forma. Lo comparan experimentalmente con varias definiciones alternativas relacionadas de un espacio de escala para formas, y encuentran que el espacio de escala de curvatura remuestreado es menos intensivo en computación, más robusto contra ruido no uniforme y menos influenciado por diferencias de forma a pequeña escala.

Reacción-difusión

En los sistemas de reacción-difusión modelados por la ecuación de Allen-Cahn , el comportamiento limitante para la reacción rápida, la difusión lenta y dos o más mínimos locales de energía con el mismo nivel de energía es que el sistema se asiente en regiones de diferentes niveles locales. mínimos, con los frentes que delimitan los límites entre estas regiones evolucionando de acuerdo con el flujo de acortamiento de la curva.

Autómata celular

El autómata celular Anneal, 1600 pasos después de un inicio aleatorio

En un autómata celular , cada celda en una cuadrícula infinita de celdas puede tener uno de un conjunto finito de estados, y todas las celdas actualizan sus estados simultáneamente basándose solo en la configuración de un pequeño conjunto de celdas vecinas. Una regla de autómata celular realista es aquella en la que la cuadrícula es la celosía cuadrada infinita, hay exactamente dos estados de celda, el conjunto de vecinos de cada celda son los ocho vecinos de la vecindad de Moore y la regla de actualización depende solo de la número de vecinos con cada uno de los dos estados en lugar de una función más complicada de esos estados. En una regla particular, similar a la vida, introducida por Gerard Vichniac y llamada regla de mayoría retorcida o regla de recocido, la regla de actualización establece el nuevo valor para que cada celda sea la mayoría entre las nueve celdas dadas por ella y sus ocho vecinas, excepto cuando estas celdas se dividen en cuatro con un estado y cinco con el otro estado, en cuyo caso el nuevo valor de la celda es la minoría en lugar de la mayoría. La dinámica detallada de esta regla es complicada, incluida la existencia de pequeñas estructuras estables. Sin embargo, en conjunto (cuando se inicia con todas las celdas en estados aleatorios) tiende a formar grandes regiones de celdas que están todas en el mismo estado entre sí, y los límites entre estas regiones evolucionan de acuerdo con el flujo de acortamiento de la curva.

Construcción de geodésicas cerradas

El flujo de acortamiento de la curva se puede utilizar para probar una desigualdad isoperimétrica para superficies cuya curvatura gaussiana es una función no creciente de la distancia desde el origen , como el paraboloide . En tal superficie, el conjunto compacto y liso que tiene un área determinada y un perímetro mínimo para esa área es necesariamente un círculo centrado en el origen. La prueba aplica el flujo de acortamiento de curvas a dos curvas, un círculo métrico y el límite de cualquier otro conjunto compacto, y compara el cambio en el perímetro de las dos curvas, ya que ambas se reducen a un punto por el flujo. El flujo de acortamiento de curvas también se puede utilizar para demostrar el teorema de las tres geodésicas , que cada variedad de Riemanniana suave topológicamente equivalente a una esfera tiene tres geodésicas que forman curvas cerradas simples .

Flujos relacionados

Otros flujos geométricos relacionados con el flujo de acortamiento de curvas incluyen los siguientes.

  • Para simular el comportamiento de cristales u otros materiales anisotrópicos , es importante tener variantes del flujo de acortamiento de curvas para las cuales la velocidad del flujo depende de la orientación de una curva así como de su curvatura. Una forma de hacer esto es definir la energía de una curva como la integral de una función suave γ de sus vectores normales, y formar el flujo gradiente de esta energía, según el cual la velocidad normal a la que fluye la curva es proporcional a un análogo anisotrópico de la curvatura. Este flujo se puede simular discretizando la curva como un polígono. En experimentos numéricos, las curvas iniciales parecen converger a la forma de Wulff para γ antes de reducirse a un punto. Alternativamente, se puede dejar que el flujo de curva con la velocidad de un ( θ ) κ + b ( θ ) , donde κ es el (usual) de curvatura y un y b son funciones suaves de la orientación θ . Cuando un ( θ + π ) = un ( θ ) y b ( θ + π ) = - b ( θ ) (de modo que el flujo es invariante bajo punto de reflexión ), el flujo resultante se puede demostrar que obedecen el principio de evitación y una análogo del teorema de Gage-Hamilton-Grayson.
  • El flujo de acortamiento de la curva afín fue investigado por primera vez por Alvarez et al. (1993) y Sapiro y Tannenbaum (1993) . En este flujo, la velocidad normal de la curva es proporcional a la raíz cúbica de la curvatura. El flujo resultante es invariante (con una escala de tiempo correspondiente) bajo las transformaciones afines del plano euclidiano, un grupo de simetría más grande que las transformaciones de similitud bajo las cuales el flujo de acortamiento de la curva es invariante. Bajo este flujo, se aplica un análogo del teorema de Gage-Hamilton-Grayson, bajo el cual cualquier curva cerrada simple eventualmente se vuelve convexa y luego converge en una elipse cuando colapsa en un punto.
  • La transformación de una curva con velocidades normales iguales en todos los puntos se denomina transformada grassfire . Las curvas evolucionadas de esta manera desarrollarán en general esquinas afiladas, cuyo trazo forma el eje medial de la curva. Una evolución de curva estrechamente relacionada que mueve segmentos rectos de una curva poligonal a velocidades iguales pero permite que las esquinas cóncavas se muevan más rápidamente que la velocidad unitaria, en cambio, forma un tipo diferente de esqueleto topológico de la curva dada, su esqueleto recto .
  • Para superficies de mayores dimensiones, hay más de una definición de curvatura, incluidas medidas extrínsecas (dependientes de la incrustación) como la curvatura media y medidas intrínsecas como la curvatura escalar y la curvatura de Ricci . En consecuencia, hay varias formas de definir los flujos geométricos basados ​​en la curvatura, incluido el flujo de curvatura media (en el que la velocidad normal de una superficie incrustada es su curvatura media), el flujo de Ricci (un flujo intrínseco en la métrica de un espacio basado en su curvatura de Ricci), el flujo de curvatura de Gauss y el flujo de Willmore (el flujo de gradiente para una energía funcional que combina la curvatura media y la curvatura de Gauss). El flujo de acortamiento de curvas es un caso especial del flujo de curvatura media y del flujo de curvatura de Gauss para curvas unidimensionales.
  • Inspirados por el flujo de acortamiento de curvas en curvas suaves, los investigadores han estudiado métodos para fluir polígonos para que permanezcan poligonales, con aplicaciones que incluyen la formación de patrones y la sincronización en sistemas distribuidos de robots. Se pueden usar flujos poligonales que preservan la longitud para resolver el problema de la regla de carpintero .
  • En la visión por computadora , el modelo de contorno activo para la detección de bordes y la segmentación de la imagen se basa en el acortamiento de la curva y evoluciona las curvas en función de una combinación de su curvatura y las características de una imagen.

Notas

Referencias