Ecuaciones de Navier-Stokes

Les **ecuaciones de Navier-Stokes** reciben el so nome de Claude-Louis Navier y George Gabriel Stokes. Trátase d'un conxuntu d'ecuaciones en derivaes parciales non lliniales que describen el movimientu d'un fluyíu. Estes ecuaciones gobiernen l'atmósfera terrestre, les corrientes oceániques y el fluxu alredor de vehículos o proyeutiles y, polo xeneral, cualquier fenómenu nel que s'arreyen fluyíos newtonianos.

Estes ecuaciones llógrense aplicando los principios de caltenimientu de la mecánica y la termodinámica a un volume fluyíu. Faciendo esto llógrase la llamada *formulación integral* de les ecuaciones. Pa llegar a la so formulación diferencial manipóliense aplicando ciertes considerancies, principalmente aquella na que los esfuercios tanxenciales guarden una rellación llinial col gradiente de velocidá (llei de mafa de Newton), llogrando d'esta manera la formulación diferencial que xeneralmente ye más útil pal resolución de los problemes que se planteguen na mecánica de fluyíos.

Como yá se dixo, les ecuaciones de Navier-Stokes son un conxuntu d'ecuaciones en derivaes parciales non lliniales. Nun se dispón d'una solución xeneral pa esti conxuntu d'ecuaciones, y salvo ciertos tipos de fluxu y situaciones bien concretes nun ye posible topar una solución analítica; polo qu'en munches ocasiones ye precisu recurrir al analís numbéricu pa determinar una solución averada. A la caña de la mecánica de fluyíos que s'ocupa del llogru d'estes soluciones por aciu métodos numbéricos denominar dinámica de fluyíos computacional (CFD, de la so acrónimu anglosaxón *Computational Fluyíi Dynamics*).

## Conceutos previos

### Derivada sustancial o material

Por cuenta de que xeneralmente adoptamos la descripción euleriana la derivada ordinaria {\partial\phi}/{\partial t} yá nun representa tola variación por unidá de tiempu d'una determinada propiedá del fluyíu (o magnitú fluyida) {\phi} siguiendo a la partícula fluyida. Esto debe al movimientu del fluyíu. Pa reflexar esta variación vamos usar la derivada sustancial (o derivada siguiendo a la partícula fluyida). La **derivada sustancial** o **derivada material** defínese como l'operador:
{=}\ \frac{\partial(\star )}{\partial t} + \mathbf{v}\cdot\nabla (\star )
||left}}
Onde \mathbf{v} ye la velocidá del fluyíu. El primer términu representa la variación de la propiedá nun puntu fixu del espaciu y por ello denominar derivada local, ente que'l segundu representa la variación de la propiedá acomuñáu al cambéu de posición de la partícula fluyida, y denominar derivada convectiva. Este ye'l procedimientu que sigue José Echegaray pa demostrar la derivada material. Vease una demostración de cómo llegar a una derivada material. Tomando les coordenaes de Euler como:
+v_y(x,y,z,t)\hat{\mathbf{j}}+v_z(x,y,z,t)\hat{\mathbf{k}}.
||left}}
Vamos Calcular l'aceleración pa estes coordenaes:
{dt}=\frac{dv_x}{dt}\hat{\mathbf{i}}+\frac{dv_y}{dt}\hat{\mathbf{j}}+\frac{dv_z}{dt}\hat{\mathbf{k}}
||left}}
Desenvolvemos cada derivada total de cada componente, asina vamos poder siguir un desenvolvimientu fácil de recordar:

Si súmase términu a términu y sácase factor común, puede llograse:
{Dt}=\frac{\partial (v_xi+v_yj+v_zk)}{\partial t}+v_x\frac{\partial (v_xi+v_yj+v_zk)}{\partial x}+v_y\frac{\partial (v_xi+v_yj+v_zk)}{\partial y}+v_z\frac{\partial (v_xi+v_yj+v_zk)}{\partial z}
:\frac{D\mathbf{v}}{Dt} = \frac{\partial \mathbf{v}}{\partial t} + 
[v_x\frac{\partial}{\partial x} + v_y\frac{\partial}{\partial y} + v_z\frac{\partial}{\partial z}]\mathbf{v} =
\frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot\nabla) \mathbf{v}
||left}}
Vemos que la parte de les derivaes parciales espaciales pueden escribise como: \mathbf{v}\cdot\nabla

Si agora sustituyimos velocidá por (\star) llogramos formalmente la espresión de la derivada material:

### Teorema del tresporte de Reynolds

Si la derivada sustancial dexa calcular la variación d'una magnitú fluyida amestada a una partícula fluyida, el teorema del tresporte de Reynolds va dexar calcular la variación d'una magnitú fluyida estensiva amestada a un volume fluyíu. Esiste por tanto una analoxía ente dambos conceutos, pos una partícula fluyida nun ye más qu'un volume fluyíu infinitesimal. Na so forma xeneral el teorema del tresporte de Reynolds esprésase como:

onde \phi ye la magnitú fluyida estensiva definida por unidá de volume (una magnitú estensiva por unidá de volume ye una magnitú intensiva), V_f ye un volume fluyíu, V_c ye un volume de control que coincide con V_f nel intre t, S_c la superficie de dichu volume de control, \mathbf{v} la velocidá del fluyíu y \mathbf{v_c} la velocidá de la superficie de control.

El segundu términu del miembru derechu representa'l fluxu convectivo de la magnitú fluyida estensiva al traviés de la superficie de control que llinda'l volume de control. Defínese'l fluxu convectivo d'una magnitú fluyida estensiva al traviés d'una superficie de control como la cantidá de dicha magnitú que, tresportada pol fluyíu, traviesa la superficie de control na unidá de tiempu.

Espresáu en términos coloquiales puede dicise que'l teorema del tresporte de Reynolds vien dicir que la variación d'una propiedá estensiva nun volume fluyíu, ye igual a la variación de dicha propiedá nel interior d'esi volume más la cantidá de dicha propiedá que traviesa la superficie del volume.

### Teorema de la diverxencia

El teorema de la diverxencia (o teorema de Gauss) dexa, so ciertes hipótesis, tresformar integrales de superficie n'integrales de volume ( y viceversa). Nel casu particular de tres dimensiones podemos espresalo como:

\iiint\limits_V\left(\nabla\cdot\mathbf{F}\right)dV=\iint\limits_{\partial V}\mathbf{F\cdot n}\; dS

## Les ecuaciones de Navier-Stokes

Esta espresión representa'l principiu de caltenimientu del momentu llinial aplicada a un fluyíu xeneral:

:\rho\frac{Du_i}{Dt}=\rho F_i-\frac{\partial P}{\partial x_i}+\frac{\partial}{\partial x_j}\left[
2\mu\left(y_{ij}-\Delta\delta_{ij}/3\right)\right].

La llei de caltenimientu de la masa escríbese:

:\frac{\partial\rho}{\partial t}+\frac{\partial\rho o_i}{\partial x_i}=0

Nestes ecuaciones ρ representa la densidá, *o<sub>i</sub>* (i = 1,2,3) les componentes cartesianes de la velocidá, *F<sub>i</sub>* el campu d'aceleraciones creáu poles fuercies aplicaes sobre'l cuerpu, como la gravedá, *P* la presión del fluyíu, y μ la mafa dinámica.

:y_{ij}=\frac{1}{2} \left(\frac{\partial o_i}{\partial x_j}+\frac{\partial o_j}{\partial x_i}\right)

onde Δ = *y<sub>ii</sub>* ye la diverxencia del fluyíu y δij la delta de Kronecker. D / Dt ye la derivada total o derivada material temporal siguiendo'l fluyíu:

:\frac{D}{Dt}(\cdot) \equiv \frac{\partial(\cdot)}{\partial t} + (\mathbf{v}\cdot\nabla)(\cdot)

La **non linealidad** de les ecuaciones débese precisamente al términu rellacionáu cola derivada total. Cuando μ ye uniforme sobremanera'l fluyíu les ecuaciones de fluyíu simplificar de la manera siguiente:

O en forma vectorial:
{Dt}=
\rho \mathbf{k} - \boldsymbol{\nabla}P +\mu \left( \frac{1}{3} \boldsymbol{\nabla}(\boldsymbol{\nabla} \cdot \mathbf{o}) + \boldsymbol{\nabla}^{2}\mathbf{o} \right)
||left}}

## Casos particulares
Para fluyíos de mafa nula, ye dicir cuando μ = 0, les ecuaciones resultantes denominar ecuaciones de Euler que s'utilicen nel estudiu de fluyíos compresibles y n'ondes de choque:

Per otra parte si considérase un fluyíu mafosu pero incompresible, entós puede ser considerada constante (como nun líquidu) y les ecuaciones resulten ser:
:\rho \left({\partial v_x \over \partial t}+ v_x {\partial v_x \over \partial x}+ v_y {\partial v_x \over \partial y}+ v_z {\partial v_x \over \partial z}\right)= \mu \left[{\partial^2 v_x \over \partial x^2}+{\partial^2 v_x \over \partial y^2}+{\partial^2 v_x \over \partial z^2}\right]-{\partial P \over \partial x} +\rho g_x 

:\rho \left({\partial v_y \over \partial t}+ v_x {\partial v_y \over \partial x}+ v_y {\partial v_y \over \partial y}+ v_z {\partial v_y \over \partial z}\right)= \mu \left[{\partial^2 v_y \over \partial x^2}+{\partial^2 v_y \over \partial y^2}+{\partial^2 v_y \over \partial z^2}\right]-{\partial P \over \partial y} +\rho g_y 

:\rho \left({\partial v_z \over \partial t}+ v_x {\partial v_z \over \partial x}+ v_y {\partial v_z \over \partial y}+ v_z {\partial v_z \over \partial z}\right)= \mu \left[{\partial^2 v_z \over \partial x^2}+{\partial^2 v_z \over \partial y^2}+{\partial^2 v_z \over \partial z^2}\right]-{\partial P \over \partial z} +\rho g_z

y la ecuación de continuidá adquier la forma siguiente:

:{\partial v_x \over \partial x}+{\partial v_y \over \partial y}+{\partial v_z \over \partial z}=0

 

## Ver tamién
* Númberu de Reynolds
* Númberu de Mach
* Ecuaciones de Saint-Venant en 1D

## Referencies

### Bibliografía
*
*
*
*
*
*
* V. Girault and P.A. Raviart. *Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms.* Springer Series in Computational Mathematics. Springer-Verlag, 1986.
*
## Enllaces esternos
* 
*