MateWiki - Contribuciones del usuario [es]

Archivo:Campo escalar1.png

2015-12-04T11:07:28Z

Isabel Roselló:

Campo escalar de presiones del fluido.

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T10:49:14Z

Isabel Roselló: /* Estudio de la presión media de los puntos del fluido */

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

[[Archivo:Campo_escalar1.png|marco|centro]]

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites x[0,4] e y[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división entre la integral y el área del canal.

[[Archivo:Integral_campo.png|marco|centro]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T10:47:19Z

Isabel Roselló: /* Estudio de la presión media de los puntos del fluido */

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

[[Archivo:Campo_escalar1.png|marco|centro]]

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites x[0,4] e y[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división de la integral y el área del canal.

[[Archivo:Integral_campo.png|marco|centro]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Archivo:Campo escalar1.png

2015-12-04T10:40:29Z

Isabel Roselló:

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T10:31:59Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites x[0,4] e y[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división de la integral y el área del canal.

[[Archivo:Integral_campo.png|marco|centro]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Archivo:Integral campo.png

2015-12-04T10:30:58Z

Isabel Roselló:

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T10:23:52Z

Isabel Roselló: /* Estudio de la presión media de los puntos del fluido */

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites x[0,4] e y[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división de la integral y el área del canal.

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T10:04:03Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites x[0,4] e y[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división de la integral y el área del canal.

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T09:59:04Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

[[Archivo:Campo_escalar.jpg|sinmarco|centro]]

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites x[0,4] e y[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división de la integral y el área del canal.

[[Archivo:Integral_presion_media1.png|sinmarco|centro]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T09:56:56Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

[[Archivo:Campo_escalar.jpg|sinmarco|centro]]

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites \[x\in[0,4]\] e y\in[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división de la integral y el área del canal.

[[Archivo:Integral_presion_media1.png|sinmarco|centro]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Archivo:Integral presion media1.png

2015-12-04T09:55:35Z

Isabel Roselló:

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T09:54:10Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

[[Archivo:Campo_escalar.jpg|sinmarco|centro]]

Una vez tenemos la expresión del campo escalar de las presiones calculamos la integral definida en los limites x\in[0,4] e y\in[0,1] y el área de nuestro canal. La presión media de los puntos del fluido será la división de la integral y el área del canal.

[[Archivo:

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T09:45:28Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

[[Archivo:Campo_escalar.jpg|200px|sinmarco|centro]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Archivo:Integral presion media.jpg

2015-12-04T09:44:07Z

Isabel Roselló:

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-04T09:42:43Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Mallado ocupado por el fluido==

==Comprobación ecuación Navier-Stokes==
Partimos de las ecuaciones de la velocidad y la presión de las partículas del fluido, respectivamente:
[[Archivo:Camp1.jpg|marco|centro]]
Con ello se pide verificar la ecuación de Navier-Stokes:
[[Archivo:Camp3.jpg|marco|centro]]
Para ello vamos a ir parte por parte. En primer lugar, tenemos:
[[Archivo:Camp4.jpg|marco|centro]]
Desarrollamos este primer término:
[[Archivo:Camp5.jpg|marco|centro]] donde todos los términos son nulos, excepto la derivada parcial de u1 con respecto a y, y esto se debe a que:
-Se ha de cumplir la expresión que nos dan como resultado del primer sumando de la ecuación, lo que implica que la primera y la tercera derivada parcial del desarrollo anterior sean nulos.

-u2,u3 son 0, ya que la ecuación velocidad que se proporciona está representada en i.
De esta forma nos queda:
[[Archivo:Camp6.jpg|marco|centro]]
[[Archivo:Camp10.jpg|marco|centro]]
Por otro lado, vamos con otro término de la ecuación:
[[Archivo:Camp7.jpg|marco|centro]]
Desarrollando, se tiene:
[[Archivo:Camp8.jpg|marco|centro]]
Por último, vamos con el segundo sumando de la ecuación a demostrar:
[[Archivo:Camp9.jpg|marco|centro]]
Finalmente, basta con sustituir cada uno de estos términos en la ecuación de Navier-Stokes:
[[Archivo:Camp11.jpg|marco|centro]]
Por lo tanto, queda así demostrado.

Por otro lado, vamos a comprobar la condición de incompresibilidad de un fluido, para lo cual partimos de la expresión facilitada y desarrollamos:
[[Archivo:Camp12.jpg|marco|centro]]
Para terminar, comprobamos que la velocidad del fluido en las paredes y=0, y=1 es nula, para lo cual simplemente sustituimos en la expresión de la velocidad dada:
[[Archivo:Camp0.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de presiones==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

En el campo de presiones podemos apreciar que la presion es más alta en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja. Esto significa que a medida que el fluido va avanzando en el canal la presión va bajando.

==Representacion del campo de velocidades==

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

En el campo de velocidades podemos observar que todos los vectores apuntan hacia la misma dirección, en sentido paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Lineas de corriente==

Ahora vamos a dibujar las líneas de corriente des campo u, es decir, las líneas que son tangentes a este en cada punto. Para esto calculamos un campo v, que es ortogonal a u.

[[Archivo:1s.jpg|marco|centro]]

Vamos a calcular el potencial escalar ψ, que se conoce como función de corriente de u :

[[Archivo:2s.jpg|marco|centro]]

[[Archivo:3s.jpg|marco|centro]]

[[Archivo:4s.jpg|marco|centro]]

[[Archivo:5s.jpg|marco|centro]]

[[Archivo:6s.jpg|marco|centro]]

Calculamos el rotacional:

[[Archivo:7s.jpg|marco|centro]]

Se comprueba que es irrotacional (por ser u de divergencia nula)

[[Archivo:8s.jpg|marco|centro]]

Se mostrara a continuacion la representacion de las lineas de corriente de u (ψ=cte) que el enunciado pide.

Usando el siguiente programa de matlab :

[[Archivo:Helloprogr.png|marco|centro]]

El grafico resultante es el siguiente:

[[Archivo:Helloo.png|marco|centro]]

Las dos conclusiones fundamentales que se pueden apreciar en la representación gráfica son que las líneas de corriente son efectivamente tangentes al campo de velocidades y, además que el campo \( \vec v \) es irrotacional debido a que su divergencia es nula.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Estudio de la presión media de los puntos del fluido==

La presión de los puntos del fluido viene dada por el campo escalar. Particularizando en los valores que nos dan para p1 y p2:

[[Archivo:Campo_escalar.jpg|sinmarco|centro]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Archivo:Campo escalar.jpg

2015-12-04T09:38:53Z

Isabel Roselló:

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-03T18:01:38Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de velocidades==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Repres3D.png|marco|centro]]

==Conclusiones de las representaciones==
En el campo de presiones podemos ver que en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), la presión es más alta y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja, esto significa que la presión baja a medida que el fluido va avanzando en el canal.

En el campo de velocidades, apreciamos que todos los vectores apuntan a la misma dirección, paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Archivo:Repres3D.png

2015-12-03T18:00:41Z

Isabel Roselló:

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-03T17:53:31Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de velocidades==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:
[[Archivo:Porgramapresion.png|marco|centro]]

La representacion en 2D:
[[Archivo:Reprpresion.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3D.png|marco|centro]]

El codigo para el campo de velocidades es el siguiente:
[[Archivo:Programavelocidad.png|marco|centro]]

La representacion en 2D:
[[Archivo:Representacion2Dvel.png|marco|centro]]

La representacion en 3D:
[[Archivo:Representacion3Dvel|marco|centro]]

==Conclusiones de las representaciones==
En el campo de presiones podemos ver que en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), la presión es más alta y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja, esto significa que la presión baja a medida que el fluido va avanzando en el canal.

En el campo de velocidades, apreciamos que todos los vectores apuntan a la misma dirección, paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Archivo:Representacion3Dvel.png

2015-12-03T17:30:12Z

Isabel Roselló:

Archivo:Reprpresion.png

2015-12-03T17:27:51Z

Isabel Roselló:

Archivo:Representacion3D.png

2015-12-03T17:27:39Z

Isabel Roselló:

Archivo:Representacion2Dvel.png

2015-12-03T17:27:28Z

Isabel Roselló:

Archivo:Programavelocidad.png

2015-12-03T17:27:18Z

Isabel Roselló:

Archivo:Porgramapresion.png

2015-12-03T17:26:59Z

Isabel Roselló:

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-03T17:17:38Z

Isabel Roselló:

{{ TrabajoED | Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC15/16|2015-16]] | Sarah Boufounas, Laura García Pedraz, Isabel Roselló Colom, Ignacio Riesco Fernández , Larisa Elena Suta, Álvaro Valbuena Pámpanas }}

El trabajo que nos ocupa es el del estudio de la visualización de campos escalares y vectoriales en fluidos. La consideración a tomar como fluido ha sido la de un fluido incompresible a traves de un canal con pareces rectas.
La variación del volumen del fluido va a ser por tanto nula.
El canal está compuesto por dos paneles rectos.
Las coordenadas en la que refiramos los apartados serán cartesianas.

==Breve introducción a la incompresibilidad de un fluido==
Un fluido incompresible es cualquier fluido cuya densidad siempre permanece constante con el tiempo, y tiene la capacidad de oponerse a la compresión del mismo bajo cualquier condición. Esto quiere decir que ni la masa ni el volumen del fluido puede cambiar. El agua es un fluido casi incompresible, es decir, la cantidad de volumen y la cantidad de masa permanecerán prácticamente iguales, aún bajo presión. De hecho, todos los fluidos son compresibles, algunos más que otros. La compresión de un fluido mide el cambio en el volumen de una cierta cantidad de líquido cuando se somete a una presión exterior.

Por esta razón, para simplificar las ecuaciones de la mecánica de fluidos, se considera que los líquidos son incompresibles. En términos matemáticos, esto significa que la densidad de tal fluido se supone constante.
Se debe prestar atención a todas las propiedades del fluido (aire, agua) para definir las condiciones de flujo. Esto se debe a que todas las propiedades están conectadas entre sí. Si la presión o la temperatura de un fluido cambia, su densidad generalmente también cambia (a menos que se trate de un fluido incompresible, como en el ejercicio que ahora estudiamos).
[[Archivo:Awa.jpg.jpg|marco|centro]]

==Estudio de la velocidad del fluido==
===Método de estudio===
El campo de velocidades realizado en el apartado anterior ha dejado de forma gráfica patente, dónde se desarrollan las velocidades máximas de nuestro fluido; sin embargo y a nivel matemático, para desarrollarlo de forma analítica bastaría con realizar la derivada del campo de velocidades con respecto a Y llegando a la conclusión de que las mayores velocidades se dan en la céntrica línea de Y=1/2
En la siguiente imagen podemos observar las operaciones realizadas para llegar a nuestra conclusión.
[[Archivo:Roti.png|marco|centro]]
===¿Esta solución tiene un sentido físico? ¿Es razonable?===
Sí, en efecto podríamos hablar de un sentido físico dado que en la circulación de fluidos por conductos o tuberías la velocidad máxima siempre tiende a localizarse en el centro. Los extremos siempre conllevan perdidas de carga continuas y fuerzas de rozamiento. Es por el centro donde la circulación es más libre y por tanto la velocidad se incrementa.
[[Archivo:Centri.jpg|marco|centro]]

==Representacion del campo de velocidades==
Para representar el campo de velocidades y presiones vamos a suponer que p1=2 y p2=1 , con mu=1.
Para ello usaremos el programa de Matlab.

El codigo para el campo de presiones es el siguiente:

El codigo para el campo de velocidades es el siguiente:

==Conclusiones de las representaciones==
En el campo de presiones podemos ver que en las zonas más cálidas (colores como rojo, naranja, amarillo y verde), la presión es más alta y en las zonas más frías (colores como verde,celeste y azul), la presión es más baja, esto significa que la presión baja a medida que el fluido va avanzando en el canal.

En el campo de velocidades, apreciamos que todos los vectores apuntan a la misma dirección, paralelo al eje x. También que la velocidad es mayor en el centro del canal que en las paredes porque los vectores tienen mayor módulo ahí, esto ocurre porque no hay barreras en el centro que pueda provocar una disminución en la velocidad.

==Estudio del rotacional de campo==
Usamos el concepto de rotacional de campo en el cálculo vectorial, el rotacional o también denominado a veces rotor es un operador vectorial para mostrar la tendencia de un campo vectorial a inducir rotación alrededor de un punto concreto.
Matemáticamente, esta idea se expresa como el límite de la circulación del campo vectorial, cuando la curva sobre la que se integra se reduce a un punto, es decir:
[[Archivo:Plapo.png|marco|centro]]
Sin embargo se puede expresar de forma más concisa con ayuda del operador nabla como un producto vectorial, calculable mediante un determinante:
[[Archivo:Deteret.png|marco|centro]]
y relacionando el operador nabla con nuestro campo,
[[Archivo:Rotixx.png|marco|centro]]

{{matlab|codigo=

x=0:0.05:4;
y=0:0.05:1;
[X,Y]=meshgrid(x,y); % Mallado

rotabs=abs(Y-(1/2)); %Rotacional

surf(X,Y,rotabs)

axis([0,4,-1,2])
}}
[[Archivo:16C-rotabs.png|marco|centro|Representación del rotacional en valor absoluto]]
El dibujo del mallado en Matlab deja patente que el rotacional crece en valor absoluto a medida que nos acercamos a las paredes de la sección. Eso quiere decir que a medida que el rotacional crece, decrecen la velocidad de las partículas.
El rotacional es mayor en los puntos cercanos a las paredes porque estos son más propensos a generar rotación alrededor de puntos en concreto al tener una velocidad más pequeña y una trayectoria menos libre.

== Estudio de la temperatura del fluido ==
=== Campo de temperaturas ===
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,0.2450,3);
[RO,THETA]=meshgrid(ro,theta); % Mallado

T= cos(THETA).*exp(-(RO-(1/2)).^2); %Temperatura
%Superficie
figure(1)
surface(RO,THETA,T)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel
figure(2)
contour(RO,THETA,T,30)
axis equal
hold on
x = [0,4.1231];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4.1231];
y = [0.245,0.245];
plot(x,y,'-k','linewidth',1.5)

%Curvas de nivel y superficie 3D
figure (3)
surf(RO,THETA,T)
hold on
contour(RO,THETA,T,30)
grid on
axis equal
xlabel('RO')
ylabel('THETA')
zlabel('Temperatura')
colorbar

}}

[[Archivo:Superficie16cf.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cf.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cf.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[ro,theta]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=(0.245/3)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.0817

=== Gradiente de temperatura ===
La temperatura es una magnitud intensiva que adopta un valor distinto para cada punto del espacio, por lo que puede ser descrito como un campo escalar.
En cálculo vectorial, el gradiente de un campo escalar '''f''' es un campo vectorial. El vector gradiente de '''f''' evaluado en un punto genérico '''x''' del dominio de '''f''' , '''f(x)''', indica la dirección en la cual el campo '''f''' varía más rápidamente y su módulo representa el ritmo de variación de '''f''' en la dirección de dicho vector gradiente.

Si suponemos que '''c : R → R''' n es una curva parametrizada con la propiedad de que '''c(t)''' pertenece al k conjunto de nivel de '''f''' para cualquier '''t ∈ R'''. Esto significa simplemente que '''f ◦ c''' es constantemente igual a k. Por lo tanto: '''d(f ◦ c) dt = 0''', porque es la derivada de una constante. Pero como, por otra parte, tenemos: '''d(f ◦ c) dt = grad f(p) · c’(t)''' , concluimos que el gradiente es perpendicular a las curvas de nivel.

Por tanto, el gradiente térmico o gradiente de temperatura es la variación de temperatura por unidad de distancia.
Siendo esto así, para cada punto del fluido, el gradiente en ese punto nos dará la dirección en la cual la temperatura aumenta más rápido. La magnitud del gradiente nos dirá cuán rápido aumenta la temperatura en esa dirección.

Mediante Matlab obtenemos la representación gráfica de las curvas de nivel de la temperatura y el gradiente térmico, y podemos comprobar gráficamente que el gradiente es perpendicular a las curvas de nivel.

{{matlab|codigo=

x=linspace(0,4,40);
y=linspace(0,1,10);
[X,Y]=meshgrid(x,y); % Mallado

T= cos(Y).*exp(-(X-(1/2)).^2); %Temperatura
Tro=(-2*X+1).*cos(Y).*exp(-(X-(1/2)).^2); %Componente ρ
Ttheta=-sin(Y).*exp(-(X-(1/2)).^2); %Componente θ
hold on
quiver(X,Y,Tro,Ttheta); %Representación del gradiente
contour(X,Y,T,30); %Representación de la temperatura

grid on
axis([0,4,-0.5,1.5])
axis equal
xlabel('X')
ylabel('Y')

hold on
x = [0,4];
y = [0,0];
plot(x,y,'-k','linewidth',1.5)

hold on
x = [0,4];
y = [1,1];
plot(x,y,'-k','linewidth',1.5)
}}

[[Archivo:gradientef16c.jpg|marco|centro|Representación gráfica de las curvas de nivel y el gradiente de temperatura.]]

[[Archivo:gradientef16czoom.jpg|marco|miniaturadeimagen|centro|Ampliación del gradiente y las curvas de nivel, donde se puede observar con más detalle que son perpendiculares entre sí. ]]

==Caudal==

El caudal es la cantidad flujo volumétrico o volumen que pasa por un área dada en la unidad de tiempo.
Para calcular el caudal, lo hacemos mediante una integral de línea:

[[Archivo:16C-10-1.png|sinmarco|centro]]
En nuestro caso,

[[Archivo:16C-10-2.png|marco|centro]]

donde p1=2, p2=1 y μ=1.
Por tanto

[[Archivo:16C-10-3.png|marco|centro]]

El ancho del canal es 1 y por tanto los límites de nuestra integral serán 0 y 1.

[[Archivo:16C-10-4.png|marco|centro]]

El caudal que lleva el canal es de 0.083 m3/s
=== Porcentaje del caudal que pasa por la mitad central del canal===

Para calcular el caudal que pasa por la mitad central del canal lo hacemos con la misma integral empleada para el caudal total, pero en el intervalo 0.25 y 0.75, teniendo el valor medio de y, en el centro de dicho intervalo.

[[Archivo:16C-10.2-1.png|marco|centro]]

Y el porcentaje del caudal que pasa por la mitad central del canal es:

[[Archivo:16C-10.2-2.png|marco|centro]]

[[Categoría:Teoría de Campos]]
[[Categoría:TC15/16]]

Visualización de campos escalares y vectoriales en fluidos. Grupo 16-C

2015-12-03T17:13:55Z