MateWiki - Contribuciones del usuario [es]

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:57:34Z

Mariorodriguezrodriguez: /* Líneas de corriente */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

Su representación gráfica en matlab es:

{{matlab|codigo=
x=0:0.25:4;
y=0:0.2:1;
[xx,yy]=meshgrid(x,y); % Mallado en [0,4]x[0,1] con equidistancia 0.1
p=(1/2).*((yy.^2/2)-(yy.^3/3)); % Potencial escalar
contour(xx,yy,p) %Líneas de nivel
axis([0,4,-1,2]) % Ejes en [0,4]x[-1,2]
}}

[[Archivo:potencial1.jpg|400px|sinmarco|centro|Potencial del campo de velocidades]]

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:
Veamos su representación gráfica:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p=((xx-1).^2)-yy.^2;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=((xx-1).^2)-yy.^2
surf(xx,yy,p)
axis([0,4,-1,2])
}}
[[Archivo:campodetemperaturas1.jpg|400px|sinmarco|centro|Campo de Temperaturas]]
[[Archivo:campodetemperaturas2.jpg|400px|sinmarco|centro|Campo de Temperaturas sobre el plano z=0]]

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

[[Archivo:fugure8.jpg|400px|sinmarco|centro|Campo de Temperaturas y lineas equipotenciales]]

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:57:05Z

Mariorodriguezrodriguez: /* Líneas de corriente */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>
Su representación gráfica es:

{{matlab|codigo=
x=0:0.25:4;
y=0:0.2:1;
[xx,yy]=meshgrid(x,y); % Mallado en [0,4]x[0,1] con equidistancia 0.1
p=(1/2).*((yy.^2/2)-(yy.^3/3)); % Potencial escalar
contour(xx,yy,p) %Líneas de nivel
axis([0,4,-1,2]) % Ejes en [0,4]x[-1,2]
}}

[[Archivo:potencial1.jpg|400px|sinmarco|centro|Potencial del campo de velocidades]]

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:
Veamos su representación gráfica:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p=((xx-1).^2)-yy.^2;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=((xx-1).^2)-yy.^2
surf(xx,yy,p)
axis([0,4,-1,2])
}}
[[Archivo:campodetemperaturas1.jpg|400px|sinmarco|centro|Campo de Temperaturas]]
[[Archivo:campodetemperaturas2.jpg|400px|sinmarco|centro|Campo de Temperaturas sobre el plano z=0]]

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

[[Archivo:fugure8.jpg|400px|sinmarco|centro|Campo de Temperaturas y lineas equipotenciales]]

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:55:29Z

Mariorodriguezrodriguez: /* Líneas de corriente */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>
Su representación gráfica es:

x=0:0.25:4;
y=0:0.2:1;
[xx,yy]=meshgrid(x,y); % Mallado en [0,4]x[0,1] con equidistancia 0.1
p=(1/2).*((yy.^2/2)-(yy.^3/3)); % Potencial escalar
contour(xx,yy,p) %Líneas de nivel
axis([0,4,-1,2]) % Ejes en [0,4]x[-1,2]

[[Archivo:potencial1.jpg|400px|sinmarco|centro|Potencial del campo de velocidades]]

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:
Veamos su representación gráfica:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p=((xx-1).^2)-yy.^2;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=((xx-1).^2)-yy.^2
surf(xx,yy,p)
axis([0,4,-1,2])
}}
[[Archivo:campodetemperaturas1.jpg|400px|sinmarco|centro|Campo de Temperaturas]]
[[Archivo:campodetemperaturas2.jpg|400px|sinmarco|centro|Campo de Temperaturas sobre el plano z=0]]

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

[[Archivo:fugure8.jpg|400px|sinmarco|centro|Campo de Temperaturas y lineas equipotenciales]]

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:54:54Z

Mariorodriguezrodriguez: /* Líneas de corriente */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>
Su representación gráfica es:
[[Archivo:potencial1.jpg|400px|sinmarco|centro|Potencial del campo de velocidades]]

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:
Veamos su representación gráfica:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p=((xx-1).^2)-yy.^2;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=((xx-1).^2)-yy.^2
surf(xx,yy,p)
axis([0,4,-1,2])
}}
[[Archivo:campodetemperaturas1.jpg|400px|sinmarco|centro|Campo de Temperaturas]]
[[Archivo:campodetemperaturas2.jpg|400px|sinmarco|centro|Campo de Temperaturas sobre el plano z=0]]

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

[[Archivo:fugure8.jpg|400px|sinmarco|centro|Campo de Temperaturas y lineas equipotenciales]]

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Archivo:Potencial1.jpg

2013-12-09T17:53:01Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:35:14Z

Mariorodriguezrodriguez: /* Campo de temperatura */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:
Veamos su representación gráfica:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p=((xx-1).^2)-yy.^2;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=((xx-1).^2)-yy.^2
surf(xx,yy,p)
axis([0,4,-1,2])
}}
[[Archivo:campodetemperaturas1.jpg|400px|sinmarco|centro|Campo de Temperaturas]]
[[Archivo:campodetemperaturas2.jpg|400px|sinmarco|centro|Campo de Temperaturas sobre el plano z=0]]

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

[[Archivo:fugure8.jpg|400px|sinmarco|centro|Campo de Temperaturas y lineas equipotenciales]]

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Archivo:Fugure8.jpg

2013-12-09T17:32:38Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:31:02Z

Mariorodriguezrodriguez: /* Campo de temperatura */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:
Veamos su representación gráfica:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p=((xx-1).^2)-yy.^2;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=((xx-1).^2)-yy.^2
surf(xx,yy,p)
axis([0,4,-1,2])
}}
[[Archivo:campodetemperaturas1.jpg|400px|sinmarco|centro|Campo de Temperaturas]]
[[Archivo:campodetemperaturas2.jpg|400px|sinmarco|centro|Campo de Temperaturas sobre el plano z=0]]

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

(falta subir gráfica)

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:29:47Z

Mariorodriguezrodriguez: /* Campo de temperatura */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:
Veamos su representación gráfica:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p=((xx-1).^2)-yy.^2;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=((xx-1).^2)-yy.^2
surf(xx,yy,p)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|campodetemperaturas1]]
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|campodetemperaturas2]]

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

(falta subir gráfica)

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Archivo:Campodetemperaturas2.jpg

2013-12-09T17:26:23Z

Mariorodriguezrodriguez:

Archivo:Campodetemperaturas1.jpg

2013-12-09T17:26:02Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:25:27Z

Mariorodriguezrodriguez: /* Campo de temperatura */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:

Veamos su representación gráfica.
(FALTA GRÁFICA)

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

(falta subir gráfica)
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:19:07Z

Mariorodriguezrodriguez: /* Puntos de velocidad nula y velocidad máxima */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y. Esto es lo que queda representado en la segunda gráfica.
La tercera gráfica muestra todo lo anterior pero en una sección x.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:

Veamos su representación gráfica.
(FALTA GRÁFICA)

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

(falta subir gráfica)

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:17:39Z

Mariorodriguezrodriguez: /* Puntos de velocidad nula y velocidad máxima */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:

Veamos su representación gráfica.
(FALTA GRÁFICA)

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

(falta subir gráfica)

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T17:15:58Z

Mariorodriguezrodriguez: /* Puntos de velocidad nula y velocidad máxima */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:velocidadampliada.jpg|400px|sinmarco|izquierda|Sólido antes y después de desplazarse]]
[[Archivo:seccionvelocidad.jpg|400px|sinmarco|derecha|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:

Veamos su representación gráfica.
(FALTA GRÁFICA)

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}

(falta subir gráfica)

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

.

Archivo:Velocidadampliada.jpg

2013-12-09T17:13:29Z

Mariorodriguezrodriguez:

Archivo:Seccionvelocidad.jpg

2013-12-09T17:12:59Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:58:42Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

[[Archivo:Imagencampodevelocidades.jpg|200px|sinmarco|centro|Campo de velocidades]]

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

[[Archivo:Imagenvalorabsolutorotacional.jpg|200px|sinmarco|centro|Valor absoluto del rotacional]]

== Campo de temperatura ==
Supongamos que la temperatura del fluido viene dado por el campo:

Veamos su representación gráfica.
(FALTA GRÁFICA)

Sea '''$\vec{w}$''' el gradiente de la temperatura. Comprobemos gráficamente que son ortogonales.

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}
{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

(falta subir gráfica)

OBSERVACIÓN: Esta propiedad, aplicada al campo de velocidades, es la utilizada para definir las líneas de corriente del mismo.

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:54:04Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

(poner Imagencampodevelocidades)

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=v_1(x,y)$\vec{i}+v_2(x,y)$\vec{j}=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\exists[0,4]x[0,1]: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

(Explicación analítica)

¿es razonable?

Entendamos el rotacional como un vector perpendicular al plano de curvatura del campo que caracteriza el efecto de giro o rotación del fluido. Para entender esto podemos verlo con el ejemplo práctico de "la paleta":

Para ver la rotación en un punto, cogemos una sección a la que pertenezca el punto y dibujamos su campo de velocidades. Como éste es constante en x, es lógico coger una sección cualquiera "x". (ver figura adjunta). Nos imaginamos al lado del campo de velocidades una "paleta" que tiene como eje de giro el punto en cuestión. El campo de velocidades "empujará" la paleta de una forma u otra en función de dónde se encuentre dicho eje. Veamos algunos ejemplos:
(EJEMPLO EN MEDIO, EN LA MITAD SUPERIOR Y EN EXTREMO)
Si el eje está en medio, el campo por arriba y por abajo son iguales, pero empujan en sentido contrario por tanto, el rotacional será nulo.

Si el eje está en el medio de la parte superior, girará en sentido antihorario porque el empuje del campo por abajo es mayor que el de arriba.

Por último, si el eje está en el extremo superior, todo el campo empujará la paleta en sentido antihorario.

En la parte inferior sería totalmente análogo al ser simétrico.

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

(poner Imagenvalorabsolutorotacional)

== Campo de temperatura ==

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}
{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:52:05Z

Mariorodriguezrodriguez: /* Líneas de corriente */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

(poner Imagencampodevelocidades)

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: [0,4]x[0,1]\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=v_1(x,y)*$\vec{i}+v_2(x,y)*$\vec{j}=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Existe
<math>\existsΨ: \mathbb{R^2}\longrightarrow\mathbb{R}</math> tal que
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

(poner Imagenvalorabsolutorotacional)

== Campo de temperatura ==

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}
{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:47:28Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Ecuaciones de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente, salvo ciertos tipos de flujo y situaciones muy concretas, y se tiene que llevar a cabo un análisis numérico para obtener una solución aproximada.

Las ecuaciones de Navier-Stokes para un fluido son la estacionaria y la condición de incompresibilidad, que se comprueban a continuación

<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center> (Estacionaria)
<center><math>\nabla \cdot \vec{u}=0 </math></center> (Incompresibilidad)

Siendo '''$\vec{u}$''' la velocidad del fluido y '''p''' la presión del mismo.

Sea un fluido con velocidad $\vec{u}$(x,y)=u_1(x,y)$\vec{i}$+u_2(x,y)$\vec{j}$=y*(1-y)\frac{p_1-p_2}{2μ}

=== Ecuación estacionaria ===
Veamos el fluido satisface la ecuación de Navier-Stokes estacionaria
<center><math>\vec{u} \cdot \nabla \vec{u} + \nabla p=μ \Delta\vec{u} </math> </center>

Para poder operar analíticamente pasamos la expresión a lenguaje indicial:

OBSERVACIÓN: Al trabajar en coordenadas cartesianas las componentes contravariantes y covariantes coinciden.
Ahora pasamos la expresión a matrices:

=== Condición de incompresibilidad ===
Continuemos la comprobación de las ecuaciones de Navier-Stokes viendo que se trata de un fluido incompresible, es decir, que ni se expande ni se contrae. Para ello debemos comprobar, como ya antes se ha mencionado:

La divergencia de '''$\vec{u}$''' se define como:

== Campo de presiones y campo de velocidades ==
Para el estudio de los campos de presiones y velocidades, consideramos el flujo del fluido a través de un canal con paredes rectas de dimensiones 1mx4m.

Trabajaremos en el plano con coordenadas cartesianas. Empezamos definiendo un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1], simulando el canal.
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|400px|sinmarco|centro|Mallado de puntos]]

=== Campo de presiones ===
Veamos el campo de presiones del fluido:

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|sinmarco|centro|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|sinmarco|centro|Campo de presiones proyectado sobre z=0]]

El campo de presiones viene representado por una ecuación, la cual sólo depende de la variable x, lo que nos conduce a afirmar que para todo valor de y arbitrario, la presión depende de x exclusivamente. Esto explica y da sentido a nuestras gráficas ya que la zona de color rojo nos indica altas presiones al ser los valores de x próximos a cero, y alcanzar la función máximos valores para cualquier valor de y. A medida que aumentamos los valores de x la presión disminuye y alcanza ese tono azul.

Para obtener la gráfica en 2D simplemente proyectamos.

=== Campo de velocidades ===

{{matlab|codigo=
%visualización del campo de velocidades
x=0:0.5:4; %Mallado en x
y=0:0.1:1; %Mallado en y
[xx,yy]=meshgrid (x,y);
ux=1/2.*(yy-yy.^2);
uy=0.*xx
figure(1)
quiver(xx,yy,ux,uy)
}}

(poner Imagencampodevelocidades)

==== Puntos de velocidad nula y velocidad máxima ====

Corroboramos lo visto analíticamente mediante la siguiente gráfica:

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|400px|sinmarco|centro|Sólido antes y después de desplazarse]]

En la primera gráfica representamos el sólido antes de desplazarse en dicho mallado de puntos, y proyectado según los ejes.
Para representar el desplazamiento longitudinal, usamos el campo de velocidades del fluido, el cual depende de la variable y (altura del punto). Sabiendo que s=vt y para tiempo igual a 1 segundo, s=v, o lo que es lo mismo, el valor del desplazamiento coincide con el de la velocidad. Por lo tanto, sumamos al punto su velocidad y vemos el desplazamiento que experimenta. De esta forma, en las paredes del canal, no hay desplazamiento por tanto la velocidad es nula, y en el centro es máxima, al tener el máximo valor de y.

==== Líneas de corriente ====

Vamos a dibujar ahora las líneas de corriente del campo de velocidades $\vec{u}$, es decir, las líneas de corriente son tangentes a $\vec{u}$en cada punto.
De forma general, siendo $\vec{u}$ un campo cualquiera que tenga Ψ como potencial escalar, las superficies equipotenciales (superficies con el mismo potencial) se caracterizan por ser ortogonales al campo en cada punto.
En el caso tratado, al trabajar en $\mathbb{R^2}$ , las superficies equipotenciales pasarán a ser líneas equipotenciales y se definirán como la intersección entre la superficie equipotencial y el plano en el que se trabaja (en nuestro caso el plano XOY)
Para definir las líneas de corriente $\vec{u}$ (tangentes al campo $\vec{u}$ en cada punto) vamos a definir un campo $\vec{v}$ ortogonal al mismo en cada punto y definir éstas como las líneas equipotenciales de $\vec{v}$.

<math> \vec{v}=\vec{k}\times \vec{u} = \displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

Considerando

<math>\vec{v}: Ω\subset\mathbb{R^2}\longrightarrow\mathbb{R^2}</math> 
<math>(x,y)\longrightarrow\vec{v}(x,y)=\displaystyle\frac{y(1-y)(p_1-p_2)}{n}\vec{j}</math>

<math>\existsΨ: \mathbb{R^2}\longrightarrow\mathbb{R}</math> 
<math>\nablaΨ\longrightarrow\vec{v}(x,y)</math>

<math>
\left .
\begin{matrix}
\frac{\partial Ψ}{\partial x}=u_1(x,y) \\
\frac{\partial Ψ}{\partial y}=u_2(x,y)
\end{matrix}
\right \}</math>

Del sistema de ecuaciones diferenciales se deduce el potencial Ψ

<math>Ψ=(1-x)(y-y^2)\frac{p_1-p_2}{2μ}+C</math>

==== Rotacional ====

{{matlab|codigo=
%Valor absoluto del rotacional
y=0:0.1:1;
x=abs(1/2*(1-2*y));
plot(x,y)
}}

(poner Imagenvalorabsolutorotacional)

== Campo de temperatura ==

{{matlab|codigo=
f=inline('((x-1).^2)-y.^2','x','y'); %Campo escalar de temperatura
fx=inline('2.*(x-1)','x','y'); %derivadas parciales del campo
fy=inline('2.*y','x','y');
x=0:0.1:4;
y=0:0.1:1;
[X,Y]=meshgrid(x,y); %hacemos el mallado
Z=f(X,Y);
contour(X,Y,Z,20,'k'); %Curvas de nivel
hold on
xx=0:0.25:4;
yy=0:0.1:1;
[XX,YY]=meshgrid(xx,yy) %Representación de las flechas
U=fx(XX,YY);
V=fy(XX,YY);
quiver(XX,YY,U,V)
axis equal
}}
{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Archivo:Imagenvalorabsolutorotacional.jpg

2013-12-09T16:27:20Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:26:07Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:17:30Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:12:22Z

Mariorodriguezrodriguez:

Archivo:Imagencampodevelocidades.jpg

2013-12-09T16:09:37Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:03:57Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-09T16:00:00Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:57:33Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente y se tiene que llevar a cabo un análisis numérico para obtener su solución.

== Campos de velocidades, presiones y temperatura ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|Sólido antes y después de desplazarse]]

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|left|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:54:11Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente y se tiene que llevar a cabo un análisis numérico para obtener su solución.

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|Sólido antes y después de desplazarse]]

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|left|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:52:42Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente y se tiene que llevar a cabo un análisis numérico para obtener su solución.

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%Sólido después de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|left|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:49:26Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente y se tiene que llevar a cabo un análisis numérico para obtener su solución.

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

Tomando <math> p_{1}=2</math>, <math> p_{2}=1</math> y <math> \mu=1</math> los campos de presiones y velocidades que resultan:
{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|left|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:45:52Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de Campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente y se tiene que llevar a cabo un análisis numérico para obtener su solución.

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

Tomando p1=2 y p2=1

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|left|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:38:34Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Las ecuaciones de Navier-Stokes son un conjunto de ecuaciones que describen el movimiento de fluidos. Con ellas se pueden analizar corrientes oceánicas, flujo alrededor de aviones o vehículos o cualquier otro fluido newtoniano presente en la atmósfera de la Tierra.
Estas ecuaciones son en derivadas parciales no lineales, y por ello no se tiene solución para las mismas. La mayoría de las veces no se pueden usar analíticamente y se tiene que llevar a cabo un análisis numérico para obtener su solución.

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|left|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:29:54Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|left|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:28:51Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb||Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb||Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:28:11Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|left|Campo de presiones]]
[[Archivo:C-12-321.jpg|200px|thumb|center|Campo de presiones proyectado sobre z=0]]
[[Archivo:C-12-332.jpg|200px|thumb|right|Campo de velocidades]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:25:41Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%Sólido antes de desplazarse
figure(1);
subplot(1,2,1)
p= 3 - xx;
surf(xx,yy,p)
%axis([0,4,-1,2])
figure(2)
subplot(1,2,1)
p=3-xx;
surf(xx,yy,p)
axis([0,4,-1,2])
figure(3)
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy))
uy=0*xx;
quiver(xx,yy,ux,uy)
axis([0,4,-1,2])
}}

[[Archivo:C-12-311.jpg|200px|thumb|right|.......]]
[[Archivo:C-12-321.jpg|200px|thumb|right|.......]]
[[Archivo:C-12-332.jpg|200px|thumb|right|.......]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Archivo:C-12-332.jpg

2013-12-06T16:24:13Z

Mariorodriguezrodriguez:

Archivo:C-12-331.jpg

2013-12-06T16:21:53Z

Mariorodriguezrodriguez:

Archivo:C-12-321.jpg

2013-12-06T16:21:27Z

Mariorodriguezrodriguez:

Archivo:C-12-311.jpg

2013-12-06T16:21:08Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-06T16:12:09Z

Mariorodriguezrodriguez: /* Introducción. Ecuación de Navier-Stokes */

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{matlab|codigo=
%Dibujo del campo
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y);
%SÛlido antes de desplazarse
figure(1);
subplot(1,2,1)
mesh(xx,yy,0*xx)
axis([0,4,-1,2])
%solido despuÈs de desplazarse
subplot(1,2,2)
ux=(1/2)*(yy.*(1-yy));
uy=0*xx;
mesh(xx+ux,yy+uy,0*xx)
axis([0,4,-1,2])
}}
[[Archivo:C-12-21.jpg|200px|thumb|right|.......]]

{{matlab|codigo=}}

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Archivo:C-12-21.jpg

2013-12-06T16:09:15Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-06T15:59:38Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

[[Archivo:C-12-2.jpg|200px|thumb|right|Mallado de puntos]]
[[Archivo:C-12-31.jpg|200px|thumb|right|Mallado de puntos]]
[[Archivo:C-12-32.jpg|200px|thumb|right|Mallado de puntos]]
== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Archivo:C-12-32.jpg

2013-12-06T15:57:14Z

Mariorodriguezrodriguez:

Archivo:C-12-31.jpg

2013-12-06T15:56:51Z

Mariorodriguezrodriguez:

Archivo:C-12-2.jpg

2013-12-06T15:56:09Z

Mariorodriguezrodriguez:

Visualización de campos en fluidos (Grupo C12)

2013-12-06T15:54:57Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

== Introducción. Ecuación de Navier-Stokes ==

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

== Campos de velocidades, presiones y temperatura ==

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]

Visualización de campos en fluidos (Grupo C12)

2013-12-05T10:11:06Z

Mariorodriguezrodriguez:

{{Trabajo|Visualización de campos en fluidos(Grupo C12)|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

Vamos a considerar el flujo de un fluido incompresible a través de un canal con paredes rectas. Trabajaremos en el plano con coordenadas cartesianas.
Para empezar definimos un mallado de ejes [0,4]x[-1,2] y dentro de ellos el fluido que ocupa el rectángulo [0,4]x[0,1]
El codigo matlab para realizarlo es el siguiente
{{matlab|codigo=
x=0:0.1:4; %Mallado en el intervalo [0,4]x[0,1] con longitud de paso de 0,1
y=0:0.1:1;
[xx,yy]=meshgrid(x,y); %Hacemos el mallado
mesh(xx,yy,xx*0)
axis([0,4,-1,2]) %Definimos los ejes }}
[[Archivo:C-12-1.jpg|200px|thumb|left|Mallado de puntos]]

{{ Beta }}
[[Categoría:Teoría de Campos]]
[[Categoría:Trabajos 2013-14]]