Articulo de referencia

interacción fluido-estructura

La interacción fluido-estructura ( IFE ) es la interacción de una estructura móvil o deformable con un flujo de fluido interno o circundante. [ 1 ] Las interacciones fluido-estr...

La interacción fluido-estructura ( IFE ) es la interacción de una 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 deformación inducida en la estructura sólida provoca su movimiento, reduciendo así la fuente de deformación, y la estructura regresa a su estado original 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 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 la interacción fluido-estructura. Una lengüeta produce sonido porque el sistema de ecuaciones que rige su dinámica tiene soluciones oscilatorias. La dinámica de las válvulas de lengüeta utilizadas en motores de dos tiempos y compresores está regida por la interacción fluido-estructura. El acto de " soplar una frambuesa " es otro ejemplo de esto. La interacción entre componentes de máquinas tribológicas, como cojinetes y engranajes , y el lubricante también es un ejemplo de interacción fluido-estructura. [ 2 ] El lubricante fluye entre los componentes sólidos en contacto y causa 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 debidas 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 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 generar cargas laterales inestables considerables en la estructura de la tobera. [ 7 ] Además de los efectos impulsados ​​por la presión, la FSI también puede tener una gran influencia en las temperaturas superficiales en vehículos supersónicos e hipersónicos. [ 8 ]

Las interacciones fluido-estructura también desempeñan un papel fundamental en la modelización adecuada del flujo sanguíneo . Los vasos sanguíneos actúan como tubos elásticos que cambian de tamaño dinámicamente cuando hay variaciones 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 cizallamiento de la pared (WSS) resultante. Este efecto es especialmente importante a la hora de analizar aneurismas. Se ha convertido en práctica habitual el uso de dinámica de fluidos computacional para analizar modelos específicos de pacientes. El cuello de un aneurisma es la parte más susceptible a los cambios en la WSS. Si la pared del aneurisma se debilita lo suficiente, corre el riesgo de romperse cuando la WSS es demasiado alta. Los modelos FSI contienen una WSS globalmente menor en comparación con los modelos no elásticos. Esto es importante porque una modelización incorrecta de los aneurismas podría llevar a los médicos a decidir realizar cirugía invasiva en pacientes que no presentaban un alto riesgo de ruptura. Si bien FSI ofrece un mejor análisis, conlleva un coste computacional mucho mayor. Los modelos no conformes tienen un tiempo de cálculo de unas pocas horas, mientras que los modelos FSI pueden tardar hasta 7 días en completarse. Esto hace que los modelos FSI sean más útiles para medidas preventivas en aneurismas detectados precozmente, pero inutilizables en situaciones de emergencia donde el aneurisma ya se haya 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, 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:

  • Enfoque monolítico: las ecuaciones que rigen el flujo y el desplazamiento de la estructura se resuelven simultáneamente con un único solucionador.
  • Enfoque particionado: las ecuaciones que rigen el flujo y el desplazamiento de la estructura se resuelven por separado, con dos solucionadores distintos.

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

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 pueden utilizarse para resolver problemas de interacción fluido-estructura (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 fluido-sólido con el método de Newton-Raphson. El sistema de ecuaciones lineales dentro de la iteración de Newton-Raphson puede resolverse sin conocer el 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 estructura para el estado en todo el dominio fluido-sólido, 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 de 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 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 cuasi-Newton de interfaz 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 ] mediante la información recopilada durante las iteraciones de acoplamiento. Esta técnica se basa en la técnica cuasi-Newton de 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 la posición de la interfaz y la distribución de tensiones en la interfaz como incógnitas. Este sistema se resuelve con iteraciones cuasi-Newton de 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 de Gauss-Seidel (por bloques), [ 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, si es que convergen, especialmente cuando la interacción entre el fluido y la estructura es fuerte debido a una alta relación de densidad fluido/estructura o a 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, solo 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 de código abierto

  • solids4Foam , una caja de herramientas para OpenFOAM con capacidades para mecánica de sólidos e interacciones fluido-sólido.
  • oomph-lib
  • Página de Elmer FSI
  • Solucionadores biomédicos CBC.solve
  • preCICE (software)
  • La biblioteca multifísica SPHinXsys proporciona API de C++ para una simulación físicamente precisa y tiene como objetivo modelar sistemas dinámicos industriales acoplados, incluyendo dinámica de fluidos, sólidos, multicuerpos y más allá, con SPH (hidrodinámica de partículas suavizadas), un método computacional sin malla que utiliza la discretización de partículas.

Códigos académicos

  • Métodos estocásticos de contorno inmerso en 3D, P. Atzberger, UCSB
  • Método de contorno inmerso para mallas adaptativas en 3D, B. Griffith, NYU.
  • Método de contorno inmerso para mallas uniformes en 2D, A. Fogelson, Utah
  • IFLS, IFL, TU Braunschweig

Códigos comerciales

  • Acoplamiento multifísico de Abaqus
  • Aplicaciones de AcuSolve FSI
  • Página principal de ADINA FSI archivada el 28/04/2021 en Wayback Machine .
  • Página principal de Ansys sobre FSI
  • Altair RADIOSS
  • Simulación CFD de Autodesk
  • Simcenter STAR-CCM+ de Siemens Digital Industries Software
  • CoLyX - FSI y transformación de mallas de EVEN - Evolutionary Engineering AG
  • Acoplamiento multifísico Fluidyn-MP FSI Archivado el 17/05/2014 en Wayback Machine
  • Página principal de COMSOL FSI
  • Página principal de MpCCI
  • MSC Software MD Nastran
  • Software MSC Dytran
  • FINE/Oofelie FSI: Totalmente integrado y fuertemente acoplado para una mejor convergencia. Archivado el 1 de abril de 2012 en Wayback Machine.
  • Página principal de LS-DYNA
  • Fluidyn-MP FSI: Interacción fluido-estructura
  • CompassFEM Tdyn
  • CompassFEM SeaFEM
  • Software CFD Cradle SC/Tetra
  • Página principal de paracaídas FSI

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 (2021-09-01). "Modelado de interacción fluido-estructura de contactos lineales lubricados elastohidrodinámicamente" . Journal of Tribology . 143 (9) 091602. doi : 10.1115/1.4049260 . ISSN 0742-4787 . S2CID 230619508 .  
  3. Kolaei, Amir; Rakheja, Subhash; Richard, Marc J. (25-01-2016). "Una metodología eficiente para simular la dinámica de balanceo de un vehículo tanque acoplado 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. (2015-09-01). "Oscilación dinámica tridimensional de líquidos en tanques horizontales parcialmente llenos sometidos a excitaciones longitudinales y laterales simultáneas". European Journal of Mechanics B . 53 : 251– 263. Bibcode : 2015EuJMB..53..251K . doi : 10.1016/j.euromechflu.2015.06.001 .
  5. Kolaei, Amir; Rakheja, Subhash; Richard, Marc J. (2014-01-06). "Rango de aplicabilidad de la teoría lineal de oscilación de fluidos para predecir la oscilación lateral transitoria y la estabilidad de balanceo de vehículos tanque". Journal of Sound and Vibration . 333 (1): 263– 282. Bibcode : 2014JSV...333..263K . doi : 10.1016/j.jsv.2013.09.002 .
  6. Kolaei, Amir; Rakheja, Subhash; Richard, Marc J. (2014-07-01). "Efectos de la sección transversal del tanque en las cargas dinámicas de chapoteo de fluido y la estabilidad de balanceo de un camión cisterna parcialmente lleno". European Journal of Mechanics B . 46 : 46– 58. Bibcode : 2014EuJMB..46...46K . doi : 10.1016/j.euromechflu.2014.01.008 .
  7. ^ Manuel, Frey (2001). Behandlung von Strömungsproblemen in Raketendüsen bei Überexpansion (Tesis) (en alemán). doi : 10.18419/opus-3650 .
  8. Daub, Dennis; Esser, Burkard; Gülhan, Ali (abril de 2020). "Experimentos sobre la 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". The American Journal of Nursing . 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. (2009-09-15). "Modelado de interacción fluido-estructura del flujo sanguíneo y aneurisma cerebral: Importancia de las formas de la arteria y el aneurisma". Métodos computacionales 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 (2019-07-01). "Riesgo de ruptura del aneurisma cerebral en relación con la lesión cerebral traumática utilizando un modelo de interacción fluido-estructura específico del paciente". Computer Methods and Programs in Biomedicine . 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: una introducción al acoplamiento de elementos finitos . Wiley ( ISBN) 978-1-119-95227-5)
  15. "Copia archivada" (PDF) . Archivado del original (PDF) el 31-10-2014 . Recuperado el 28-11-2014 .{{cite web}}: CS1 mantenimiento: copia archivada como título ( enlace )
  16. Zhang, Chi; Rezavand, Massoud; Hu, Xiangyu (2021-03-15). "Un método SPH multirresolución para interacciones fluido-estructura". Journal of Computational Physics . 429 110028. arXiv : 1911.13255 . Bibcode : 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 computacionales 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". International Journal for Numerical Methods in Engineering . 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 . Lecture Notes in Computational Science and Engineering. Vol. Interacción fluido-estructura: modelado, simulación, optimización. Springer-Verlag . pp. 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". Computers and Structures . 81 ( 8–11 ): 805–812 . CiteSeerX 10.1.1.487.5577 . doi : 10.1016/S0045-7949(02)00409-1 . 
  21. 1 2 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 precondicionado por subiteración para la interacción fluido-estructura". Computer Methods in Applied Mechanics and Engineering . 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: Mathematical Modelling and Numerical Analysis . 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". 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. 1 2 P. Causin; J.-F. Gerbeau; F. Nobile (2005). "Efecto de masa añadida en el diseño de algoritmos particionados para problemas 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 una interfaz fluido-sólido curva 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 añadida artificial en algoritmos de interacción fluido-estructura escalonada secuencial . 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 añadida artificial en el acoplamiento secuencial escalonado de estructuras no lineales y flujos viscosos incompresibles" . Computer Methods in Applied Mechanics and Engineering . 196 (7): 1278– 1293. Bibcode : 2007CMAME.196.1278F . doi : 10.1016/j.cma.2006.09.002 .

Lecturas adicionales

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 Introduce el tema de las 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.