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 esta se mueva de tal manera que la fuente de tensión se reduce y la estructura vuelve 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, aeronaves, 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 de Tacoma Narrows, es probablemente uno de los ejemplos más infames de falla a gran escala. Las alas de las aeronaves y las palas de las turbinas pueden romperse debido a las oscilaciones de FSI. Una lengüeta en realidad produce sonido porque el sistema de ecuaciones que gobierna su dinámica tiene soluciones oscilatorias. La dinámica de las válvulas de lengüeta utilizadas en los motores de dos tiempos y compresores está gobernada por FSI. El acto de " soplar una frambuesa " es otro ejemplo de este tipo. La interacción entre los componentes de la máquina tribológica, como los cojinetes y los engranajes , y el lubricante también es un ejemplo de FSI. [2] El lubricante fluye entre los componentes sólidos en contacto y causa una deformación elástica en ellos durante este proceso. 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 considerables cargas laterales inestables en la estructura de la boquilla. [7] Además de los efectos impulsados ​​por la presión, la 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 sanguínea 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. Este efecto es especialmente imperativo de tener en cuenta 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 del paciente. El cuello de un aneurisma es el más susceptible a los cambios en WSS. Si la pared aneurismática se debilita lo suficiente, corre el riesgo de romperse cuando la WSS se vuelve demasiado alta. Los modelos FSI contienen una WSS general más baja en comparación con los modelos no flexibles. Esto es significativo porque el modelado incorrecto de los aneurismas podría llevar a los médicos a decidir realizar cirugía invasiva en pacientes que no tenían un alto riesgo de ruptura. Si bien la FSI ofrece un mejor análisis, tiene el costo de un tiempo computacional muy aumentado. Los modelos que no cumplen con las normas tienen un tiempo de cálculo de unas pocas horas, mientras que los modelos FSI pueden tardar hasta siete días en completarse. Esto hace que los modelos FSI sean más útiles para las medidas preventivas de los aneurismas detectados a tiempo, pero inutilizables para situaciones de emergencia en las que el aneurisma puede haberse roto ya. [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, por lo que 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 un solucionador de flujo existente y un solucionador estructural están acoplados. 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 que es una ventaja atractiva. Sin embargo, se debe 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 del fluido que se desplaza 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 conforme y métodos de malla no conforme. [15] Otras clasificaciones pueden ser métodos basados ​​en malla y métodos sin malla. [16]

Simulación numérica

El método de Newton-Raphson o una iteración de punto fijo diferente se pueden utilizar 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 lineal 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 matriz , utilizando una aproximación de diferencias finitas del producto jacobiano-vector.

Mientras que los métodos de Newton-Raphson resuelven el problema de flujo y estructural 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 se puede escribir como un problema de búsqueda de raíces o 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 de física reducida lineal. [23] [24] El método de interfaz cuasi-Newton con aproximación para el inverso del jacobiano a partir de un modelo de mínimos cuadrados acopla un solucionador de flujo de caja negra y un solucionador estructural [25] por medio de la información que se ha recopilado durante las iteraciones de acoplamiento. Esta técnica se basa en la técnica de interfaz de bloque cuasi-Newton 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 la posición de la interfaz y la distribución de tensiones en la interfaz como incógnitas. Este sistema se resuelve con iteraciones de bloque cuasi-Newton del tipo Gauss-Seidel y los jacobianos del solucionador de flujo y el solucionador estructural se aproximan por medio de modelos de mínimos cuadrados. [26]

El problema de punto fijo se puede resolver con iteraciones de punto fijo, también llamadas iteraciones de Gauss-Seidel (en bloque), [21] lo que significa que el problema de flujo y el problema estructural se resuelven sucesivamente hasta que el cambio es menor que el criterio de convergencia. Sin embargo, las iteraciones convergen lentamente, si es que lo hacen, 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 se puede estabilizar y acelerar 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 denominados escalonados o débilmente acoplados no imponen el equilibrio en la interfaz fluido-estructura dentro de un paso de tiempo, pero son adecuados para la simulación de 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]

Véase 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 acoplada con el chapoteo transitorio del fluido". Journal of Vibration and 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 dinámico tridimensional de líquidos en tanques horizontales parcialmente llenos sujetos a excitaciones longitudinales y laterales simultáneas". European Journal of Mechanics B . 53 : 251–263. Bibcode :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 de balanceo de los vehículos cisterna". Revista de sonido y vibración . 333 (1): 263–282. Código Bibliográfico :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 en las cargas dinámicas de chapoteo del fluido y la estabilidad de balanceo de un camión cisterna parcialmente lleno". Revista Europea de Mecánica B . 46 : 46–58. Código Bibliográfico :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}}: Requiere citar revista |journal=( ayuda )
  8. ^ Daub, Dennis; Esser, Burkard; Gülhan, Ali (abril de 2020). "Experimentos sobre interacción fluido-estructura hipersónica de alta temperatura con deformación plástica". AIAA Journal . 58 (4): 1423–1431. Bibcode :2020AIAAJ..58.1423D. doi : 10.2514/1.J059150 . ISSN  0001-1452.
  9. ^ Taren, James A. (1965). "Aneurisma cerebral". 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, AK; 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". Journal of Physics: Conference Series . 722 (1): 012042. Bibcode :2016JPhCS.722a2042K. doi : 10.1088/1742-6596/722/1/012042 . ISSN  1742-6588.
  12. ^ Torii, Ryo; Oshima, Marie; Kobayashi, Toshio; Takagi, Kiyoshi; Tezduyar, Tayfun E. (15 de septiembre de 2009). "Modelado de la interacción fluido-estructura del flujo sanguíneo y aneurisma cerebral: importancia de las formas de las arterias y los aneurismas". Métodos informáticos en mecánica aplicada e ingeniería . Modelos y métodos en mecánica vascular y cardiovascular computacional. 198 (45): 3613–3621. Bibcode :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 ruptura 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 2014-10-31 . Consultado el 2014-11-28 .{{cite web}}: CS1 maint: 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 Bibliográfico :2021JCoPh.42910028Z. doi :10.1016/j.jcp.2020.110028. S2CID  208513116.
  17. ^ M. Heil (2004). "Un solucionador eficiente para la solución completamente acoplada de problemas de interacción fluido-estructura de gran desplazamiento". Métodos informáticos en mecánica aplicada e ingeniería . 193 (1–2): 1–23. Bibcode :2004CMAME.193....1H. doi :10.1016/j.cma.2003.09.006.
  18. ^ K.-J. Bathe ; 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. Bibcode :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 la interacción fluido-estructura con aplicación en biomecánica . Apuntes de clase en ciencia computacional e ingeniería. 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 particionados para la interacción fluido-estructura". Computers and Structures . 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 aplicada e ingeniería . 195 (17–18): 2028–2049. Bibcode :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 preacondicionado por subiteración para interacción fluido-estructura". Métodos informáticos en mecánica aplicada e ingeniería . 195 (17–18): 2124–2148. Bibcode :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 flujos sanguíneos" (PDF) . ESAIM: Modelado matemático 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 en geometrías basadas en imágenes médicas". Computers and Structures . 83 (2–3): 155–165. doi :10.1016/j.compstruc.2004.03.083.
  25. ^ J. Degroote; K.-J. Bathe; J. Vierendeels (2009). "Rendimiento de un nuevo procedimiento particionado frente a un procedimiento monolítico en la interacción fluido-estructura". Computers and Structures . 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 añadida en el diseño de algoritmos particionados para problemas de fluido-estructura" (PDF) . Métodos informáticos en mecánica aplicada e ingeniería . 194 (42–44): 4506–4527. Bibcode :2005CMAME.194.4506C. doi :10.1016/j.cma.2004.12.005. S2CID  122528121.
  28. ^ U. Küttler; W. Wall (2008). "Solucionadores de interacción fluido-estructura de punto fijo con relajación dinámica". Mecánica Computacional . 43 (1): 61–72. Bibcode :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. Jaiman; X. Jiao; P. Geubelle; E. Loth (2006). "Transferencia de carga conservativa a lo largo de la interfaz fluido-sólido curvada con mallas no coincidentes". Journal of Computational Physics . 218 (1): 372–397. Bibcode :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". AIAA Journal . 43 (12): 2549–2557. Bibcode :2005AIAAJ..43.2549V. doi :10.2514/1.3660.
  32. ^ Christiane Förster; Wolfgang A. Wall; Ekkehard Ramm (2006). P. Wesseling; E. Oñate; J. Périaux (eds.). El efecto de masa artificial añadida en algoritmos secuenciales de interacción fluido-estructura escalonados . Conferencia Europea sobre Dinámica de Fluidos Computacional ECCOMAS CFD 2006. Países Bajos.
  33. ^ Christiane Förster; Wolfgang A. Wall; Ekkehard Ramm (2007). "Inestabilidades de masa artificial añadida en acoplamiento escalonado secuencial de estructuras no lineales y flujos viscosos incompresibles". Métodos informáticos en mecánica aplicada e ingeniería . 196 (7): 1278–1293. Bibcode :2007CMAME.196.1278F. doi :10.1016/j.cma.2006.09.002.

Lectura adicional

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 Presenta el tema de las interacciones fluido-estructura (FSI) a estudiantes y profesionales y analiza las principales ideas en FSI con el objetivo de proporcionar una comprensión fundamental a los lectores que poseen una comprensión limitada o nula del tema.