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.
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 smearing 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 emparejan 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 delgada) 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
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
PACE3D – Algoritmos paralelos para la evolución de cristales en 3D es un paquete de simulación de campo de fase paralelizado que incluye transformaciones multifase y multicomponente, estructuras de grano a gran escala y acoplamiento con interacciones elásticas, plásticas, magnéticas y de flujo de fluidos. Se desarrolló en la Universidad de Ciencias Aplicadas de Karlsruhe y el Instituto Tecnológico de Karlsruhe .
El Proyecto de Simulación de Microestructura a Mesoescala (MMSP) es una colección de clases de C++ para simulación de microestructura basada en cuadrícula.
El software de simulación de evolución de microestructuras (MICRESS) es un paquete de simulación de campo multifase y multicomponente acoplado a bases de datos termodinámicas y cinéticas. Lo desarrolla y mantiene ACCESS eV.
MOOSE, marco de elementos finitos multifísicos C++ de código abierto y masivamente paralelo con soporte para simulaciones de campo de fase desarrollado en el Laboratorio Nacional de Idaho.
PhasePot es una herramienta de simulación de microestructura basada en Windows, que utiliza una combinación de modelos de campo de fase y Monte Carlo Potts.
OpenPhase es un software de código abierto para la simulación de la formación de microestructuras en sistemas que experimentan una transformación de fase de primer orden basada en el modelo de campo multifásico.
mef90/vDef es un simulador de fractura de campo de fase variacional de código abierto basado en la teoría desarrollada en. [3] [4] [5]
MicroSim es una pila de software que consta de códigos de campo de fase que ofrecen flexibilidad con discretización, modelos y hardware de computación de alto rendimiento (CPU/GPU) en el que pueden ejecutarse.
PRISMS-PF es un código de elementos finitos masivamente paralelo para realizar simulaciones de campo de fase y otras simulaciones relacionadas con la evolución de la microestructura. [27] Se basa en la biblioteca de elementos finitos deal.II y es desarrollado y mantenido por el Centro PRISMS de la Universidad de Michigan.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ GJ Fix, en Problemas de contorno libre: teoría y aplicaciones, Ed. A. Fasano y M. Primicerio, pág. 580, Pitman (Boston, 1983).
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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.
^ 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. Bibcode :2018JMPSo.110...80T. doi :10.1016/j.jmps.2017.09.006. S2CID 20139734.
^ 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". Revista de mecánica y física de sólidos . 142 : 104027. Bibcode :2020JMPSo.14204027K. doi : 10.1016/j.jmps.2020.104027 .
^ 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.
^ "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 .
^ 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
Chen, Long-Qing (2002). "Modelos de campo de fase para la evolución de la microestructura". Revisión anual de investigación de materiales . 32 : 113–140. doi :10.1146/annurev.matsci.32.112001.132041.
Moelans, Nele; Blanpain, Bart; Wollants, Patrick (2008). "Introducción al modelado de campo de fase de la evolución de la microestructura". Calphad . 32 (2): 268. doi :10.1016/j.calphad.2007.11.003.
Steinbach, Ingo (2009). "Modelos de campo de fase en la ciencia de los materiales". Modelado y simulación en ciencia e ingeniería de materiales . 17 (7): 073001. Bibcode :2009MSMSE..17g3001S. doi :10.1088/0965-0393/17/7/073001. S2CID 3383625.
Fries, Suzana G.; Boettger, Bernd; Eiken, Janin; Steinbach, Ingo (2009). "Actualización de CALPHAD para simulación de microestructuras: el método de campo de fases". Revista Internacional de Investigación de Materiales . 100 (2): 128. Bibcode :2009IJMR..100..128F. doi :10.3139/146.110013. S2CID 138203262.
Qin, RS; Bhadeshia, HK (2010). "Método de campo de fase" (PDF) . Ciencia y tecnología de materiales . 26 (7): 803. Bibcode :2010MatST..26..803Q. doi :10.1179/174328409X453190. S2CID 136124682.
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.
Gonzalez-Cinca, R.; Folch, R.; Benitez, R.; Ramirez-Piscina, L.; Casademunt, J.; Hernandez-Machado, A. (2003). "Modelos de campo de fase en la formación de patrones interfaciales fuera de equilibrio". En Advances in Condensed Matter and Statistical Mechanics, ed. Por E. Korutcheva y R. Cuerno, Nova Science Publishers (Nueva York, ), Pp . 2004 : 203–236. arXiv : cond-mat/0305058 . Código Bibliográfico :2003cond.mat..5058G.Una revisión de los modelos de campo de fase.
Provatas, Nikolas; Elder, Ken (2010). Métodos de campo de fase en ciencia e ingeniería de materiales. Weinheim, Alemania: Wiley-VCH Verlag GmbH & Co. KGaA. doi :10.1002/9783527631520. ISBN 9783527631520
Steinbach, I.: "Concepto de materia de campo de fase cuántica: gravedad emergente en el universo dinámico", Zeitschrift für Naturforschung A 72 1 (2017) doi :10.1515/zna-2016-0270
Schmitz, GJ: "Un enfoque combinado de entropía y campo de fases para la gravedad", Entropy 2017 , 19 (4) 151; doi :10.3390/e19040151