4 votos

La no dimensionalización de Navier-Stokes incompresible

Tengo una pregunta sobre la adimensionalización de las ecuaciones incompresibles de Navier-Stokes (NS). Tengo entendido que el propósito de la adimensionalización es "colapsar" las soluciones en una curva para que el espacio de soluciones pueda con menos parámetros. También puede dar una idea de la física. Consideremos la ecuación NS dimensional incompresible:

$\partial_t u_i + u_j \partial_j u_i = -\rho^{-1}\partial_i p + \nu \partial_{jj} u_i, \qquad \partial_j u_j = 0.$

La no-dimensión utilizando las escalas $ [x,y,z] = L, [u] = U, [p] = \rho U^2$ . Más

$[t] = L/U, \qquad \qquad \text{convective time}$

$[t] = L^2/\nu, \qquad \qquad \text{diffusive time}$

Resultados en dos ecuaciones diferentes:

$\partial_t u_i + u_j \partial_j u_i = -\partial_i p + Re^{-1} \partial_{jj} u_i, \qquad \qquad \qquad \text{using convective time scale}$

$Re^{-1}\partial_t u_i + u_j \partial_j u_i = - \partial_i p + Re^{-1}\partial_{jj} u_i, \qquad \qquad \text{using diffusive time scale}$

Aquí $Re = \frac{UL}{\nu}$ es el número de Reynolds.

Preguntas

1) ¿Por qué estas dos ecuaciones no son iguales dado que $Re$ se define exactamente igual entre las dos ecuaciones?

2) ¿Se pueden relacionar de alguna manera las soluciones de cada ecuación?

3) ¿Es la extensión de ambos espacios de parámetros la misma para ambas ecuaciones?

Busco una visión metódica para responder a estas preguntas, ya que me interesa que la comprensión de las respuestas se extienda a un problema más complicado.

Agradezco cualquier ayuda.

1voto

Liza Puntos 11

Las dos ecuaciones sólo se diferencian por un reescalado trivial de la coordenada temporal, por lo que son equivalentes. Si $u(x,t)$ es una solución de la primera ecuación, entonces $u(x,{\it Re}\, t)$ es una solución de la segunda. Re es físicamente relevante porque gobierna la importancia relativa de los términos advectivos y disipativos. Multiplicando $\partial_t u$ por una constante es sólo un reescalado, y no tiene ningún significado físico.

0 votos

Estoy de acuerdo con todo lo que has dicho, pero me gustaría preguntar sobre tu 2ª y última frase. Creo que tu 2ª frase no es muy general. Supongamos que hay más términos en la ecuación NS. No me queda claro cómo tu notación de comparar las soluciones $u(x,t)$ y $u(x,Re t)$ funcionaría comparando, por ejemplo, el término de advección y una fuerza electromagnética de Lorentz. La última frase que has escrito parece un enfoque más general, pero ¿cómo se puede reescalar un solo término de una ecuación manteniendo los mismos parámetros (Re, por ejemplo)?

0 votos

@Charlie Si sólo hay uno $\partial_t$ en la EDP, entonces multiplicar esos términos por una constante sólo cambia la unidad de tiempo, y por lo tanto no tiene ningún significado físico.

0 votos

Si se añaden otros términos a la ecuación de Navier-Stokes, relacionados con el flujo de calor, términos MHD, etc., entonces, por supuesto, aparecen relaciones adimensionales adicionales (la wikipedia tiene toda una tabla, es.wikipedia.org/wiki/Números_sin_dimensiones_en_mecánica_de_fluidos ) pero eso no tiene nada que ver con su capacidad de reescalar t.

-1voto

akostadinov Puntos 380

... el propósito de la adimensionalización es "colapsar" las soluciones en una curva para que el espacio de soluciones pueda ser explorado con menos parámetros.

Ese es uno de los "propósitos" de la adimensionalización; otros dos son la identificación de las escalas características de longitud, velocidad, presión, etc. y el análisis de la ecuación bajo diferentes regímenes, en este caso $\mathrm{Re}\ll1$ (dominado por la viscosidad) y $\mathrm{Re}\gg1$ (dominado por la convección).

Las escalas de longitudes características son las escalas del sistema que caracterizan de forma única la dinámica, de forma que todos los términos de las ecuaciones (incluidas las condiciones iniciales/de frontera) son $O(1)$ o más pequeño. Ha elegido un conjunto de escalas con las que adimensionaliza la ecuación, pero ¿son escalas características? Pues vamos a ponerlo a prueba.

En primer lugar, considere el régimen dominado por la convección para el que se determina:

$$\partial_t u_i + u_j \partial_j u_i = -\partial_i p + \mathrm{Re}^{-1} \partial_{jj} u_i, \qquad \qquad \qquad \text{using convective time scale}$$

En este régimen, $\mathrm{Re}\gg1$ y suponiendo que el escalado se haya hecho correctamente, es decir, que todas las variables no dimensionales sean $O(1)$ entonces vemos que todos los términos son $O(1)$ o menor (ya que $O(\mathrm{Re}^{-1})\ll O(1)$ ). Esto indica que el escalado se hizo correctamente y que las escalas son las características en este régimen. Dado que el término viscoso es despreciable en comparación con el término convectivo, se nos permite prescindir completamente de él. Nótese que el término de presión es del mismo orden que el término convectivo; el gradiente de presión debe ser siempre del mismo orden que el término dominante en la ecuación para equilibrar los términos.

Consideremos ahora el régimen viscoso donde $\mathrm{Re}\ll1$ se encuentra (después de dividir por $\mathrm{Re}^{-1}$ ):

$$\partial_t u_i + \mathrm{Re} u_j \partial_j u_i = - \mathrm{Re} \partial_i p + \partial_{jj} u_i, \qquad \qquad \text{using diffusive time scale}$$

Vemos aquí que de nuevo todos los términos son $O(1)$ o menor (ya que $O(\mathrm{Re})\ll O(1)$ ). En este régimen, podemos despreciar los términos convectivos ya que los términos viscosos son dominantes. Sin embargo, aparentemente también es despreciable el gradiente de presión que, de hecho, debería ser del mismo orden que los términos viscosos. Esto indica que el escalamiento propuesto en este régimen no es correcto.

Para resolver este problema es necesario reescalar la escala de presión. Reescalemos $[p]^* = \mathrm{Re}^{\alpha}[p]$ donde $\alpha$ es una constante a determinar. Sustituyendo en las ecuaciones se obtiene:

$$\partial_t u_i + \mathrm{Re} u_j \partial_j u_i = - \mathrm{Re}^{1+\alpha} \partial_i p^* + \partial_{jj} u_i, \qquad \qquad \text{using diffusive time scale}$$

Para que el gradiente de presión sea del mismo orden que los términos viscosos, requerimos $1+\alpha=0$ o $\alpha=-1$ . Entonces obtenemos:

$$\partial_t u_i + \mathrm{Re} u_j \partial_j u_i = - \partial_i p^* + \partial_{jj} u_i, \qquad \qquad \text{using diffusive time scale}$$

y la escala de presión se reescaló a $[p]^* = \mathrm{Re}^{-1}[p] = \frac{\mu U}{L}$ que es claramente una escala de presión viscosa como se desprende de la presencia de la viscosidad $\mu$ . Esto contrasta con la escala de presión original $[p] = \rho U^2$ que claramente no contiene una viscosidad y sólo era apropiado en el régimen convectivo $\mathrm{Re}\gg1$ . Por lo tanto, podemos referirnos a ella como una escala de presión convectiva.

Para terminar de responder a sus preguntas:

  1. Porque las ecuaciones sólo son válidas para dos regímenes diferentes según su elección de escalas, $\mathrm{Re}\ll1$ y $\mathrm{Re}\gg1$ .
  2. La única solución "relacionada" se encontrará para $\mathrm{Re}=1$ De lo contrario, las soluciones serán completamente diferentes. De hecho, las soluciones analíticas generalmente sólo son posibles para $\mathrm{Re}\ll1$ porque las ecuaciones se vuelven lineales. Para $\mathrm{Re}\gg1$ Las ecuaciones son altamente no lineales, por lo que sólo es posible obtener soluciones numéricas.
  3. No, véase la respuesta a la pregunta 2.

0 votos

Aprecio el enfoque semi-metódico (utilizando una ecuación $1+\alpha = 0$ para determinar una escala adecuada). Sin embargo, esta afirmación: "... el gradiente de presión debe ser siempre del mismo orden que el término dominante en la ecuación para equilibrar los términos" no es, en general, cierta. Un ejemplo contrario es el primer problema de Stokes, no hay gradiente de presión. El otro problema que tengo es que pareces sugerir que una solución es incorrecta si se utilizan las escalas equivocadas. Esto no me parece bien. Como ha señalado @Thomas, ambas soluciones pueden estar relacionadas para este caso sencillo.

0 votos

@Charlie - Si el gradiente de presión es lo que impulsa el flujo (ya sea viscoso o inercial), entonces siempre el gradiente de presión tiene que ser necesariamente del mismo orden que los términos viscosos o inerciales respectivamente. Evidentemente, esto no se aplica al primer problema de Stokes. Independientemente de la elección de las escalas, las dos ecuaciones darán las mismas soluciones en términos dimensionales, ya que sólo se trata de un reescalado de la coordenada temporal. Sin embargo, si tu objetivo es determinar las escalas que dominan, es evidente que no obtendrás la solución correcta.

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