Usuario:Grupo 4 - Historial de revisiones

Grupo 4: /* Representación gráfica e interpretación del borde de la capa limite */

2016-04-30T01:12:23Z

‎Representación gráfica e interpretación del borde de la capa limite

Grupo 4: /* Resolución numérica y representación gráfica de la componente U1 del campo de velocidades */

2016-04-30T01:06:45Z

‎Resolución numérica y representación gráfica de la componente U1 del campo de velocidades

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434

Grupo 4: /* Interpretación de la derivada de la función de forma f(η) */

2016-04-30T00:54:03Z

‎Interpretación de la derivada de la función de forma f(η)

Grupo 4: /* Resolución mediante Euler. */

2016-04-30T00:46:28Z

‎Resolución mediante Euler.

Grupo 4: /* Resolución mediante Runge-Kutta de orden 4. */

2016-04-30T00:42:33Z

‎Resolución mediante Runge-Kutta de orden 4.

Grupo 4: /* Resolución mediante Euler modificado. */

2016-04-30T00:33:54Z

‎Resolución mediante Euler modificado.

Grupo 4: /* Resolución mediante Euler modificado. */

2016-04-30T00:33:22Z

‎Resolución mediante Euler modificado.

Grupo 4: /* Representación gráfica e interpretación del borde de la capa limite */

2016-04-30T00:10:58Z

‎Representación gráfica e interpretación del borde de la capa limite

Grupo 4: /* Representación gráfica e interpretación del borde de la capa limite */

2016-04-28T12:57:40Z

‎Representación gráfica e interpretación del borde de la capa limite

Grupo 4: /* Representación gráfica e interpretación del borde de la capa limite */

2016-04-28T12:57:22Z

‎Representación gráfica e interpretación del borde de la capa limite

@@ Línea 367: / Línea 367: @@
 [[Archivo:Ejemploreal1.png|sinmarco|centro|460px]]
-<br />
-[[Archivo:APARTADO5.png|sin marco|derecha|460px]]
 {{matlab|codigo=
-clc
+%Se definen los datos del problema
-clear all
+%Vector x (valor de la longitud de la chapa)
 x0=0;
 xn=10;
 h=0.05;
 x=x0:h:xn;
-eta=5.225;
+%Valor de la variable de similaridad en el borde de la capa
 for i=1:length(x)
-     y(i)=eta*sqrt(x(i)/2);
+     y(i)=n*sqrt(x(i)/2);
 end
 hold on
 plot(x,y);

@@ Línea 319: / Línea 319: @@
 [[Archivo:APARTADO4.png|marco|derecha]]
 {{matlab|codigo=
-clc
+%Se definen los datos del problema
-clear all
+%En este caso la variable "t" es la variable de similaridad
- %Datos iniciales
+%Numero de ecuacionesne=3;
-ne=3;
 h=0.05;
 x=[0.05 0.2 0.4 0.6 0.8];
 y0=0;
@@ Línea 330: / Línea 330: @@
 P=length(y);
 N=round((yn-y0)/h);
 w=inline('[f(2);f(3);-1/2*f(1)*f(3)]','t','f');
 f=zeros(ne,N+1);
 for j=1:length(x)
     H2=h*sqrt(2/x(j));
@@ Línea 343: / Línea 345: @@
     end
     Z(j,:)=f(2,:);
     hold on
     plot(y,2*Z(j,:))

@@ Línea 242: / Línea 242: @@
 {{matlab|codigo=
-clc
+%Se definen los datos del problema
-clear all
+%En este caso la variable "t" es la variable de similaridad
-%Datos iniciales
 %Numero de ecuaciones
 ne=3;
-%Tiempo inicial
+%Se define el intervalo temporal
 t0=0;
 tn=20;
 %Tamaño del salto
@@ Línea 257: / Línea 255: @@
 %Determinacion del vector t
 t=linspace(t0,tn,N+1);
-%Vector k, valor inicial de f(3)
+%Determinacion del vector t
 k0=0.1;
 kf=1;
 k=k0:0.01:kf;
 G=length(k);
-%Sistema de ecuaciones
+%Definicion de la función w siendo esta la parametrizacion del sistema en forma general
 w=inline('[f(2);f(3);-1/2*f(1)*f(3)]','t','f');
 f=zeros(ne,N+1);
-%Resolucion
+%Se resuelve el problema mediante el método de Runge-Kutta de orden 4 aplicado a sistemas para k=0.33
-   f0=[0 0 k(24)];
+f0=[0 0 k(24)];
-   f(:,1)=f0;
+f(:,1)=f0;
-   for i=1:N
+for i=1:N
-  K1=w(t(i),f(:,i));
+   K1=w(t(i),f(:,i));
-  t(i)=t(i)+h/2;
+   t(i)=t(i)+h/2;
-  z(:,i)=f(:,i)+K1*h/2;
+   z(:,i)=f(:,i)+K1*h/2;
-  K2=w(t(i),z(:,i));
+   K2=w(t(i),z(:,i));
-  v(:,i)=f(:,i)+K2*h/2;
+   v(:,i)=f(:,i)+K2*h/2;
-  K3=w(t(i),v(:,i));
+   K3=w(t(i),v(:,i));
-  c(:,i)=f(:,i)+K3*h;
+   c(:,i)=f(:,i)+K3*h;
-  K4=w(t(i)-h/2,c(:,i));
+   K4=w(t(i)-h/2,c(:,i));
-  f(:,i+1)=f(:,i)+(h/6)*(K1+2*K2+2*K3+K4);
+   f(:,i+1)=f(:,i)+(h/6)*(K1+2*K2+2*K3+K4);
-   end
+end
- O=ones(1,N+1);
+O=ones(1,N+1);
-    plot(t,f(2,:))
+%Representación de  f´ para k=0,33
-  for j=1:N+1
+hold on
     if abs(f(2,j)-O)<0.01
         t(j)
@@ Línea 290: / Línea 296: @@
      b(m,:)=abs(f(2,m)-O);
 end
 figure(2)
 hold on

@@ Línea 174: / Línea 174: @@
 [[Archivo:EULER33.png|marco|derecha]]
 {{matlab|codigo=
-clc
+%Se definen los datos del problema
-clear all
+%En este caso la variable "t" es la variable de similaridad
-%Datos iniciales
 %Numero de ecuaciones
 ne=3;
-%Tiempo inicial
+%Se define el intervalo temporal
 t0=0;
 tn=20;
 %Tamaño del salto
@@ Línea 189: / Línea 187: @@
 %Determinacion del vector t
 t=linspace(t0,tn,N+1);
-%Vector k, valor inicial de f(3)
+%Se establece el vector K
 k0=0.1;
 kf=1;
 k=k0:0.01:kf;
 G=length(k);
-%Sistema de ecuaciones
+%Definicion de la función w siendo esta la parametrizacion del sistema en forma general
 w=inline('[f(2);f(3);-1/2*f(1)*f(3)]','t','f');
 f=zeros(ne,N+1);
-%Resolucion
+%Se resuelve el problema mediante el método de Euler aplicado a sistemas
 for j=1:G
     f0=[0 0 k(j)];
@@ Línea 206: / Línea 204: @@
     Y(j)=f(2,end);
 end
 hold on
 plot(k,Y,'g')
@@ Línea 213: / Línea 212: @@
 title('Representación de f(20)´ para los posibles k')
 hold off
 for q=1:G
     a(q)=abs(Y(q)-1);

@@ Línea 106: / Línea 106: @@
 [[Archivo:RUNGEKUTTA2.png|marco|derecha]]
 {{matlab|codigo=
-clc
+%Se definen los datos del problema
-clear all
+%En este caso la variable "t" es la variable de similaridad
-%Datos iniciales
 %Numero de ecuaciones
 ne=3;
-%Tiempo inicial
+%Se define el intervalo temporal
 t0=0;
 tn=20;
 %Tamaño del salto
@@ Línea 121: / Línea 119: @@
 %Determinacion del vector t
 t=linspace(t0,tn,N+1);
-%Vector k, valor inicial de f(3)
+%Se establece el vector K
 k0=0.1;
 kf=1;
 k=k0:0.01:kf;
 G=length(k);
-%Sistema de ecuaciones
+%Definicion de la función w siendo esta la parametrizacion del sistema en forma general
 w=inline('[f(2);f(3);-1/2*f(1)*f(3)]','t','f');
 f=zeros(ne,N+1);
-%Resolucion
+%Se resuelve el problema mediante el método de Runge-Kutta de orden 4 aplicado a sistemas
 for j=1:G
     f0=[0 0 k(j)];
     f(:,1)=f0;
 for i=1:N
        K1=w(t(i),f(:,i));
        t(i)=t(i)+h/2;
@@ Línea 146: / Línea 145: @@
     Y(j)=f(2,end);
 end
 hold on
 plot(k,Y,'r')
@@ Línea 153: / Línea 153: @@
 title('Representación de f(20)´ para los posibles k')
 hold off
 for q=1:G
      a(q)=abs(Y(q)-1);
@@ Línea 165: / Línea 166: @@
 }}
-Parametro de k para el cual f'(20) se acerca mas a 1 es k=0.33
+El parámetro de k para el cual f'(20) se acerca mas a 1 es k=0.33
 ===Resolución mediante Euler.===

@@ Línea 355: / Línea 355: @@
 <br /><br />
-[[Archivo:121212.jpg|sinmarco|centro|460px]]
+[[Archivo:Ejemploreal1.png|sinmarco|centro|460px]]
 <br />
 [[Archivo:APARTADO5.png|sin marco|derecha|460px]]

@@ Línea 354: / Línea 354: @@
 En esta última ilustración se observa gráficamente el borde de la capa límite exactamente igual como si se viera en un laboratorio realizando el experimento físicamente con los datos indicados, dando a entender el comportamiento real del fluido en esta situación. Por debajo del borde el fluido comienza a reducir su velocidad hasta llegar a ser nula justo en la superficie de la placa, análogamente, a medida que nos alejamos del origen de ella, el fluido tiende a volver a ganar la velocidad constante anterior al contacto.
 <br />
-[[Archivo:121212.jpg|sinmarco|centro]]
+[[Archivo:121212.jpg|sinmarco|centro|460px]]
 <br />
 [[Archivo:APARTADO5.png|sin marco|derecha|460px]]