Diferencia entre revisiones de «Flujo de acortamiento de una curva»

De Wikipedia, la enciclopedia libre
Contenido eliminado Contenido añadido
Creación de «Contracción de una curva»
(Sin diferencias)

Revisión del 20:09 28 jun 2023

Plantilla:Short description Plantilla:Good article

Convergence of a convex curve to a circle under the curve-shortening flow. Inner curves (lighter color) are flowed versions of the outer curves. Time steps between curves are not uniform.

En matemáticas, el flujo de acortamiento de la curva es un proceso que modifica un curva en el Plano (geometría) moviendo sus puntos perpendicularmente a la curva a una velocidad proporcional al curvatura. El flujo de acortamiento de la curva es un ejemplo de geometric flow y es el caso unidimensional de flujo de curvatura promedio. Otros nombres para el mismo proceso incluyen flujo de acortamiento euclidiano, flujo de calor geométrico,[1]​ y evolución de la longitud del arco.

Como los puntos de cualquier teorema de la curva de Jordan suave se mueven de esta manera, la curva sigue siendo simple y suave. Pierde área a un ritmo constante, y su perímetro disminuye lo más rápido posible para cualquier evolución de curva continua. Si la curva no es convexa, su total absolute curvature decrece monótonamente, hasta volverse convexa. Una vez convexo, el isoperimetric ratio de la curva disminuye a medida que la curva converge en una forma circular, antes de colapsar en un solo punto de singularidad. Si dos curvas cerradas suaves 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 entre sí o tienen una longitud infinita mantienen su forma, incluida la curva de la parca, una curva infinita que se traslada hacia arriba, y espiral que giran. manteniendo 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 polígono y usando método de las diferencias finitas para calcular el movimiento de cada vértice del polígono. Los métodos alternativos incluyen calcular un convolución de vértices de polígonos y luego volver a muestrear los vértices en la curva resultante, o aplicar repetidamente un median filter a un imagen digital cuyos píxeles en blanco y negro representan el interior y el exterior de la curva.

El flujo de acortamiento de la curva se estudió originalmente como modelo para annealing de láminas de metal. Más tarde, 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 cellular automata. El flujo de acortamiento de curvas se puede usar para encontrar closed geodesic en Variedad de Riemann y como modelo para el comportamiento de flujos de dimensiones superiores.

Definiciones

Un flow es un proceso en el que los puntos de un espacio cambian continuamente sus ubicaciones o propiedades en el tiempo. Más específicamente, en un geometric flow unidimensional como el flujo de acortamiento de la curva, los puntos que experimentan el flujo pertenecen a un curva, y lo que cambia es la forma de la curva, su encaje (matemática) en el plano euclidiano determinado por las ubicaciones de cada uno de sus puntos[2]​ En el flujo de acortamiento de curvas, cada punto de una curva se mueve en la dirección de normal vector hacia la curva, a una velocidad proporcional a curvatura. Para una curva en evolución representada por una función de dos parámetros C(s,t) donde s parametriza la longitud del arco en la curva y t parametriza un tiempo en la evolución de la curva, el flujo de acortamiento de la curva puede ser descrito por el ecuación parabólica en derivadas parciales

una forma de ecuación del calor, donde κ es la curvatura y n es el vector unitario normal. [3]

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, equivariant) . Si el plano se escala por un factor de dilatación constante, el flujo permanece esencialmente sin cambios, pero se ralentiza o acelera por el mismo factor. [4]

Curvas irregulares

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 convierte en analytic, y permanece así hasta llegar a una singularidad en la que la curvatura explota. Para una curva suave sin cruces, la única singularidad posible ocurre cuando la curva colapsa en un punto, pero immersed curves puede tener otros tipos de singularidad.[5]​ En tales casos, con algo de cuidado es posible continuar el flujo más allá de estas singularidades hasta que toda la curva se reduzca a un solo punto.[6]

Para una curva cerrada simple, usando una extensión del flujo a curvas no suaves basadas en método del conjunto de nivel, solo hay dos posibilidades. Las curvas con cero Medida de Lebesgue (incluidos todos los polígono 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, los Curva de Osgood con una medida distinta de cero evolucionan inmediatamente a un annulus topológico con un área distinta de cero y límites uniformes. [7]​ El seno del topólogo es un ejemplo que instantáneamente se vuelve suave, a pesar de no ser siquiera locally connected; 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 finita de tiempo. [8]

Superficies no euclidianas

El flujo de acortamiento de curvas, y muchos de los resultados sobre el flujo de acortamiento de curvas, se pueden generalizar desde el plano euclidiano a cualquier Variedad de Riemann bidimensional. Para evitar tipos adicionales de singularidad, es importante que la variedad sea "convexa en el infinito"; esto es definitivamente


ined para significar que cada espacio compacto tiene un envolvente convexa compacto, como se define usando geodesic convexity. El flujo que acorta 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 de la variedad. [9]

Curvas espaciales

El flujo de acortamiento de curvas también se ha estudiado para curvas en Espacio euclídeo 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 Frenet–Serret frame. 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; toda singularidad es asintótica a un plano. [10]​ Sin embargo, se sabe que las curvas esféricas y las curvas que se pueden proyectar ortogonalmente en una curva plana convexa regular siguen siendo simples. [11]​ El flujo de acortamiento de curvas para curvas espaciales se ha utilizado como una forma de definir el flujo más allá de singularidades en curvas planas. [12]

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 rectifiable varifolds o método del conjunto de nivel. Sin embargo, estas definiciones extendidas pueden permitir que partes de las curvas desaparezcan instantáneamente o se conviertan en conjuntos de áreas distintas de cero.[13]

For networks of curves, extending the curve-shortening flow past a singularity may result in ambiguity or fattening.

Una variación del problema comúnmente estudiada implica redes de curvas suaves disjuntas en el interior, 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 2Plantilla:Pi/3 (las mismas condiciones que se observan en un Problema del árbol de Steiner óptimo o en un espuma bidimensional de pompa de jabón), el flujo está bien definido a corto plazo. Sin embargo, eventualmente puede llegar a 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 tal singularidad. [14]

Comportamiento

Principio de evitación, radio y factor de estiramiento

Si dos simple closed curves suaves disjuntos experimentan el flujo de acortamiento de la curva simultáneamente, permanecen disjuntos a medida que avanza el flujo. La razón es que, si dos curvas suaves se mueven de una manera que crean un cruce, entonces, en el momento del primer cruce, las curvas serían necesariamente 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 el principio de evitación. [15]

El principio de evasión implica que cualquier curva cerrada suave debe eventualmente alcanzar una singularidad, como un punto de curvatura infinita. Porque, si una curva suave dada C está rodeada por un círculo, ambas permanecerán disjuntas mientras existan. Pero el círculo envolvente se contrae bajo el flujo de curvatura, permaneciendo circular, hasta que colapsa, y por el principio de evasió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 de smallest circle that encloses 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. [15]

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 llamado 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 del diámetro de un círculo, en cuyo caso el factor de estiramiento es constante en Plantilla:Pi. Esta propiedad de monotonicidad implica el principio de evitación, ya que si la curva alguna vez se tocara a sí misma, el factor de estiramiento se volvería infinito en los dos puntos de contacto.

Longitud

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 la integral se toma sobre la curva, κ es la curvatura y s es la longitud del arco en 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. Más generalmente, para cualquier evolución de curvas cuya velocidad normal sea f, la tasa de cambio en la longitud es

que puede interpretarse como un espacio prehilbertiano negado entre la evolución dada y el flujo de acortamiento de la curva. Por lo tanto, el flujo que acorta la curva se puede describir como el campo vectorial de longitud, el flujo que (localmente) disminuye la longitud de la curva lo más rápido posible en relación con el L2 norm del flujo. Esta propiedad es la que le da su nombre al flujo de acortamiento de la curva.[16]

Área

Para

curva cerrada simple, el área encerrado por la curva se contrae, a una tasa constante de 2Plantilla:Pi unidades de área por unidad de tiempo, independientemente de la curva. Por lo tanto, el tiempo total para que una curva se reduzca a un punto es proporcional a su área, independientemente de su forma inicial.[17]​ Debido a que el área de una curva se reduce a una tasa constante, y (por isoperimetría) un círculo tiene el área más grande posible entre las curvas cerradas simples de una longitud dada, se deduce que los círculos son las curvas más lentas para colapsar hasta un punto bajo el flujo de acortamiento de la curva. Todas las demás curvas tardan menos en contraerse que un círculo de la misma longitud. [18]

La tasa constante de reducción del área es la única leyes de conservación (física) satisfecha por el flujo de acortamiento de la curva. Esto implica que no es posible expresar el "punto de fuga" donde la curva finalmente colapsa como una integral sobre la curva de cualquier función de sus puntos y sus derivados, porque tal expresión conduciría a una segunda ley de conservación prohibida. [19]​ Sin embargo, al combinar la tasa constante de pérdida de área con el principio de evitación, es posible demostrar que el punto de fuga siempre se encuentra dentro de un círculo, concéntrico con el círculo envolvente mínimo, cuya área es la diferencia de áreas entre el círculo envolvente y la curva dada. [20]

Curvatura absoluta total

El total absolute curvature de una curva suave es la integral del valor absoluto de la curvatura en la longitud del arco de la curva,

También se puede expresar como la suma de los ángulos entre los vectores normales en pares consecutivos de punto de inflexión. Es 2Plantilla:Pi para curvas convexas y mayor para curvas no convexas, sirviendo como medida de no convexidad de una curva. [21]

El flujo de acortamiento de la curva no puede crear nuevos puntos de inflexión.[22]​ Cada uno de los ángulos en la representación de la curvatura absoluta total como suma decrece monótonamente, excepto en los instantes en que dos puntos de inflexión consecutivos alcanzan el mismo ángulo o posición entre sí 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 2Plantilla:Pi y para curvas no convexas decrece monótonamente.[23]

Teorema de Gage-Hamilton-Grayson

Si una curva cerrada simple y suave se somete al flujo de acortamiento de la curva, permanece incrustada suavemente sin autointersecciones. Eventualmente se convertirá en convex, y una vez que lo haga, permanecerá 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 la curva completa se reduce a un solo punto. Este comportamiento a veces se resume diciendo que cada curva cerrada simple se reduce a un "punto redondo".[24]

Este resultado se debe a Michael Gage, Richard Hamilton (profesor) y Matthew Grayson. Plantilla:Harvs demostró la convergencia a un círculo para curvas convexas que se contraen en un punto. Más específicamente, Gage demostró que isoperimetric ratio (la relación entre la longitud de la curva al cuadrado y el área, un número que es 4Plantilla:Pi para un círculo y mayor para cualquier otra curva convexa) disminuye monótona y rápidamente. Gage y Hamilton (1986) demostró que todas las curvas convexas suaves finalmente se contraen hasta un punto sin formar ninguna otra singularidad, y Grayson (1987) demostró que todas las curvas no convexas eventualmente se volverán convexas.[25]Andrews y Bryan (2011) proporciona una prueba más simple del resultado de Grayson, basada en la monotonicidad del factor de estiramiento.

The limiting shape for all networks of two collinear rays and two curves connecting the endpoints of the two rays. The central lens has the shape of a vesica piscis.

Se pueden extender resultados similares de curvas cerradas a curvas ilimitadas que satisfacen un Lipschitz condition local. Para tales curvas, si ambos lados de la curva tienen un área infinita, entonces la curva evolucionada permanece uniforme y libre de singularidades todo el tiempo. 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. . [26]​ 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. [27]​ Para redes formadas por dos rayos disjuntos en la misma línea, junto con dos curvas suaves que conectan los extremos 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 vesica piscis.[28]

Singularidades de las curvas autocruzadas

Las curvas que tienen autocruces pueden alcanzar singularidades antes de contraerse en un punto. Por ejemplo, si un lemniscata (cualquier immersed curve suave con un solo cruce, parecido a una figura 8 o símbolo de infinito) tiene áreas desiguales en sus dos lóbulos, 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 en la evolución de la curva, y la relación isoperimétrica divergirá a medida que la curva colapsa en una singularidad. [4]

Cuando una curva localmente convexa que se cruza a sí misma se aproxima a una singularidad a medida que uno de sus bucles se contrae, se contrae en una autosimilar

se acerca o asintóticamente a la curva de la parca (descrita a continuación) a medida que se contrae. Cuando un bucle colapsa en una singularidad, la cantidad de curvatura absoluta total que se pierde es al menos 2Plantilla:Pi o exactamente Plantilla:Pi. [29]

Sobre variedades de Riemann

En una variedad de Riemann, cualquier curva cerrada simple suave seguirá siendo suave y simple a medida que evoluciona, al igual que en el caso euclidiano. O se derrumbará hasta un punto en un tiempo finito, o permanecerá suave y simple para siempre. En este último caso, la curva necesariamente converge a un closed geodesic de la superficie.[30]

Las curvas sumergidas en las variedades de Riemann, con un número finito de autocruces, se vuelven autotangentes solo en un conjunto discreto de tiempos, en cada uno de los cuales pierden un cruce. Como consecuencia, el número de puntos de autocruce no aumenta. [31]

A pelota de tenis

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

Fórmula de monotonicidad de Huisken

De acuerdo con Fórmula de monotonicidad de Huisken, la convolución de una curva en evolución con un kernel de calor invertido en el tiempo no es creciente. Este resultado se puede utilizar para analizar las singularidades de la evolución. [33]

Curvas específicas

Curvas con evolución autosimilar

The grim reaper curve and translated copies of it produced by the curve-shortening flow

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 cruces automáticos) o no están cerradas (se extienden hasta el infinito) y mantienen su forma. En particular,[34]

  • Cada line 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 la curva,[34]​, aunque existen redes de curvas estables más complejas, como la teselado hexagonal del avión.
  • El grim reaper curve y= − log cos x se mueve hacia arriba sin cambiar su forma. De la misma manera, cualquier curva similar a la parca es translated por el flujo de acortamiento de la curva, desplazada en la dirección de la symmetry axis de la curva sin cambiar su forma u orientación. La parca es la única curva con esta propiedad.[34]​ También se le llama modelo de horquilla en la literatura de física.[35]
  • Una familia de curvas cerradas autocruzadas, derivadas de proyecciones de torus knot, encogen homothetically pero siguen siendo autosimilares bajo el flujo de acortamiento de la curva.[34]​ Estos han llegado a ser conocidos como Abresch–Langer curves, después del trabajo de Abresch y Langer (1986),[36]​, aunque Mullins (1956) los mencionó anteriormente y Epstein y Weinstein (1987) los redescubrió de forma independiente. Estas curvas son localmente convexas y, por lo tanto, pueden describirse mediante sus support function. Las versiones adecuadamente escaladas de estas funciones de soporte obedecen a ecuación diferencial
que tiene soluciones periódicas positivas (correspondientes a curvas con evolución autosimilar) para cualquier periodo que se encuentre estrictamente entre Plantilla:Pi y . [36]
  • Otras curvas, incluidas algunas espiral infinitas, siguen siendo similares a sí mismas con movimientos más complicados que incluyen rotación o combinaciones de rotación, contracción o expansión y traslación.[34]
  • Para redes de curvas suaves, reunidas de a tres en uniones con ángulos de 2Plantilla:Pi/3, las soluciones de contracción autosimilares incluyen un double bubble que rodea dos áreas iguales, una forma lens (vesica piscis) delimitada por dos arcos de círculos congruentes junto con dos colineales rayos 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 autosimilar implica un mayor número de curvas. [37]​ Otra familia de redes crece homotéticamente y sigue siendo autosimilar; estas son redes de curvas en forma de árbol, que se encuentran en ángulos de 2Plantilla:Pi/3 en uniones triples, asymptotic a un abanico de dos o más rays 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 se puede definir utilizando variantes generalizadas del flujo de acortamiento de curvas, como el de los pliegues múltiples. 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 la curva a partir de un abanico de rayos.[38]

Soluciones antiguas

Un ancient solution para un problema de flujo es una curva cuya evolución puede extrapolarse hacia atrás para todos los tiempos, sin singularidades. Todas las soluciones autosimilares que se encogen o mantienen el mismo tamaño en lugar de crecer son soluciones antiguas en este sentido; se pueden extrapolar hacia atrás invirtiendo la transformación autosimilitud que t

sufrirían por el flujo de acortamiento de la curva hacia adelante. Así, por ejemplo, el círculo, la parca y las curvas de Abresch-Langer son todas soluciones antiguas. [39]

También hay ejemplos que no son autosimilares. Un ejemplo explícito es la solución Angenent oval 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 límite bajo la evolución inversa un par de curvas de la muerte que se acercan entre sí desde direcciones opuestas. [40]​ En el Coordenadas cartesianas, pueden estar dadas por la ecuación implicit curve[41]

En la literatura de física, las mismas formas se conocen como "modelo de sujetapapeles".[35]

Las soluciones del óvalo de Angenent y del círculo que se contrae son las únicas soluciones antiguas cuyos intervalos de tiempo unían conjuntos convexos acotados. [39]​ The Grim Reaper, el medio espacio estacionario y las soluciones de franjas estacionarias son los únicos ejemplos cuyos intervalos de tiempo unen conjuntos convexos ilimitados. [42]​ Existen muchos ejemplos convexos "localmente" adicionales (no incrustados), así como muchos ejemplos incrustados (no convexos) adicionales. [43][44]

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 Front tracking se han utilizado durante mucho tiempo en fluidodinámica para modelar y rastrear el movimiento de los límites entre diferentes materiales, de gradientes pronunciados en las propiedades del material, como frente (meteorología), o de 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.[45]​ Los mismos métodos también se pueden usar para simular el flujo de acortamiento de la curva, incluso cuando la curva que atraviesa 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 las 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.[46]​ Aunque el flujo de acortamiento de la curva se define por el movimiento de una curva perpendicular 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 en 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 en la curva en situaciones en las que el movimiento perpendicular haría que se agruparan.[47]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 simples curvas cerradas, 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 paso de tiempo debe elegirse ad hoc".[48]​ Otro método de diferenciación finita por Crandall y Lions (1996) modifica la fórmula para la curvatura en cada vértice al agregarle un término pequeño basado en Operador laplaciano. Esta modificación se llama elliptic regularization y puede usarse para ayudar a probar la existencia de flujos generalizados, así como en su simulación numérica.[49]​ Usá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.[50]​ Para obtener una comparación empírica de los métodos de diferencias finitas forward Euler, backward Euler y Crank–Nicolson más precisos, consulte Balažovjech y Mikula (2009).

Convolución remuestreada

Mokhtarian y Mackworth (1992) sugiere 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 en un espacio uniforme, medido por la longitud de arco normalizada.
  • Convolve las ubicaciones de los puntos con un Función gaussiana con desviación estándar pequeña, en efecto reemplazando la ubicación de cada punto con un weighted average de las ubicaciones de los puntos cercanos en la curva, con pesos gaussianos. La desviación estándar de la Gaussiana debe elegirse para que sea lo suficientemente pequeña como para que, después de este paso, los puntos de muestra aún tengan un espacio casi uniforme.

Como muestran, este método converge a la distribución de acortamiento de la curva en el límite a medida que crece el número de puntos de muestra y se reduce la longitud de arco normalizada del radio de convolución.[51]

Filtrado mediano

Merriman, Bence y Osher (1992) describe un esquema que opera en una cuadrícula cuadrada bidimensional, efectivamente una matriz de píxel. 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 el función indicatriz para el interior de la curva. Esta representación se actualiza alternando dos pasos:

  • Convolucionar la imagen pixelada con un kernel de calor para simular su evolución bajo el ecuación del calor por un paso de tiempo corto. El resultado es un Desenfoque gaussiano de la imagen, o equivalentemente el Weierstrass transform de la función indicadora, con radio proporcional a la raíz cuadrada del paso de tiempo.
  • Establezca cada píxel con valor numérico inferior a 1/2 a 0, y cada píxel con valor numérico superior a 1/2 a 1, thresholding la imagen vuelva a sus valores originales en nuevas posiciones.

Para que este esquema sea preciso, el paso de tiempo debe ser lo suficientemente grande como para hacer que la curva se mueva al menos un píxel incluso en puntos de baja curvatura, pero lo suficientemente pequeño como para hacer 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 elegir un paso de tiempo intermedio adecuado.

El método se puede generalizar a la evolución de redes de curvas que se encuentran en los cruces y dividen el plano en más de tres regiones, aplicando el mismo método simultáneamente a cada región. [52]​ En lugar de desenfoque y umbralización, este método puede describirse alternativamente como aplicar un median filter con Gaussian weights a cada píxel. Es posible usar kernels distintos al kernel de calor, o refinar adaptativamente la cuadrícula para que tenga alta resolución cerca de la curva pero no pierda tiempo y memoria en píxeles lejos de la curva que no contribuyen al resultado.[53]​ 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íxel representan la distancia firmada a la curva puede lograr una precisión de subpíxel y requiere una resolución más baja. [54]

Aplicaciones

Recocido de chapas

Una referencia temprana al flujo de acortamiento de la curva por Plantilla:Harvs lo motiva como un modelo para el proceso físico de annealing, en el que el tratamiento térmico hace que los límites entre los granos de metal cristalizado se desplacen. A diferencia de los soap film, que se ven obligados por las diferencias en presión atmosférica a convertirse en superficies de curvatura media constante, los límites de grano en el recocido están sujetos solo a efectos locales, lo que hace que se muevan de acuerdo con el flujo de curvatura media. El caso unidimensional de este flujo, el flujo de acortamiento de curvas, corresponde al recocido de láminas de metal que son lo suficientemente delgadas para que los granos se vuelvan efectivamente bidimensionales y sus límites se vuelvan unidimensionales.[55]

Análisis de forma

En procesamiento digital de imágenes y visión artificial, Mokhtarian y Mackworth (1992) sugiere 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 scale space que brinde 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 punto de inflexión de la curva a medida que avanzan a través del flujo y dibujar un gráfico que represente las posiciones de los puntos de inflexión alrededor de la curva contra el parámetro de tiempo. Los puntos de inflexión generalmente se eliminarán de la curva en pares como el

La curva se vuelve convexa (según el teorema de Gage-Hamilton-Grayson) y el tiempo de vida de un par de puntos corresponde a la prominencia de una característica de la forma. Debido al método de convolución remuestreado que describen para calcular una aproximación numérica del flujo de acortamiento de la curva, llaman a su método "espacio de escala de curvatura remuestreado". Observan que este espacio de escala es invariable bajo las transformaciones euclidianas de la forma dada y afirman que determina de forma única la forma y es robusto frente a pequeñas variaciones en la forma. Lo comparan experimentalmente con varias definiciones alternativas relacionadas de un espacio de escala para formas, y descubren que el espacio de escala de curvatura remuestreado es menos computacionalmente intensivo, más robusto contra el 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 Ecuación de Allen-Cahn, el comportamiento límite para reacción rápida, difusión lenta y dos o más mínimos locales de energía con el mismo nivel de energía entre sí es que el sistema se asiente en regiones de diferentes mínimos locales, con los frentes delimitando los límites entre estas regiones evolucionan de acuerdo con el flujo de acortamiento de la curva. [56]

Autómatas celulares

The Anneal cellular automaton, 1600 steps after a random start

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 Life-like cellular automaton es aquella en la que la cuadrícula es la retícula cuadrada infinita, hay exactamente dos estados de celda, el conjunto de vecinos de cada celda son los ocho vecinos de Vecindad de Moore, y la regla de actualización depende solo del número de vecinos con cada uno de los dos estados en lugar de cualquier función más complicada de esos estados. En una regla similar a la vida en particular, presentada por Gerard Vichniac y llamada regla de la mayoría torcida o regla de recocido, la regla de actualización establece que el nuevo valor para cada celda sea la mayoría entre las nueve celdas dadas por ella y sus ocho vecinas, excepto cuando estas celdas se dividen entre 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. [57]​ 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, con los límites entre estas regiones evolucionando de acuerdo con el flujo de acortamiento de la curva.[58]

Construcción de geodésicas cerradas

El flujo de acortamiento de la curva se puede usar para probar un isoperimetría para superficies cuyo Curvatura de Gauss es una función no creciente de la distancia desde el origin, como el paraboloide. En tal superficie, el conjunto compacto liso que tiene un área determinada y un perímetro mínimo para esa área es necesariamente un círculo con centro 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 cuando el flujo las reduce a un punto.[59]​ El flujo de acortamiento de la curva también se puede usar para probar el theorem of the three geodesics, que cada variedad riemanniana suave topológicamente equivalente a una esfera tiene tres geodésicas que forman teorema de la curva de Jordan. [60]

Flujos relacionados

Otros geometric flow relacionados con el flujo de acortamiento de la curva incluyen los siguientes.

  • Para simular el comportamiento de cristal u otros materiales anisotropic, es importante contar con 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 hacerlo es definir la energía de una curva como la integral de un función infinitamente diferenciable γ 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 en Wulff shape para γ antes de reducirse a un punto.[61]​ Alternativamente, se puede dejar que la curva fluya con la velocidad a(θ)κ + b(θ), donde κ es la curvatura (habitual) y a y b son funciones uniformes de la orientación θ. Cuando a(θ + Plantilla:Pi)= a(θ) y b(θ + Plantilla:Pi)= −b(θ) (de modo que el flujo es invariable bajo simetría central), se puede demostrar que el flujo resultante obedece al principio de evitación y un análogo del teorema de Gage-Hamilton-Grayson.[62]
  • El affine curve-shortening flow 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.[63]​ El flujo resultante es invariante (con una escala de tiempo correspondiente) bajo los transformación afín del plano euclidiano, un grupo de simetría más grande que el similarity transformations bajo el cual 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 a

un elipse a medida que colapsa en un punto.[64]

  • La transformación de una curva con velocidades normales iguales en todos los puntos se ha denominado grassfire transform. Las curvas desarrolladas de esta manera en general desarrollarán esquinas agudas, cuya traza forma el eje medio de la curva. [65]​ Una evolución de curva estrechamente relacionada que mueve segmentos rectos de una curva poligonal a la misma velocidad pero permite que las esquinas cóncavas se muevan más rápido que la velocidad unitaria en su lugar forma un tipo diferente de cálculo del esqueleto de la curva dada, su straight skeleton. [66]
  • Para superficies en dimensiones más altas, hay más de una definición de curvatura, incluidas medidas extrínsecas (dependientes de la incrustación) como curvatura media y medidas intrínsecas como curvatura escalar de Ricci y Tensor de Ricci. En consecuencia, hay varias formas de definir flujos geométricos basados en la curvatura, incluido el flujo de curvatura promedio (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 Gauss curvature flow y el Energía de Willmore (el flujo de gradiente para un funcional de energía que combina la curvatura media y la curvatura gaussiana). 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. [18]
  • En real-time path planning para robot móvil, se ha utilizado una versión modificada del flujo de acortamiento de curvas con fuerzas adicionales para encontrar caminos que logren un equilibrio entre ser cortos y mantenerse libres de obstáculos. [67]
  • Inspirándose en el flujo de acortamiento de curvas en curvas suaves, los investigadores han estudiado métodos para hacer fluir polígono de modo que permanezcan poligonales, con aplicaciones que incluyen la formación de patrones y la sincronización en sistemas distribuidos de robots. [68]​ Los flujos poligonales que conservan la longitud se pueden utilizar para resolver el carpenter's rule problem. [69]
  • En visión artificial, active contour model para detector de bordes y segmentación (procesamiento de imágenes) se basa en el acortamiento de curvas y desarrolla curvas en función de una combinación de su curvatura y las características de una imagen. [70]

Notas

  1. The phrase "geometric heat flow" has also been used for flows on other kinds of object than curves, such as forma diferencials.
  2. Devadoss y O'Rourke (2011), p.140: "a geometric flow [is] an evolution of the geometry of C over time t."
  3. Devadoss y O'Rourke, 2011.
  4. a b Grayson, 1989a.
  5. Grayson (1989a); White (2002).
  6. Angenent (1991a); Altschuler y Grayson (1992).
  7. Lauer, 2013.
  8. Lam y Lauer, 2016.
  9. Ritoré y Sinestrari, 2010, p. 72.
  10. Altschuler, 1991.
  11. Minarčík y Beneš, 2020.
  12. Altschuler y Grayson, 1992.
  13. Brakke (1978); White (1989); Cao (2003), "4.7.1 Brakke's varifold solution", p. 100. Lauer (2013).
  14. Ilmanen, Neves y Schulze, 2014.
  15. a b White, 2002.
  16. Chou y Zhu (2001), p.  vii; White (2002), p. 526.
  17. Brakke (1978), Appendix B, Proposition 1, p. 230; Chou y Zhu (2001), p.  vii; White (2002), Theorem 1, p. 527.
  18. a b White, 1989.
  19. Bryant y Griffiths, 1995.
  20. Kimmel, 2004.
  21. Brook, Bruckstein y Kimmel, 2005.
  22. Cao (2003), p. 143.
  23. Brakke (1978), Appendix B, Proposition 2, p. 230; Chou y Zhu (2001), Lemma 5.5, p. 130; "6.1 The decrease in total absolute curvature", pp. 144–147.
  24. Chou y Zhu (2001), p.  vii; White (2002), Theorems 2 and 3, pp. 527–528; Cao (2003), Theorem 3.26, p. 47; Devadoss y O'Rourke (2011), p. 141.
  25. Chou y Zhu (2001), p.  vii; Cao (2003), p. 47; Devadoss y O'Rourke (2011), p. 141.
  26. Chou y Zhu, 1998.
  27. Ishimura, 1995.
  28. Schnürer et al. (2011); Bellettini y Novaga (2011).
  29. Angenent, 1991b.
  30. Grayson (1989b); White (2002), p. 528; Ritoré y Sinestrari (2010), Theorem 2.2.1, p. 73. This result was already stated as a conjetura by Gage y Hamilton (1986).
  31. Angenent, 1991a.
  32. Angenent, 1999.
  33. Huisken, 1990.
  34. a b c d e Mullins (1956); Abresch y Langer (1986); Epstein y Weinstein (1987); Chou y Zhu (2001), "2. Invariant solutions for the curve-shortening flow", pp. 27–44; Halldórsson (2012); Altschuler et al. (2013).
  35. a b Lukyanov, Vitchev y Zamolodchikov (2004); Huisken y Sinestrari (2015).
  36. a b Au, 2010.
  37. Schnürer et al., 2011.
  38. The two-ray case was already described by Mullins (1956). For the generalization to two or more rays and issues of non-uniqueness see Brakke (1978), Appendix C, pp. 235–237 and Ilmanen, Neves y Schulze (2014).
  39. a b Daskalopoulos, Hamilton y Sesum, 2010.
  40. Angenent, 1992.
  41. Broadbridge y Vassiliou, 2011.
  42. Bourni, Langford y Tinaglia, 2020.
  43. Angenent y You, 2021.
  44. You, 2014.
  45. See, e.g., Scriven (1960); Holden y Risebro (2015).
  46. Merriman, Bence y Osher (1992); Mikula y Ševčovič (1999); Cao (2003), "5.1.1 Finite difference methods", pp. 107–108.
  47. Kimura (1994); Deckelnick y Dziuk (1995); Mikula y Ševčovič (2001); Barrett, Garcke y Nürnberg (2011); Elliott y Fritz (2017).
  48. Cao (2003), "5.1.1 Finite difference methods", pp. 107–108.
  49. Ilmanen (1994), p. 1.
  50. Crandall y Lions (1996); Deckelnick (2000); Cao (2003), "5.2.3 A monotone and convergent finite difference schemes", p. 109.
  51. Mokhtarian y Mackworth (1992), pp. 796–797; Cao (2003), pp. 10–11.
  52. Merriman, Bence y Osher, 1992.
  53. Cao (2003), "5.2.4 Bence, Merriman and Osher scheme for mean curvature motion", pp. 109–110. For the correctness of median filtering with other isotropic kernels, see section 4.4.1, pp. 90–92.
  54. Esedoḡlu, Ruuth y Tsai, 2010.
  55. Mullins (1956); Rhines, Craig y DeHoff (1974); Brakke (1978), Appendix A, pp. 224–228.
  56. Rubinstein, Sternberg y Keller, 1989.
  57. Pickover, 1993.
  58. Vichniac (1986); Chopard y Droz (1998).
  59. Benjamini y Cao (1996); Ritoré y Sinestrari (2010), Theorem 2.3.1, p. 75.
  60. Grayson, 1989b.
  61. Dziuk (1999); Haußer y Voigt (2006).
  62. Chou y Zhu (2001), Chapter 6: A Class of Non-convex Anisotropic Flows, pp. 143–177.
  63. Cao (2003), "3.2.3 The affine invariant flow: the simplest affine invariant curve flow", pp. 42–46.
  64. Angenent, Sapiro y Tannenbaum (1998); Cao (2003), Theorem 3.28, p. 47.
  65. Sapiro y Tannenbaum, 1993.
  66. Aichholzer et al., 1995.
  67. Huptych y Röck, 2021.
  68. Smith, Broucke y Francis, 2007.
  69. Cantarella et al., 2004.
  70. Kichenassamy et al., 1995.

Referencias