stringtranslate.com

Modelo de campo de fase

Un modelo de campo de fases es un modelo matemático para resolver problemas interfaciales. Se ha aplicado principalmente a la dinámica de solidificación, [1] pero también se ha aplicado a otras situaciones como la digitación viscosa , [2] la mecánica de fracturas , [ 3] [4] [5] [6] la fragilización por hidrógeno , [7] y la dinámica de vesículas . [8] [9] [10] [11]

El método sustituye las condiciones de contorno en la interfaz por una ecuación diferencial parcial para la evolución de un campo auxiliar (el campo de fase) que asume el papel de un parámetro de orden . Este campo de fase toma dos valores distintos (por ejemplo, +1 y −1) en cada una de las fases, con un cambio suave entre ambos valores en la zona alrededor de la interfaz, que entonces es difusa con un ancho finito. Una ubicación discreta de la interfaz puede definirse como la colección de todos los puntos donde el campo de fase toma un cierto valor (por ejemplo, 0).

Un modelo de campo de fases se construye generalmente de tal manera que en el límite de una anchura de interfaz infinitesimal (el llamado límite de interfaz agudo) se recupera la dinámica interfacial correcta. Este enfoque permite resolver el problema mediante la integración de un conjunto de ecuaciones diferenciales parciales para todo el sistema, evitando así el tratamiento explícito de las condiciones de contorno en la interfaz.

Los modelos de campo de fase fueron introducidos por primera vez por Fix [12] y Langer, [13] y han experimentado un creciente interés en la solidificación y otras áreas.

Ecuaciones del modelo de campo de fase

Los modelos de campo de fases se construyen generalmente para reproducir una dinámica de interfaz dada. Por ejemplo, en problemas de solidificación, la dinámica del frente está dada por una ecuación de difusión para la concentración o la temperatura en la masa y algunas condiciones de contorno en la interfaz (una condición de equilibrio local y una ley de conservación), [14] que constituye el modelo de interfaz nítido.

Se muestra una microestructura de dos fases y el perfil de parámetros de orden en una línea que atraviesa el dominio. El cambio gradual de los parámetros de orden de una fase a otra muestra la naturaleza difusa de la interfaz.

Varias formulaciones del modelo de campo de fase se basan en una función de energía libre que depende de un parámetro de orden (el campo de fase) y de un campo difusivo (formulaciones variacionales). Las ecuaciones del modelo se obtienen entonces utilizando relaciones generales de física estadística . Dicha función se construye a partir de consideraciones físicas, pero contiene un parámetro o una combinación de parámetros relacionados con el ancho de la interfaz. Los parámetros del modelo se eligen entonces estudiando el límite del modelo con este ancho que tiende a cero, de tal manera que se pueda identificar este límite con el modelo de interfaz nítido deseado.

Otras formulaciones parten de la redacción directa de las ecuaciones de campo de fases, sin hacer referencia a ningún funcional termodinámico (formulaciones no variacionales). En este caso la única referencia es el modelo de interfase nítida, en el sentido de que debe recuperarse al realizar el límite de ancho de interfase pequeño del modelo de campo de fases.

Las ecuaciones de campo de fase reproducen en principio la dinámica interfacial cuando el ancho de la interfaz es pequeño en comparación con la escala de longitud más pequeña del problema. En la solidificación, esta escala es la longitud capilar , que es una escala microscópica. Desde un punto de vista computacional, la integración de ecuaciones diferenciales parciales que resuelven una escala tan pequeña es prohibitiva. Sin embargo, Karma y Rappel introdujeron el límite de interfaz delgada, [15] que permitió relajar esta condición y abrió el camino a simulaciones cuantitativas prácticas con modelos de campo de fase. Con el aumento de la potencia de las computadoras y el progreso teórico en el modelado de campo de fase, los modelos de campo de fase se han convertido en una herramienta útil para la simulación numérica de problemas interfaciales.

Formulaciones variacionales

Se puede construir un modelo para un campo de fases mediante argumentos físicos si se tiene una expresión explícita para la energía libre del sistema. Un ejemplo simple para problemas de solidificación es el siguiente:

donde es el campo de fases, , es la entalpía local por unidad de volumen, es una cierta función polinómica de , y (donde es el calor latente , es la temperatura de fusión y es el calor específico). El término con corresponde a la energía interfacial. La función se toma generalmente como un potencial de doble pozo que describe la densidad de energía libre de la masa de cada fase, que a su vez corresponden a los dos mínimos de la función . Las constantes y tienen dimensiones respectivamente de energía por unidad de longitud y energía por unidad de volumen. El ancho de la interfaz se da entonces por . El modelo de campo de fases se puede obtener entonces a partir de las siguientes relaciones variacionales: [16]

donde D es un coeficiente de difusión para la variable , y y son términos estocásticos que dan cuenta de las fluctuaciones térmicas (y cuyas propiedades estadísticas se pueden obtener a partir del teorema de disipación de fluctuaciones ). La primera ecuación da una ecuación para la evolución del campo de fases, mientras que la segunda es una ecuación de difusión, que normalmente se reescribe para la temperatura o para la concentración (en el caso de una aleación). Estas ecuaciones son, escalando el espacio con y los tiempos con :

donde es el ancho de la interfaz adimensional, , y , son ruidos no dimensionalizados.

Funciones de densidad de energía alternativas

La elección de la función de energía libre, , puede tener un efecto significativo en el comportamiento físico de la interfaz, y debe seleccionarse con cuidado. La función de doble pozo representa una aproximación de la ecuación de estado de Van der Waals cerca del punto crítico, y se ha utilizado históricamente por su simplicidad de implementación cuando el modelo de campo de fase se emplea únicamente para fines de seguimiento de la interfaz. Pero esto ha llevado al fenómeno de contracción espontánea de la gota observado con frecuencia, por el cual la alta miscibilidad de fase predicha por una ecuación de estado cerca del punto crítico permite una interpenetración significativa de las fases y eventualmente puede llevar a la desaparición completa de una gota cuyo radio está por debajo de algún valor crítico. [17] Minimizar las pérdidas de continuidad percibidas durante la duración de una simulación requiere límites en el parámetro de movilidad, lo que resulta en un delicado equilibrio entre el borrón interfacial debido a la convección, la reconstrucción interfacial debido a la minimización de la energía libre (es decir, la difusión basada en la movilidad) y la interpenetración de fases, también dependiente de la movilidad. Una revisión reciente de funciones de densidad de energía alternativas para aplicaciones de seguimiento de interfaz ha propuesto una forma modificada de la función de doble obstáculo que evita los fenómenos de contracción espontánea de la gota y los límites de movilidad [18] , con resultados comparativos que proporcionan una serie de simulaciones de referencia utilizando la función de doble pozo y la técnica de interfaz nítida de volumen de fluido . La implementación propuesta tiene una complejidad computacional solo ligeramente mayor que la de la función de doble pozo, y puede resultar útil para aplicaciones de seguimiento de interfaz del modelo de campo de fase donde la duración/naturaleza de los fenómenos simulados introduce preocupaciones de continuidad de fase (es decir, gotas pequeñas, simulaciones extendidas, interfaces múltiples, etc.).

Límite de interfaz agudo de las ecuaciones de fase-campo

Se puede construir un modelo de campo de fase para reproducir deliberadamente una dinámica interfacial dada tal como se representa mediante un modelo de interfaz nítido. En tal caso, se debe realizar el límite de interfaz nítido (es decir, el límite cuando el ancho de la interfaz tiende a cero) del conjunto propuesto de ecuaciones de campo de fase. Este límite se toma generalmente mediante expansiones asintóticas de los campos del modelo en potencias del ancho de la interfaz . Estas expansiones se realizan tanto en la región interfacial (expansión interna) como en el volumen (expansión externa), y luego se combinan asintóticamente orden por orden. El resultado proporciona una ecuación diferencial parcial para el campo difusivo y una serie de condiciones de contorno en la interfaz, que deben corresponder al modelo de interfaz nítido y cuya comparación con él proporciona los valores de los parámetros del modelo de campo de fase.

Mientras que en los primeros modelos de campo de fases se realizaban estas expansiones hasta el orden inferior , los modelos más recientes utilizan asintóticas de orden superior (límites de interfaz delgados) para cancelar efectos espurios no deseados o para incluir nueva física en el modelo. Por ejemplo, esta técnica ha permitido cancelar efectos cinéticos, [15] tratar casos con difusividades desiguales en las fases, [19] modelar digitación viscosa [2] y flujos de Navier-Stokes bifásicos, [20] incluir fluctuaciones en el modelo, [21] etc.

Modelos de campos multifásicos

Los parámetros de orden múltiple describen una microestructura de material policristalino.

En los modelos de campo multifásico, la microestructura se describe mediante un conjunto de parámetros de orden, cada uno de los cuales está relacionado con una fase específica u orientación cristalográfica. Este modelo se utiliza principalmente para transformaciones de fase de estado sólido donde evolucionan múltiples granos (por ejemplo, crecimiento de grano , recristalización o transformación de primer orden como austenita a ferrita en aleaciones ferrosas). Además de permitir la descripción de múltiples granos en una microestructura, los modelos de campo multifásico permiten especialmente la consideración de múltiples fases termodinámicas que ocurren, por ejemplo, en grados de aleaciones técnicas. [22]

Modelos de campo de fase en gráficos

Muchos de los resultados de los modelos de campo de fase continuo tienen análogos discretos para los gráficos, simplemente reemplazando el cálculo con cálculo sobre gráficos .

Modelado de campos de fases en mecánica de fracturas

Las fracturas en sólidos se analizan a menudo numéricamente dentro de un contexto de elementos finitos utilizando representaciones de grietas discretas o difusas. Los enfoques que utilizan una representación de elementos finitos a menudo hacen uso de fuertes discontinuidades incrustadas en el nivel intraelemento y a menudo requieren criterios adicionales basados ​​en, por ejemplo, tensiones, densidades de energía de deformación o tasas de liberación de energía u otros tratamientos especiales como técnicas de cierre de grietas virtuales y remallado para determinar las trayectorias de las grietas. Por el contrario, los enfoques que utilizan una representación de grietas difusas conservan la continuidad del campo de desplazamiento, como los modelos de daño continuo y las teorías de fractura de campo de fase. Este último se remonta a la reformulación del principio de Griffith en una forma variacional y tiene similitudes con los modelos de tipo de daño mejorado por gradiente. Quizás la característica más atractiva de los enfoques de campo de fase para la fractura es que la iniciación de grietas y las trayectorias de grietas se obtienen automáticamente a partir de un problema de minimización que acopla las energías elásticas y de fractura. En muchas situaciones, la nucleación de grietas se puede explicar adecuadamente siguiendo las ramas de los puntos críticos asociados con las soluciones elásticas hasta que pierden estabilidad. En particular, los modelos de fractura de campo de fase pueden permitir la nucleación incluso cuando la densidad de energía de deformación elástica es constante en el espacio. [23] Una limitación de este enfoque es que la nucleación se basa en la densidad de energía de deformación y no en la tensión. Una visión alternativa basada en la introducción de una fuerza impulsora de la nucleación busca abordar esta cuestión. [24]

Modelos de campo de fases para la migración celular colectiva

Un grupo de células biológicas puede autopropulsarse de una manera compleja debido al consumo de trifosfato de adenosina . Las interacciones entre células como la cohesión o varias señales químicas pueden producir movimiento de manera coordinada, este fenómeno se llama "migración celular colectiva". Un modelo teórico para estos fenómenos es el modelo de campo de fase [25] [26] e incorpora un campo de fase para cada especie celular y variables de campo adicionales como la concentración de agente quimiotáctico . Tal modelo puede usarse para fenómenos como el cáncer, la cicatrización de heridas, la morfogénesis y los fenómenos de ectoplasma .

Software

Referencias

  1. ^ Boettinger, WJ; Warren, JA; Beckermann, C .; Karma, A. (2002). "Simulación de solidificación por campo de fase". Revisión anual de investigación de materiales . 32 : 163–194. doi :10.1146/annurev.matsci.32.101901.155803.
  2. ^ ab Folch, R.; Casademunt, J.; Hernández-Machado, A.; Ramírez-Piscina, L. (1999). "Modelo de campo de fases para flujos de Hele-Shaw con contraste de viscosidad arbitrario. II. Estudio numérico". Physical Review E . 60 (2): 1734–40. arXiv : cond-mat/9903173 . Código Bibliográfico :1999PhRvE..60.1734F. doi :10.1103/PhysRevE.60.1734. PMID  11969955. S2CID  8488585.
  3. ^ ab Bourdin, B.; Francfort, GA; Marigo, JJ. (abril de 2000). "Experimentos numéricos en fractura frágil revisada". Revista de mecánica y física de sólidos . 48 (4): 797–826. Bibcode :2000JMPSo..48..797B. doi :10.1016/S0022-5096(99)00028-9.
  4. ^ ab Bourdin, Blaise (2007). "Implementación numérica de la formulación variacional para fractura frágil cuasiestática". Interfaces y límites libres . 9 (3): 411–430. doi : 10.4171/IFB/171 . ISSN  1463-9963.
  5. ^ ab Bourdin, Blaise; Francfort, Gilles A.; Marigo, Jean-Jacques (abril de 2008). "El enfoque variacional de la fractura". Journal of Elasticity . 91 (1–3): 5–148. doi :10.1007/s10659-007-9107-3. ISSN  0374-3535. S2CID  120498253.
  6. ^ Karma, Alain; Kessler, David; Levine, Herbert (2001). "Modelo de campo de fase de fractura dinámica de modo III". Physical Review Letters . 87 (4): 045501. arXiv : cond-mat/0105034 . Código Bibliográfico :2001PhRvL..87d5501K. doi :10.1103/PhysRevLett.87.045501. PMID  11461627. S2CID  42931658.
  7. ^ Martinez-Paneda, Emilio; Golahmar, Alireza; Niordson, Christian (2018). "Una formulación de campo de fase para el craqueo asistido por hidrógeno". Métodos informáticos en mecánica aplicada e ingeniería . 342 : 742–761. arXiv : 1808.03264 . Código Bibliográfico : 2018CMAME.342..742M. doi : 10.1016/j.cma.2018.07.021. S2CID  : 52360579.
  8. ^ Biben, Thierry; Kassner, Klaus; Misbah, Chaouqi (2005). "Enfoque de campo de fase para dinámica de vesículas tridimensionales". Physical Review E . 72 (4): 041921. Bibcode :2005PhRvE..72d1921B. doi :10.1103/PhysRevE.72.041921. PMID  16383434.
  9. ^ Ashour, Mohammed; Valizadeh, Navid; Rabczuk, Timon (2021). "Análisis isogeométrico para un problema de optimización restringido por campo de fase de la evolución morfológica de vesículas en campos eléctricos". Métodos informáticos en mecánica aplicada e ingeniería . 377 . Elsevier BV: 113669. Bibcode :2021CMAME.377k3669A. doi :10.1016/j.cma.2021.113669. ISSN  0045-7825. S2CID  233580102.
  10. ^ Valizadeh, Navid; Rabczuk, Timon (2022). "Análisis isogeométrico de la hidrodinámica de vesículas utilizando un enfoque monolítico de campo de fase". Métodos informáticos en mecánica aplicada e ingeniería . 388 . Elsevier BV: 114191. Bibcode :2022CMAME.388k4191V. doi :10.1016/j.cma.2021.114191. ISSN  0045-7825. S2CID  240657318.
  11. ^ Valizadeh, Navid; Rabczuk, Timon (2019). "Análisis isogeométrico para modelos de campo de fase de ecuaciones diferenciales parciales geométricas y ecuaciones diferenciales parciales de orden superior en superficies estacionarias y evolutivas". Métodos informáticos en mecánica aplicada e ingeniería . 351 . Elsevier BV: 599–642. Bibcode :2019CMAME.351..599V. doi :10.1016/j.cma.2019.03.043. ISSN  0045-7825. S2CID  145903238.
  12. ^ GJ Fix, en Problemas de contorno libre: teoría y aplicaciones, Ed. A. Fasano y M. Primicerio, pág. 580, Pitman (Boston, 1983).
  13. ^ Langer, JS (1986). "Modelos de formación de patrones en transiciones de fase de primer orden". Directions in Condensed Matter Physics . Serie sobre Directions in Condensed Matter Physics. Vol. 1. Singapur: World Scientific. págs. 165–186. Código Bibliográfico : 1986dcmp.book..165L. doi : 10.1142/9789814415309_0005. ISBN: 978-9971-978-42-6.
  14. ^ Langer, JS (1980). "Inestabilidades y formación de patrones en el crecimiento de cristales". Reseñas de Física Moderna . 52 (1): 1–28. Bibcode :1980RvMP...52....1L. doi :10.1103/RevModPhys.52.1.
  15. ^ ab Karma, Alain; Rappel, Wouter-Jan (1998). "Modelado cuantitativo de campo de fase del crecimiento dendrítico en dos y tres dimensiones". Physical Review E . 57 (4): 4323. Bibcode :1998PhRvE..57.4323K. doi :10.1103/PhysRevE.57.4323.
  16. ^ Hohenberg, P.; Halperin, B. (1977). "Teoría de los fenómenos críticos dinámicos". Reseñas de Física Moderna . 49 (3): 435. Bibcode :1977RvMP...49..435H. doi :10.1103/RevModPhys.49.435. S2CID  122636335.
  17. ^ Yue, Pengtao; Zhou, Chunfeng; Feng, James J. (2007). "Contracción espontánea de gotas y conservación de masa en simulaciones de campo de fase". Journal of Computational Physics . 223 (1): 1–9. Bibcode :2007JCoPh.223....1Y. CiteSeerX 10.1.1.583.2109 . doi :10.1016/j.jcp.2006.11.020. 
  18. ^ Donaldson, AA; Kirpalani, DM; MacChi, A. (2011). "Seguimiento de la interfaz difusa de fluidos inmiscibles: mejora de la continuidad de fase mediante la selección de densidad de energía libre". Revista internacional de flujo multifásico . 37 (7): 777. Bibcode :2011IJMF...37..777D. doi :10.1016/j.ijmultiphaseflow.2011.02.002.
  19. ^ McFadden, GB; Wheeler, AA; Anderson, DM (2000). "Asintótica de interfaz delgada para un enfoque de energía/entropía para modelos de campo de fase con conductividades desiguales". Physica D: Nonlinear Phenomena . 144 (1–2): 154–168. Bibcode :2000PhyD..144..154M. doi :10.1016/S0167-2789(00)00064-6. hdl : 2060/20000014455 . S2CID  119641692.
  20. ^ Jacqmin, David (1999). "Cálculo de flujos de Navier-Stokes bifásicos mediante modelado de campo de fase". Journal of Computational Physics . 155 (1): 96–127. Bibcode :1999JCoPh.155...96J. doi :10.1006/jcph.1999.6332.
  21. ^ Benítez, R.; Ramírez-Piscina, L. (2005). "Proyección de interfaz nítida de un modelo de campo de fase fluctuante". Physical Review E . 71 (6): 061603. arXiv : cond-mat/0409707 . Bibcode :2005PhRvE..71f1603B. doi :10.1103/PhysRevE.71.061603. PMID  16089744. S2CID  28956874.
  22. ^ Schmitz, GJ; Böttger, B.; Eiken, J.; Apel, M.; Viardin, A.; Carré, A.; Laschet, G. (2011). "Simulación basada en campo de fase de la evolución de la microestructura en grados de aleación técnica". Revista internacional de avances en ciencias de la ingeniería y matemáticas aplicadas . 2 (4): 126. doi :10.1007/s12572-011-0026-y. S2CID  121915897.
  23. ^ Tanné, E.; Li, T.; Bourdin, B.; Marigo, J.-J.; Maurini, C. (2018). "Nuclearización de grietas en modelos de campo de fase variacional de fractura frágil" (PDF) . Revista de mecánica y física de sólidos . 110 : 80–99. Código Bibliográfico :2018JMPSo.110...80T. doi :10.1016/j.jmps.2017.09.006. S2CID  20139734.
  24. ^ Kumar, A.; Bourdin, B.; Francfort, GA; Lopez-Pamies, O. (2020). "Revisitando la nucleación en el enfoque de campo de fase para la fractura frágil". Journal of the Mechanics and Physics of Solids . 142 : 104027. Bibcode :2020JMPSo.14204027K. doi : 10.1016/j.jmps.2020.104027 .
  25. ^ Najem, Sara; Grant, Martin (9 de mayo de 2016). "Modelo de campo de fase para la migración celular colectiva". Physical Review E . 93 (5): 052405. Bibcode :2016PhRvE..93e2405N. doi :10.1103/PhysRevE.93.052405. PMID  27300922.
  26. ^ "Modelo de campo de fase para monocapas celulares: un estudio de migración de células cancerosasautores: Benoit Palmieri y Martin Grant | Perimeter Institute". www2.perimeterinstitute.ca . Consultado el 5 de noviembre de 2021 .
  27. ^ DeWitt, S.; Rudraraju, S.; Montiel, D.; Montiel, D.; Andrews, WB; Thornton, K. (2020). "PRISMS-PF: Un marco general para el modelado de campo de fase con un método de elementos finitos sin matriz". npj Comput Mater . 6 : 29. Bibcode :2020npjCM...6...29D. doi : 10.1038/s41524-020-0298-5 .

Lectura adicional