stringtranslate.com

Dinámica de fluidos computacional

La dinámica de fluidos computacional ( CFD ) es una rama de la mecánica de fluidos que utiliza análisis numéricos y estructuras de datos para analizar y resolver problemas que involucran flujos de fluidos . Se utilizan computadoras para realizar los cálculos necesarios para simular el flujo libre del fluido y la interacción del fluido ( líquidos y gases ) con superficies definidas por condiciones de contorno . Con supercomputadoras de alta velocidad , se pueden lograr mejores soluciones, y a menudo se requieren para resolver los problemas más grandes y complejos. La investigación en curso produce software que mejora la precisión y la velocidad de escenarios de simulación complejos, como flujos transónicos o turbulentos . La validación inicial de dicho software generalmente se realiza utilizando aparatos experimentales como túneles de viento . Además, el análisis analítico o empírico realizado previamente de un problema particular se puede utilizar para la comparación. A menudo se realiza una validación final utilizando pruebas a escala real, como pruebas de vuelo .

La CFD se aplica a una amplia gama de problemas de investigación e ingeniería en muchos campos de estudio e industrias, incluidos la aerodinámica y el análisis aeroespacial, la hipersónica , la simulación meteorológica , las ciencias naturales y la ingeniería ambiental , el diseño y análisis de sistemas industriales, la ingeniería biológica , los flujos de fluidos y la transferencia de calor , el análisis de motores y combustión y los efectos visuales para películas y juegos.

Antecedentes e historia

Una simulación por computadora del flujo de aire de alta velocidad alrededor del transbordador espacial durante el reingreso
Una simulación del vehículo estatorreactor Hyper-X en funcionamiento a Mach -7

La base fundamental de casi todos los problemas de CFD son las ecuaciones de Navier-Stokes , que definen muchos flujos de fluidos monofásicos (gas o líquido, pero no ambos). Estas ecuaciones se pueden simplificar eliminando los términos que describen acciones viscosas para obtener las ecuaciones de Euler . Si se simplifica aún más, al eliminar los términos que describen la vorticidad se obtienen las ecuaciones de potencial completas . Finalmente, para pequeñas perturbaciones en flujos subsónicos y supersónicos (no transónicos o hipersónicos ), estas ecuaciones se pueden linealizar para obtener las ecuaciones de potencial linealizadas.

Históricamente, los primeros métodos se desarrollaron para resolver las ecuaciones de potencial linealizadas. Los métodos bidimensionales (2D), que utilizan transformaciones conformes del flujo alrededor de un cilindro al flujo alrededor de un perfil aerodinámico , se desarrollaron en la década de 1930. [1] [2]

Uno de los primeros tipos de cálculos que se asemejan a la CFD moderna son los de Lewis Fry Richardson , en el sentido de que estos cálculos utilizaban diferencias finitas y dividían el espacio físico en celdas. Aunque fracasaron estrepitosamente, estos cálculos, junto con el libro de Richardson Weather Prediction by Numerical Process [3] , sentaron las bases para la CFD moderna y la meteorología numérica. De hecho, los primeros cálculos de CFD durante la década de 1940 utilizando ENIAC utilizaban métodos cercanos a los del libro de Richardson de 1922. [4]

La potencia informática disponible impulsó el desarrollo de métodos tridimensionales . Probablemente el primer trabajo que utilizó ordenadores para modelar el flujo de fluidos, tal como lo rigen las ecuaciones de Navier-Stokes, se realizó en el Laboratorio Nacional de Los Álamos , en el grupo T3. [5] [6] Este grupo fue dirigido por Francis H. Harlow , considerado ampliamente uno de los pioneros de la CFD. Desde 1957 hasta finales de los años 1960, este grupo desarrolló una variedad de métodos numéricos para simular flujos de fluidos bidimensionales transitorios, como el método de partículas en celda , [7] el método de fluido en celda, [8] el método de función de corriente de vorticidad, [9] y el método de marcador y celda . [10] El método de función de corriente de vorticidad de Fromm para flujo incompresible transitorio 2D fue el primer tratamiento de flujos incompresibles fuertemente contorsionados en el mundo.

El primer artículo con modelo tridimensional fue publicado por John Hess y AMO Smith de Douglas Aircraft en 1967. [11] Este método discretizó la superficie de la geometría con paneles, dando lugar a que esta clase de programas se llamaran Métodos de Panel. Su método en sí mismo fue simplificado, en el sentido de que no incluía flujos de sustentación y, por lo tanto, se aplicó principalmente a cascos de barcos y fuselajes de aeronaves. El primer Código de Panel de sustentación (A230) fue descrito en un artículo escrito por Paul Rubbert y Gary Saaris de Boeing Aircraft en 1968. [12] Con el tiempo, se desarrollaron Códigos de Panel tridimensionales más avanzados en Boeing (PANAIR, A502), [13] Lockheed (Quadpan), [14] Douglas (HESS), [15] McDonnell Aircraft (MACAERO), [16] NASA (PMARC) [17] y Analytical Methods (WBAERO, [18] USAERO [19] y VSAERO [20] [21] ). Algunos (PANAIR, HESS y MACAERO) eran códigos de orden superior, que utilizaban distribuciones de orden superior de singularidades de superficie, mientras que otros (Quadpan, PMARC, USAERO y VSAERO) utilizaban singularidades individuales en cada panel de superficie. La ventaja de los códigos de orden inferior era que se ejecutaban mucho más rápido en las computadoras de la época. Hoy, VSAERO ha crecido hasta convertirse en un código de orden múltiple y es el programa más utilizado de esta clase. Se ha utilizado en el desarrollo de muchos submarinos , barcos de superficie , automóviles , helicópteros , aviones y, más recientemente, turbinas eólicas . Su código hermano, USAERO, es un método de panel inestable que también se ha utilizado para modelar cosas como trenes de alta velocidad y yates de carreras . El código PMARC de la NASA de una versión temprana de VSAERO y un derivado de PMARC, llamado CMARC, [22] también está disponible comercialmente.

En el ámbito bidimensional, se han desarrollado varios códigos de panel para el análisis y diseño de perfiles aerodinámicos. Los códigos suelen incluir un análisis de capa límite , de modo que se puedan modelar los efectos viscosos. Richard Eppler  [de] desarrolló el código PROFILE, en parte con financiación de la NASA, que estuvo disponible a principios de los años 1980. [23] A éste le siguió pronto el código XFOIL de Mark Drela . [24] Tanto PROFILE como XFOIL incorporan códigos de panel bidimensionales, con códigos de capa límite acoplados para el trabajo de análisis de perfiles aerodinámicos. PROFILE utiliza un método de transformación conforme para el diseño de perfiles aerodinámicos inversos, mientras que XFOIL tiene tanto una transformación conforme como un método de panel inverso para el diseño de perfiles aerodinámicos.

Un paso intermedio entre los códigos de panel y los códigos de potencial completo fueron los códigos que utilizaban las ecuaciones de pequeñas perturbaciones transónicas. En particular, el código tridimensional WIBCO [25] , desarrollado por Charlie Boppe de Grumman Aircraft a principios de los años 1980, ha sido ampliamente utilizado.

Una simulación de la nave espacial SpaceX durante el reingreso

Los desarrolladores recurrieron a los códigos Full Potential, ya que los métodos de panel no podían calcular el flujo no lineal presente a velocidades transónicas . La primera descripción de un medio para utilizar las ecuaciones Full Potential fue publicada por Earll Murman y Julian Cole de Boeing en 1970. [26] Frances Bauer, Paul Garabedian y David Korn del Courant Institute de la Universidad de Nueva York (NYU) escribieron una serie de códigos de perfil aerodinámico Full Potential bidimensionales que se utilizaron ampliamente, siendo el más importante el denominado Programa H. ​​[27] Bob Melnik y su grupo en Grumman Aerospace desarrollaron una versión posterior del Programa H con el nombre de Grumfoil. [28] Antony Jameson , originalmente en Grumman Aircraft y el Courant Institute de NYU, trabajó con David Caughey para desarrollar el importante código Full Potential tridimensional FLO22 [29] en 1975. Muchos códigos Full Potential surgieron después de esto, culminando en el código Tranair (A633) de Boeing, [30] que todavía se utiliza mucho.

El siguiente paso fueron las ecuaciones de Euler, que prometían proporcionar soluciones más precisas de los flujos transónicos. La metodología utilizada por Jameson en su código tridimensional FLO57 [31] (1981) fue utilizada por otros para producir programas como el programa TEAM de Lockheed [32] y el programa MGAERO de IAI/Analytical Methods. [33] MGAERO es único por ser un código de malla cartesiana estructurada , mientras que la mayoría de los demás códigos de este tipo utilizan cuadrículas estructuradas ajustadas al cuerpo (con la excepción del exitoso código CART3D de la NASA, [34] el código SPLITFLOW de Lockheed [35] y el NASCART-GT de Georgia Tech ). [36] Antony Jameson también desarrolló el código tridimensional AIRPLANE [37] que hacía uso de cuadrículas tetraédricas no estructuradas.

En el ámbito bidimensional, Mark Drela y Michael Giles, entonces estudiantes de posgrado en el MIT, desarrollaron el programa Euler ISES [38] (en realidad, un conjunto de programas) para el diseño y análisis de perfiles aerodinámicos. Este código estuvo disponible por primera vez en 1986 y se ha desarrollado aún más para diseñar, analizar y optimizar perfiles aerodinámicos de uno o varios elementos, como el programa MSES. [39] MSES se utiliza ampliamente en todo el mundo. Un derivado de MSES, para el diseño y análisis de perfiles aerodinámicos en cascada, es MISES, [40] desarrollado por Harold Youngren mientras era estudiante de posgrado en el MIT.

Las ecuaciones de Navier-Stokes fueron el objetivo final del desarrollo. Surgieron los primeros códigos bidimensionales, como el código ARC2D de la NASA Ames. Se desarrollaron varios códigos tridimensionales (ARC3D, OVERFLOW y CFL3D son tres contribuciones exitosas de la NASA), lo que dio lugar a numerosos paquetes comerciales.

Recientemente, los métodos de CFD han ganado terreno para modelar el comportamiento del flujo de materiales granulares en diversos procesos químicos de ingeniería. Este enfoque ha surgido como una alternativa rentable, que ofrece una comprensión matizada de fenómenos de flujo complejos y, al mismo tiempo, minimiza los gastos asociados con los métodos experimentales tradicionales. [41] [42]

Jerarquía de ecuaciones de flujo de fluidos

La CFD puede considerarse un grupo de metodologías computacionales (que se analizan a continuación) que se utilizan para resolver ecuaciones que rigen el flujo de fluidos. En la aplicación de la CFD, un paso fundamental es decidir qué conjunto de suposiciones físicas y ecuaciones relacionadas se deben utilizar para el problema en cuestión. [43] Para ilustrar este paso, a continuación se resumen las suposiciones físicas/simplificaciones adoptadas en las ecuaciones de un flujo monofásico (véase flujo multifásico y flujo bifásico ), monoespecie (es decir, que consta de una especie química), no reactivo y (a menos que se indique lo contrario) compresible. Se descuida la radiación térmica y se consideran las fuerzas del cuerpo debidas a la gravedad (a menos que se indique lo contrario). Además, para este tipo de flujo, la siguiente discusión destaca la jerarquía de ecuaciones de flujo resueltas con CFD. Tenga en cuenta que algunas de las siguientes ecuaciones podrían derivarse de más de una forma.

Metodología

En todos estos enfoques se sigue el mismo procedimiento básico.

Métodos de discretización

La estabilidad de la discretización seleccionada se establece generalmente de forma numérica en lugar de analítica, como en el caso de los problemas lineales simples. También se debe tener especial cuidado para garantizar que la discretización maneje sin problemas las soluciones discontinuas. Las ecuaciones de Euler y las ecuaciones de Navier-Stokes admiten choques y superficies de contacto.

Algunos de los métodos de discretización que se utilizan son:

Método de volumen finito

El método de volumen finito (FVM) es un enfoque común utilizado en códigos CFD, ya que tiene una ventaja en el uso de memoria y la velocidad de solución, especialmente para problemas grandes, flujos turbulentos con alto número de Reynolds y flujos dominados por términos fuente (como la combustión). [55]

En el método de volumen finito, las ecuaciones diferenciales parciales que rigen (normalmente las ecuaciones de Navier-Stokes, las ecuaciones de conservación de masa y energía y las ecuaciones de turbulencia) se reformulan en una forma conservativa y luego se resuelven sobre volúmenes de control discretos. Esta discretización garantiza la conservación de flujos a través de un volumen de control particular. La ecuación de volumen finito produce ecuaciones que rigen en la forma:

donde es el vector de variables conservadas, es el vector de flujos (ver ecuaciones de Euler o ecuaciones de Navier-Stokes ), es el volumen del elemento de volumen de control, y es el área de superficie del elemento de volumen de control.

Método de elementos finitos

El método de elementos finitos (MEF) se utiliza en el análisis estructural de sólidos, pero también es aplicable a fluidos. Sin embargo, la formulación MEF requiere un cuidado especial para garantizar una solución conservadora. La formulación MEF se ha adaptado para su uso con ecuaciones que rigen la dinámica de fluidos. [56] [57] Aunque el MEF debe formularse con cuidado para que sea conservador, es mucho más estable que el enfoque de volumen finito. [58] Sin embargo, el MEF puede requerir más memoria y tiene tiempos de solución más lentos que el MFF. [59]

En este método se forma una ecuación residual ponderada:

donde es el residuo de la ecuación en el vértice de un elemento , es la ecuación de conservación expresada sobre una base de elementos, es el factor de peso y es el volumen del elemento.

Método de diferencias finitas

El método de diferencias finitas (FDM) tiene importancia histórica [57] y es fácil de programar. Actualmente, solo se utiliza en unos pocos códigos especializados que manejan geometría compleja con alta precisión y eficiencia mediante el uso de límites integrados o cuadrículas superpuestas (con la solución interpolada en cada cuadrícula). [ cita requerida ]

donde es el vector de variables conservadas, y , , y son los flujos en las direcciones , , y respectivamente.

Método de elementos espectrales

El método de elementos espectrales es un método de tipo de elementos finitos. Requiere que el problema matemático (la ecuación diferencial parcial) se formule en una formulación débil. Esto se hace típicamente multiplicando la ecuación diferencial por una función de prueba arbitraria e integrando sobre todo el dominio. Puramente matemáticamente, las funciones de prueba son completamente arbitrarias: pertenecen a un espacio de funciones de dimensión infinita. Claramente, un espacio de funciones de dimensión infinita no se puede representar en una malla de elementos espectrales discretos; aquí es donde comienza la discretización de elementos espectrales. Lo más crucial es la elección de las funciones de interpolación y prueba. En un método FEM estándar de orden bajo en 2D, para elementos cuadriláteros la elección más típica es la función de prueba o interpolación bilineal de la forma . Sin embargo, en un método de elementos espectrales, las funciones de interpolación y prueba se eligen para que sean polinomios de un orden muy alto (típicamente, por ejemplo, del décimo orden en aplicaciones CFD). Esto garantiza la convergencia rápida del método. Además, se deben utilizar procedimientos de integración muy eficientes, ya que el número de integraciones que se deben realizar en los códigos numéricos es grande. Por lo tanto, se emplean cuadraturas de integración de Gauss de alto orden, ya que logran la mayor precisión con el menor número de cálculos a realizar. En este momento, existen algunos códigos académicos de CFD basados ​​en el método de elementos espectrales y algunos más están actualmente en desarrollo, ya que surgen nuevos esquemas de pasos de tiempo en el mundo científico.

Método de Boltzmann en red

El método de Boltzmann en red (LBM), con su representación cinética simplificada en una red, proporciona una descripción computacionalmente eficiente de la hidrodinámica. A diferencia de los métodos CFD tradicionales, que resuelven numéricamente las ecuaciones de conservación de propiedades macroscópicas (es decir, masa, momento y energía), el LBM modela el fluido que consiste en partículas ficticias, y dichas partículas realizan procesos consecutivos de propagación y colisión sobre una malla de red discreta. En este método, se trabaja con la versión discreta en el espacio y el tiempo de la ecuación de evolución cinética en la forma Boltzmann Bhatnagar-Gross-Krook (BGK) .

Método del vórtice

El método de vórtices, también llamado método de partículas de vórtices lagrangianos, es una técnica sin malla para la simulación de flujos turbulentos incompresibles. En él, la vorticidad se discretiza sobre partículas lagrangianas , a las que se denominan vórtices, vórtices o partículas de vórtice. [60] Los métodos de vórtices se desarrollaron como una metodología sin malla que no estaría limitada por los efectos de suavizado fundamentales asociados con los métodos basados ​​en mallas. Sin embargo, para ser prácticos, los métodos de vórtices requieren medios para calcular rápidamente las velocidades a partir de los elementos de vórtice; en otras palabras, requieren la solución a una forma particular del problema de N cuerpos (en el que el movimiento de N objetos está ligado a sus influencias mutuas). Este avance se produjo en la década de 1980 con el desarrollo de los algoritmos de Barnes-Hut y del método multipolar rápido (FMM). Estos allanaron el camino para el cálculo práctico de las velocidades a partir de los elementos de vórtice.

El software basado en el método de vórtices ofrece un nuevo medio para resolver problemas complejos de dinámica de fluidos con una mínima intervención del usuario. [ cita requerida ] Todo lo que se requiere es la especificación de la geometría del problema y el establecimiento de condiciones iniciales y de contorno. Entre las ventajas significativas de esta tecnología moderna;

Método de elementos de contorno

En el método de elementos límite, el límite ocupado por el fluido se divide en una malla de superficie.

Esquemas de discretización de alta resolución

Los esquemas de alta resolución se utilizan cuando hay choques o discontinuidades. Para captar cambios bruscos en la solución se requiere el uso de esquemas numéricos de segundo orden o de orden superior que no introduzcan oscilaciones espurias. Esto suele requerir la aplicación de limitadores de flujo para garantizar que la solución disminuya la variación total . [ cita requerida ]

Modelos de turbulencia

En el modelado computacional de flujos turbulentos, un objetivo común es obtener un modelo que pueda predecir magnitudes de interés, como la velocidad del fluido, para su uso en diseños de ingeniería del sistema que se está modelando. En el caso de los flujos turbulentos, la variedad de escalas de longitud y la complejidad de los fenómenos involucrados en la turbulencia hacen que la mayoría de los enfoques de modelado sean prohibitivamente costosos; la resolución requerida para resolver todas las escalas involucradas en la turbulencia está más allá de lo que es computacionalmente posible. El enfoque principal en tales casos es crear modelos numéricos para aproximarse a los fenómenos no resueltos. Esta sección enumera algunos modelos computacionales de uso común para flujos turbulentos.

Los modelos de turbulencia se pueden clasificar en función del gasto computacional, que corresponde al rango de escalas que se modelan en comparación con las que se resuelven (cuanto más escalas turbulentas se resuelvan, más precisa será la resolución de la simulación y, por lo tanto, mayor será el costo computacional). Si no se modela la mayoría o la totalidad de las escalas turbulentas, el costo computacional es muy bajo, pero la contrapartida se presenta en forma de una menor precisión.

Además de la amplia gama de escalas de longitud y tiempo y el costo computacional asociado, las ecuaciones que rigen la dinámica de fluidos contienen un término de convección no lineal y un término de gradiente de presión no lineal y no local. Estas ecuaciones no lineales deben resolverse numéricamente con las condiciones iniciales y de contorno adecuadas.

Navier-Stokes promediado por Reynolds

Aerodinámica externa del modelo DrivAer, calculada utilizando URANS (arriba) y DDES (abajo)
Simulación del paquete aerodinámico de un Porsche Cayman (987.2)

Las ecuaciones de Navier-Stokes promediadas por Reynolds (RANS) son el enfoque más antiguo para el modelado de turbulencia. Se resuelve una versión de conjunto de las ecuaciones gobernantes, que introduce nuevas tensiones aparentes conocidas como tensiones de Reynolds . Esto agrega un tensor de incógnitas de segundo orden para el cual varios modelos pueden proporcionar diferentes niveles de cierre. Es un error común pensar que las ecuaciones RANS no se aplican a flujos con un flujo medio variable en el tiempo porque estas ecuaciones están "promediadas en el tiempo". De hecho, los flujos estadísticamente inestables (o no estacionarios) también pueden tratarse. A esto a veces se lo denomina URANS. No hay nada inherente al promedio de Reynolds que lo impida, pero los modelos de turbulencia utilizados para cerrar las ecuaciones son válidos solo mientras el tiempo durante el cual ocurren estos cambios en la media sea grande en comparación con las escalas de tiempo del movimiento turbulento que contiene la mayor parte de la energía.

Los modelos RANS se pueden dividir en dos enfoques amplios:

Hipótesis de Boussinesq
Este método implica el uso de una ecuación algebraica para las tensiones de Reynolds que incluyen la determinación de la viscosidad turbulenta y, dependiendo del nivel de sofisticación del modelo, la resolución de ecuaciones de transporte para determinar la energía cinética turbulenta y la disipación. Los modelos incluyen k-ε ( Launder y Spalding ), [61] el modelo de longitud de mezcla ( Prandtl ), [62] y el modelo de ecuación cero (Cebeci y Smith ). [62] Los modelos disponibles en este enfoque a menudo se mencionan por el número de ecuaciones de transporte asociadas con el método. Por ejemplo, el modelo de longitud de mezcla es un modelo de "ecuación cero" porque no se resuelven ecuaciones de transporte; el es un modelo de "dos ecuaciones" porque se resuelven dos ecuaciones de transporte (una para y una para ).
Modelo de estrés de Reynolds (MSR)
Este enfoque intenta resolver ecuaciones de transporte para las tensiones de Reynolds. Esto implica la introducción de varias ecuaciones de transporte para todas las tensiones de Reynolds y, por lo tanto, este enfoque es mucho más costoso en cuanto a esfuerzo de CPU. [ cita requerida ]

Simulación de grandes remolinos

Representación volumétrica de una llama en espiral no premezclada simulada por LES

La simulación de grandes remolinos (LES, por sus siglas en inglés) es una técnica en la que se eliminan las escalas más pequeñas del flujo mediante una operación de filtrado y se modela su efecto utilizando modelos de escala de submalla. Esto permite resolver las escalas más grandes e importantes de la turbulencia, al tiempo que reduce en gran medida el costo computacional que implican las escalas más pequeñas. Este método requiere mayores recursos computacionales que los métodos RANS, pero es mucho más económico que el DNS.

Simulación de remolino desprendido

Las simulaciones de remolinos desprendidos (DES) son una modificación de un modelo RANS en el que el modelo cambia a una formulación de escala de subcuadrícula en regiones lo suficientemente finas para los cálculos LES. Las regiones cercanas a los límites sólidos y donde la escala de longitud turbulenta es menor que la dimensión máxima de la cuadrícula se asignan al modo de solución RANS. Como la escala de longitud turbulenta excede la dimensión de la cuadrícula, las regiones se resuelven utilizando el modo LES. Por lo tanto, la resolución de cuadrícula para DES no es tan exigente como LES puro, lo que reduce considerablemente el costo del cálculo. Aunque DES se formuló inicialmente para el modelo Spalart-Allmaras (Philippe R. Spalart et al., 1997), se puede implementar con otros modelos RANS (Strelets, 2001), modificando adecuadamente la escala de longitud que está explícita o implícitamente involucrada en el modelo RANS. Entonces, mientras que DES basado en el modelo Spalart-Allmaras actúa como LES con un modelo de pared, DES basado en otros modelos (como modelos de dos ecuaciones) se comporta como un modelo híbrido RANS-LES. La generación de cuadrículas es más complicada que en el caso de un RANS o LES simple debido al cambio entre RANS y LES. DES es un enfoque no zonal y proporciona un único campo de velocidad uniforme en las regiones RANS y LES de las soluciones.

Simulación IDDES del BMW de Karel Motorsports. Se trata de un tipo de simulación DES realizada en OpenFOAM. El gráfico representa el coeficiente de presión.

Simulación numérica directa

La simulación numérica directa (DNS) resuelve todo el rango de escalas de longitud turbulentas. Esto margina el efecto de los modelos, pero es extremadamente costoso. El costo computacional es proporcional a . [63] La DNS es intratable para flujos con geometrías complejas o configuraciones de flujo.

Simulación de vórtice coherente

El enfoque de simulación de vórtices coherentes descompone el campo de flujo turbulento en una parte coherente, que consiste en un movimiento vorticial organizado, y la parte incoherente, que es el flujo de fondo aleatorio. [64] Esta descomposición se realiza utilizando filtrado wavelet . El enfoque tiene mucho en común con LES, ya que utiliza la descomposición y resuelve solo la parte filtrada, pero se diferencia en que no utiliza un filtro lineal de paso bajo. En cambio, la operación de filtrado se basa en wavelets, y el filtro se puede adaptar a medida que evoluciona el campo de flujo. Farge y Schneider probaron el método CVS con dos configuraciones de flujo y demostraron que la parte coherente del flujo exhibía el espectro de energía exhibido por el flujo total, y correspondía a estructuras coherentes ( tubos de vórtice ), mientras que las partes incoherentes del flujo componían ruido de fondo homogéneo, que no exhibía estructuras organizadas. Goldstein y Vasilyev [65] aplicaron el modelo FDV a una simulación de grandes remolinos, pero no asumieron que el filtro wavelet eliminara todos los movimientos coherentes de las escalas del subfiltro. Al emplear tanto el filtrado LES como el CVS, demostraron que la disipación de SFS estaba dominada por la porción coherente del campo de flujo de SFS.

Métodos PDF

Los métodos de función de densidad de probabilidad (PDF) para la turbulencia, introducidos por primera vez por Lundgren , [66] se basan en el seguimiento de la PDF de un punto de la velocidad, , que da la probabilidad de que la velocidad en el punto esté entre y . Este enfoque es análogo a la teoría cinética de los gases , en la que las propiedades macroscópicas de un gas se describen mediante una gran cantidad de partículas. Los métodos PDF son únicos en el sentido de que se pueden aplicar en el marco de varios modelos de turbulencia diferentes; las principales diferencias ocurren en la forma de la ecuación de transporte PDF. Por ejemplo, en el contexto de la simulación de grandes remolinos , la PDF se convierte en la PDF filtrada. [67] Los métodos PDF también se pueden utilizar para describir reacciones químicas, [68] [69] y son particularmente útiles para simular flujos que reaccionan químicamente porque el término de fuente química está cerrado y no requiere un modelo. La PDF se rastrea comúnmente utilizando métodos de partículas lagrangianas; cuando se combina con la simulación de grandes remolinos, esto conduce a una ecuación de Langevin para la evolución de partículas de subfiltro.

Método de confinamiento de vorticidad

El método de confinamiento de vorticidad (VC) es una técnica euleriana utilizada en la simulación de estelas turbulentas. Utiliza un enfoque similar al de las ondas solitarias para producir una solución estable sin dispersión numérica. El VC puede capturar las características de pequeña escala en tan solo 2 celdas de la cuadrícula. Dentro de estas características, se resuelve una ecuación diferencial no lineal en lugar de una ecuación diferencial finita . El VC es similar a los métodos de captura de choques , donde se satisfacen las leyes de conservación, de modo que las cantidades integrales esenciales se calculan con precisión.

Modelo de remolino lineal

El modelo de remolino lineal es una técnica utilizada para simular la mezcla convectiva que tiene lugar en el flujo turbulento. [70] Específicamente, proporciona una forma matemática de describir las interacciones de una variable escalar dentro del campo de flujo vectorial. Se utiliza principalmente en representaciones unidimensionales de flujo turbulento, ya que se puede aplicar en una amplia gama de escalas de longitud y números de Reynolds. Este modelo se utiliza generalmente como un bloque de construcción para representaciones de flujo más complicadas, ya que proporciona predicciones de alta resolución que se mantienen en una amplia gama de condiciones de flujo.

Flujo de dos fases

Simulación de una horda de burbujas utilizando el método del volumen de fluido

El modelado del flujo bifásico aún se encuentra en desarrollo. Se han propuesto diferentes métodos, incluido el método del volumen del fluido , el método de ajuste de nivel y el seguimiento frontal. [71] [72] Estos métodos a menudo implican un equilibrio entre mantener una interfaz definida o conservar la masa [ ¿según quién? ] . Esto es crucial ya que la evaluación de la densidad, la viscosidad y la tensión superficial se basa en los valores promediados sobre la interfaz. [ cita requerida ]

Algoritmos de solución

La discretización en el espacio produce un sistema de ecuaciones diferenciales ordinarias para problemas no estacionarios y ecuaciones algebraicas para problemas estacionarios. Generalmente se utilizan métodos implícitos o semi-implícitos para integrar las ecuaciones diferenciales ordinarias, produciendo un sistema de ecuaciones algebraicas (generalmente) no lineales. La aplicación de una iteración de Newton [ ancla rota ] o de Picard produce un sistema de ecuaciones lineales que es asimétrico en presencia de advección e indefinido en presencia de incompresibilidad. Dichos sistemas, particularmente en 3D, son con frecuencia demasiado grandes para los solucionadores directos, por lo que se utilizan métodos iterativos, ya sean métodos estacionarios como la sobrerelajación sucesiva o los métodos de subespacio de Krylov . Los métodos de Krylov como GMRES , que normalmente se utilizan con preacondicionamiento , funcionan minimizando el residuo sobre subespacios sucesivos generados por el operador preacondicionado.

La multimalla tiene la ventaja de un rendimiento asintóticamente óptimo en muchos problemas. Los solucionadores y preacondicionadores tradicionales [ ¿según quién? ] son ​​eficaces para reducir los componentes de alta frecuencia del residuo, pero los componentes de baja frecuencia suelen requerir muchas iteraciones para su reducción. Al operar en múltiples escalas, la multimalla reduce todos los componentes del residuo en factores similares, lo que da lugar a un número de iteraciones independiente de la malla. [ cita requerida ]

En el caso de sistemas indefinidos, los preacondicionadores como la factorización LU incompleta , el Schwarz aditivo y la multimalla tienen un rendimiento deficiente o fallan por completo, por lo que se debe utilizar la estructura del problema para un preacondicionamiento efectivo. [73] Los métodos comúnmente utilizados en CFD son los algoritmos SIMPLE y Uzawa que exhiben tasas de convergencia dependientes de la malla, pero los avances recientes basados ​​en la factorización LU en bloque combinada con la multimalla para los sistemas definidos resultantes han llevado a preacondicionadores que ofrecen tasas de convergencia independientes de la malla. [74]

Aerodinámica inestable

La CFD tuvo un gran avance a finales de los años 70 con la introducción de LTRAN2, un código 2-D para modelar perfiles aerodinámicos oscilantes basado en la teoría de pequeñas perturbaciones transónicas de Ballhaus y asociados. [75] Utiliza un algoritmo de conmutación Murman-Cole para modelar las ondas de choque en movimiento. [26] Más tarde se amplió a 3-D con el uso de un esquema de diferencia rotada por AFWAL/Boeing que dio como resultado LTRAN3. [76] [77]

Ingeniería biomédica

Simulación del flujo sanguíneo en una aorta humana

Las investigaciones de CFD se utilizan para aclarar las características del flujo aórtico en detalles que están más allá de las capacidades de las mediciones experimentales. Para analizar estas condiciones, se extraen modelos CAD del sistema vascular humano empleando técnicas de imágenes modernas como la resonancia magnética o la tomografía computarizada . Se reconstruye un modelo 3D a partir de estos datos y se puede calcular el flujo del fluido. Se deben tener en cuenta las propiedades de la sangre, como la densidad y la viscosidad, y las condiciones límite realistas (por ejemplo, la presión sistémica). Por lo tanto, es posible analizar y optimizar el flujo en el sistema cardiovascular para diferentes aplicaciones. [78]

CPU versus GPU

Tradicionalmente, las simulaciones CFD se realizan en CPU. [79]

En una tendencia más reciente, las simulaciones también se realizan en GPU. Estas suelen contener procesadores más lentos pero con más procesadores. Para los algoritmos CFD que presentan un buen rendimiento de paralelismo (es decir, una buena aceleración al agregar más núcleos), esto puede reducir en gran medida los tiempos de simulación. Los métodos de partículas implícitas en fluidos [80] y de Boltzmann en red [81] son ​​ejemplos típicos de códigos que escalan bien en GPU.

Véase también

Referencias

  1. ^ Milne-Thomson, Louis Melville (1973). Aerodinámica teórica . Courier Corporation. ISBN 978-0-486-61980-4.[ página necesaria ]
  2. ^ McMurtry, Patrick A.; Gansauge, Todd C.; Kerstein, Alan R.; Krueger, Steven K. (abril de 1993). "Simulaciones de remolinos lineales de mezcla en un flujo turbulento homogéneo". Física de fluidos A: dinámica de fluidos . 5 (4): 1023–1034. Código Bibliográfico :1993PhFlA...5.1023M. doi :10.1063/1.858667.
  3. ^ Richardson, LF; Chapman, S. (1965). Predicción meteorológica mediante procesos numéricos . Publicaciones de Dover.
  4. ^ Hunt, JCR (enero de 1998). "Lewis Fry Richardson y sus contribuciones a las matemáticas, la meteorología y los modelos de conflicto". Revista anual de mecánica de fluidos . 30 (1): xiii–xxxvi. Código Bibliográfico :1998AnRFM..30D..13H. doi :10.1146/annurev.fluid.30.1.0.
  5. ^ "El legado del Grupo T-3" . Consultado el 13 de marzo de 2013 .
  6. ^ Harlow, Francis H. (abril de 2004). "Dinámica de fluidos en el Laboratorio Nacional de Los Álamos Grupo T-3". Revista de Física Computacional . 195 (2): 414–433. Bibcode :2004JCoPh.195..414H. doi :10.1016/j.jcp.2003.09.031.
  7. ^ Harlow, Francis Harvey; Evans, Martha; Richtmyer, Robert D. (1955). Un método de cálculo de máquina para problemas hidrodinámicos . Laboratorio Científico de Los Álamos de la Universidad de California. hdl :2027/mdp.39015095283399. OCLC  1288309947.[ página necesaria ]
  8. ^ Gentry, Richard A; Martin, Robert E; Daly, Bart J (agosto de 1966). "Un método de diferenciación euleriano para problemas de flujo compresible no estacionario". Journal of Computational Physics . 1 (1): 87–118. Bibcode :1966JCoPh...1...87G. doi :10.1016/0021-9991(66)90014-3.
  9. ^ Fromm, Jacob E.; Harlow, Francis H. (julio de 1963). "Solución numérica del problema del desarrollo de calles de vórtices". Física de fluidos . 6 (7): 975–982. Bibcode :1963PhFl....6..975F. doi :10.1063/1.1706854.
  10. ^ Harlow, Francis H.; Welch, J. Eddie (diciembre de 1965). "Cálculo numérico del flujo viscoso incompresible dependiente del tiempo de un fluido con superficie libre". Física de fluidos . 8 (12): 2182–2189. Bibcode :1965PhFl....8.2182H. doi :10.1063/1.1761178.
  11. ^ Hess, JL; Smith, AMO (1967). "Cálculo del flujo potencial sobre cuerpos arbitrarios". Progreso en las ciencias aeroespaciales . 8 : 1–138. Bibcode :1967PrAeS...8....1H. doi :10.1016/0376-0421(67)90003-6.
  12. ^ Rubbert, P.; Saaris, G. (1972). "Revisión y evaluación de un método computacional de flujo de potencial de sustentación tridimensional para configuraciones arbitrarias". 10.ª Reunión de Ciencias Aeroespaciales . doi :10.2514/6.1972-188.
  13. ^ Carmichael, R.; Erickson, L. (1981). "PAN AIR - Un método de panel de orden superior para predecir flujos potenciales lineales subsónicos o supersónicos sobre configuraciones arbitrarias". 14.ª Conferencia sobre dinámica de fluidos y plasma . doi :10.2514/6.1981-1255.
  14. ^ Youngren, H.; Bouchard, E.; Coopersmith, R.; Miranda, L. (1983). "Comparación de formulaciones de métodos de panel y su influencia en el desarrollo de QUADPAN, un método avanzado de bajo orden". Conferencia de Aerodinámica Aplicada . doi :10.2514/6.1983-1827.
  15. ^ Hess, J.; Friedman, D. (1983). "Análisis de configuraciones de entrada complejas utilizando un método de panel de orden superior". Applied Aerodynamics Conference . doi :10.2514/6.1983-1828.
  16. ^ Bristow, DR, "Desarrollo de métodos de panel para análisis y diseño subsónico", NASA CR-3234, 1980.
  17. ^ Ashby, Dale L.; Dudley, Michael R.; Iguchi, Steve K.; Browne, Lindsey y Katz, Joseph, "Teoría de flujo potencial y guía de operación para el código de panel PMARC", NASA NASA-TM-102851 1991.
  18. ^ Woodward, FA, Dvorak, FA y Geller, EW, "Un programa informático para cuerpos de elevación tridimensionales en flujo subsónico no viscoso", Informe técnico USAAMRDL, TR 74-18, Ft. Eustis, Virginia, abril de 1974.
  19. ^ Katz, Joseph; Maskew, Brian (abril de 1988). "Modelo aerodinámico inestable de baja velocidad para configuraciones completas de aeronaves". Journal of Aircraft . 25 (4): 302–310. doi :10.2514/3.45564.
  20. ^ Maskew, Brian (febrero de 1982). "Predicción de características aerodinámicas subsónicas: un caso para métodos de panel de orden bajo". Journal of Aircraft . 19 (2): 157–163. doi :10.2514/3.57369.
  21. ^ Maskew, Brian, "Documento teórico del programa VSAERO: un programa informático para calcular características aerodinámicas no lineales de configuraciones arbitrarias", NASA CR-4023, 1987.
  22. ^ Pinella, David y Garrison, Peter, "Túnel de viento digital CMARC; códigos de panel tridimensionales de orden bajo", Aerologic, 2009.
  23. ^ Eppler, R.; Somers, DM, "Un programa informático para el diseño y análisis de perfiles aerodinámicos de baja velocidad", NASA TM-80210, 1980.
  24. ^ Drela, Mark, "XFOIL: Un sistema de análisis y diseño para perfiles aerodinámicos de bajo número de Reynolds", en Springer-Verlag Lecture Notes in Engineering, No. 54, 1989.
  25. ^ Boppe, C. (1977). "Cálculo de flujos transónicos en alas mediante incrustación en rejilla". 15.ª Reunión de Ciencias Aeroespaciales . doi :10.2514/6.1977-207.
  26. ^ ab Murman, Earll M.; Cole, Julian D. (enero de 1971). "Cálculo de flujos transónicos planos estacionarios". AIAA Journal . 9 (1): 114–121. Bibcode :1971AIAAJ...9..114C. doi :10.2514/3.6131.
  27. ^ Una teoría de secciones de alas supercríticas, con programas informáticos y ejemplos . Apuntes de clase sobre economía y sistemas matemáticos. Vol. 66. 1972. doi :10.1007/978-3-642-80678-0. ISBN 978-3-540-05807-6.[ página necesaria ]
  28. ^ Mead, HR; Melnik, RE, "GRUMFOIL: Un código informático para el flujo transónico viscoso sobre perfiles aerodinámicos", NASA CR-3806, 1985.
  29. ^ Jameson, A.; Caughey, D. (1977). "Un método de volumen finito para cálculos de flujo potencial transónico". Tercera Conferencia sobre Dinámica de Fluidos Computacional . doi :10.2514/6.1977-635.
  30. ^ Samant, S.; Bussoletti, J.; Johnson, F.; Burkhart, R.; Everson, B.; Melvin, R.; Young, D.; Erickson, L.; Madson, M. (1987). "TRANAIR - Un código informático para análisis transónicos de configuraciones arbitrarias". 25.ª Reunión de Ciencias Aeroespaciales de la AIAA . doi :10.2514/6.1987-34.
  31. ^ Jameson, A.; Schmidt, Wolfgang; Turkel, ELI (1981). "Solución numérica de las ecuaciones de Euler mediante métodos de volumen finito utilizando esquemas de pasos de tiempo de Runge Kutta". 14.ª Conferencia sobre dinámica de fluidos y plasma . doi :10.2514/6.1981-1259.
  32. ^ Raj, Pradeep; Brennan, James E. (1989). "Mejoras de un método aerodinámico de Euler para el análisis de flujo transónico". Journal of Aircraft . 26 : 13–20. doi :10.2514/3.45717.
  33. ^ Tidd, D.; Strash, D.; Epstein, B.; Luntz, A.; Nachshon, A.; Rubin, T. (1991). "Aplicación de un método de Euler multimalla 3-D eficiente (MGAERO) para completar configuraciones de aeronaves". 9.ª Conferencia de Aerodinámica Aplicada . doi :10.2514/6.1991-3236.
  34. ^ Melton, John; Berger, Marsha; Aftosmis, Michael; Wong, Michael (1995). "Aplicaciones 3D de un método de Euler de cuadrícula cartesiana". 33.ª Reunión y exposición de ciencias aeroespaciales . doi :10.2514/6.1995-853.
  35. ^ Karman, l, Jr, Steve (1995). "SPLITFLOW - Un código CFD de cuadrícula cartesiana/prismática no estructurada en 3D para geometrías complejas". 33.ª Reunión y exposición de ciencias aeroespaciales . doi :10.2514/6.1995-343.{{cite book}}: CS1 maint: multiple names: authors list (link)
  36. ^ Marshall, David; Ruffin, Stephen (2004). "Un esquema de cuadrícula cartesiana de límites integrados para flujos viscosos utilizando un nuevo tratamiento de condición de límite de pared viscosa". 42.ª Reunión y exposición de ciencias aeroespaciales de la AIAA. doi :10.2514/6.2004-581. ISBN 978-1-62410-078-9.
  37. ^ Jameson, A.; Baker, T.; Weatherill, N. (1986). "Cálculo del flujo transónico no viscoso sobre una aeronave completa". 24.ª Reunión de Ciencias Aeroespaciales . doi :10.2514/6.1986-103.
  38. ^ Giles, M.; Drela, M.; Thompkins, Jr., W. (1985). "Solución de Newton de ecuaciones de Euler transónicas directas e inversas". 7.ª Conferencia de Física Computacional . doi :10.2514/6.1985-1530.
  39. ^ Drela, Mark (1990). "Solución de Newton de flujos aerodinámicos acoplados viscosos/no viscosos de múltiples elementos". XXI Conferencia sobre dinámica de fluidos, dinámica de plasma y láseres . doi :10.2514/6.1990-1470.
  40. ^ Drela, M. y Youngren H., "Guía del usuario para MISES 2.53", Laboratorio de Ciencias Computacionales del MIT, diciembre de 1998.
  41. ^ Jop, Pierre; Forterre, Yoël; Pouliquen, Olivier (junio de 2006). "Una ley constitutiva para flujos granulares densos". Nature . 441 (7094): 727–730. arXiv : cond-mat/0612110 . Código Bibliográfico :2006Natur.441..727J. doi :10.1038/nature04801. ISSN  1476-4687.
  42. ^ Biroun, Mehdi H.; Mazzei, Luca (junio de 2024). "Flujos granulares no canalizados: efecto de la geometría inicial de la columna granular en la dinámica de fluidos". Chemical Engineering Science . 292 : 119997. doi :10.1016/j.ces.2024.119997. ISSN  0009-2509.
  43. ^ Ferziger, JH y Peric, M. (2002). Métodos computacionales para dinámica de fluidos . Springer-Verlag.{{cite book}}: CS1 maint: multiple names: authors list (link)
  44. ^ "Ecuaciones de Navier-Stokes" . Consultado el 7 de enero de 2020 .
  45. ^ abcdefghij Panton, RL (1996). Flujo incompresible . John Wiley and Sons.
  46. ^ abcd Landau, LD y Lifshitz, EM (2007). Mecánica de fluidos . Elsevier.{{cite book}}: CS1 maint: multiple names: authors list (link)
  47. ^ ab Fox, RW y McDonald, AT (1992). Introducción a la mecánica de fluidos . John Wiley and Sons.{{cite book}}: CS1 maint: multiple names: authors list (link)
  48. ^ ab Poinsot, T. y Veynante, D. (2005). Combustión teórica y numérica . RT Edwards.{{cite book}}: CS1 maint: multiple names: authors list (link)
  49. ^ abcd Kundu, P. (1990). Mecánica de fluidos . Prensa académica.
  50. ^ ab "Favre promedió las ecuaciones de Navier-Stokes" . Consultado el 7 de enero de 2020 .
  51. ^ Bailly, C. y Daniel J. (2000). "Solución numérica de problemas de propagación acústica utilizando ecuaciones de Euler linealizadas". AIAA Journal . 38 (1): 22–29. Bibcode :2000AIAAJ..38...22B. doi :10.2514/2.949.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  52. ^ Harley, JC y Huang, Y. y Bau, HH y Zemel, JN (1995). "Flujo de gas en microcanales". Journal of Fluid Mechanics . 284 : 257–274. Bibcode :1995JFM...284..257H. doi :10.1017/S0022112095000358. S2CID  122833857.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  53. ^ "Ecuaciones de Euler unidimensionales" . Consultado el 12 de enero de 2020 .
  54. ^ Cavazzuti, M. y Corticelli, MA y Karayiannis, TG (2019). "Flujos de Fanno compresibles en microcanales: un modelo numérico cuasi-2D mejorado para flujos laminares". Progreso en ciencia e ingeniería térmica . 10 : 10–26. doi : 10.1016/j.tsep.2019.01.003 . hdl : 11392/2414220 .{{cite journal}}: CS1 maint: multiple names: authors list (link)
  55. ^ Patankar, Suhas V. (1980). Transferencia numérica de calor y flujo de fluidos . Hemisphere Publishing Corporation. ISBN 978-0891165224.
  56. ^ "Explicación detallada del método de elementos finitos (MEF)". www.comsol.com . Consultado el 15 de julio de 2022 .
  57. ^ ab Anderson, John David (1995). Dinámica de fluidos computacional: conceptos básicos con aplicaciones. McGraw-Hill. ISBN 978-0-07-113210-7.
  58. ^ Surana, KA; Allu, S.; Tenpas, PW; Reddy, JN (febrero de 2007). "versión k del método de elementos finitos en dinámica de gases: soluciones numéricas de diferenciabilidad global de orden superior". Revista internacional de métodos numéricos en ingeniería . 69 (6): 1109–1157. Bibcode :2007IJNME..69.1109S. doi :10.1002/nme.1801. S2CID  122551159.
  59. ^ Huebner, KH; Thornton, EA; y Byron, TD (1995). El método de elementos finitos para ingenieros (tercera edición). Wiley Interscience.
  60. ^ Cottet, Georges-Henri; Koumoutsakos, Petros D. (2000). Métodos Vortex: teoría y práctica . Cambridge, Reino Unido: Cambridge Univ. Press. ISBN 0-521-62186-0.
  61. ^ Launder, BE; DB Spalding (1974). "El cálculo numérico de flujos turbulentos". Métodos informáticos en mecánica aplicada e ingeniería . 3 (2): 269–289. Bibcode :1974CMAME...3..269L. doi :10.1016/0045-7825(74)90029-2.
  62. ^ de Wilcox, David C. (2006). Modelado de turbulencia para CFD (3.ª edición). DCW Industries, Inc. ISBN 978-1-928729-08-2.
  63. ^ Pope, SB (2000). Flujos turbulentos . Cambridge University Press. ISBN 978-0-521-59886-6.
  64. ^ Farge, Marie ; Schneider, Kai (2001). "Simulación de vórtices coherentes (CVS), un modelo de turbulencia semideterminista utilizando wavelets". Flujo, turbulencia y combustión . 66 (4): 393–426. doi :10.1023/A:1013512726409. S2CID  53464243.
  65. ^ Goldstein, Daniel; Vasilyev, Oleg (1995). "Método de simulación de grandes remolinos adaptativo coherente estocástico". Física de fluidos A . 24 (7): 2497. Bibcode :2004PhFl...16.2497G. CiteSeerX 10.1.1.415.6540 . doi :10.1063/1.1736671. 
  66. ^ Lundgren, TS (1969). "Ecuación modelo para turbulencia no homogénea". Física de fluidos A . 12 (3): 485–497. Código Bibliográfico :1969PhFl...12..485L. doi :10.1063/1.1692511.
  67. ^ Colucci, PJ; Jaberi, FA; Givi, P.; Pope, SB (1998). "Función de densidad filtrada para simulación de grandes remolinos de flujos reactivos turbulentos". Física de fluidos A . 10 (2): 499–515. Código Bibliográfico :1998PhFl...10..499C. doi :10.1063/1.869537.
  68. ^ Fox, Rodney (2003). Modelos computacionales para flujos reactivos turbulentos . Cambridge University Press. ISBN 978-0-521-65049-6.
  69. ^ Pope, SB (1985). "Métodos PDF para flujos reactivos turbulentos". Progreso en la ciencia de la energía y la combustión . 11 (2): 119–192. Bibcode :1985PrECS..11..119P. doi :10.1016/0360-1285(85)90002-4.
  70. ^ Krueger, Steven K. (1993). "Simulaciones de remolinos lineales de mezcla en un flujo turbulento homogéneo". Física de fluidos . 5 (4): 1023–1034. Código Bibliográfico :1993PhFlA...5.1023M. doi :10.1063/1.858667.
  71. ^ Hirt, CW; Nichols, BD (enero de 1981). "Método de volumen de fluido (VOF) para la dinámica de límites libres". Journal of Computational Physics . 39 (1): 201–225. Bibcode :1981JCoPh..39..201H. doi :10.1016/0021-9991(81)90145-5.
  72. ^ Unverdi, Salih Ozen; Tryggvason, Grétar (mayo de 1992). "Un método de seguimiento frontal para flujos multifluídicos viscosos e incompresibles". Journal of Computational Physics . 100 (1): 25–37. Bibcode :1992JCoPh.100...25U. doi :10.1016/0021-9991(92)90307-K. hdl : 2027.42/30059 .
  73. ^ Benzi, Michele; Golub, Gene H.; Liesen, Jörg (mayo de 2005). "Solución numérica de problemas de punto de silla". Acta Numerica . 14 : 1–137. Bibcode :2005AcNum..14....1B. CiteSeerX 10.1.1.409.4160 . doi :10.1017/S0962492904000212. S2CID  122717775. 
  74. ^ Elman, Howard; Howle, VE; Shadid, John; Shuttleworth, Robert; Tuminaro, Ray (enero de 2008). "Una taxonomía y comparación de preacondicionadores multinivel de bloques paralelos para las ecuaciones incompresibles de Navier-Stokes". Journal of Computational Physics . 227 (3): 1790–1808. Bibcode :2008JCoPh.227.1790E. doi :10.1016/j.jcp.2007.09.026. S2CID  16365489.
  75. ^ Adamson, MR (enero de 2006). "Biografías". IEEE Annals of the History of Computing . 28 (1): 99–103. doi :10.1109/MAHC.2006.5.
  76. ^ Jameson, Antony (mayo de 1974). "Solución iterativa de flujos transónicos sobre perfiles aerodinámicos y alas, incluidos flujos a Mach 1". Communications on Pure and Applied Mathematics . 27 (3): 283–309. doi :10.1002/cpa.3160270302.
  77. ^ Borland, CJ, "XTRAN3S - Aerodinámica transónica estable y no estable para aplicaciones aeroelásticas", AFWAL-TR-85-3214, Laboratorios aeronáuticos Wright de la Fuerza Aérea, Base de la Fuerza Aérea Wright-Patterson, Ohio, enero de 1986
  78. ^ Kaufmann, TAS, Graefe, R., Hormes, M., Schmitz-Rode, T. y Steinseiferand, U., "Dinámica de fluidos computacional en ingeniería biomédica", Dinámica de fluidos computacional: teoría, análisis y aplicaciones, págs. 109-136
  79. ^ Lao, Shandong; Holt, Aaron; Vaidhynathan, Deepthi; Sitaraman, Hariswaran; Hrenya, Christine M.; Hauser, Thomas (2021). "Comparación de rendimiento del solucionador CFD-DEM MFiX-Exa en GPU y CPU". arXiv : 2108.08821 [cs.DC].
  80. ^ Wu, Kui; Truong, Nghia; Yuksel, Cem; Hoetzlein, Rama (mayo de 2018). "Simulaciones rápidas de fluidos con volúmenes dispersos en la GPU". Computer Graphics Forum . 37 (2): 157–167. doi :10.1111/cgf.13350. S2CID  43945038.
  81. ^ "Soporte de la aplicación Intersect 360 HPC" (PDF) .

Notas

Enlaces externos