La hidrodinámica de partículas suavizadas ( SPH ) es un método computacional utilizado para simular la mecánica de medios continuos, como la mecánica de sólidos y los flujos de fluidos . Fue desarrollado por Gingold y Monaghan [2] y Lucy [3] en 1977, inicialmente para problemas astrofísicos. Se ha utilizado en muchos campos de investigación, incluidos la astrofísica , la balística , la vulcanología y la oceanografía . Es un método lagrangiano sin malla (donde las coordenadas se mueven con el fluido) y la resolución del método se puede ajustar fácilmente con respecto a variables como la densidad .
Método
Ventajas
Por construcción, SPH es un método sin malla , lo que lo hace ideal para simular problemas dominados por dinámicas de límites complejas, como flujos de superficie libre o grandes desplazamientos de límites.
La falta de una malla simplifica significativamente la implementación del modelo y su paralelización, incluso para arquitecturas de múltiples núcleos . [4] [5]
El modelo SPH se puede ampliar fácilmente a una amplia variedad de campos e hibridar con algunos otros modelos, como se analiza en Modelado de la física.
Como se discutió en la sección sobre SPH débilmente compresible, el método tiene excelentes características de conservación.
El costo computacional de las simulaciones SPH por número de partículas es significativamente menor que el costo de las simulaciones basadas en cuadrículas por número de celdas cuando la métrica de interés está relacionada con la densidad del fluido (por ejemplo, la función de densidad de probabilidad de las fluctuaciones de densidad). [6] Esto es así porque en SPH la resolución se pone donde está la materia.
Limitaciones
El establecimiento de condiciones de contorno en SPH, como entradas y salidas [7] y paredes [8], es más difícil que con los métodos basados en cuadrículas. De hecho, se ha afirmado que "el tratamiento de las condiciones de contorno es sin duda uno de los puntos técnicos más difíciles del método SPH". [9] Este desafío se debe en parte a que en SPH las partículas cercanas al límite cambian con el tiempo. [10] No obstante, existen condiciones de contorno de pared para SPH. [8] [10] [11]
El costo computacional de las simulaciones SPH por número de partículas es significativamente mayor que el costo de las simulaciones basadas en cuadrículas por número de celdas cuando la métrica de interés no está relacionada (directamente) con la densidad (por ejemplo, el espectro de energía cinética). [6] Por lo tanto, pasando por alto los problemas de aceleración paralela , la simulación de flujos de densidad constante (por ejemplo, aerodinámica externa ) es más eficiente con métodos basados en cuadrículas que con SPH.
Ejemplos
Dinámica de fluidos
La hidrodinámica de partículas suavizadas también se utiliza cada vez más para modelar el movimiento de fluidos . Esto se debe a varias ventajas sobre las técnicas tradicionales basadas en cuadrículas. En primer lugar, la hidrodinámica de partículas suavizadas garantiza la conservación de la masa sin cálculos adicionales, ya que las partículas en sí mismas representan la masa. En segundo lugar, la hidrodinámica de partículas suavizadas calcula la presión a partir de las contribuciones ponderadas de las partículas vecinas en lugar de resolver sistemas lineales de ecuaciones. Por último, a diferencia de las técnicas basadas en cuadrículas, que deben rastrear los límites de los fluidos, la hidrodinámica de partículas suavizadas crea una superficie libre para fluidos que interactúan en dos fases directamente, ya que las partículas representan el fluido más denso (normalmente agua) y el espacio vacío representa el fluido más ligero (normalmente aire). Por estas razones, es posible simular el movimiento de fluidos utilizando la hidrodinámica de partículas suavizadas en tiempo real. Sin embargo, tanto las técnicas basadas en cuadrículas como las de la hidrodinámica de partículas suavizadas aún requieren la generación de una geometría de superficie libre renderizable utilizando una técnica de poligonización como metaballs y marching cubes , point splatting o visualización de "alfombra". Para la dinámica de gases, es más apropiado utilizar la propia función kernel para producir una representación de la densidad de la columna de gas (por ejemplo, como se hace en el paquete de visualización SPLASH).
Una desventaja de las técnicas basadas en cuadrículas es la necesidad de un gran número de partículas para producir simulaciones de resolución equivalente. En la implementación típica de las cuadrículas uniformes y las técnicas de partículas SPH, se utilizarán muchos vóxeles o partículas para llenar volúmenes de agua que nunca se renderizan. Sin embargo, la precisión puede ser significativamente mayor con técnicas sofisticadas basadas en cuadrículas, especialmente aquellas acopladas con métodos de partículas (como conjuntos de niveles de partículas), ya que es más fácil hacer cumplir la condición de incompresibilidad en estos sistemas. La SPH para simulación de fluidos se está utilizando cada vez más en animación y juegos en tiempo real donde la precisión no es tan crítica como la interactividad.
Los trabajos recientes en SPH para simulación de fluidos han aumentado el rendimiento, la precisión y las áreas de aplicación:
B. Solenthaler, 2009, desarrolla SPH predictivo-correctivo (PCISPH) para permitir mejores restricciones de incompresibilidad [12]
M. Ihmsen et al., 2010, introducen el manejo de límites y el paso de tiempo adaptativo para PCISPH para interacciones precisas de cuerpos rígidos [13]
K. Bodin et al., 2011, reemplazan la ecuación estándar de presión de estado con una restricción de densidad y aplican un integrador de tiempo variacional [14]
R. Hoetzlein, 2012, desarrolla un SPH basado en GPU eficiente para escenas grandes en Fluids v.3 [15]
N. Akinci et al., 2012, presentan una técnica versátil de manejo de límites y acoplamiento SPH-rígido bidireccional que se basa completamente en fuerzas hidrodinámicas; el enfoque es aplicable a diferentes tipos de solucionadores SPH [16].
M. Macklin et al., 2013 simula flujos incompresibles dentro del marco de dinámica basada en la posición, para pasos de tiempo más grandes [17]
N. Akinci et al., 2013, presentan una técnica versátil de tensión superficial y adhesión fluido-sólido bidireccional que permite simular una variedad de efectos físicos interesantes que se observan en la realidad [18].
J. Kyle y E. Terrell, 2013, aplican SPH a la lubricación de película completa [19]
A. Mahdavi y N. Talebbeydokhti, 2015, proponen un algoritmo híbrido para la implementación de la condición de límite sólido y simulan el flujo sobre un vertedero de cresta aguda [20]
S. Tavakkol et al., 2016, desarrollan curvSPH, que hace que el tamaño horizontal y vertical de las partículas sean independientes y genera una distribución de masa uniforme a lo largo de los límites curvos [21].
W. Kostorz y A. Esmail-Yakas, 2020, proponen un método general, eficiente y simple para evaluar los factores de normalización cerca de los límites planos por partes [11]
Colagrossi et al., 2019, estudian el flujo alrededor de un cilindro cerca de una superficie libre y lo comparan con otras técnicas [1]
Astrofísica
La resolución adaptativa de la hidrodinámica de partículas suavizadas, la conservación numérica de cantidades físicamente conservadas y la capacidad de simular fenómenos que abarcan muchos órdenes de magnitud la hacen ideal para los cálculos en astrofísica teórica . [22]
La SPH se utiliza para modelar flujos hidrodinámicos, incluidos los posibles efectos de la gravedad . La incorporación de otros procesos astrofísicos que pueden ser importantes, como la transferencia radiativa y los campos magnéticos , es un área activa de investigación en la comunidad astronómica y ha tenido un éxito limitado. [25] [26]
Mecánica de sólidos
Libersky y Petschek [27] [28]
extendieron la SPH a la mecánica de sólidos. La principal ventaja de la SPH en esta aplicación es la posibilidad de tratar con una distorsión local mayor que los métodos basados en cuadrículas. Esta característica se ha explotado en muchas aplicaciones en mecánica de sólidos: conformado de metales, impacto, crecimiento de grietas, fractura, fragmentación, etc.
Otra ventaja importante de los métodos sin malla en general, y de SPH en particular, es que los problemas de dependencia de la malla se evitan naturalmente dada la naturaleza sin malla del método. En particular, la alineación de la malla está relacionada con problemas que involucran grietas y se evita en SPH debido al soporte isotrópico de las funciones kernel. Sin embargo, las formulaciones clásicas de SPH sufren inestabilidades de tracción [29]
y falta de consistencia. [30]
En los últimos años, se han introducido diferentes correcciones para mejorar la precisión de la solución SPH, lo que llevó al RKPM de Liu et al. [31]
Randles y Libersky [32]
y Johnson y Beissel [33]
intentaron resolver el problema de la consistencia en su estudio de los fenómenos de impacto.
Dyka et al. [34] [35]
y Randles y Libersky [36]
introdujeron la integración de puntos de tensión en SPH y Ted Belytschko et al. [37]
demostraron que la técnica de puntos de tensión elimina la inestabilidad debida a modos singulares espurios, mientras que las inestabilidades de tracción se pueden evitar utilizando un núcleo lagrangiano. Se pueden encontrar muchos otros estudios recientes en la literatura dedicados a mejorar la convergencia del método SPH.
Las recientes mejoras en la comprensión de la convergencia y la estabilidad de los SPH han permitido aplicaciones más generalizadas en la mecánica de sólidos. Otros ejemplos de aplicaciones y desarrollos del método incluyen:
Simulaciones de conformado de metales. [38]
Método basado en SPH SPAM (Smoothed Particle Applied Mechanics) para fractura por impacto en sólidos por William G. Hoover . [39]
SPH modificado (SPH/MLSPH) para fractura y fragmentación. [40]
Taylor-SPH (TSPH) para la propagación de ondas de choque en sólidos. [41]
La SPH de coordenadas generalizadas (GSPH) asigna partículas de manera no homogénea en el sistema de coordenadas cartesianas y las organiza mediante el mapeo en un sistema de coordenadas generalizado en el que las partículas están alineadas con un espaciado uniforme. [42]
Herramientas numéricas
Interpolaciones
El método de hidrodinámica de partículas suavizadas (SPH) funciona dividiendo el fluido en un conjunto de elementos móviles discretos , denominados partículas. Su naturaleza lagrangiana permite fijar su posición mediante la integración de su velocidad como:
Estas partículas interactúan a través de una función kernel con un radio característico conocido como "longitud de suavizado", que se representa típicamente en ecuaciones por . Esto significa que la cantidad física de cualquier partícula se puede obtener sumando las propiedades relevantes de todas las partículas que se encuentran dentro del rango del kernel, utilizándose este último como una función de ponderación . Esto se puede entender en dos pasos. Primero, un campo arbitrario se escribe como una convolución con :
El error al realizar la aproximación anterior es de orden . En segundo lugar, la integral se aproxima utilizando una suma de Riemann sobre las partículas:
donde la sumatoria incluye todas las partículas en la simulación. es el volumen de la partícula , es el valor de la cantidad de la partícula y denota la posición. Por ejemplo, la densidad de la partícula se puede expresar como:
donde denota la masa de la partícula y la densidad de la partícula, mientras que es una notación corta para . El error cometido al aproximar la integral mediante una suma discreta depende de , del tamaño de la partícula (es decir , , siendo la dimensión espacial) y de la disposición de las partículas en el espacio. Este último efecto aún es poco conocido. [43]
Las funciones kernel que se usan comúnmente incluyen la función gaussiana , la spline quintica y el kernel de Wendland. [44] Los dos últimos kernels tienen un soporte compacto (a diferencia del gaussiano, donde hay una pequeña contribución a cualquier distancia finita), con un soporte proporcional a . Esto tiene la ventaja de ahorrar esfuerzo computacional al no incluir las contribuciones relativamente menores de partículas distantes.
Aunque el tamaño de la longitud de suavizado se puede fijar tanto en el espacio como en el tiempo , esto no aprovecha toda la potencia de SPH. Al asignar a cada partícula su propia longitud de suavizado y permitir que varíe con el tiempo, se puede hacer que la resolución de una simulación se adapte automáticamente en función de las condiciones locales. Por ejemplo, en una región muy densa donde muchas partículas están juntas, la longitud de suavizado se puede hacer relativamente corta, lo que produce una alta resolución espacial. Por el contrario, en regiones de baja densidad donde las partículas individuales están muy separadas y la resolución es baja, se puede aumentar la longitud de suavizado, optimizando el cálculo para las regiones de interés.
Discretización de ecuaciones de gobierno
Para partículas de masa constante, diferenciar la densidad interpolada con respecto al tiempo da como resultado
es evidente que su lado derecho es una aproximación de ; por lo tanto, se define un operador de divergencia discreto de la siguiente manera:
Este operador proporciona una aproximación SPH de en la partícula para un conjunto dado de partículas con masas , posiciones y velocidades dadas .
La otra ecuación importante para un fluido no viscoso compresible es la ecuación de Euler para el equilibrio del momento:
De manera similar a la continuidad, la tarea es definir un operador de gradiente discreto para escribir
Una opción es
que tiene la propiedad de ser asfálticamente adjunta con el operador de divergencia anterior, en el sentido de que
Esta es una versión discreta de la identidad continua.
Esta propiedad conduce a buenas propiedades de conservación. [45]
Nótese también que esta elección conduce a un operador de divergencia simétrico y a un gradiente antisimétrico. Aunque existen varias formas de discretizar el gradiente de presión en las ecuaciones de Euler, la forma antisimétrica anterior es la más reconocida. Respalda la conservación estricta del momento lineal y angular. Esto significa que una fuerza que se ejerce sobre partícula por partícula es igual a la que se ejerce sobre partícula por partícula , incluido el cambio de signo de la dirección efectiva, gracias a la propiedad de antisimetría .
Sin embargo, se han propuesto otros operadores que pueden tener un mejor rendimiento numérico o físico. Por ejemplo, una desventaja de estos operadores es que, si bien la divergencia es consistente en orden cero (es decir, produce cero cuando se aplica a un campo vectorial constante), se puede observar que el gradiente no lo es. Se han propuesto varias técnicas para evitar este problema, lo que conduce a operadores renormalizados (véase, por ejemplo, [46] ).
Principio variacional
Las ecuaciones de gobierno de SPH anteriores se pueden derivar de un principio de mínima acción , comenzando desde el Lagrangiano de un sistema de partículas:
Cuando se aplica al Lagrangiano anterior, se obtiene la siguiente ecuación de momento:
donde se ha utilizado la regla de la cadena, ya que depende de , y esta última, de la posición de las partículas. Utilizando la propiedad termodinámica podemos escribir
Al conectar la interpolación de densidad SPH y diferenciar explícitamente se llega a
que es la ecuación de momento SPH ya mencionada, donde reconocemos el operador. Esto explica por qué se conserva el momento lineal, y permite que se conserve también la conservación del momento angular y la energía. [47]
Integración temporal
A partir del trabajo realizado en los años 80 y 90 sobre la integración numérica de partículas puntuales en grandes aceleradores, se han desarrollado integradores temporales adecuados con propiedades de conservación precisas a largo plazo; se denominan integradores simplécticos . El más popular en la literatura sobre SPH es el esquema de salto de rana , que dice para cada partícula :
donde es el paso de tiempo, los superíndices representan iteraciones de tiempo mientras que es la aceleración de la partícula, dada por el lado derecho de la ecuación del momento.
Existen otros integradores simplécticos (véase el libro de referencia [48] ). Se recomienda utilizar un esquema simpléctico (incluso de orden bajo) en lugar de un esquema no simpléctico de orden alto, para evitar la acumulación de errores después de muchas iteraciones.
La integración de la densidad no se ha estudiado en profundidad (consulte a continuación para obtener más detalles).
Los esquemas simplécticos son conservadores pero explícitos, por lo que su estabilidad numérica requiere condiciones de estabilidad, análogas a la condición de Courant-Friedrichs-Lewy (ver más abajo).
Técnicas de límites
En caso de que la convolución SPH se practique cerca de un límite, es decir, más cerca que s · h , entonces el soporte integral se trunca. De hecho, cuando la convolución se ve afectada por un límite, la convolución se dividirá en 2 integrales,
donde B( r ) es la bola de soporte compacta centrada en r , con radio s · h , y Ω( r ) denota la parte del soporte compacto dentro del dominio computacional, Ω ∩ B( r ) . Por lo tanto, la imposición de condiciones de contorno en SPH se basa completamente en la aproximación de la segunda integral en el lado derecho. Por supuesto, lo mismo se puede aplicar al cálculo de operadores diferenciales,
Se han introducido varias técnicas en el pasado para modelar límites en SPH.
Descuido integral
El modelo de límite más sencillo es descuidar la integral,
de modo que sólo se tengan en cuenta las interacciones en masa,
Este es un enfoque popular cuando se considera la superficie libre en simulaciones monofásicas. [49]
La principal ventaja de esta condición de contorno es su evidente simplicidad. Sin embargo, se deben tener en cuenta varios problemas de coherencia cuando se aplica esta técnica de contorno. [49] De hecho, esto supone una importante limitación para sus posibles aplicaciones.
Extensión de fluido
Probablemente la metodología más popular, o al menos la más tradicional, para imponer condiciones de contorno en SPH, es la técnica de Extensión de Fluidos. Dicha técnica se basa en poblar el soporte compacto a través del contorno con las llamadas partículas fantasma, imponiendo convenientemente sus valores de campo. [50]
En esta línea, la metodología de negligencia integral puede considerarse como un caso particular de extensiones de fluidos, donde el campo, A , se desvanece fuera del dominio computacional.
El principal beneficio de esta metodología es la simplicidad, ya que la contribución de los límites se calcula como parte de las interacciones en masa. Además, esta metodología ha sido analizada en profundidad en la literatura. [51] [50] [52]
Por otra parte, la implementación de partículas fantasma en el dominio truncado no es una tarea trivial, por lo que modelar formas de contorno complejas se vuelve engorroso. Los dos enfoques más populares para poblar el dominio vacío con partículas fantasma son las partículas reflejadas [53] y las partículas fijas [50] .
Integral de frontera
La técnica de límites más nueva es la metodología de integral de límites. [54] En esta metodología, la integral de volumen vacío se reemplaza por una integral de superficie y una renormalización:
donde n j es la normal del elemento de contorno genérico j -ésimo. El término de superficie también se puede resolver considerando una expresión semianalítica. [54]
Modelado de la física
Hidrodinámica
Enfoque débilmente compresible
Otra forma de determinar la densidad se basa en el propio operador de suavizado SPH. Por lo tanto, la densidad se estima a partir de la distribución de partículas utilizando la interpolación SPH. Para superar los errores no deseados en la superficie libre mediante el truncamiento del núcleo, la formulación de densidad se puede integrar nuevamente en el tiempo. [54]
El SPH débilmente compresible en dinámica de fluidos se basa en la discretización de las ecuaciones de Navier-Stokes o ecuaciones de Euler para fluidos compresibles. Para cerrar el sistema, se utiliza una ecuación de estado adecuada para vincular la presión y la densidad . Generalmente, en SPH se utiliza la denominada ecuación de Cole [55]
(a veces denominada erróneamente como la " ecuación de Tait "). Su texto es el siguiente:
donde es la densidad de referencia y la velocidad del sonido . Para el agua, se utiliza comúnmente. Se agrega la presión de fondo para evitar valores de presión negativos.
Los fluidos reales casi incompresibles, como el agua, se caracterizan por velocidades del sonido muy altas, del orden de . Por lo tanto, la información sobre la presión viaja más rápido en comparación con el flujo real, lo que da lugar a números de Mach muy pequeños . La ecuación del momento da como resultado la siguiente relación:
donde es el cambio de densidad y el vector de velocidad. En la práctica, se adopta un valor de c menor que el real para evitar pasos de tiempo demasiado pequeños en el esquema de integración temporal. Generalmente, se adopta una velocidad numérica del sonido tal que se permiten variaciones de densidad menores al 1%. Este es el llamado supuesto de compresibilidad débil. Esto corresponde a un número de Mach menor que 0,1, lo que implica:
donde es necesario estimar la velocidad máxima , por ejemplo mediante la ley de Torricelli o una estimación aproximada. Como solo se producen pequeñas variaciones de densidad, se puede adoptar una ecuación de estado lineal: [56]
Generalmente los esquemas débilmente compresibles se ven afectados por un ruido espurio de alta frecuencia en los campos de presión y densidad. [57]
Este fenómeno es causado por la interacción no lineal de las ondas acústicas y por el hecho de que el esquema es explícito en el tiempo y centrado en el espacio. [58]
A lo largo de los años se han propuesto diversas técnicas para solucionar este problema, que pueden clasificarse en tres grupos diferentes:
los esquemas que adoptan filtros de densidad,
los modelos que añaden un término difusivo en la ecuación de continuidad,
los esquemas que emplean solucionadores de Riemann para modelar la interacción de partículas.
Técnica de filtrado de densidad
Los esquemas del primer grupo aplican un filtro directamente sobre el campo de densidad para eliminar el ruido numérico espurio. Los filtros más utilizados son el MLS (mínimos cuadrados móviles) y el filtro Shepard [57]
que se pueden aplicar en cada paso de tiempo o cada n pasos de tiempo. Cuanto más frecuente sea el uso del procedimiento de filtrado, más regulares se obtendrán los campos de densidad y presión. Por otro lado, esto conduce a un aumento de los costos computacionales. En simulaciones de largo plazo, el uso del procedimiento de filtrado puede conducir a la interrupción del componente de presión hidrostática y a una inconsistencia entre el volumen global del fluido y el campo de densidad. Además, no garantiza el cumplimiento de la condición de contorno de superficie libre dinámica.
Técnica de difusión del término
Una forma diferente de suavizar el campo de densidad y presión es agregar un término difusivo dentro de la ecuación de continuidad (grupo 2):
Los primeros esquemas que adoptaron este enfoque fueron descritos en Ferrari [59]
y en Molteni [56],
donde el término difusivo se modeló como un laplaciano del campo de densidad. Un enfoque similar también fue utilizado en Fatehi y Manzari [60] .
En Antuono et al. [61] se propuso
una corrección al término difusivo de Molteni [56] para eliminar algunas inconsistencias cerca de la superficie libre. En este caso, el término difusivo adoptado es equivalente a un operador diferencial de orden superior en el campo de densidad. [62]
El esquema se denomina δ-SPH y conserva todas las propiedades de conservación del SPH sin difusión (por ejemplo, momentos lineales y angulares, energía total, véase [63]
) junto con una representación suave y regular de los campos de densidad y presión.
En el tercer grupo se encuentran aquellos esquemas SPH que emplean flujos numéricos obtenidos a través de solucionadores de Riemann para modelar las interacciones de partículas. [64] [65] [66]
Técnica de solución de Riemann
Para un método SPH basado en solucionadores de Riemann, se construye un problema de Riemann entre partículas a lo largo de un vector unitario que apunta de partícula a partícula . En este problema de Riemann, los estados iniciales izquierdo y derecho están en las partículas y , respectivamente. Los estados y son
La solución del problema de Riemann da como resultado tres ondas que emanan de la discontinuidad. Dos ondas, que pueden ser ondas de choque o de rarefacción, que viajan con la velocidad de onda más pequeña o más grande. La onda del medio es siempre una discontinuidad de contacto y separa dos estados intermedios, denotados por y . Al suponer que el estado intermedio satisface y , un solucionador de Riemann linealizado para flujos suaves o con choques moderadamente fuertes se puede escribir como
donde y son promedios entre partículas. Con la solución del problema de Riemann, es decir y , la discretización del método SPH es
donde . Esto indica que la velocidad y la presión promedio entre partículas se reemplazan simplemente por la solución del problema de Riemann. Al comparar ambos, se puede ver que la velocidad y la presión intermedias de los promedios entre partículas equivalen a disipación implícita, es decir, regularización de la densidad y viscosidad numérica, respectivamente.
Dado que la discretización anterior es muy disipativa, una modificación sencilla es aplicar un limitador para disminuir las disipaciones numéricas implícitas introducidas al limitar la presión intermedia en [67].
donde el limitador se define como
Nótese que garantiza que no haya disipación cuando el fluido está bajo la acción de una onda de expansión, es decir , y que el parámetro , se utiliza para modular la disipación cuando el fluido está bajo la acción de una onda de compresión, es decir , . Los experimentos numéricos encontraron que el es generalmente efectivo. Nótese también que la disipación introducida por la velocidad intermedia no está limitada.
Enfoque incompresible
Modelado de viscosidad
En general, la descripción de flujos hidrodinámicos requiere un tratamiento conveniente de los procesos difusivos para modelar la viscosidad en las ecuaciones de Navier-Stokes . Esto requiere una consideración especial porque involucra al operador diferencial laplaciano . Dado que el cálculo directo no proporciona resultados satisfactorios, se han propuesto varios enfoques para modelar la difusión.
Viscosidad artificial
Introducida por Monaghan y Gingold [68],
la viscosidad artificial se utilizó para tratar flujos de fluidos con números de Mach altos .
Aquí, se controla una viscosidad volumétrica mientras actúa de manera similar a la viscosidad artificial de Neumann Richtmeyr. Se define por
donde η h es una pequeña fracción de h (por ejemplo, 0,01 h ) para evitar posibles infinitos numéricos a distancias cercanas.
También se ha demostrado que la viscosidad artificial mejora la estabilidad general de las simulaciones de flujo general. Por lo tanto, se aplica a problemas no viscosos de la siguiente forma:
Es posible no sólo estabilizar simulaciones no viscosas sino también modelar la viscosidad física mediante este enfoque.
se sustituye en la ecuación anterior, donde es el número de dimensiones espaciales del modelo. Este enfoque introduce la viscosidad volumétrica .
Morris
Para números de Reynolds bajos se propuso el modelo de viscosidad de Morris [69]
.
LoShao
Física adicional
Tensión superficial
Transferencia de calor
Turbulencia
Extensiones multifase
Astrofísica
En astrofísica, a menudo se desea modelar la autogravedad además de la hidrodinámica pura. La naturaleza basada en partículas de SPH hace que sea ideal para combinarlo con un solucionador de gravedad basado en partículas, por ejemplo, tree gravity code , [70] particle mesh o particle-particle-mesh .
Mecánica de sólidos e interacción fluido-estructura (FSI)
Formulación lagrangiana total para mecánica de sólidos
Para discretizar las ecuaciones que rigen la dinámica de sólidos, primero se introduce una matriz de corrección [71] [72]
para reproducir la rotación del cuerpo rígido como
dónde
representa el gradiente de la función kernel evaluada en la configuración de referencia inicial. Nótese que los subíndices y se utilizan para indicar partículas sólidas, y la longitud de suavizado es idéntica a la de la discretización de ecuaciones de fluidos.
Utilizando la configuración inicial como referencia, la densidad del sólido se evalúa directamente como
donde es el determinante jacobiano del tensor de deformación .
Ahora podemos discretizar la ecuación del momento en la siguiente forma
donde la primera tensión de Piola-Kirchhoff promediada entre partículas
se define como
También y corresponden a la presión del fluido y a las fuerzas viscosas que actúan sobre la partícula sólida , respectivamente.
Acoplamiento fluido-estructural
En el acoplamiento fluido-estructura, la estructura sólida circundante se comporta como un límite móvil para el fluido, y la condición de límite sin deslizamiento se impone en la interfaz fluido-estructura. Las fuerzas de interacción y que actúan sobre una partícula de fluido , debido a la presencia de la partícula sólida vecina , se pueden obtener como [73]
y
Aquí, la presión y la velocidad imaginarias se definen mediante
donde denota la dirección normal de la superficie de la estructura sólida, y la densidad de partículas imaginarias se calcula a través de la ecuación de estado.
En consecuencia, las fuerzas de interacción y que actúan sobre una partícula sólida están dadas por
y
La propiedad antisimétrica de la derivada de la función kernel asegurará la conservación del momento para cada par de partículas interactuantes y .
^ ab Colagrossi (2019). "Flujo viscoso que pasa por un cilindro cerca de una superficie libre: puntos de referencia con respuestas estables, periódicas y metaestables, resueltos mediante esquemas basados en malla y sin malla". Computers and Fluids . 181 : 345–363. doi :10.1016/j.compfluid.2019.01.007. S2CID 128143912.
^ Gingold, Robert A.; Monaghan, Joseph J. (1977). "Hidrodinámica de partículas suavizadas: teoría y aplicación a estrellas no esféricas". Monthly Notices of the Royal Astronomical Society . 181 (3): 375–89. Bibcode :1977MNRAS.181..375G. doi : 10.1093/mnras/181.3.375 .
^ LB Lucy (1977). "Un enfoque numérico para la prueba de la hipótesis de la fisión". Astron. J . 82 : 1013–1024. Bibcode :1977AJ.....82.1013L. doi :10.1086/112164.
^ Takahiro Harada; Seiichi Koshizuka; Yoichiro Kawaguchi (2007). Hidrodinámica de partículas suavizadas en GPU . Computer Graphics International . págs. 63–70.
^ Alejandro Crespo; Jose M. Dominguez; Anxo Barreiro; Moncho Gomez-Gesteira; Benedict D. Rogers (2011). "GPUs, una nueva herramienta de aceleración en CFD: eficiencia y confiabilidad en métodos de hidrodinámica de partículas suavizadas". PLOS ONE . 6 (6): e20685. Bibcode :2011PLoSO...620685C. doi : 10.1371/journal.pone.0020685 . PMC 3113801 . PMID 21695185.
^ ab Price, DJ (2011). "Hidrodinámica de partículas suavizadas: cosas que desearía que mi madre me hubiera enseñado". Avances en astrofísica computacional: métodos . 453 : 249. arXiv : 1111.1259 . Código Bibliográfico : 2012ASPC..453..249P.
^ "El método de hidrodinámica de partículas suavizadas frente a los métodos numéricos de volumen finito". 2018-03-21 . Consultado el 2018-08-30 .
^ ab Adami, S. y Hu, XY y Adams, NA. (2012). "Una condición de contorno de pared generalizada para la hidrodinámica de partículas suavizadas". Journal of Computational Physics . 231 (21): 7057–7075. Bibcode :2012JCoPh.231.7057A. doi :10.1016/j.jcp.2012.05.005.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ Shadloo, MS y Oger, G. y Touze, DL. (2016). "Método de hidrodinámica de partículas suavizadas para flujos de fluidos, hacia aplicaciones industriales: motivaciones, estado actual y desafíos". Computadoras y fluidos . 136 : 11–34. doi :10.1016/j.compfluid.2016.05.029.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ ab Fraser, K. y Kiss, LI y St-George, L. (2016). "Una condición de contorno de pared generalizada para la hidrodinámica de partículas suavizadas". 14.ª Conferencia internacional LS-DYNA .{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ ab Kostorz (2020). "Un método de integral de contorno semianalítico para funciones radiales con aplicación a la hidrodinámica de partículas suavizadas". Journal of Computational Physics . 417 : 109565. Bibcode :2020JCoPh.41709565K. doi :10.1016/j.jcp.2020.109565. S2CID 219438340.
^ Imhsen (2010). "Manejo de límites y pasos de tiempo adaptativos para PCISPH". Taller sobre interacción de realidad virtual y simulación física VRIPHYS .
^ Bodin (2011). "Fluidos de restricción". IEEE Transactions on Visualization and Computer Graphics . 18 (3): 516–26. doi :10.1109/TVCG.2011.29. PMID 22241284. S2CID 14023161.
^ Hoetzlein (2012). "Fluids v.3, un simulador de fluidos de código abierto y a gran escala".{{cite journal}}: Requiere citar revista |journal=( ayuda )
^ Akinci (2012). "Acoplamiento rígido-fluídico versátil para SPH incompresible". ACM Transactions on Graphics . 31 (4): 1–8. doi :10.1145/2185520.2185558. S2CID 5669154.
^ Macklin (2013). "Fluidos basados en la posición". ACM Transactions on Graphics . 32 (4): 1–12. doi :10.1145/2461912.2461984. S2CID 611962.
^ Akinci (2013). "Tensión superficial versátil y adhesión para fluidos SPH". ACM Transactions on Graphics . 32 (6): 1–8. CiteSeerX 10.1.1.462.8293 . doi :10.1145/2508363.2508395. S2CID 12550964.
^ Revista de tribología (2013). "Aplicación de la hidrodinámica de partículas suavizadas a la lubricación de película completa".{{cite journal}}: Requiere citar revista |journal=( ayuda )
^ Mahdavi y Talebbeydokhti (2015). "Un algoritmo híbrido de tratamiento de límites sólidos para la hidrodinámica de partículas suavizadas". Scientia Iranica, Transacción A, Ingeniería civil . 22 (4): 1457–1469.
^ Revista internacional de métodos numéricos en fluidos (2016). "Hidrodinámica de partículas suavizadas curvilíneas". Revista internacional de métodos numéricos en fluidos . 83 (2): 115–131. Bibcode :2017IJNMF..83..115T. doi :10.1002/fld.4261. S2CID 124425780.
^ Price, Daniel J (2009). "Hidrodinámica astrofísica de partículas lisas". New Astron.rev . 53 (4–6): 78–104. arXiv : 0903.5075 . Código Bibliográfico :2009NewAR..53...78R. doi :10.1016/j.newar.2009.08.007. S2CID 129246.
^ Rosswog, Stephan (2015). "Métodos SPH en el modelado de objetos compactos". Living Rev Comput Astrophys . 1 (1): 1. arXiv : 1406.4224 . Código Bibliográfico :2015LRCA....1....1R. doi :10.1007/lrca-2015-1. S2CID 119119783.
^ Price, Daniel J; Rockefeller, Gabriel; Warren, Michael S (2006). "SNSPH: Un código hidrodinámico de radiación de partículas suavizadas en 3D paralelas". Astrophys. J. 643 : 292–305. arXiv : astro-ph/0512532 . doi :10.1086/501493. S2CID : 16733573.
^ "Formación de estrellas con transferencia radiativa".
^ Price, Daniel J. (1 de febrero de 2012). "Hidrodinámica y magnetohidrodinámica de partículas suavizadas". Journal of Computational Physics . Número especial: Física computacional del plasma. 231 (3): 759–794. arXiv : 1012.1885 . doi :10.1016/j.jcp.2010.12.011. ISSN 0021-9991.
^ Libersky, LD; Petschek, AG (1990). Hidrodinámica de partículas lisas con resistencia de materiales, avances en el método de Lagrange libre . Apuntes de clase en física. Vol. 395. págs. 248–257. doi :10.1007/3-540-54960-9_58. ISBN978-3-540-54960-4.
^ LD Libersky; AG Petschek; AG Carney; TC Hipp; JR Allahdadi; FA High (1993). "Hidrodinámica lagrangiana de deformación: un código SPH tridimensional para la respuesta dinámica del material". J. Comput. Phys . 109 (1): 67–75. Bibcode :1993JCoPh.109...67L. doi :10.1006/jcph.1993.1199.
^ JW Swegle; DA Hicks; SW Attaway (1995). "Análisis de estabilidad hidrodinámica de partículas lisas". J. Comput. Phys . 116 (1): 123–134. Código Bibliográfico :1995JCoPh.116..123S. doi :10.1006/jcph.1995.1010.
^ T. Belytschko; Y. Krongauz; J. Dolbow; C. Gerlach (1998). "Sobre la completitud de los métodos de partículas sin malla". Int. J. Numer. Methods Eng . 43 (5): 785–819. Bibcode :1998IJNME..43..785B. CiteSeerX 10.1.1.28.491 . doi :10.1002/(sici)1097-0207(19981115)43:5<785::aid-nme420>3.0.co;2-9.
^ WK Liu; S. Jun; YF Zhang (1995). "Reproducción de métodos de partículas de núcleo". Int. J. Numer. Methods Eng . 20 (8–9): 1081–1106. Código Bibliográfico :1995IJNMF..20.1081L. doi :10.1002/fld.1650200824.
^ PW Randles; LD Libersky (1997). "Mejoras recientes en el modelado SPH del impacto a hipervelocidad". Int. J. Impact Eng . 20 (6–10): 525–532. doi :10.1016/s0734-743x(97)87441-6.
^ GR Johnson; SR Beissel (1996). "Funciones de suavizado normalizadas para cálculos de impacto de SPH". Int. J. Numer. Methods Eng . 39 (16): 2725–2741. Bibcode :1996IJNME..39.2725J. doi :10.1002/(sici)1097-0207(19960830)39:16<2725::aid-nme973>3.0.co;2-9.
^ CT Dyka; RP Ingel (1995). "Un enfoque para la inestabilidad de la tensión en la hidrodinámica de partículas suavizadas". Comput. Struct . 57 (4): 573–580. doi :10.1016/0045-7949(95)00059-p.
^ CT Dyka; PW Randles; RP Ingel (1997). "Puntos de tensión para la inestabilidad de la tensión en SPH". Int. J. Numer. Methods Eng . 40 (13): 2325–2341. Bibcode :1997IJNME..40.2325D. doi :10.1002/(sici)1097-0207(19970715)40:13<2325::aid-nme161>3.0.co;2-8.
^ PW Randles; LD Libersky (2000). "SPH normalizado con puntos de estrés". Int. J. Numer. Methods Eng . 48 (10): 1445–1462. Bibcode :2000IJNME..48.1445R. doi :10.1002/1097-0207(20000810)48:10<1445::aid-nme831>3.0.co;2-9.
^ T. Belytschko; Y. Guo; WK Liu; SP Xiao (2000). "Un análisis de estabilidad unificado de métodos de partículas sin malla". Int. J. Numer. Methods Eng . 48 (9): 1359–1400. Bibcode :2000IJNME..48.1359B. doi :10.1002/1097-0207(20000730)48:9<1359::aid-nme829>3.0.co;2-u.
^ J. Bonet; S. Kulasegaram (2000). "Corrección y estabilización de métodos de hidrodinámica de partículas lisas con aplicaciones en simulaciones de conformado de metales". Int. J. Numer. Methods Eng . 47 (6): 1189–1214. Bibcode :2000IJNME..47.1189B. doi :10.1002/(sici)1097-0207(20000228)47:6<1189::aid-nme830>3.0.co;2-i.
^ WG Hoover; CG Hoover (2001). "Recetas basadas en SPAM para simulaciones de continuo". Computing in Science and Engineering . 3 (2): 78–85. Bibcode :2001CSE.....3b..78H. doi : 10.1109/5992.909007 .
^ T. Rabczuk; J. Eibl; L. Stempniewski (2003). "Simulación de la fragmentación del hormigón a alta velocidad utilizando SPH/MLSPH". Int. J. Numer. Methods Eng . 56 (10): 1421–1444. Bibcode :2003IJNME..56.1421R. doi :10.1002/nme.617. S2CID 119799557.
^ MI Herreros; M. Mabssout (2011). "Un esquema de discretización temporal en dos pasos utilizando el método SPH para la propagación de ondas de choque". Comput. Methods Appl. Mech. Engrg . 200 (21–22): 1833–1845. Bibcode :2011CMAME.200.1833H. doi :10.1016/j.cma.2011.02.006.
^ S. Yashiro; T. Okabe (2015). "Hidrodinámica de partículas suavizadas en un sistema de coordenadas generalizado con un modelo constitutivo de deformación finita". Int. J. Numer. Methods Eng . 103 (11): 781–797. Bibcode :2015IJNME.103..781Y. doi :10.1002/nme.4906. hdl : 2324/4476056 . S2CID 121935150.
^ NJ Quinlan; M. Basa; M. Lastiwka (2006). "Error de truncamiento en métodos de partículas sin malla" (PDF) . Revista internacional de métodos numéricos en ingeniería . 66 (13): 2064–2085. Bibcode :2006IJNME..66.2064Q. doi :10.1002/nme.1617. hdl : 10379/1170 . S2CID 53701216.
^ H. Wendland (1995). "Funciones radiales definidas positivas, polinómicas por partes y con soporte compacto de grado mínimo". Avances en Matemática Computacional . 4 (4): 389–396. doi :10.1007/BF02123482. S2CID 36452865.
^ A. Mayrhofer; BD Rogers; D. Violeau; M. Ferrand (2013). "Investigación de flujos delimitados por paredes usando SPH y las condiciones de contorno de pared semianalíticas unificadas". Computer Physics Communications . 184 (11): 2515–2527. arXiv : 1304.3692 . Bibcode :2013CoPhC.184.2515M. CiteSeerX 10.1.1.770.4985 . doi :10.1016/j.cpc.2013.07.004. S2CID 35008128.
^ J. Bonet; TS Lok (1999). "Aspectos variacionales y de conservación del momento en formulaciones de hidrodinámica de partículas suavizadas". Métodos informáticos en ingeniería mecánica aplicada . 180 (1–2): 97–115. Bibcode :1999CMAME.180...97B. doi :10.1016/S0045-7825(99)00051-1.
^ JJ Monaghan (2005). "Hidrodinámica de partículas suavizadas". Informes sobre el progreso en física . 68 (8): 1703–1759. Bibcode :2005RPPh...68.1703M. doi :10.1088/0034-4885/68/8/R01. S2CID 5987481.
^ E. Hairer; C. Lubich; G. Wanner (2006). Integración numérica geométrica. Saltador. ISBN978-3-540-30666-5.
^ ab Andrea Colagrossi; Matteo Antuono; David Le Touzè (2009). "Consideraciones teóricas sobre el papel de la superficie libre en el modelo de hidrodinámica de partículas suavizadas". Physical Review E . 79 (5): 056701. Bibcode :2009PhRvE..79e6701C. doi :10.1103/PhysRevE.79.056701. PMID 19518587.
^ abc Bejamin Bouscasse; Andrea Colagrossi; Salvatore Marrone; Matteo Antuono (2013). "Interacción no lineal de ondas de agua con cuerpos flotantes en SPH". Revista de fluidos y estructuras . 42 : 112–129. Código Bibliográfico :2013JFS....42..112B. doi :10.1016/j.jfluidstructs.2013.05.010.
^ Fabricio Macià; Matteo Antuono; Leo M González; Andrea Colagrossi (2011). "Análisis teórico de la aplicación de la condición de contorno de no deslizamiento en métodos SPH". Progress of Theoretical Physics . 125 (6): 1091–1121. Bibcode :2011PThPh.125.1091M. doi : 10.1143/PTP.125.1091 .
^ Jose Luis Cercos-Pita; Matteo Antuono; Andrea Colagrossi; Antonio Souto (2017). "Conservación de energía SPH para interacciones fluido--sólido". Métodos informáticos en mecánica aplicada e ingeniería . 317 : 771–791. Bibcode :2017CMAME.317..771C. doi :10.1016/j.cma.2016.12.037.
^ J. Campbell; R. Vignjevic; L. Libersky (2000). "Un algoritmo de contacto para la hidrodinámica de partículas suavizadas". Métodos informáticos en mecánica aplicada e ingeniería . 184 (1): 49–65. Bibcode :2000CMAME.184...49C. doi :10.1016/S0045-7825(99)00442-9.
^ abc M. Ferrand, DR Laurence, BD Rogers, D. Violeau, C. Kassiotis (2013). "Condiciones de contorno de pared semianalíticas unificadas para flujos no viscosos, laminares o turbulentos en el método SPH sin malla". Revista internacional de métodos numéricos en fluidos . 71 (4). Int. J. Numer. Meth. Fluids: 446–472. Bibcode :2013IJNMF..71..446F. doi :10.1002/fld.3666. S2CID 124465492.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ HR Cole (1948). Explosiones submarinas. Princeton, Nueva Jersey: Princeton University Press.
^ abc D. Molteni, A. Colagrossi (2009). "Un procedimiento simple para mejorar la evaluación de la presión en el contexto hidrodinámico utilizando el SPH". Computer Physics Communications . 180 (6): 861–872. Bibcode :2009CoPhC.180..861M. doi :10.1016/j.cpc.2008.12.004.
^ ab Colagrossi, Andrea; Landrini, Maurizio (2003). "Simulación numérica de flujos interfaciales mediante hidrodinámica de partículas suavizadas". Journal of Computational Physics . 191 (2): 448–475. Bibcode :2003JCoPh.191..448C. doi :10.1016/S0021-9991(03)00324-3.
^ Randall J. LeVeque (2007). Métodos de diferencias finitas para ecuaciones diferenciales ordinarias y parciales: problemas de estado estacionario y dependientes del tiempo . Siam.
^ A. Ferrari, M. Dumbser, E. Toro, A. Armanini (2009). "Un nuevo esquema SPH paralelo 3D para flujos de superficie libre". Computadoras y fluidos . 38 (6). Elsevier: 1203–1217. doi :10.1016/j.compfluid.2008.11.012.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ Fatehi, R y Manzari, MT (2011). "Un remedio para las oscilaciones numéricas en la hidrodinámica de partículas suavizadas débilmente compresibles". Revista internacional de métodos numéricos en fluidos . 67 (9). Biblioteca Wiley Online: 1100–1114. Bibcode :2011IJNMF..67.1100F. doi :10.1002/fld.2406. S2CID 121381641.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ M. Antuono, A. Colagrossi, S. Marrone, D. Molteni (2010). "Flujos de superficie libre resueltos mediante esquemas SPH con términos numéricos difusivos". Computer Physics Communications . 181 (3). Elsevier: 532–549. Bibcode :2010CoPhC.181..532A. doi :10.1016/j.cpc.2009.11.002.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ M. Antuono, A. Colagrossi, S. Marrone (2012). "Términos difusivos numéricos en esquemas SPH débilmente compresibles". Computer Physics Communications . 183 (12). Elsevier: 2570–2580. Bibcode :2012CoPhC.183.2570A. doi :10.1016/j.cpc.2012.07.006.{{cite journal}}: CS1 maint: multiple names: authors list (link)
^ Antuono Matteo y Marrone S y Colagrossi A y Bouscasse B (2015). "Balance de energía en el esquema δ-SPH". Métodos informáticos en mecánica aplicada e ingeniería . 289 . Elsevier: 209–226. Bibcode :2015CMAME.289..209A. doi :10.1016/j.cma.2015.02.004.
^ JP Vila (1999). "Sobre métodos ponderados de partículas e hidrodinámica de partículas suaves". Modelos y métodos matemáticos en ciencias aplicadas . 9 (2). World Scientific: 161–209. doi :10.1142/S0218202599000117.
^ Marongiu Jean-Christophe y Leboeuf Francis y Caro Joëlle y Parkinson Etienne (2010). "Simulaciones de flujos de superficie libre en turbinas Pelton utilizando un método híbrido SPH-ALE" (PDF) . Journal of Hydraulic Research . 48 (S1). Taylor & Francis: 40–49. doi :10.1080/00221686.2010.9641244. S2CID 121493014.
^ De Leffe, Matthieu (2011). Modelización de gases viscosos por el método SPH en vista de aplicación de hidrodinámica naval . Tesis doctoral, Escuela Central de Nantes.
^ Chi Zhang y Xiangyu Hu y Nikolaus Adams (2017). "Un método SPH débilmente compresible basado en un solucionador de Riemann de baja disipación". Journal of Computational Physics . 335 : 605–620. Bibcode :2017JCoPh.335..605Z. doi :10.1016/j.jcp.2017.01.027.
^ Monaghan, Joseph J.; Gingold, Robert A. (1983). "Simulación de choque mediante el método de partículas". Journal of Computational Physics . 52 (2): 347–389. Bibcode :1983JCoPh..52..374M. doi :10.1016/0021-9991(83)90036-0.
^ JP Morris; PJ Fox; Y. Zhu (1997). "Modelado de flujos incompresibles con bajo número de Reynolds utilizando SPH". Revista de Física Computacional . 136 (1): 214–226. Código Bibliográfico :1997JCoPh.136..214M. doi :10.1006/jcph.1997.5776.
^ Marios D. Dikaiakos; Joachim Stadel, PKDGRAV The Parallel kD Tree Gravity Code , consultado el 1 de febrero de 2017
^ Vignjevic, Rade; Reveles, Juan R.; Campbell, James (2006). "SPH en un formalismo lagrangiano total". Modelado computacional en ingeniería y ciencias . 44 : 181–198.
^ Han, Luhui; Hu, Xiangyu (2018). "Modelado SPH de la interacción fluido-estructura". Revista de hidrodinámica . 30 (1): 62–69. Código Bibliográfico :2018JHyDy..30...62H. doi :10.1007/s42241-018-0006-9. S2CID 125369012.
^ Chi Zhang; Massoud Rezavand; Xiangyu Hu (2020). "Un método SPH de resolución múltiple para interacciones fluido-estructura". Revista de Física Computacional . 429 : 110028. arXiv : 1911.13255 . doi :10.1016/j.jcp.2020.110028. ISSN 0021-9991. S2CID 208513116.
Lectura adicional
Hoover, WG (2006); Mecánica aplicada de partículas suaves: el estado del arte, World Scientific.
Stellingwerf, RF; Wingate, California; "Modelado de impacto con SPH", Memorie della Societa Astronomia Italiana, vol. 65, pág. 1117 (1994).
Amada, T.; Imura, M.; Yasumuro, Y.; Manabe, Y.; y Chihara, K. (2004); "Simulación de fluidos basada en partículas en GPU", en Actas del taller ACM sobre informática de propósito general en procesadores gráficos (agosto de 2004, Los Ángeles, California).
Desbrun, M.; y Cani, M.-P. (1996). "Partículas suavizadas: un nuevo paradigma para animar cuerpos altamente deformables" en Actas del Taller Eurographics sobre Animación y Simulación por Computadora (agosto de 1996, Poitiers, Francia).
Hegeman, K.; Carr, NA; y Miller, GSP; "Simulación de fluidos basada en partículas en la GPU", en Actas de la Conferencia Internacional sobre Ciencias Computacionales (Reading, Reino Unido, mayo de 2006), Lecture Notes in Computer Science v. 3994/2006 (Springer-Verlag).
Kelager, M. (2006) Dinámica de fluidos lagrangiana utilizando hidrodinámica de partículas suavizadas , tesis de maestría, Univ. Copenhague.
Kolb, A.; y Cuntz, N. (2005); "Acoplamiento dinámico de partículas para simulación de fluidos basada en GPU", en Actas del 18º Simposio sobre Técnicas de Simulación (2005), págs. 722–727.
Liu, GR; y Liu, MB; Hidrodinámica de partículas suavizadas: un método de partículas sin malla , Singapur: World Scientific (2003).
Monaghan, Joseph J. (1992). "Hidrodinámica de partículas suavizadas", Revista anual de astronomía y astrofísica (1992). 30: 543–74.
Muller, M.; Charypar, D.; y Gross, M.; "Simulación de fluidos basada en partículas para aplicaciones interactivas", en Breen, D; y Lin, M. (eds.), Actas del Simposio Eurographics/SIGGRAPH sobre animación por computadora (2003).
Vesterlund, M.; Simulación y representación de un fluido viscoso utilizando hidrodinámica de partículas suavizadas , tesis de maestría, Universidad de Umea, Suecia.
Violeau, D.; Mecánica de fluidos y el método SPH , Oxford University Press (2012).
Enlaces externos
Primera gran simulación de formación estelar utilizando SPH
SPHERIC (Comunidad internacional de investigación e ingeniería de SPH)
ITVO es el sitio web del Observatorio Virtual Teórico Italiano creado para consultar una base de datos de archivos de simulación numérica.
La galería de imágenes SPHC muestra una amplia variedad de casos de prueba, validaciones experimentales y aplicaciones comerciales del código SPH SPHC.
Una derivación del modelo SPH a partir de las ecuaciones de Navier-Stokes
Software
Algodoo es un marco de simulación 2D para educación que utiliza SPH
AQUAgpusph es la SPH gratuita (GPLv3) de los investigadores, por los investigadores, para los investigadores
Dive Solutions es un software de ingeniería SPH comercial basado en la web para fines de CFD.
DualSPHysics es un código SPH en su mayor parte de código abierto basado en SPHysics y que utiliza computación por GPU. Los componentes de código abierto están disponibles bajo la licencia LGPL.
FLUIDS v.1 es una implementación de SPH 3D en tiempo real, simple y de código abierto (Zlib) en C++ para líquidos para CPU y GPU.
Fluidix es una API de simulación de partículas basada en GPU disponible en OneZero Software
GADGET [1] es un código disponible libremente ( GPL ) para simulaciones cosmológicas de N cuerpos/SPH
Simulador GPUSPH SPH con viscosidad (GPLv3)
Pasimodo es un paquete de programas para métodos de simulación basados en partículas, por ejemplo, SPH
LAMMPS es un código de dinámica molecular clásica, de código abierto y masivamente paralelo que puede realizar simulaciones SPH
Physics Abstraction Layer es un sistema de abstracción de código abierto que admite motores de física en tiempo real con soporte SPH
PreonLab es un software de ingeniería comercial desarrollado por FIFTY2 Technology que implementa un método SPH implícito
Punto es una herramienta de visualización disponible gratuitamente para simulaciones de partículas.
pysph Marco de código abierto para la hidrodinámica de partículas suavizadas en Python (nueva licencia BSD)
Py-SPHViewer Herramienta de visualización de Python de código abierto para simulaciones de hidrodinámica de partículas suavizadas. [1]
Solucionador SPH comercial RealFlow para la industria cinematográfica.
RheoCube es un producto SaaS comercial de Lorenz Research para el estudio y la predicción de la reología y la estabilidad de fluidos complejos.
SPHinXsys es una biblioteca SPH multifísica y multiresolución de código abierto. Proporciona API de C++ para simulación física precisa y tiene como objetivo modelar sistemas dinámicos industriales acoplados, incluidas dinámicas de fluidos, sólidos y multicuerpo, entre otros.
SPHysics es una implementación de SPH de código abierto en Fortran
SPLASH es una herramienta de visualización de código abierto (GPL) para simulaciones SPH
SYMPLER: un simulador de partículas simbólicas gratuito de la Universidad de Friburgo.
Nauticle es una herramienta computacional de propósito general para métodos numéricos basados en partículas.
NDYNAMICS es un software de simulación de fluidos comercial basado en SPH implícito desarrollado por CENTROID LAB, actualmente utilizado para aplicaciones de ingeniería nuclear, química y de inundaciones internas y externas.
^ Benítez-Llambay, Alejandro (28 de julio de 2015), "Py-Sphviewer: Py-Sphviewer V1.0.0", Zenodo , Bibcode :2015zndo.....21703B, doi :10.5281/zenodo.21703 , recuperado en 2022- 03-30