Método de contorno inmerso

En dinámica de fluidos computacional, el método de los límites inmersos se refería originalmente a un enfoque desarrollado por Charles Peskin en 1972 para simular las interacciones entre fluidos y estructuras (fibras).[1]​ El tratamiento del acoplamiento entre las deformaciones de la estructura y el flujo de fluidos plantea una serie de problemas difíciles para las simulaciones numéricas (el límite elástico cambia el flujo del fluido y el fluido mueve el límite elástico simultáneamente). En el método de los límites inmersos, el fluido se representa en un sistema de coordenadas de Euler y la estructura se representa en coordenadas de Lagrange. Para los fluidos newtonianos regidos por las ecuaciones de Navier-Stokes, las ecuaciones del fluido son:

y si el flujo es incompresible, tenemos la condición adicional de que

Las estructuras sumergidas se representan normalmente como un conjunto de fibras unidimensionales, denotadas por . Cada fibra puede considerarse como una curva paramétrica , donde es la coordenada lagrangiana a lo largo de la fibra y es el tiempo. La física de la fibra se representa mediante una función de distribución de la fuerza de la fibra . En este término se pueden incorporar fuerzas elásticas, resistencia a la flexión o cualquier otro tipo de comportamiento. A continuación, la fuerza ejercida por la estructura sobre el fluido se interpola como un término fuente en la ecuación del momento utilizando

donde es la función δ de Dirac. El forzamiento puede extenderse a múltiples dimensiones para modelar superficies elásticas o sólidos tridimensionales. Suponiendo una estructura sin masa, la fibra elástica se mueve con la velocidad local del fluido y puede interpolarse mediante la función delta.

donde denota todo el dominio del fluido. La discretización de estas ecuaciones se puede realizar suponiendo una malla euleriana sobre el fluido y una malla lagrangiana separada sobre la fibra. Las aproximaciones de la distribución Delta mediante funciones más suaves nos permitirán interpolar entre las dos mallas. Cualquier solucionador de fluidos existente se puede acoplar a un solucionador para las ecuaciones de la fibra con el fin de resolver las ecuaciones de contorno inmerso. Se han aplicado variantes de este enfoque básico para simular una amplia variedad de sistemas mecánicos que implican estructuras elásticas que interactúan con fluidos.

Desde el desarrollo original de este método por Peskin, se han desarrollado diversos enfoques. Entre ellos se incluyen formulaciones estocásticas para sistemas microscópicos, materiales blandos viscoelásticos y fluidos complejos, como los métodos estocásticos de límites inmersos de Atzberger, Kramer y Peskin,[2][3]​ métodos para simular flujos sobre cuerpos sólidos inmersos complicados en rejillas que no se ajustan a la superficie del cuerpo Mittal e Iaccarino,[4]​ y otros enfoques que incorporan masa y grados de libertad rotacionales Olson, Lim, Cortez.[5]​ Los métodos para formas corporales complicadas incluyen el método de interfaz inmersa, el método de cuadrícula cartesiana, el método de fluido fantasma y los métodos de celda cortada, que clasifican los métodos de frontera inmersa en métodos de «fuerza continua» y «fuerza discreta». Se han desarrollado métodos para simulaciones de fluidos viscoelásticos, interfaces de fluidos curvos, sistemas biofísicos microscópicos (proteínas en membranas de bicapa lipídica, nadadores) y dispositivos de ingeniería, como los métodos estocásticos de límites inmersos de Atzberger, Kramer y Peskin,[6][7]​ Métodos lagrangianos eulerianos estocásticos de Atzberger,[8][9][10]​ Métodos de límites inmersos masivos de Mori,[11]​ y los métodos de límites inmersos rotacionales de Olson, Lim y Cortez.[5]

En general, para los métodos de fronteras inmersas y variantes relacionadas, existe una comunidad de investigación activa que sigue desarrollando nuevas técnicas e implementaciones de software relacionadas, e incorporando técnicas relacionadas en paquetes de simulación y software de ingeniería CAD. Para más detalles, véase más abajo.

Véase también

Referencias

  1. Peskin, Charles S (1 de octubre de 1972). «Patrones de flujo alrededor de las válvulas cardíacas: A numerical method». Journal of Computational Physics 10 (2): 252-271. Bibcode:1972JCoPh..10..252P. ISSN 0021-9991. doi:10.1016/0021-9991(72)90065-4. 
  2. Atzberger, Paul J. (2011). «Métodos estocásticos de Euler-Lagrange para interacciones fluido-estructura con fluctuaciones térmicas». Journal of Computational Physics 230 (8): 2821-2837. Bibcode:2011JCoPh.230.2821A. S2CID 6067032. arXiv:1009.5648. doi:10.1016/j.jcp.2010.12.028. 
  3. Atzberger, Paul (2013), «Incorporating Shear into Stochastic Eulerian Lagrangian Methods for Rheological Studies of Complex Fluids and Soft Materials», Physica D 265: 57-70, arXiv:2212.10651, doi:10.1016/j.physd.2013.09.002 .
  4. Mittal y Iaccarino, 2005.
  5. a b Olson, S.; Lim, S.; Cortez, R. (2013). «Modelización de la dinámica de una varilla elástica con curvatura intrínseca y torsión utilizando una formulación regularizada de Stokes». Journal of Computational Physics 238: 169-187. doi:10.1016/j.jcp.2012.12.026. 
  6. Atzberger, Paul J. (2011). «Métodos lagrangianos eulerianos estocásticos para interacciones fluido-estructura con fluctuaciones térmicas». Journal of Computational Physics 230 (8): 2821-2837. Bibcode:2011JCoPh.230.2821A. S2CID 6067032. arXiv:1009.5648. doi:10.1016/j.jcp.2010.12.028. 
  7. Rower, David A.; Padidar, Misha; Atzberger, Paul J. (abril de 2022). «Métodos de hidrodinámica fluctuante superficial para la dinámica de difusión-difusión de partículas y microestructuras dentro de interfaces fluidas curvas». Journal of Computational Physics 455: 110994. arXiv:1906.01146. doi:10.1016/j.jcp.2022.110994. 
  8. Atzberger, Paul J. (2011). «Métodos lagrangianos eulerianos estocásticos para interacciones entre estructuras fluidas con fluctuaciones térmicas». Journal of Computational Physics 230 (8): 2821-2837. Bibcode:2011JCoPh.230.2821A. S2CID 6067032. arXiv:1009.5648. doi:10.1016/j.jcp.2010.12.028. 
  9. Atzberger, Paul (2013), «Incorporating Shear into Stochastic Eulerian Lagrangian Methods for Rheological Studies of Complex Fluids and Soft Materials», Physica D 265: 57-70, arXiv:2212.10651, doi:10.1016/j.physd.2013.09.002 .
  10. Atzberger, Paul (2016). «Hydrodynamic Coupling of Particle Inclusions Embedded in Curved Lipid Bilayer Membranes». Soft Matter 12 (32): 6685-6707. PMID 27373277. arXiv:1601.06461. doi:10.1039/C6SM00194G. 
  11. Mori, Yoichiro; Peskin, Charles S. (2008). «Métodos implícitos de frontera inmersa de segundo orden con masa en la frontera». Métodos informáticos en mecánica aplicada e ingeniería 197 (25–28): 2049-2067. Bibcode:2008CMAME.197.2049M. doi:10.1016/j.cma.2007.05.028. 

Software: Códigos numéricos

Bibliografía