Buenas tardes, nos surgió una duda con la discretización de Euler en la primera parte del tp. Cuándo vamos a aplicarla en H1 y H2 nos queda que f1=-2Q/omega1 y f2=2Q/omega2 respectivamente. Entonces no dependen de H y de t. Tenemos que poner el caudal en función de H1 y H2? O simplemente no depende y lo único que va a variar es Hn y siempre vamos a sumarle el mismo hf. Si es así cómo sacaríamos el paso h? porque nos queda En+1=En.
Gracias! Saludos
O deberíamos poner f1=-2Qn/omega1 y entonces variaría también ese término. Lo mismo para f2. Pero igualmente no sabemos cómo calcular el paso h, y si hay que agregarle además tn a f
Guido, ante todo el paso h se elige. Deberán luego chequear si hay problemas de estabilidad con el paso elegido.
No hay que hallar ningun caudal en funcion de la altura. El ejemplo para H1 y metodo de Euler es el siguiente:
H1(n+1) = H1(n) + h * f(H1(n), H2(n), Q(n),tn)
Sucede en este caso particular que f no depende de H1 ni de H2 ni de t. Solo es función del caudal que en la condicion inicial es nulo: Q(n=0)=0.
Asi se debe seguir discretizando las otras 2 variables.
Saludos