MateWiki - Contribuciones del usuario [es]

Visualización de campos escalares y vectoriales en fluidos. (Grupo 23-C)

2013-12-10T02:44:47Z

Vicente:

{{beta}}
{{Trabajo|Visualización de campos escalares y vectoriales en fluidos. Grupo 23-C|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

==Introducción==
Vamos a estudiar el flujo de un fluido incompresible alrededor de un obstáculo con forma circular de radio igual a 1, ósea un círculo unidad. Para la realización de este trabajo usaremos tanto coordendas cartesianas como con coordenadas Cilíndricas (polares).
== Mallado ==
La región ocupada por el fluido será el exterior del círculo unidad tal y como indica a continuación el mallado que hemos dibujado.
Hemos tomado un mallado del anillo comprendido entre los radio 1 y 5 y centro el origen, con unos ejes dentro del intervalo [-4,4]x[-4,4].

El código matlab con el que lo hemos realizado es el siguiente:
{{matlab|codigo=
u=linspace(1,5,100); % Creamos el vector u en el intervalo[1,5]
v=linspace(0,2*pi,100); % Creamos el vector v en el intervalo [0,2*pi]
[uu,vv]=meshgrid(u,v); % matrices de las coordenadas u y v
xx=uu.*cos(vv); % Parametrizamos en coordendas cartesianas y así definimos el anillo
yy=uu.*sin(vv);
mesh(xx,yy,0*xx) % Dibujamos el mallado
hold on
x=cos(v);
y=sin(v);
plot(x,y,'k','Linewidth',1) % Dibujamos el obstáculo: círculo unidad
axis([-4,4,-4,4]) % Definimos los ejes en [-4,4]x[-4,4]
view(2) % Visualizamos lo dibujado en el plano x e y
}}
[[Archivo:Malladoc.png|marco|centro|Mallado de la región]]
== Velocidad de las partículas del fluido ==
La velocidad de las partículas del fluido viene dada por un campo vectorial <math>\vec u </math> cuya función potencial es φ=(ρ+1/p)*cos⁡θ. Luego <math>\vec u</math> es el gradiente de la función potencial: : <math>\vec u=∇φ</math>
Realizaremos los cálculos tanto en cartesianas como en cilíndricas: : <math>φ(ρ,θ)=(ρ+1/p)*cos⁡θ = φ(x,y)=x+\frac{x}{x^2+y^2}</math>
En el caso de Cartesianas:: <math>φ(x,y)=x+\frac{x}{x^2+y^2}</math>

: <math>\vec u=∇φ(x,y)=\frac{∂φ}{∂x} \vec i +\frac{∂φ}{∂y} \vec j =\frac{1+(y^2-x^2)}{(x^2+y^2 )^2 } \vec i -\frac{(2xy}{(x^2+y^2 )^2 )} \vec j</math>

En el caso de Cilíndricas calculamos tanto las coordenadas contravariantes como las covariantes: : <math>φ(ρ,θ)=(ρ+1/p)*cos⁡θ</math>
: <math>\vec u=∇φ(ρ,θ)=\frac{∂φ}{∂ρ} \vec g^ρ +\frac{∂φ}{∂θ} \vec g^θ =cos⁡θ (1-\frac{1}{ρ^2 }) \vec g^ρ -sin⁡θ \frac{ρ+1}{ρ} \vec g^θ =cos⁡θ(1-\frac {1}{ρ^2} ) \vec g_ρ -\frac{sin⁡θ}{ρ^2} \frac{ρ+1}{ρ} \vec g_θ</math>

Las representaciones de estos dos campos, el campo vectorial <math>\vec u </math> y el campo escalar φ, función potencial de <math>\vec u </math>, son las siguientes:

{{matlab|codigo=
%Representación del Campo escalar φ
u=linspace(1,5,100);
v=linspace(0,2*pi,100);
[uu,vv]=meshgrid(u,v); %mallado sobre el que representar
xx=uu.*cos(vv);
yy=uu.*sin(vv);
fi=inline('x+x./(x.^2+y.^2)','x','y'); %Definimos la función potencial fi
f=fi(xx,yy);
contour(xx,yy,f,30) %Dibujamos sus líneas de nivel
axis([-4,4,-4,4])
hold on
x=cos(v);
y=sin(v);
plot(x,y,'k','Linewidth',1) % Dibujamos el campo escalar fi
axis([-4,4,-4,4]) % Seleccionamos la región donde lo dibuja
axis equal

%Representación del Campo vectorial u
u=linspace(1,5,15);
v=linspace(0,2*pi,15);
[uu,vv]=meshgrid(u,v); %mallado sobre el que representar
xx=uu.*cos(vv);
yy=uu.*sin(vv);
fx=inline('1+(y.^2-x.^2)./((y.^2+x.^2).^2)','x','y'); %Definimos la función del campo u de velocidades en x
fy=inline('-2*(x.*y)./((y.^2+x.^2).^2)','x','y'); %Definimos la función del campo u de velocidades en y
u1=fx(xx,yy);
u2=fy(xx,yy);
quiver(xx,yy,u1,u2) %dibujamos el campo vectorial u
axis([-4,4,-4,4])
view(2) %Para visualizar lo dibujado

hold off
}}
[[Archivo:Campo ortogonal a.png|marco|centro|Campo vectorial <math>\vec u </math> ortogonal a las curvas de nivel de φ ]]

Como podemos observar, el campo vectorial <math>\vec u </math> de las velocidades de las partículas es ortogonal a las curvas de nivel de φ .
[[Archivo:Campo ortogonal a más de cerca.png|marco|centro|Campo <math>\vec u </math> ortogonal a las líneas de nivel de φ zoom]]
Esto se debe a que
== Rotacional de <math>\vec u</math> nulo==

El rotacional de un campo nos indica la tendencia de este a inducir rotación alrededor de un punto.

El rotacional de nuestro campo <math>\vec u</math> es nulo ya que todo campo potencial, es decir, expresable como el gradiente de un potencial escalar es irrotacional.

Vamos a comprobarlo entonces analíticamente realizando los correspondientes cálculo en coordenadas cilíndricas:

<math>\nabla\times\vec{u}= \frac {1}{g^\frac{1}{2}}\begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\ \frac{ \partial}{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z}\\ u_1 & u_2 & u_3 \end{pmatrix}= \frac {1}{ρ} \begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\ \frac{ \partial}{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z}\\ cos(θ)(1 - \frac {1}{ρ^2})& - sen (θ) (ρ + \frac {1}{ρ}) & 0 \end{pmatrix}=</math>
<math>=\frac {1}{ρ}(\frac{\partial}{\partial z}((cos(θ)(1 - \frac {1}{ρ^2}))\vec{g_θ} ) - (\frac{\partial}{\partial ρ}(sen (θ) (ρ + \frac {1}{ρ}))\vec{g_z})-(\frac{\partial}{\partial θ}((cos(θ)(1 - \frac {1}{ρ^2}))\vec{g_z} + (\frac{\partial}{\partial z}(sen (θ) (ρ + \frac {1}{ρ})\vec{g_ρ}) )=</math>
<math>\frac {1}{ρ}(-sen(θ) (1- \frac{1}{ρ^2}) \vec{g_z} + sen(θ) (1- \frac{1}{ρ^2}) \vec{g_z})= 0 </math>

== Divergencia de <math>\vec u</math> nula==

Que la divergencia de <math>\vec u</math> sea nula se debe a que se trata de un fluido incompresible que implica que este tiene densidad constante en todo el conjunto y a lo largo del tiempo.

La divergencia mide la diferencia de flujo saliente y flujo entrante de un campo vectorial sobre la superficie que rodea un volumen control. Si la divergencia en un punto es positiva, se dice que el campo posee fuentes. Si la divergencia es negativa, se dice que tiene sumideros.

Vamos a comprobarlo analíticamente realizando los orrespondientes cálculo en coordenadas cilíndricas::

<math>\nabla·\vec{u}= \frac {1}{g^\frac{1}{2}}\frac{\partial }{\partial x^i}(g^\frac{1}{2} u^i)= </math>
<math>= \frac {1}{ρ}(\frac{\partial }{\partial ρ}(ρ u^ρ) + \frac{\partial }{\partial θ}(ρ u^θ) + \frac{\partial }{\partial z}(ρ u^z))=</math>
<math>\frac {1}{ρ}(\frac{\partial }{\partial ρ}((ρ - \frac{1}{ρ})cos(θ)) - \frac{\partial }{\partial θ}(1 + \frac{1}{ρ^2})sen(θ) + 0 )=</math>
<math>((1 + \frac{1}{ρ^2})cos(θ) - (1 + \frac{1}{ρ^2})cos(θ))= 0</math>

==Líneas de corriente==
Las líneas de corriente del campo <math>\vec u</math> son las líneas que son tangentes a <math>\vec u</math> en cada punto.
Para calcularlas primeramente hallamos el campo <math>\vec u</math> que en cada punto es ortogonal a <math>\vec u</math> , luego hacemos:
<math>\vec{v}=\vec{k}\times\vec{u}</math>
: <math>\vec{v}=\vec{k}\times\vec{u}= {g^\frac{1}{2}}\begin{pmatrix} \vec{i} & \vec{j} & \vec{k} \\ 0 & 0 & 1 \\ u^1 & u^2 & u^3 \end{pmatrix}= 1 (\vec{u^1} \vec{j} - \vec{u^2} \vec{i}) = (\frac{-2xy}{(x^2+y^2 )^2})\vec{i} + (1 + \frac{y^2-x^2}{(x^2+y^2 )^2})\vec{j}</math>

Este nuevo campo vectorial <math>\vec v</math> es irrotacional, es decir que su rotaciones es 0; como demostramos a continuación ya que <math>\vec u</math> es de divergencia nula.

<math>\nabla\times\vec{v}= \frac {1}{g^\frac{1}{2}}\begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\\frac{\partial }{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z} \\ v_1 & v_2 & v_3 \end{pmatrix}=</math>
<math>=\frac{1}{ρ}\begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\ \frac{\partial }{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z} \\ (ρ+\frac{1}{ρ}) \frac{sen(θ)}{ρ} & (ρ-\frac{1}{ρ}) cos(θ) & 0 \end{pmatrix}=</math>
<math>=\frac{1}{ρ}((1+\frac{1}{ρ^2})cos(θ)- (ρ+\frac{1}{ρ})\frac{cos(θ)}{ρ})\vec{g_z}= \frac{1}{ρ}((1+\frac{1}{ρ^2}) - (1+\frac{1}{ρ^2}))cos(θ)\vec{g_z} = 0</math>

Y tiene un potencial escalar <math>ψ</math> que se conoce como función de corriente de <math>\vec u</math> y que vamos a calcular:

: <math>\vec{v}=\vec{grad}ψ= \frac{∂ψ}{∂x}\vec{i} + \frac{∂ψ}{∂y}\vec{j} + \frac{∂ψ}{∂z}\vec{k} = (\frac{-2xy}{(x^2+y^2 )^2})\vec{i} + (1 + \frac{y^2-x^2}{(x^2+y^2 )^2})\vec{j}</math>

De esta fórmula anterior podemos sacar que:
: <math>\frac{∂ψ}{∂x}= \frac{-2xy}{(x^2+y^2 )^2}</math>
: <math>\frac{∂ψ}{∂y}= -1- \frac{y^2-x^2}{(x^2+y^2 )^2}</math>
Despejando <math>ψ</math> de <math>\frac{∂ψ}{∂x}</math> e Integrando obtenemos la función de <math>ψ(x,y)</math>:
<math>ψ=\int \frac{-2xy}{(x^2+y^2 )^2} dx= -\frac{y}{x^2+y^2} + h(y)</math>

A continuación derivamos la <math>ψ</math> obtenida respecto "y" y la igualamos a la expresión que teníamos anterior de <math>\frac{∂ψ}{∂y}</math> de manera que así podemos sacar <math>h(y)</math> y completar la función <math>ψ</math>
: <math>\frac{∂ψ}{∂y}=-\frac{y^2-x^2}{(x^2+y^2 )^2}+h'(y)= -1- \frac{y^2-x^2}{(x^2+y^2 )^2}</math>
Despejamos <math>h'(y)</math> y sacamos integrándola <math>h(y)</math>: <math>h'(y)=-1</math>
: <math>h(y)=-y</math>
Por lo tanto el potencial escalar <math>ψ</math> es igual a: <math>ψ=-\frac{y}{x^2+y^2}-y</math>

== Estudio de la velocidad del fluido ==

A continuación se realiza el estudio de la velocidad del fluido en los puntos de la frontera S, que son aquellos en los que se particulariza la función ϒ para ρ=1, definida por la función determinada por el gradiente de la función ϒ.

Para determinar los valores máximo y mínimo de la velocidad consideramos su módulo:

<math>\vec{v}= cos(θ)(1-\frac{1}{ρ^2})\vec{g_ρ} - sen (θ) (\frac{1}{ρ} + \frac{1}{ρ^3})\vec{g_θ}</math>

: <math>|\vec{v}|=\sqrt(((1- \frac{1}{ρ^2})^2 (cos(θ))^2 + (\frac{1}{ρ}+\frac{1}{ρ^3})^2 (sen(θ)^2)))=f(ρ,θ)</math>

Los valores máximo y mínimos se obtienen sustituyendo en la f(ρ,θ) para ρ=1:

<math>f(1,θ)=\sqrt((4(sen(θ))^2))=2|sen(θ)|</math>

Por tanto, la velocidad será máxima cuando la función seno tome su máximo, es decir, cuando toma el valor 1:
<math>|\vec u|(\frac{ \pi}{ 2})= 2</math>
: <math>|\vec u| (\frac{ 3\pi}{ 2})= 2</math>

En los puntos en coordenadas polares (1,π⁄2) y (1,3π⁄2), la velocidad máxima será igual a 2.
Por su parte, la velocidad será mínima cuando la función seno tome su valor mínimo, es decir, cuando toma el valor 0:
<math>|\vec u| (0)= 0 </math>
: <math>|\vec u| (\pi)= 0 </math>

En los puntos en coordenadas polares (1,0) y (1,π), la velocidad mínima será igual a 0. Dichos puntos serán puntos de remanso, pues en ellos su velocidad es nula.

== Estudio de la presión del fluído ==

Basándonos en la ecuación de Bernouilli:

<math> \frac{1}{2}ρ (|\vec{u}|)^2 + p = cte </math>

con densidad <math>ρ = 2</math>.
La mayor presión se alcanzará al tener la menor velocidad y la menor presión se alcanzará al tener la mayor velocidad,y en el gráfico de la presión, se pueden distinguir el valor de la presión en cualquier punto del fluido.
[[Archivo:Grafica Bernouilli.png|marco|centro|Gráfica de la ecuación de Bernouilli]]
Si se consideran los puntos de mínima velocidad (v = 0), el primer término se anula y entonces la presión alcanza su valor máximo, que se corresponde con el valor de la constante. Por el contrario, si se toman los puntos de velocidad máxima se tendrá el menor valor de presión. Sabiendo que la velocidad máxima es 2 y la densidad del fluido es 2, el valor del primer término es 4 y se puede determinar que el valor de la constante tiene que ser mayor que 4 para que se obtenga un valor de presión positivo, en este caso usamos un valor de 50 para poder observar bien la gráfica.

== Trayectoria de una partícula del fluído ==

Supongamos que nuestra partícula va por una trayectoria horizontal paralela y próxima al eje del diámetro horizontal del objeto circular, por el hecho de estar en el fluido llevara una determinada velocidad y tendrá una determinada presión, al aproximarse al objeto empezara a curvar su trayectoria con el objetivo de rodear el objeto, como hemos supuesto que lleva una trayectoria próxima al diámetro horizontal, mas o menos recorrerá el objeto por su superficie <math> ρ=1 </math>, a medida que empiece a rodear el objeto empezara a disminuir su presión y a aumentar su velocidad, hasta el punto de llegar a ser tangente a la circunferencia (corte con el diámetro vertical) dónde la velocidad será máxima, una vez alcanzado este punto empezará a disminuir, además por el hecho de ser un objeto simétrico, tomará los mismos valores en su puntos diametralmente opuesto, solo que en este caso su velocidad irá disminuyendo y su presión aumentando, hasta volver a una trayectoria horizontal paralela y próxima al eje del diámetro horizontal del objeto circular.

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

Visualización de campos escalares y vectoriales en fluidos. (Grupo 23-C)

2013-12-10T02:31:53Z

Vicente:

{{beta}}
{{Trabajo|Visualización de campos escalares y vectoriales en fluidos. Grupo 23-C|[[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:Trabajos 2013-14|2013-14]]}}

==Introducción==
Vamos a estudiar el flujo de un fluido incompresible alrededor de un obstáculo con forma circular de radio igual a 1, ósea un círculo unidad. Para la realización de este trabajo usaremos tanto coordendas cartesianas como con coordenadas Cilíndricas (polares).
== Mallado ==
La región ocupada por el fluido será el exterior del círculo unidad tal y como indica a continuación el mallado que hemos dibujado.
Hemos tomado un mallado del anillo comprendido entre los radio 1 y 5 y centro el origen, con unos ejes dentro del intervalo [-4,4]x[-4,4].

El código matlab con el que lo hemos realizado es el siguiente:
{{matlab|codigo=
u=linspace(1,5,100); % Creamos el vector u en el intervalo[1,5]
v=linspace(0,2*pi,100); % Creamos el vector v en el intervalo [0,2*pi]
[uu,vv]=meshgrid(u,v); % matrices de las coordenadas u y v
xx=uu.*cos(vv); % Parametrizamos en coordendas cartesianas y así definimos el anillo
yy=uu.*sin(vv);
mesh(xx,yy,0*xx) % Dibujamos el mallado
hold on
x=cos(v);
y=sin(v);
plot(x,y,'k','Linewidth',1) % Dibujamos el obstáculo: círculo unidad
axis([-4,4,-4,4]) % Definimos los ejes en [-4,4]x[-4,4]
view(2) % Visualizamos lo dibujado en el plano x e y
}}
[[Archivo:Malladoc.png|marco|centro|Mallado de la región]]
== Velocidad de las partículas del fluido ==
La velocidad de las partículas del fluido viene dada por un campo vectorial <math>\vec u </math> cuya función potencial es φ=(ρ+1/p)*cos⁡θ. Luego <math>\vec u</math> es el gradiente de la función potencial: : <math>\vec u=∇φ</math>
Realizaremos los cálculos tanto en cartesianas como en cilíndricas: : <math>φ(ρ,θ)=(ρ+1/p)*cos⁡θ = φ(x,y)=x+\frac{x}{x^2+y^2}</math>
En el caso de Cartesianas:: <math>φ(x,y)=x+\frac{x}{x^2+y^2}</math>

: <math>\vec u=∇φ(x,y)=\frac{∂φ}{∂x} \vec i +\frac{∂φ}{∂y} \vec j =\frac{1+(y^2-x^2)}{(x^2+y^2 )^2 } \vec i -\frac{(2xy}{(x^2+y^2 )^2 )} \vec j</math>

En el caso de Cilíndricas calculamos tanto las coordenadas contravariantes como las covariantes: : <math>φ(ρ,θ)=(ρ+1/p)*cos⁡θ</math>
: <math>\vec u=∇φ(ρ,θ)=\frac{∂φ}{∂ρ} \vec g^ρ +\frac{∂φ}{∂θ} \vec g^θ =cos⁡θ (1-\frac{1}{ρ^2 }) \vec g^ρ -sin⁡θ \frac{ρ+1}{ρ} \vec g^θ =cos⁡θ(1-\frac {1}{ρ^2} ) \vec g_ρ -\frac{sin⁡θ}{ρ^2} \frac{ρ+1}{ρ} \vec g_θ</math>

Las representaciones de estos dos campos, el campo vectorial <math>\vec u </math> y el campo escalar φ, función potencial de <math>\vec u </math>, son las siguientes:

{{matlab|codigo=
%Representación del Campo escalar φ
u=linspace(1,5,100);
v=linspace(0,2*pi,100);
[uu,vv]=meshgrid(u,v); %mallado sobre el que representar
xx=uu.*cos(vv);
yy=uu.*sin(vv);
fi=inline('x+x./(x.^2+y.^2)','x','y'); %Definimos la función potencial fi
f=fi(xx,yy);
contour(xx,yy,f,30) %Dibujamos sus líneas de nivel
axis([-4,4,-4,4])
hold on
x=cos(v);
y=sin(v);
plot(x,y,'k','Linewidth',1) % Dibujamos el campo escalar fi
axis([-4,4,-4,4]) % Seleccionamos la región donde lo dibuja
axis equal

%Representación del Campo vectorial u
u=linspace(1,5,15);
v=linspace(0,2*pi,15);
[uu,vv]=meshgrid(u,v); %mallado sobre el que representar
xx=uu.*cos(vv);
yy=uu.*sin(vv);
fx=inline('1+(y.^2-x.^2)./((y.^2+x.^2).^2)','x','y'); %Definimos la función del campo u de velocidades en x
fy=inline('-2*(x.*y)./((y.^2+x.^2).^2)','x','y'); %Definimos la función del campo u de velocidades en y
u1=fx(xx,yy);
u2=fy(xx,yy);
quiver(xx,yy,u1,u2) %dibujamos el campo vectorial u
axis([-4,4,-4,4])
view(2) %Para visualizar lo dibujado

hold off
}}
[[Archivo:Campo ortogonal a.png|marco|centro|Campo vectorial <math>\vec u </math> ortogonal a las curvas de nivel de φ ]]

Como podemos observar, el campo vectorial <math>\vec u </math> de las velocidades de las partículas es ortogonal a las curvas de nivel de φ .
[[Archivo:Campo ortogonal a más de cerca.png|marco|centro|Campo <math>\vec u </math> ortogonal a las líneas de nivel de φ zoom]]
Esto se debe a que
== Rotacional de <math>\vec u</math> nulo==

El rotacional de un campo nos indica la tendencia de este a inducir rotación alrededor de un punto.

El rotacional de nuestro campo <math>\vec u</math> es nulo ya que todo campo potencial, es decir, expresable como el gradiente de un potencial escalar es irrotacional.

Vamos a comprobarlo entonces analíticamente realizando los correspondientes cálculo en coordenadas cilíndricas:

<math>\nabla\times\vec{u}= \frac {1}{g^\frac{1}{2}}\begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\ \frac{ \partial}{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z}\\ u_1 & u_2 & u_3 \end{pmatrix}= \frac {1}{ρ} \begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\ \frac{ \partial}{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z}\\ cos(θ)(1 - \frac {1}{ρ^2})& - sen (θ) (ρ + \frac {1}{ρ}) & 0 \end{pmatrix}=</math>
<math>=\frac {1}{ρ}(\frac{\partial}{\partial z}((cos(θ)(1 - \frac {1}{ρ^2}))\vec{g_θ} ) - (\frac{\partial}{\partial ρ}(sen (θ) (ρ + \frac {1}{ρ}))\vec{g_z})-(\frac{\partial}{\partial θ}((cos(θ)(1 - \frac {1}{ρ^2}))\vec{g_z} + (\frac{\partial}{\partial z}(sen (θ) (ρ + \frac {1}{ρ})\vec{g_ρ}) )=</math>
<math>\frac {1}{ρ}(-sen(θ) (1- \frac{1}{ρ^2}) \vec{g_z} + sen(θ) (1- \frac{1}{ρ^2}) \vec{g_z})= 0 </math>

== Divergencia de <math>\vec u</math> nula==

Que la divergencia de <math>\vec u</math> sea nula se debe a que se trata de un fluido incompresible que implica que este tiene densidad constante en todo el conjunto y a lo largo del tiempo.

La divergencia mide la diferencia de flujo saliente y flujo entrante de un campo vectorial sobre la superficie que rodea un volumen control. Si la divergencia en un punto es positiva, se dice que el campo posee fuentes. Si la divergencia es negativa, se dice que tiene sumideros.

Vamos a comprobarlo analíticamente realizando los orrespondientes cálculo en coordenadas cilíndricas::

<math>\nabla·\vec{u}= \frac {1}{g^\frac{1}{2}}\frac{\partial }{\partial x^i}(g^\frac{1}{2} u^i)= </math>
<math>= \frac {1}{ρ}(\frac{\partial }{\partial ρ}(ρ u^ρ) + \frac{\partial }{\partial θ}(ρ u^θ) + \frac{\partial }{\partial z}(ρ u^z))=</math>
<math>\frac {1}{ρ}(\frac{\partial }{\partial ρ}((ρ - \frac{1}{ρ})cos(θ)) - \frac{\partial }{\partial θ}(1 + \frac{1}{ρ^2})sen(θ) + 0 )=</math>
<math>((1 + \frac{1}{ρ^2})cos(θ) - (1 + \frac{1}{ρ^2})cos(θ))= 0</math>

==Líneas de corriente==
Las líneas de corriente del campo <math>\vec u</math> son las líneas que son tangentes a <math>\vec u</math> en cada punto.
Para calcularlas primeramente hallamos el campo <math>\vec u</math> que en cada punto es ortogonal a <math>\vec u</math> , luego hacemos:
<math>\vec{v}=\vec{k}\times\vec{u}</math>
: <math>\vec{v}=\vec{k}\times\vec{u}= {g^\frac{1}{2}}\begin{pmatrix} \vec{i} & \vec{j} & \vec{k} \\ 0 & 0 & 1 \\ u^1 & u^2 & u^3 \end{pmatrix}= 1 (\vec{u^1} \vec{j} - \vec{u^2} \vec{i}) = (\frac{-2xy}{(x^2+y^2 )^2})\vec{i} + (1 + \frac{y^2-x^2}{(x^2+y^2 )^2})\vec{j}</math>

Este nuevo campo vectorial <math>\vec v</math> es irrotacional, es decir que su rotaciones es 0; como demostramos a continuación ya que <math>\vec u</math> es de divergencia nula.

<math>\nabla\times\vec{v}= \frac {1}{g^\frac{1}{2}}\begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\\frac{\partial }{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z} \\ v_1 & v_2 & v_3 \end{pmatrix}=</math>
<math>=\frac{1}{ρ}\begin{pmatrix} \vec{g_ρ} & \vec{g_θ} & \vec{g_z} \\ \frac{\partial }{\partial ρ} & \frac{\partial }{\partial θ} & \frac{\partial }{\partial z} \\ (ρ+\frac{1}{ρ}) \frac{sen(θ)}{ρ} & (ρ-\frac{1}{ρ}) cos(θ) & 0 \end{pmatrix}=</math>
<math>=\frac{1}{ρ}((1+\frac{1}{ρ^2})cos(θ)- (ρ+\frac{1}{ρ})\frac{cos(θ)}{ρ})\vec{g_z}= \frac{1}{ρ}((1+\frac{1}{ρ^2}) - (1+\frac{1}{ρ^2}))cos(θ)\vec{g_z} = 0</math>

Y tiene un potencial escalar <math>ψ</math> que se conoce como función de corriente de <math>\vec u</math> y que vamos a calcular:

: <math>\vec{v}=\vec{grad}ψ= \frac{∂ψ}{∂x}\vec{i} + \frac{∂ψ}{∂y}\vec{j} + \frac{∂ψ}{∂z}\vec{k} = (\frac{-2xy}{(x^2+y^2 )^2})\vec{i} + (1 + \frac{y^2-x^2}{(x^2+y^2 )^2})\vec{j}</math>

De esta fórmula anterior podemos sacar que:
: <math>\frac{∂ψ}{∂x}= \frac{-2xy}{(x^2+y^2 )^2}</math>
: <math>\frac{∂ψ}{∂y}= -1- \frac{y^2-x^2}{(x^2+y^2 )^2}</math>
Despejando <math>ψ</math> de <math>\frac{∂ψ}{∂x}</math> e Integrando obtenemos la función de <math>ψ(x,y)</math>:
<math>ψ=\int \frac{-2xy}{(x^2+y^2 )^2} dx= -\frac{y}{x^2+y^2} + h(y)</math>

A continuación derivamos la <math>ψ</math> obtenida respecto "y" y la igualamos a la expresión que teníamos anterior de <math>\frac{∂ψ}{∂y}</math> de manera que así podemos sacar <math>h(y)</math> y completar la función <math>ψ</math>
: <math>\frac{∂ψ}{∂y}=-\frac{y^2-x^2}{(x^2+y^2 )^2}+h'(y)= -1- \frac{y^2-x^2}{(x^2+y^2 )^2}</math>
Despejamos <math>h'(y)</math> y sacamos integrándola <math>h(y)</math>: <math>h'(y)=-1</math>
: <math>h(y)=-y</math>
Por lo tanto el potencial escalar <math>ψ</math> es igual a: <math>ψ=-\frac{y}{x^2+y^2}-y</math>

== Estudio de la velocidad del fluido ==

A continuación se realiza el estudio de la velocidad del fluido en los puntos de la frontera S, que son aquellos en los que se particulariza la función ϒ para ρ=1, definida por la función determinada por el gradiente de la función ϒ.

Para determinar los valores máximo y mínimo de la velocidad consideramos su módulo:

<math>\vec{v}= cos(θ)(1-\frac{1}{ρ^2})\vec{g_ρ} - sen (θ) (\frac{1}{ρ} + \frac{1}{ρ^3})\vec{g_θ}</math>

: <math>|\vec{v}|=((1- \frac{1}{ρ^2})^2 (cos(θ))^2 + (\frac{1}{ρ}+\frac{1}{ρ^3})^2 (sen(θ)^2))^(\frac{1}{2})=f(ρ,θ)</math>

Los valores máximo y mínimos se obtienen sustituyendo en la f(ρ,θ) para ρ=1:

<math>f(1,θ)=(4(sen(θ))^2)^(1/2)=2|sen(θ)|</math>

Por tanto, la velocidad será máxima cuando la función seno tome su máximo, es decir, cuando toma el valor 1:
<math>|\vec u|(\frac{ \pi}{ 2})= 2</math>
: <math>|\vec u| (\frac{ 3\pi}{ 2})= 2</math>

En los puntos en coordenadas polares (1,π⁄2) y (1,3π⁄2), la velocidad máxima será igual a 2.
Por su parte, la velocidad será mínima cuando la función seno tome su valor mínimo, es decir, cuando toma el valor 0:
<math>|\vec u| (0)= 0 </math>
: <math>|\vec u| (\pi)= 0 </math>

En los puntos en coordenadas polares (1,0) y (1,π), la velocidad mínima será igual a 0. Dichos puntos serán puntos de remanso, pues en ellos su velocidad es nula.

== Estudio de la presión del fluído ==

Basándonos en la ecuación de Bernouilli:

<math> \frac{1}{2}ρ (|\vec{u}|)^2 + p = cte </math>

con densidad <math>ρ = 2</math>.
La mayor presión se alcanzará al tener la menor velocidad y la menor presión se alcanzará al tener la mayor velocidad,y en el gráfico de la presión, se pueden distinguir el valor de la presión en cualquier punto del fluido.
[[Archivo:Grafica Bernouilli.png|marco|centro|Gráfica de la ecuación de Bernouilli]]
Si se consideran los puntos de mínima velocidad (v = 0), el primer término se anula y entonces la presión alcanza su valor máximo, que se corresponde con el valor de la constante. Por el contrario, si se toman los puntos de velocidad máxima se tendrá el menor valor de presión. Sabiendo que la velocidad máxima es 2 y la densidad del fluido es 2, el valor del primer término es 4 y se puede determinar que el valor de la constante tiene que ser mayor que 4 para que se obtenga un valor de presión positivo, en este caso usamos un valor de 50 para poder observar bien la gráfica.

== Trayectoria de una partícula del fluído ==

Supongamos que nuestra partícula va por una trayectoria horizontal paralela y próxima al eje del diámetro horizontal del objeto circular, por el hecho de estar en el fluido llevara una determinada velocidad y tendrá una determinada presión, al aproximarse al objeto empezara a curvar su trayectoria con el objetivo de rodear el objeto, como hemos supuesto que lleva una trayectoria próxima al diámetro horizontal, mas o menos recorrerá el objeto por su superficie <math> ρ=1 </math>, a medida que empiece a rodear el objeto empezara a disminuir su presión y a aumentar su velocidad, hasta el punto de llegar a ser tangente a la circunferencia (corte con el diámetro vertical) dónde la velocidad será máxima, una vez alcanzado este punto empezará a disminuir, además por el hecho de ser un objeto simétrico, tomará los mismos valores que cuando estaba haciendo la "subida", solo que en este caso su velocidad irá disminuyendo y su presión aumentando, hasta volver a una trayectoria horizontal paralela y próxima al eje del diámetro horizontal del objeto circular.

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

Visualización de campos escalares y vectoriales en fluidos. (Grupo 23-C)

2013-12-10T01:30:53Z

Vicente:

Visualización de campos escalares y vectoriales en fluidos. (Grupo 23-C)

2013-12-10T00:14:08Z

Vicente:

Usuario discusión:Vicente

2013-12-09T23:57:49Z

Vicente: Nueva sección: /* Visualización de campos escalares y vectoriales en fluidos (Grupo C-23) */

'''Bienvenido a ''MateWiki''!'''
Esperamos que contribuyas mucho y bien.
Probablemente desearás leer las [[Help:Contenidos|páginas de ayuda]].
Nuevamente, bienvenido y diviértete! [[Usuario:Herraiz|Herraiz]] ([[Usuario discusión:Herraiz|discusión]]) 10:49 30 nov 2013 (CET)

== Visualización de campos escalares y vectoriales en fluidos (Grupo C-23) ==

{{Beta}}