Exactamente ESTO es lo que tenia que hacer...

Estoy pensando que tengo que ir a Campus y la sensación térmica es de 31 grados...

(05-11-2012 14:04)Brich escribió:  Exactamente ESTO es lo que tenia que hacer...

Seguramente se puede hacer con Matlab o Octave.

PD: Mas difícil de hacer que arrastrando la formula en celdas, pero a fin de cuentas termina siendo más simple de manejar

Con el mathematica creo que encima podes sacar "reportes" de los datos (ya sea que plotees una funcion y queres quedarte con el grafico, quieras sacar un conjunto de valores,etc...).
Una vez trate de instalarlo.
No andó,pero seguro que era una masa.

Yo lo tengo instalado...al 8, pero no la tengo tan clara como para hacer tantas ecuaciones. Como tendria que hacer para arrastras las ecuaciones onda el excel?...osea...no tengo ganas de hacer 500 ecuaciones a mano

Básicamente: programando. La ecuación la "especificas" una ves, y después "loopeas" o "ciclas" con "for".

ejemplo para tu caso (todo lo que empieza con % es un comentario):

 



%%% Inicialización de variables



dz = .25; % cada paso es 1/4 metros

Nz = 400; % número de pasos en profundidad (debería ser mínimo 100 m)

Nt = 5000; % número de pasos en tiempo

dt = (365*24*60*60)/Nt; %Longitud de cada paso en tiempo en segundos (~ 6.3*10^3 segundos, or ~105 minutos)

K =  2*10^-6; % K canonica es 10e-6 m^2/s

T = 15*ones(Nz+1,Nt+1); % Creo matriz de temperatura con Nz+1 filas y Nt+1 columnas. Inicialmente es 15 en todos lados.

time = [0:12/Nt:12];

T(1, = 15-10*sin(2*pi*time/12);  %Determino la temperatura de la superficie.



%%% Código del método



maxiter = 500 %cantidad de iteraciones máximas, corta antes si converge



for iter = 1:maxiter

    Tlast = T; %Guardo la última suposición

    T(:,1) = Tlast(:,end); %Inicializo la temperatura a t=0  a la última temperatura



    for i=2:Nt+1,

        depth_2D = (T(1:end-2,i-1)-2*T(2:end-1,i-1)+T(3:end,i-1))/dz^2;

        time_1D = K*depth_2D;

        T(2:end-1,i) = time_1D*dt + T(2:end-1,i-1);

        T(end,i) = T(end-1,i); % Fuezo el límite inferior

    end



    err(iter) = max(abs(T(-Tlast());  %Calculo el error



    if err(iter)<1E-4

        break; % Corto, si las soluciones son muy similares, hay convergencia, con un error < 1E-4

    end



end



if iter==maxiter; % si realicé el maximo de iteraciones especificadas...

    warning('Convergencia no alcanzada') %...no llegué a convergencia

end

rrhh: donde preguntarte cuando podes ir a una entrevista con el cliente sin decirte donde queda el cliente no solo es moneda corriente, es el standard

(05-11-2012 17:50).py escribió:  rrhh: donde preguntarte cuando podes ir a una entrevista con el cliente sin decirte donde queda el cliente no solo es moneda corriente, es el standard

Y que te dan 2 dias para elegir cuando vas.... en uno de esos dos dias.

Mas competencia (?

http://www.lanacion.com.ar/1523924-la-pr...ingenieros

(05-11-2012 17:19)Dem0 escribió:  Básicamente: programando. La ecuación la "especificas" una ves, y después "loopeas" o "ciclas" con "for".

ejemplo para tu caso (todo lo que empieza con % es un comentario):

 



%%% Inicialización de variables



dz = .25; % cada paso es 1/4 metros

Nz = 400; % número de pasos en profundidad (debería ser mínimo 100 m)

Nt = 5000; % número de pasos en tiempo

dt = (365*24*60*60)/Nt; %Longitud de cada paso en tiempo en segundos (~ 6.3*10^3 segundos, or ~105 minutos)

K =  2*10^-6; % K canonica es 10e-6 m^2/s

T = 15*ones(Nz+1,Nt+1); % Creo matriz de temperatura con Nz+1 filas y Nt+1 columnas. Inicialmente es 15 en todos lados.

time = [0:12/Nt:12];

T(1, = 15-10*sin(2*pi*time/12);  %Determino la temperatura de la superficie.



%%% Código del método



maxiter = 500 %cantidad de iteraciones máximas, corta antes si converge



for iter = 1:maxiter

    Tlast = T; %Guardo la última suposición

    T(:,1) = Tlast(:,end); %Inicializo la temperatura a t=0  a la última temperatura



    for i=2:Nt+1,

        depth_2D = (T(1:end-2,i-1)-2*T(2:end-1,i-1)+T(3:end,i-1))/dz^2;

        time_1D = K*depth_2D;

        T(2:end-1,i) = time_1D*dt + T(2:end-1,i-1);

        T(end,i) = T(end-1,i); % Fuezo el límite inferior

    end



    err(iter) = max(abs(T(-Tlast());  %Calculo el error



    if err(iter)<1E-4

        break; % Corto, si las soluciones son muy similares, hay convergencia, con un error < 1E-4

    end



end



if iter==maxiter; % si realicé el maximo de iteraciones especificadas...

    warning('Convergencia no alcanzada') %...no llegué a convergencia

end

Gracias loco, todavía no la entiendo muy bien...cuando tenga un tiempo me pongo y te cuento como me fue.

No hay drama. Igual te aviso que en mi vida estudié termodinámica, así que quizás el modelo del código no es exactamente el de tu problema

(05-11-2012 21:58)Dem0 escribió:  No hay drama. Igual te aviso que en mi vida estudié termodinámica, así que quizás el modelo del código no es exactamente el de tu problema

Se nota...por que es Estabilidad II
Igual me viene bien...para ir viendo el tema,
El problema basico es que excel es genial para el tema de las ecuaciones...pero tarda banda con las iteraciones que tiene. Lo deje 4 horas con F9 apretado xD

jajajaja.

PD: Mas allá de mi penoso conocimiento de física y mecánica, el método es diferencias finitas, de eso no te preocupes.

(05-11-2012 22:06)Brich escribió:  pero tarda banda con las iteraciones que tiene. Lo deje 4 horas con F9 apretado xD

Ufff. Lo necesitas urgentemente entonces. Cualquier cosa también tenes Octave, que es muy similar y es gratis.

Tengo un recuperatorio

http://www.clarin.com/politica/servicio-...19526.html

Si, ahora se sumó el Mitre de Tigre.

Ayer fue un quilombo el tren. Tardé como 2 horas en llegar a Lugano.