Navier-Stokes incompresible no dimensionalizante

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

t tu i + tu j j tu i = ρ 1 i pag + v j j tu i , j tu j = 0.

No dimensionalizar usando las escalas [ X , y , z ] = L , [ tu ] = tu , [ pag ] = ρ tu 2 . Más

[ t ] = L / tu , tiempo convectivo

[ t ] = L 2 / v , tiempo difusivo

Da como resultado dos ecuaciones diferentes:

t tu i + tu j j tu i = i pag + R mi 1 j j tu i , utilizando la escala de tiempo convectiva

R mi 1 t tu i + tu j j tu i = i pag + R mi 1 j j tu i , utilizando la escala de tiempo difusiva

Aquí R mi = tu L v es el número de Reynolds.

Preguntas

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

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

3) ¿La amplitud de ambos espacios de parámetros es la misma para ambas ecuaciones?

Estoy buscando una visión metódica de responder estas preguntas ya que estoy interesado en cómo la comprensión de las respuestas se extiende a un problema más complicado.

Agradezco cualquier ayuda!

Respuestas (2)

Las dos ecuaciones solo difieren en un cambio de escala trivial de la coordenada de tiempo y, por lo tanto, son equivalentes. Si tu ( X , t ) es una solución de la primera ecuación, entonces tu ( X , R mi t ) es una solución de la segunda. Re es físicamente relevante porque gobierna la importancia relativa de los términos advectivo y disipativo. multiplicando t tu por una constante es solo un cambio de escala y no tiene significado físico.

Estoy de acuerdo con todo lo que ha dicho, pero me gustaría preguntarle acerca de su segunda y última oración. Creo que tu segunda oración no es muy general. Supongamos que hay más términos en la ecuación NS. No tengo claro cómo su notación de comparar las soluciones tu ( X , t ) y tu ( X , R mi t ) funcionaría comparando, por ejemplo, el término de advección y una fuerza electromagnética de Lorentz. La última oración que escribió parece ser un enfoque más general, pero ¿cómo puede volver a escalar un solo término en una ecuación mientras mantiene los mismos parámetros (Re, por ejemplo)?
@Charlie Si solo hay uno t término en el PDE, luego multiplicar esos términos por una constante solo cambia la unidad de tiempo y, por lo tanto, no tiene significado físico.
Si agrega 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 (wikipedia tiene una tabla completa, en.wikipedia.org/wiki/Dimensionless_numbers_in_fluid_mechanics ) pero eso no tiene nada que ver con su capacidad para cambiar la escala de t.
Al volver a escalar la coordenada de tiempo, cambiará la escala de tiempo 'convectiva' a una escala de tiempo 'difusiva', de modo que mientras las ecuaciones son matemáticamente equivalentes, físicamente no son equivalentes. Traté de explicar esto en mi respuesta, pero tal vez esto no estaba claro. @Charlie: en mi respuesta se muestra un ejemplo de cómo cambiar la escala de un solo término, donde se cambió la escala de la presión.
¿Cómo puede cambiar algo la solución cambiando la unidad de tiempo de [t]=1 segundo a [t]=17 min?
@Thomas: Hay muchas opciones de escalas para cada término adicional cuando aparecen proporciones más adimensionales. por ejemplo, considere la advección, la difusión y la fuerza EM. Un conjunto de escalas produce coeficientes 1 , R mi 1 , H a 2 / ( R mi R mi metro ) , otro juego de escalas produce 1 , R mi 1 , R mi metro . ¿Cómo se pueden relacionar las soluciones obtenidas a partir de estas diferentes escalas? Esto ya no es un simple tu ( X , t ) y tu ( X , R mi t ) relación.

... el propósito de la no dimensionalización es "colapsar" las soluciones en una curva para que el espacio de la solución se pueda explorar con menos parámetros.

Ese es un 'propósito' de la no dimensionalizació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 R mi 1 (viscosidad dominada) y R mi 1 (dominado por convección).

Las escalas de longitud característica son las escalas del sistema que caracterizan de manera única la dinámica, de modo que todos los términos de las ecuaciones (incluidas las condiciones iniciales/de contorno) son O ( 1 ) o más pequeño. Has elegido un conjunto de escalas con las que adimensionalizas la ecuación, pero ¿son escalas características? Bueno, vamos a ponerlo a prueba.

Primero, considere el régimen dominado por convección para el cual determina:

t tu i + tu j j tu i = i pag + R mi 1 j j tu i , utilizando la escala de tiempo convectiva

En este régimen, R mi 1 y asumiendo que la escala se hizo correctamente, es decir, todas las variables no dimensionales son O ( 1 ) , entonces vemos que todos los términos son O ( 1 ) o menor (ya que O ( R mi 1 ) O ( 1 ) ). Esto indica que el escalado se realizó correctamente y que las escalas son las escalas características en este régimen. Dado que el término viscoso es insignificantemente pequeño en comparación con el término convectivo, podemos ignorarlo por completo. Tenga en cuenta que el término de presión es del mismo orden que el término convectivo; el gradiente de presión siempre debe ser del mismo orden que el término dominante en la ecuación para equilibrar los términos.

Consideremos ahora el régimen viscoso donde R mi 1 , encuentras (después de dividir por R mi 1 ):

t tu i + R mi tu j j tu i = R mi i pag + j j tu i , utilizando la escala de tiempo difusiva

Vemos aquí que de nuevo todos los términos son O ( 1 ) o menor (ya que O ( R mi ) O ( 1 ) ). En este régimen, podemos despreciar los términos convectivos ya que los términos viscosos son dominantes. Sin embargo, aparentemente el gradiente de presión también es insignificante, lo 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 requerimos un reescalado de la escala de presión. vamos a reescalar [ pag ] = R mi α [ pag ] dónde α es una constante por determinar. Sustituyendo en las ecuaciones se obtiene:

t tu i + R mi tu j j tu i = R mi 1 + α i pag + j j tu i , utilizando la escala de tiempo difusiva

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

t tu i + R mi tu j j tu i = i pag + j j tu i , utilizando la escala de tiempo difusiva

y la escala de presión se reescaló a [ pag ] = R mi 1 [ pag ] = m tu L que es claramente una escala de presión viscosa como es evidente por la presencia de la viscosidad m . Esto contrasta con la escala de presión original. [ pag ] = ρ tu 2 que claramente no contiene una viscosidad y solo era apropiado en el régimen convectivo R mi 1 . Por lo tanto, podemos referirnos a ella como una escala de presión convectiva.

Para concluir respondiendo a sus preguntas:

  1. Debido a que las ecuaciones son válidas solo para dos regímenes diferentes establecidos por su elección de escalas, R mi 1 y R mi 1 .
  2. La única solución 'relacionada' se encontrará para R mi = 1 , de lo contrario, las soluciones serán completamente diferentes. De hecho, las soluciones analíticas generalmente solo son posibles para R mi 1 porque las ecuaciones se vuelven lineales. Para R mi 1 , las ecuaciones son altamente no lineales, por lo que solo es numéricamente posible obtener soluciones.
  3. No, ver respuesta a 2.
Aprecio el enfoque semi-metódico (usando una ecuación 1 + α = 0 para determinar una escala adecuada). Sin embargo, esta declaración: "... el gradiente de presión siempre debe ser del mismo orden que el término dominante en la ecuación para equilibrar los términos". no es, en general, cierto. Un contraejemplo 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 usan las escalas incorrectas. Esto no me sienta bien. Como señaló @Thomas, ambas soluciones pueden estar relacionadas para este caso simple.
@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 debe ser necesariamente del mismo orden que los términos viscosos o inerciales, respectivamente. Claramente, 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 es solo un cambio de escala de la coordenada de tiempo. Sin embargo, si su objetivo es determinar las escalas que dominan, es evidente que no obtendrá la solución correcta.