stringtranslate.com

Interacción fluido-estructura

La interacción fluido-estructura ( FSI ) es la interacción de alguna estructura móvil o deformable con un flujo de fluido interno o circundante. [1] Las interacciones fluido-estructura pueden ser estables u oscilatorias. En las interacciones oscilatorias, la tensión inducida en la estructura sólida hace que se mueva de tal manera que la fuente de tensión se reduce y la estructura regresa a su estado anterior solo para que el proceso se repita.

Propagación de una onda de presión a través de un fluido incompresible en un tubo flexible.

Ejemplos

Las interacciones fluido-estructura son una consideración crucial en el diseño de muchos sistemas de ingeniería, por ejemplo, automóviles, aviones, naves espaciales, motores y puentes. No considerar los efectos de las interacciones oscilatorias puede ser catastrófico, especialmente en estructuras que comprenden materiales susceptibles a la fatiga . El puente Tacoma Narrows (1940) , el primer puente Tacoma Narrows, es probablemente uno de los ejemplos más infames de fallas a gran escala. Las alas de los aviones y las palas de las turbinas pueden romperse debido a las oscilaciones del FSI. En realidad, una caña produce sonido porque el sistema de ecuaciones que gobierna su dinámica tiene soluciones oscilatorias. La dinámica de las válvulas de láminas utilizadas en motores y compresores de dos tiempos está gobernada por FSI. El acto de " soplar una frambuesa " es otro ejemplo de ello. La interacción entre los componentes tribológicos de la máquina, como rodamientos y engranajes , y el lubricante también es un ejemplo de FSI. [2] El lubricante fluye entre los componentes sólidos en contacto y durante este proceso provoca en ellos una deformación elástica. Las interacciones fluido-estructura también ocurren en contenedores en movimiento, donde las oscilaciones del líquido debido al movimiento del contenedor imponen magnitudes sustanciales de fuerzas y momentos a la estructura del contenedor que afectan la estabilidad del sistema de transporte de contenedores de una manera altamente adversa. [3] [4] [5] [6] Otro ejemplo destacado es el arranque de un motor de cohete, por ejemplo, el motor principal del transbordador espacial (SSME) , donde la FSI puede provocar cargas laterales inestables considerables en la estructura de la boquilla. [7] Además de los efectos impulsados ​​por la presión, el FSI también puede tener una gran influencia en las temperaturas de la superficie de los vehículos supersónicos e hipersónicos. [8]

Las interacciones fluido-estructura también desempeñan un papel importante en el modelado apropiado del flujo sanguíneo . Los vasos sanguíneos actúan como tubos flexibles que cambian de tamaño dinámicamente cuando hay cambios en la presión arterial y la velocidad del flujo. [9] No tener en cuenta esta propiedad de los vasos sanguíneos puede llevar a una sobreestimación significativa de la tensión de corte de la pared (WSS) resultante. Es especialmente imperativo tener en cuenta este efecto al analizar los aneurismas. Se ha convertido en una práctica común utilizar la dinámica de fluidos computacional para analizar modelos específicos de pacientes. El cuello de un aneurisma es el más susceptible a los cambios en el WSS. Si la pared del aneurisma se debilita lo suficiente, corre el riesgo de romperse cuando el WSS aumenta demasiado. Los modelos FSI contienen un WSS general más bajo en comparación con los modelos que no cumplen. Esto es importante porque un modelado incorrecto de los aneurismas podría llevar a los médicos a decidir realizar una cirugía invasiva en pacientes que no tenían un alto riesgo de ruptura. Si bien FSI ofrece un mejor análisis, tiene el costo de un tiempo de cálculo mucho mayor. Los modelos no conformes tienen un tiempo de cálculo de unas pocas horas, mientras que los modelos FSI podrían tardar hasta 7 días en terminar de ejecutarse. Esto lleva a que los modelos FSI sean más útiles como medidas preventivas para aneurismas detectados tempranamente, pero inutilizables en situaciones de emergencia en las que el aneurisma ya puede haberse roto. [10] [11] [12] [13]

Análisis

Los problemas de interacción fluido-estructura y los problemas multifísicos en general suelen ser demasiado complejos para resolverlos analíticamente y, por lo tanto, deben analizarse mediante experimentos o simulación numérica . La investigación en los campos de la dinámica de fluidos computacional y la dinámica estructural computacional aún está en curso, pero la madurez de estos campos permite la simulación numérica de la interacción fluido-estructura. [14] Existen dos enfoques principales para la simulación de problemas de interacción fluido-estructura:

El enfoque monolítico requiere un código desarrollado para esta combinación particular de problemas físicos, mientras que el enfoque particionado preserva la modularidad del software porque están acoplados un solucionador de flujo y un solucionador estructural existentes. Además, el enfoque particionado facilita la solución de las ecuaciones de flujo y las ecuaciones estructurales con técnicas diferentes, posiblemente más eficientes, que se han desarrollado específicamente para ecuaciones de flujo o ecuaciones estructurales. Por otro lado, se requiere el desarrollo de un algoritmo de acoplamiento estable y preciso en simulaciones particionadas. En conclusión, el enfoque particionado permite reutilizar el software existente, lo cual es una ventaja atractiva. Sin embargo, es necesario tener en cuenta la estabilidad del método de acoplamiento. Esto es especialmente difícil si la masa de la estructura en movimiento es pequeña en comparación con la masa de fluido que es desplazada por el movimiento de la estructura.

Además, el tratamiento de las mallas introduce otras clasificaciones del análisis FSI. Por ejemplo, se pueden clasificar como métodos de malla conformes y métodos de malla no conformes. [15] Otras clasificaciones pueden ser métodos basados ​​en malla y métodos sin malla. [dieciséis]

Simulación numérica

Se puede utilizar el método de Newton-Raphson o una iteración de punto fijo diferente para resolver problemas FSI. Los métodos basados ​​en la iteración de Newton-Raphson se utilizan tanto en el enfoque monolítico [17] [18] [19] como en el particionado [20] [21] . Estos métodos resuelven las ecuaciones de flujo no lineales y las ecuaciones estructurales en todo el dominio de fluidos y sólidos con el método de Newton-Raphson. El sistema de ecuaciones lineales dentro de la iteración de Newton-Raphson se puede resolver sin conocimiento del jacobiano con un método iterativo sin matrices , utilizando una aproximación en diferencias finitas del producto vectorial jacobiano.

Mientras que los métodos de Newton-Raphson resuelven el problema estructural y de flujo para el estado en todo el dominio de fluidos y sólidos, también es posible reformular un problema FSI como un sistema con solo los grados de libertad en la posición de la interfaz como incógnitas. Esta descomposición del dominio condensa el error del problema FSI en un subespacio relacionado con la interfaz. [22] Por lo tanto, el problema FSI puede escribirse como un problema de búsqueda de raíces o como un problema de punto fijo, con la posición de la interfaz como incógnita.

Los métodos de interfaz Newton-Raphson resuelven este problema de búsqueda de raíces con iteraciones de Newton-Raphson, por ejemplo, con una aproximación del jacobiano a partir de un modelo lineal de física reducida. [23] [24] El método cuasi-Newton de interfaz con aproximación para la inversa del jacobiano a partir de un modelo de mínimos cuadrados acopla un solucionador de flujo de caja negra y un solucionador estructural [25] mediante la información recopilada durante el proceso. iteraciones de acoplamiento. Esta técnica se basa en la técnica cuasi-Newton del bloque de interfaz con una aproximación para los jacobianos a partir de modelos de mínimos cuadrados que reformula el problema FSI como un sistema de ecuaciones con tanto la posición de la interfaz como la distribución de tensiones en la interfaz como incógnitas. Este sistema se resuelve con iteraciones cuasi-Newton en bloque del tipo Gauss-Seidel y los jacobianos del solucionador de flujo y del solucionador estructural se aproximan mediante modelos de mínimos cuadrados. [26]

El problema de punto fijo se puede resolver con iteraciones de punto fijo, también llamadas iteraciones (en bloque) de Gauss-Seidel, [21] lo que significa que el problema de flujo y el problema estructural se resuelven sucesivamente hasta que el cambio sea menor que el criterio de convergencia. Sin embargo, las iteraciones convergen lentamente, en todo caso, especialmente cuando la interacción entre el fluido y la estructura es fuerte debido a una alta relación de densidad fluido/estructura o la incompresibilidad del fluido. [27] La ​​convergencia de las iteraciones de punto fijo puede estabilizarse y acelerarse mediante la relajación de Aitken y la relajación de descenso más pronunciado, que adaptan el factor de relajación en cada iteración en función de las iteraciones anteriores. [28]

Si la interacción entre el fluido y la estructura es débil, sólo se requiere una iteración de punto fijo dentro de cada paso de tiempo. Estos métodos llamados escalonados o débilmente acoplados no imponen el equilibrio en la interfaz fluido-estructura dentro de un intervalo de tiempo, pero son adecuados para la simulación de la aeroelasticidad con una estructura pesada y bastante rígida. Varios estudios han analizado la estabilidad de algoritmos particionados para la simulación de la interacción fluido-estructura [27] [29] [30] . [31] [32] [33]

Ver también

Códigos fuente abiertos

Códigos Académicos

Códigos Comerciales

Referencias

  1. ^ Bungartz, Hans-Joachim; Schäfer, Michael, eds. (2006). Interacción fluido-estructura: modelado, simulación, optimización . Springer-Verlag . ISBN 978-3-540-34595-4.
  2. ^ Singh, Kushagra; Sadeghi, Farshid; Russell, Thomas; Lorenz, Steven J.; Peterson, Wyatt; Villarreal, Jaret; Jinmon, Takumi (1 de septiembre de 2021). "Modelado de interacción fluido-estructura de contactos de línea lubricados elastohidrodinámicamente". Revista de Tribología . 143 (9). doi :10.1115/1.4049260. ISSN  0742-4787. S2CID  230619508.
  3. ^ Kolaei, Amir; Rakheja, Subhash; Richard, Marc J. (25 de enero de 2016). "Una metodología eficiente para simular la dinámica de balanceo de un vehículo cisterna junto con un chapoteo transitorio de fluido". Revista de Vibración y Control . 23 (19): 3216–3232. doi :10.1177/1077546315627565. ISSN  1077-5463. S2CID  123621791.
  4. ^ Kolaei, Amir; Rakheja, Subhash; Richard, Marc J. (1 de septiembre de 2015). "Chapoteo de líquido dinámico tridimensional en tanques horizontales parcialmente llenos sujetos a excitaciones longitudinales y laterales simultáneas". Revista Europea de Mecánica B. 53 : 251–263. Código Bib : 2015EJMF...53..251K. doi :10.1016/j.euromechflu.2015.06.001.
  5. ^ Kolaei, Amir; Rakheja, Subhash; Richard, Marc J. (6 de enero de 2014). "Rango de aplicabilidad de la teoría del chapoteo lineal del fluido para predecir el chapoteo lateral transitorio y la estabilidad del balanceo de vehículos cisterna". Revista de Sonido y Vibración . 333 (1): 263–282. Código Bib : 2014JSV...333..263K. doi :10.1016/j.jsv.2013.09.002.
  6. ^ Kolaei, Amir; Rakheja, Subhash; Richard, Marc J. (1 de julio de 2014). "Efectos de la sección transversal del tanque sobre las cargas dinámicas de chapoteo de fluidos y la estabilidad de balanceo de un camión cisterna parcialmente lleno". Revista Europea de Mecánica B. 46 : 46–58. Código Bib : 2014EJMF...46...46K. doi :10.1016/j.euromechflu.2014.01.008.
  7. ^ Manuel, Frey (2001). "Behandlung von Strömungsproblemen in Raketendüsen bei Überexpansion" (en alemán). doi :10.18419/opus-3650. {{cite journal}}: Citar diario requiere |journal=( ayuda )
  8. ^ Embadurnamiento, Dennis; Esser, Burkard; Gülhan, Ali (abril de 2020). "Experimentos sobre la interacción estructura-fluido hipersónico de alta temperatura con deformación plástica". Revista AIAA . 58 (4): 1423-1431. Código Bib : 2020AIAAJ..58.1423D. doi : 10.2514/1.J059150 . ISSN  0001-1452.
  9. ^ Taren, James A. (1965). "Aneurisma cerebral". La Revista Estadounidense de Enfermería . 65 (4): 88–91. doi :10.2307/3453223. ISSN  0002-936X. JSTOR  3453223. PMID  14258014. S2CID  31190911.
  10. ^ Sforza, Daniel M.; Putman, Christopher M.; Cebral, Juan R. (junio de 2012). "Dinámica de fluidos computacional en aneurismas cerebrales". Revista internacional de métodos numéricos en ingeniería biomédica . 28 (6–7): 801–808. doi :10.1002/cnm.1481. ISSN  2040-7939. PMC 4221804 . PMID  25364852. 
  11. ^ Khe, Alaska; Cherevko, AA; Chupakhin, AP; Bobkova, MS; Krivoshapkin, AL; Orlov, K Yu (junio de 2016). "Hemodinámica del aneurisma cerebral gigante: una comparación entre los modelos FSI de pared rígida, unidireccional y bidireccional". Revista de Física: Serie de conferencias . 722 (1): 012042. Código bibliográfico : 2016JPhCS.722a2042K. doi : 10.1088/1742-6596/722/1/012042 . ISSN  1742-6588.
  12. ^ Torii, Ryo; Oshima, María; Kobayashi, Toshio; Takagi, Kiyoshi; Tezduyar, Tayfun E. (15 de septiembre de 2009). "Modelado de interacción fluido-estructura del flujo sanguíneo y aneurisma cerebral: importancia de las formas de arterias y aneurismas". Métodos Informáticos en Mecánica e Ingeniería Aplicadas . Modelos y métodos en mecánica vascular y cardiovascular computacional. 198 (45): 3613–3621. Código Bib : 2009CMAME.198.3613T. doi :10.1016/j.cma.2008.08.020. ISSN  0045-7825.
  13. ^ Razaghi, Reza; Biglari, Hasan; Karimi, Alireza (1 de julio de 2019). "Riesgo de rotura del aneurisma cerebral en relación con una lesión cerebral traumática utilizando un modelo de interacción fluido-estructura específico del paciente". Métodos y Programas Informáticos en Biomedicina . 176 : 9–16. doi :10.1016/j.cmpb.2019.04.015. ISSN  0169-2607. PMID  31200915. S2CID  155305862.
  14. ^ JF Sigrist (2015). Interacción fluido-estructura: introducción al acoplamiento de elementos finitos . Wiley ( ISBN 978-1-119-95227-5 ) 
  15. ^ "Copia archivada" (PDF) . Archivado desde el original (PDF) el 31 de octubre de 2014 . Consultado el 28 de noviembre de 2014 .{{cite web}}: Mantenimiento CS1: copia archivada como título ( enlace )
  16. ^ Zhang, Chi; Rezavand, Massoud; Hu, Xiangyu (15 de marzo de 2021). "Un método SPH de resolución múltiple para interacciones fluido-estructura". Revista de Física Computacional . 429 : 110028. arXiv : 1911.13255 . Código Bib : 2021JCoPh.42910028Z. doi : 10.1016/j.jcp.2020.110028. S2CID  208513116.
  17. ^ M.Heil (2004). "Un solucionador eficiente para la solución totalmente acoplada de problemas de interacción fluido-estructura de gran desplazamiento". Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 193 (1–2): 1–23. Código Bib : 2004CMAME.193....1H. doi :10.1016/j.cma.2003.09.006.
  18. ^ K.-J. Bañarse ; H. Zhang (2004). "Desarrollos de elementos finitos para flujos de fluidos generales con interacciones estructurales". Revista internacional de métodos numéricos en ingeniería . 60 (1): 213–232. Código Bib : 2004IJNME..60..213B. CiteSeerX 10.1.1.163.1531 . doi :10.1002/nme.959. S2CID  17143434. 
  19. ^ J. Hron, S. Turek (2006). H.-J. Bungartz; M. Schäfer (eds.). "Un solucionador monolítico FEM/multigrid para la formulación ALE de interacción fluido-estructura con aplicación en biomecánica" . Apuntes de conferencias sobre ingeniería y ciencias computacionales. vol. Interacción fluido-estructura: modelado, simulación, optimización. Springer-Verlag . págs. 146-170. ISBN 978-3-540-34595-4.
  20. ^ H. Matthies; J. Steindorf (2003). "Algoritmos de acoplamiento fuerte particionado para la interacción fluido-estructura". Computadoras y Estructuras . 81 (8–11): 805–812. CiteSeerX 10.1.1.487.5577 . doi :10.1016/S0045-7949(02)00409-1. 
  21. ^ ab H. Matthies; R. Niekamp; J. Steindorf (2006). "Algoritmos para procedimientos de acoplamiento fuerte". Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 195 (17-18): 2028-2049. Código Bib : 2006CMAME.195.2028M. doi :10.1016/j.cma.2004.11.032.
  22. ^ C. Michler; E. van Brummelen; R. de Borst (2006). "Análisis de amplificación de errores de GMRES precondicionado por subiteración para la interacción fluido-estructura". Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 195 (17–18): 2124–2148. Código Bib : 2006CMAME.195.2124M. doi :10.1016/j.cma.2005.01.018.
  23. ^ J.-F. Gerbeau; M. Vidrascu (2003). "Un algoritmo cuasi-Newton basado en un modelo reducido para problemas de interacción fluido-estructura en los flujos sanguíneos" (PDF) . ESAIM: Modelación Matemática y Análisis Numérico . 37 (4): 631–648. doi :10.1051/m2an:2003049.
  24. ^ J.-F. Gerbeau; M. Vidrascu; P. Frey (2005). "Interacción fluido-estructura en flujos sanguíneos sobre geometrías basadas en imágenes médicas". Computadoras y Estructuras . 83 (2–3): 155–165. doi :10.1016/j.compstruc.2004.03.083.
  25. ^ J. Degroote; K.-J. Bañarse; J. Vierendeels (2009). "Rendimiento de un nuevo procedimiento particionado versus un procedimiento monolítico en la interacción fluido-estructura". Computadoras y Estructuras . 87 (11–12): 793–801. CiteSeerX 10.1.1.163.827 . doi :10.1016/j.compstruc.2008.11.013. 
  26. ^ J. Vierendeels; L. Lanoye; J. Degroote; P. Verdonck (2007). "Acoplamiento implícito de problemas de interacción fluido-estructura particionada con modelos de orden reducido". Computadoras y Estructuras . 85 (11–14): 970–976. doi :10.1016/j.compstruc.2006.11.006.
  27. ^ ab P. Causin; J.-F. Gerbeau; F.Nobile (2005). "Efecto de masa agregada en el diseño de algoritmos particionados para problemas de estructura de fluidos" (PDF) . Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 194 (42–44): 4506–4527. Código Bib : 2005CMAME.194.4506C. doi :10.1016/j.cma.2004.12.005. S2CID  122528121.
  28. ^ U. Küttler; W. Muro (2008). "Solucionadores de interacción fluido-estructura de punto fijo con relajación dinámica". Mecánica Computacional . 43 (1): 61–72. Código Bib : 2008CompM..43...61K. doi :10.1007/s00466-008-0255-5. S2CID  122209351.
  29. ^ J. Degroote; P. Bruggeman; R. Haelterman; J. Vierendeels (2008). "Estabilidad de una técnica de acoplamiento para solucionadores particionados en aplicaciones FSI". Computadoras y Estructuras . 86 (23–24): 2224–2234. doi :10.1016/j.compstruc.2008.05.005. hdl : 1854/LU-533350 .
  30. ^ R. Jaimán; X. Jiao; P. Geubelle; E. Loth (2006). "Transferencia de carga conservadora a lo largo de una interfaz curva fluido-sólido con mallas que no coinciden". Revista de Física Computacional . 218 (1): 372–397. Código Bib : 2006JCoPh.218..372J. CiteSeerX 10.1.1.147.4391 . doi : 10.1016/j.jcp.2006.02.016. 
  31. ^ J. Vierendeels; K. Dumont; E. Dick; P. Verdonck (2005). "Análisis y estabilización del algoritmo de interacción fluido-estructura para el movimiento de cuerpos rígidos". Revista AIAA . 43 (12): 2549–2557. Código Bib : 2005AIAAJ..43.2549V. doi : 10.2514/1.3660.
  32. ^ Christiane Forster; Wolfgang A. Muro; Ekkehard Ramm (2006). P. Wesseling; E. Oñate; J. Périaux (eds.). "El efecto de masa agregada artificial en algoritmos secuenciales escalonados de interacción fluido-estructura" . Conferencia europea sobre dinámica de fluidos computacional ECCOMAS CFD 2006. Países Bajos.
  33. ^ Christiane Forster; Wolfgang A. Muro; Ekkehard Ramm (2007). "Inestabilidades de masa artificiales añadidas en acoplamiento secuencial escalonado de estructuras no lineales y flujos viscosos incompresibles". Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 196 (7): 1278-1293. Código Bib : 2007CMAME.196.1278F. doi :10.1016/j.cma.2006.09.002.

Otras lecturas

Modarres-Sadeghi, Yahya: Introducción a las interacciones fluido-estructura, 2021, Springer Nature, 978-3-030-85882-7, http://dx.doi.org/10.1007/978-3-030-85884-1 Introduces el tema de Interacciones fluido-estructura (FSI) a estudiantes y profesionales y analiza las ideas principales en FSI con el objetivo de proporcionar una comprensión fundamental a los lectores que poseen una comprensión limitada o nula del tema.