14 votos

Por qué no se pueden simular isotrópica de flujo de fluido en un cuadrado de celosía?

Hay métodos fáciles para las simulaciones de dispersión de gases en dos dimensiones. Si usted toma una gran plaza de celosía, cada celda de la que se supone que contiene a la mayoría de una molécula de gas, y mover las moléculas de la célula a la célula adyacente al azar, la gran escala de los resultados son en muchos sentidos, una buena simulación de gases de comportamiento. En particular, aunque son moléculas individualmente en movimiento sólo en los cuatro puntos cardinales, a gran escala de la simulación no muestran ningún prejuicio a estos cuatro direcciones. Si usted comienza con un cuadrado grande bulto de material, rápidamente se difunden en un círculo. Esta propiedad de la simulación se llama la isotropía.

Hace mucho tiempo, cuando yo estaba involucrado en la investigación en esta área me dijeron que esto funcionó muy bien para el flujo de gas de simulaciones, pero no para fluidos incompresibles. Si uno quería hacer una isotrópica de dos dimensiones de la simulación de fluido incompresible de flujo, se me dijo, que había que utilizar una red hexagonal en lugar de un cuadrado de celosía. Me dijeron que se podía deducir de las ecuaciones de Navier-Stokes que cualquier simulación de fluido incompresible en un cuadrado de celosía necesariamente debería ser anisotrópico.

Es esto correcto? Si es así, ¿cómo funciona la prueba de ir? Si el argumento es complejo, es posible obtener una idea intuitiva de por qué compresible e incompresible de los líquidos son diferentes en este sentido? ¿Qué es una referencia que incluye la totalidad de la prueba?

Anexo: a Pesar de la respuesta que he publicado a continuación, todavía no entiendo ninguno de los detalles. Yo estaría encantado de premio de la recompensa a alguien que pudiera explicarlo.

5voto

MJD Puntos 37705

Parece que esta podría ser la fuente original: Lattice gas autómata por las ecuaciones de Navier-Stokes, por Frisch et al. (PDF)

[Plaza de celosía] CH autómata es invariante bajo $\pi/2$ rotaciones. Un entramado de simetría es insuficiente para garantizar la isotropía del cuarto grado de tensor relacionadas con el impulso del flujo cuadrática en términos de la velocidad. ...
Observar que cuando el subyacente micromundo es de dos dimensiones y invariantes bajo el hexagonal de la rotación de grupo (múltiplos de $\pi/3$) el tensor de la $T$ es isotrópico...

Lamentablemente será un largo tiempo antes de que yo soy capaz de entender esto correctamente; yo preferiría ver una respuesta que explica lo que está pasando.

3voto

theog Puntos 585

Como usted dice que usted no está muy familiarizado con los tensores (y yo no me llamo a mí mismo un experto), voy a tratar de explicar todo en el lineal del lenguaje algebraico. Para citar el papel,

...para la CH modelo, el impulso de flujo tensor tiene la forma $$P_{\alpha\beta} = p\delta_{ab} + T_{\alpha\beta\gamma\epsilon}u_\gamma u_\epsilon + O(u^4).$$ ...El tensor de la $T$ es, por construcción, simétrica en ambos $(\alpha,\beta)$$(\gamma,\epsilon)$.

Aquí $T$ es "el cuarto grado de tensor relacionadas con el impulso del flujo cuadrática en términos de la velocidad". Podemos pensar en los de segundo orden tensores $P_{\alpha\beta}$ $u_\gamma u_\epsilon$ simétrica $2\times2$ matrices ($P_{\alpha\beta}$ debe ser simétrico para conservar el momento angular, mientras que$u_\gamma u_\epsilon$$uu^T$, lo que es claramente simétrica), y $T$ como una transformación lineal de asignación de la última a la primera.

Un tensor es isotrópico si lo ven de la misma bajo cualquier rotación del cuadro de coordenadas. En este caso, eso significa que si $T$ mapas matriz $X$ a una matriz de $Y$, y girar el marco de coordenadas de modo que $X$ hace $X'$ $Y$ hace $Y'$, $T$ los mapas de $X'$$Y'$. Formalmente, si una rotación transforma vectores como $v' = Rv$, entonces las matrices de transformación de la $A' = RAR^T$, y queremos que $T(RXR^T)=RT(X)R^T$.

Las cosas se vuelven más claras si hacemos una elección acertada de base para el espacio tridimensional de la simétrica $2\times2$ matrices, es decir, $$E_1 = \begin{bmatrix}1 & 0\\0 & 1\end{bmatrix},\quad E_2 = \begin{bmatrix}1 & 0\\0 & -1\end{bmatrix},\quad E_3 = \begin{bmatrix}0 & 1\\1 & 0\end{bmatrix}.$$ Entonces resulta que una rotación por un ángulo de $\theta$ transforma esta base como $$\begin{align} E_1' &= E_1, \\ E_2' &= E_2\cos2\theta + E_3\sin2\theta, \\ E_3' &= E_3\cos2\theta - E_2\sin2\theta, \end{align}$$ que también se puede ver por la visualización de ellos como funciones cuadráticas de la forma$x \mapsto x^TEx$$\mathbb R^2$. En otras palabras, una rotación transforma $E_2$ $E_3$ en cada uno de los otros, así que para $T$ a ser invariante bajo una transformación, que tiene que hacer "lo mismo" que tanto $E_2$$E_3$, en algún sentido. Pero si sólo mira las rotaciones de $\pi/2$, no se ve esta dependencia, debido a que, a continuación,$E_2' = -E_2$$E_3' = -E_3$. Así que un tensor de que el tratado $E_2$ $E_3$ diferente podría ser invariante bajo rotaciones de $\pi/2$, pero no iba a ser invariantes bajo rotaciones arbitrarias. Por ejemplo, podríamos optar $T$ tal que $$T(E_1) = T(E_2) = 0,\quad T(E_3) = E_3.$$ Puedes verificar que este tensor es invariante bajo una rotación de $\pi/2$, pero no bajo una rotación de $\pi/4$ (donde$E_2' = E_3$$E_3' = -E_2$), por lo que no es isotrópica.


De todos modos, dado que el título de la pregunta todavía se pregunta "¿por Qué no puedes simular isotrópica de flujo de fluido en un cuadrado de celosía?", permítanme reiterar mi comentario de que este es un problema con entramado de gas autómatas en particular, no con simulaciones de fluidos en la plaza de las mallas, en general. El problema con la isotropía viene cuando se intenta demostrar la convergencia del modelo de autómata celular para las ecuaciones de Navier-Stokes. Hasta donde yo soy consciente, usted puede simular isotrópica de flujo de fluido en un cuadrado de la cuadrícula a través de los métodos habituales de directamente de la discretización de las ecuaciones de Navier-Stokes y resolver numéricamente, sin pasar a través de un intermedio modelo de autómatas celulares.

0voto

ArK Puntos 5284

Tengo una existencia prueba que usted puede hacer una asintóticamente angularmente isotrópica tres dimensiones de autómatas celulares regla que funciona en un cartesiano de celosía. Crea planeadores que se pueden mover en (asintóticamente) cualquier dirección y en (asintóticamente) a cualquier velocidad. Sin embargo, la norma escrita no se manejan bien con las colisiones entre estos planeadores. Pero al menos demuestra que con un simple " regla puede producir angularmente isotrópica movimiento comportamiento.

https://docs.google.com/document/d/1ocynhGqTrHZ_-GK5EBZjSBAUzj5pPmRlFiYclf1l3-s/edit

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X