MateWiki - Contribuciones del usuario [es]

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

2015-12-10T22:27:38Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

Las líneas de corriente siguen una dirección ortogonal al campo de velocidades, se resolveran del siguiente modo:
[[Archivo:Hiii.png|marco|centro]]

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

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

[[Archivo:Hi4.png|marco|centro]]
Siendo K una constante. Entonces las líneas de corriente siguen la siguiente función:
[[Archivo:Hi5.png|marco|centro]]

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

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

Para representarlo hemos utilizado matlab con el siguiente código:

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

Cuyo grafico es el siguiente:

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

Ahora vamos a demostrar que las líneas de corriente del campo de velocidades \( \vec u \) es lo mismo que calcular la función potencial del campo

Campo irrotacional :

Como se demostró antes
vamos a ver si esto implica que el campo sea irrotacional
[[Archivo:Hi10.png|marco|centro]]

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

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

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

Queda demostrado que el campo es irrotacional por ser su rotación igual a cero.

==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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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-10T22:24:17Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

Las líneas de corriente siguen una dirección ortogonal al campo de velocidades, se resolveran del siguiente modo:
[[Archivo:Hiii.png|marco|centro]]

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

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

[[Archivo:Hi4.png|marco|centro]]
Siendo K una constante. Entonces las líneas de corriente siguen la siguiente función:
[[Archivo:Hi5.png|marco|centro]]

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

Para representarlo hemos utilizado matlab con el siguiente código:

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

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

Ahora vamos a demostrar que las líneas de corriente del campo de velocidades \( \vec u \) es lo mismo que calcular la función potencial del campo

Campo irrotacional :

Como se demostró antes
[[Archivo:Hi8.png|marco|centro]]

vamos a ver si esto implica que el campo sea irrotacional
[[Archivo:Hi10.png|marco|centro]]

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

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

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

Como queríamos demostrar el campo es irrotacional por ser su rotación igual a cero.

==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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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-10T22:23:06Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

Las líneas de corriente siguen una dirección ortogonal al campo de velocidades, se resolveran del siguiente modo:
[[Archivo:Hiii.png|marco|centro]]

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

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

[[Archivo:Hi4.png|marco|centro]]
Siendo K una constante. Entonces las líneas de corriente siguen la siguiente función:
[[Archivo:Hi5.png|marco|centro]]

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

Para representarlo hemos utilizado matlab con el siguiente código:

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

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

Ahora vamos a demostrar que las líneas de corriente del campo de velocidades \( \vec u \) es lo mismo que calcular la función potencial del campo

Campo irrotacional :

Como se demostró antes
[[Archivo:Hi9.png|marco|centro]]

vamos a ver si esto implica que el campo sea irrotacional
[[Archivo:Hi10.png|marco|centro]]

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

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

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

Como queríamos demostrar el campo es irrotacional por ser su rotación igual a cero.

==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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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-10T22:22:18Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

Las líneas de corriente siguen una dirección ortogonal al campo de velocidades, se resolveran del siguiente modo:
[[Archivo:Hiii.png|marco|centro]]

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

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

[[Archivo:Hi4.png|marco|centro]]
Siendo K una constante. Entonces las líneas de corriente siguen la siguiente función:
[[Archivo:Hi5.png|marco|centro]]

[[Archivo:Hi6.png|marco|centro]]
Para representarlo hemos utilizado matlab con el siguiente código:
[[Archivo:Hi7.png|marco|centro]]

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

Ahora vamos a demostrar que las líneas de corriente del campo de velocidades \( \vec u \) es lo mismo que calcular la función potencial del campo

Campo irrotacional :

Como se demostró antes
[[Archivo:Hi9.png|marco|centro]]

vamos a ver si esto implica que el campo sea irrotacional
[[Archivo:Hi10.png|marco|centro]]

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

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

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

Como queríamos demostrar el campo es irrotacional por ser su rotación igual a cero.

==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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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-10T22:21:02Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

Las líneas de corriente siguen una dirección ortogonal al campo de velocidades, se resolveran del siguiente modo:
[[Archivo:Hiii.png|marco|centro]]

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

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

[[Archivo:Hi4.png|marco|centro]]
Siendo K una constante. Entonces las líneas de corriente siguen la siguiente función:
[[Archivo:Hi5.png|marco|centro]]

[[Archivo:Hi6.png|marco|centro]]
Para representarlo hemos utilizado matlab con el siguiente código:
[[Archivo:Hi8.png|marco|centro]]

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

Ahora vamos a demostrar que las líneas de corriente del campo de velocidades \( \vec u \) es lo mismo que calcular la función potencial del campo

Campo irrotacional :

Como se demostró antes
[[Archivo:Hi9.png|marco|centro]]

vamos a ver si esto implica que el campo sea irrotacional
[[Archivo:Hi10.png|marco|centro]]

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

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

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

Como queríamos demostrar el campo es irrotacional por ser su rotación igual a cero.

==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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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-10T22:19:30Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

Las líneas de corriente siguen una dirección ortogonal al campo de velocidades, se resolveran del siguiente modo:
[[Archivo:Hiii.png|marco|centro]]

Ese campo que hemos hallado por ser ortogonal al anterior eso significa
[[Archivo:Hi2.png|marco|centro]]
Vamos a proceder a calcular esta función:
[[Archivo:Hi3.png|marco|centro]]

[[Archivo:Hi4.png|marco|centro]]
Siendo K una constante. Entonces las líneas de corriente siguen la siguiente función:
[[Archivo:Hi5.png|marco|centro]]

[[Archivo:Hi6.png|marco|centro]]
Para representarlo hemos utilizado matlab con el siguiente código:
[[Archivo:Hi8.png|marco|centro]]

[[Archivo:Hi7.png|marco|centro]]
Ahora vamos a demostrar que las líneas de corriente del campo de velocidades \( \vec u \) es lo mismo que calcular la función potencial del campo

Campo irrotacional :

Como se demostró antes
[[Archivo:Hi9.png|marco|centro]]
vamos a ver si esto implica que el campo sea irrotacional
[[Archivo:Hi10.png|marco|centro]]

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

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

[[Archivo:Hi13.png|marco|centro]]
Como queríamos demostrar el campo \( \vec v \) es irrotacional por ser su rotación igual a cero.

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 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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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-10T22:12:33Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

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

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

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

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

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

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

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

[[Archivo:Hi7.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.

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

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

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

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

[[Archivo:Hi13.png|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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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]]

Archivo:Hi13.png

2015-12-10T22:10:02Z

Sarah:

Archivo:Hi12.png

2015-12-10T22:09:47Z

Sarah:

Archivo:Hi11.png

2015-12-10T22:09:35Z

Sarah:

Archivo:Hi10.png

2015-12-10T22:09:16Z

Sarah:

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

2015-12-10T21:43:51Z

Sarah: /* Lineas de corriente */

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

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

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

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

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

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

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

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

[[Archivo:Hi7.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 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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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]]

Archivo:Hi9.png

2015-12-10T21:40:42Z

Sarah:

Archivo:Hi8.png

2015-12-10T21:40:29Z

Sarah:

Archivo:Hi7.png

2015-12-10T21:40:10Z

Sarah:

Archivo:Hi6.png

2015-12-10T21:39:52Z

Sarah:

Archivo:Hi5.png

2015-12-10T21:39:38Z

Sarah:

Archivo:Hi4.png

2015-12-10T21:39:23Z

Sarah:

Archivo:Hi3.png

2015-12-10T21:39:03Z

Sarah:

Archivo:Hi2.png

2015-12-10T21:38:45Z

Sarah:

Archivo:Hiii.png

2015-12-10T21:38:25Z

Sarah:

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

2015-12-10T20:33:29Z

Sarah:

{{ 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==
Se pide dibujar el mallado que represente los puntos interiores del fluido que circula por el canal rectangular dado. Dicho rectángulo ocupa la región [0,4]x[0,1], y se especifica que los ejes estén en el intervalo [0,4]x[-1,2], respectivamente.
Para ello, se ha empleado un sencillo programa en Matlab que nos muestra la malla.
{{matlab|codigo=

%Dibujamos el mallado de los puntos interiores del fluido en el interior
%del canal
[X,Y]=meshgrid([0:0.1:4],[0:0.1:1])
%Hemos dejado un espaciado mínimo de 0.1 en la malla.
mesh(X,Y,0*X)
%A continuación dibujamos la región pedida y visualizamos
grid on
axis([0,4,-1,2])
view(2)

}}

El resultado obtenido de este programa es la siguiente malla:
[[Archivo:Malladomallado.jpg|marco|centro]]

==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]]

==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==

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 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]]

==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 ==

La temperatura del fluido viene dada por la siguiente expresión:
[[Archivo:Temperatura16c.png|marco|ninguna]]

Como esta función está expresada en coordenadas cilíndricas, para el estudio de la temperatura trabajaremos en estas coordenadas. Para ello se debe establecer el dominio de ro y theta teniendo en cuenta su relación con x e y:
[[Archivo:Cilcar.png|marco|centro]]
[[Archivo:Maxmin16c.png|marco|centro]]

Por lo tanto el dominio con el que trabajaremos en coordenadas cilíndricas ser el siguiente:
[[Archivo:domcil.png|marco|ninguna]]
=== Campo de temperaturas ===

El campo de temperaturas asocia una temperatura a cada punto del fluido. Para su representación gráfica se utilizará el siguiente programa de matlab:
{{matlab|codigo=

%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[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 = [pi/2,pi/2];
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 = [pi/2,pi/2];
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:Superficie16cff.jpg|marco|centro|Representación gráfica del campo de temperaturas]]
[[Archivo:Curvas16cff.jpg||marco|centro|Curvas de nivel del campo de temperaturas]]
[[Archivo:Curvasysuperficie16cff.jpg|marco|centro|Representación del campo de temperaturas y sus curvas de nivel de manera conjunta]]

===Temperatura máxima===

Para el cálculo de la temperatura máxima del fluido, así como la deteterminación del punto donde esta se alcanza se utilizará el sguiente programa de MATLAB:
{{matlab|codigo=
%Temperatura máxima
Tmax=max(max(T))
%Coordenadas para Tmax
[theta,ro]=find(T==Tmax);
rom=(4.1231/41)*ro
thetam=((pi/2)/10)*theta
}}

La temperatura máxima es T=0.9998 y se alcanza en ro=0.6034 y theta=0.1571. Al fijarse en los gráficos anteriores, que representaban el campo de temperaturas, se observa que este resultado es completamente coherente, pues ese punto claramente está en el entorno en el que se observan las temperaturas más altas.

=== 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=

%Mallado coordenadas cilíndricas
%Mallado coordenadas cilíndricas
ro=linspace(0,4.1231,41);
theta=linspace(0,pi/2,10);
[RO,THETA]=meshgrid(ro,theta); % Mallado

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

grid on
axis equal
xlabel('RO')
ylabel('THETA')

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

hold on
x = [0,4.1231];
y = [pi/2,pi/2];
plot(x,y,'-k','linewidth',1.5)
}}

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

[[Archivo:gradienteff16czoom.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-03T23:51:17Z

Sarah:

{{ 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==

==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í. ]]

==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-03T23:48:35Z

Sarah:

{{ 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==

==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]]

==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]]

==Conclusiones de las representaciones==

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.

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í. ]]

==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-03T23:46:48Z

Sarah:

{{ 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==

==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 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.

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í. ]]

==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-03T23:39:41Z

Sarah: /* Lineas de corriente */

{{ 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==

==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 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.

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]]

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

==Conclusiones==

En base a estos cálculos, llegamos a distintas conclusiones que se pueden apreciar en la representación gráfica, donde vemos 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í. ]]

==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:Helloprogr.png

2015-12-03T23:33:31Z

Sarah:

Archivo:Helloo.png

2015-12-03T23:33:12Z

Sarah:

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

2015-12-03T23:19:04Z

Sarah:

{{ 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 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.

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 :

==Conclusiones==

En base a estos cálculos, llegamos a distintas conclusiones que se pueden apreciar en la representación gráfica, donde vemos 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í. ]]

==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-03T23:02:13Z

Sarah:

{{ 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.

==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]]

Calculamos el rotacional:

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

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

[[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 :

==Conclusiones==

En base a estos cálculos, llegamos a distintas conclusiones que se pueden apreciar en la representación gráfica, donde vemos 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í. ]]

==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-03T22:51:55Z

Sarah:

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

2015-12-03T22:50:29Z

Sarah:

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

2015-12-03T22:49:00Z

Sarah:

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

2015-12-03T22:44:11Z

Sarah:

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

2015-12-03T22:21:48Z

Sarah:

Archivo:8s.jpg

2015-12-03T22:18:04Z

Sarah:

Archivo:7s.jpg

2015-12-03T22:17:26Z

Sarah:

Archivo:6s.jpg

2015-12-03T22:17:07Z

Sarah:

Archivo:5s.jpg

2015-12-03T22:16:46Z

Sarah:

Archivo:4s.jpg

2015-12-03T22:16:22Z

Sarah:

Archivo:3s.jpg

2015-12-03T22:15:53Z

Sarah:

Archivo:2s.jpg

2015-12-03T22:15:19Z

Sarah:

Archivo:1s.jpg

2015-12-03T22:14:50Z

Sarah: