Categoría:TC22/23

</nowiki></pre>

1 INTRODUCCIÓN

Con este proyecto buscamos analizar y visualizar tanto campos escalares como vectoriales (vistos durante en grado) en fluidos. En nuestro caso tendremos un fluido incompresible cuyo caudal fluctúa a través de un canal de paredes horizontales rectas y a su vez sometido a tres campos, dos escalares: presión y temperatura, y uno vectorial, velocidad.

La cinemática de los fluidos trata del movimiento de los mismos sin considerar las causas que lo forman. Se especializa en las trayectorias, velocidades y aceleraciones. A su vez, sabemos que un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición.

Para la realización del trabajo nos apoyamos el programa informático Matlab el cual nos ayudo a visualizar el comportamiento del fluido a través de gráficos.

2 SUPERFICIE DE TRABAJO

La superficie en la que nos vamos a basar va a ser [0,8]x[0,1] en {j,k}, por lo que trabajamos en el plano x=0, entonces definimos la y en [0,8] y la z en [0,1], aunque el eje z lo definimos como [-1,2]. Para obsérvalo vamos a recurrir a octave:

y=0:0.05:8;
z=0:0.05:1;
[yy,zz]=meshgrid(y,z);
figure(1);
mesh(yy,zz,0*yy);
grid on
axis([0,8,-1,2]);
view(2);

La superficie de trabajo es la que se muestra en la siguiente imagen:

Superficie de trabajo. Mallado

3 ECUACIÓN DE NAVIER-STOKES ESTACIONARIA

En dinámica de fluidos , las ecuaciones de Navier-Stokes son ecuaciones que describen el movimiento tridimensional de sustancias fluidas viscosas.

Se pueden usar para predecir el clima, las corrientes oceánicas, el flujo de agua en una tubería o en un reactor, el estudio del flujo sanguíneo y muchas otras cosas como el diseño de submarinos.

En este epígrafe, en concreto, trataremos de demostrar que el campo de velocidad y el campo de presión al que está sometido el fluido verifica la ecuación estacionaria de Navier-Stokes, esto supondría que el fluido es incompresible pues estas ecuaciones determinan el comportamiento de los fluidos newtonianos, esto es, un fluido cuya resistencia a deformaciones puede considerarse constante en el tiempo.

La ecuación estacionaria de Navier-Stokes es la siguiente:

Trabajando en componentes tenemos que:

[math](\vec{u} \cdot \nabla) \vec{u}=u_j\frac{\partial u_i}{\partial x_j}\vec{e_i} [/math] , [math]\Delta\vec {u}=u_i\vec{e_i}[/math], [math]\vec{u}=u_i\vec{e_i}[/math]

En nuestro caso, [math]\vec{u}(y,z)=f(z)\vec{j} [/math] es el campo vectorial velocidad , [math]\ p(x,y)=p_1+(p_2-p_1)(x-1)[/math] es el campo de presiones del fluido y [math]μ[/math] es la viscosidad del fluido.

Para la demostración, se trabajará en componentes respecto de la base cartesiana [math]\{\vec{j} {,} \vec{k} \}[/math].

El primer térmido de la ecuación será:

Susutituyendo, nuestro campo y operando tenemos:

Cálculamos ahora, el gradiente del campo de presiones, [math]\nabla p[/math], siendo [math]\ p(x,y)=p_1+(p_2-p_1)(x-1)[/math]

Por último, calculamos el Laplaciano,[math]\Delta\vec{u} [/math] , del campo de velocidades, [math]\vec{u}(y,z)=f(z)\vec{j} [/math]. Para calcularlo, utilizaremos la siguiente fórmula:

1. Cálculo de la divergergencia:

Así [math]\nabla(\nabla \cdot \vec{u})=0[/math]

De este resultado, que a priori, se puede interpretar como una casualidad, se deduce que el fluido es incompresible, pues la condición de incompresibilidad es [math]\nabla \cdot \vec{u}=0[/math]. Teniendo en cuenta que la divergencia mide el cambio de volumen del fluido inducido por el campo, ésta al ser nula, implica que el movimiento de las partículas no afecta al volumen, coincidiendo con la definición dada de fluido incompresible.

2. Cálculo del rotacional:

3. Cálculo de [math]\nabla\times(\nabla\times\vec u)[/math]

Sustituimos los resultados anteriores en [math]\Delta\vec{u} = \nabla(\nabla \cdot \vec{u}) - \nabla \times (\nabla \times \vec{u})[/math], obtenemos que:

Sustituimos por último, todos los términos en la expresión de la ecuación de Navier-Stokes estacionaria [math](\vec{u} \cdot \nabla) \vec{u} + \nabla p=μ \Delta\vec{u} [/math]; resultando:

Por lo que:

[math]\vec{u}f''(z)=\frac{p_2-p_1}{-μ}[/math]

[[Categoría:Teoría de Campos]]

[[Categoría:TC22/23]]

(Editar alguno de los artículos de esta categoría para ver cómo hay que ponerlo.)

4 Referencias