Analisis numerico: Sistema de Helmholtz-Thompson

ÍNDICE

Introducción

Nuestro trabajo para la asignatura optativa de Análisis Numérico versa sobre la caracterización y el estudio del comportamiento del sistema de Helmohotz-Thompson definido con esta ecuación diferencial:

[pic 1]

Comenzaremos el estudio del sistema por la parte más fácil , intentar clasificar las partes en que se compone la ecuación diferencial. Para empezar, observamos que contiene tres parámetros que serán fijados por las condiciones iniciales del problema, es decir nosotros deberemos asignarles un valor desde el principio del problema. Estos parámetros descritos son: β, F y ω y deberemos darles un valor inicial.

Si seguimos observando la ecuación a estudiar vemos que contiene una variable x que se encuentra en función del tiempo dentro de la ecuación.

Caso lineal

Pensamos que el estudio del sistema se facilitará si suponemos que el valor de la variable x es constante en el tiempo, con lo que los diferenciales que contienen esta variable serán ceros. Esto es así porque la derivada y la segunda derivada de una constante es cero, ya que como su propio nombre indica no varía su valor en el tiempo.

Al sustituir x por un valor constante resultaría:

[pic 2]

Si observamos la ecuación anterior podemos sacar una conclusión importante, algunos de los dos parámetros que aparecen en la misma, bien  F o bien ω, deben valer cero. Este hecho al principio nos causó confusión ya que cuando nos pusimos a realizar el estudio no deparamos en este matiz fundamental, que explicaremos a continuación.

Debe ser de este modo ya que sino se llegará a una contradicción. Para explicarlo vamos a suponer el caso contrario para que se vea que llegará a algo incoherente. Supongamos que el valor de x es constante, como hemos planteado anteriormente, y que F y ω toman valores distintos a cero. Si como hemos dicho, x es constante la parte izquierda de la igualdad siempre obtendrá el mismo valor, al  estar compuesta de dos valores constantes sin embargo si nos fijamos en la parte derecha de la igualdad vemos que no es constante ya que el valor del seno variará en el tiempo, como ya sabemos, llegaríamos a un absurdo, que la igualdad de la ecuación no se cumple.

Por lo tanto, sabiendo que x es constante, hemos concluido que o bien F o bien ω son cero. Por consiguiente resulta:

[pic 3]

Resolviendo la ecuación anterior podemos observar que las soluciones al sistema de Helmohotz-Thompson cuando la variable x sabemos que es constante son las siguientes.

[pic 4]

Una vez hallada la solución al sistema cuando la variable x es constante, seguimos nuestro estudio observando  el comportamiento de la ecuación  alrededor de los puntos que son las soluciones de la ecuación lineal.

Esto es así para comprobar la estabilidad de los mismos, ya que puede que alguno de ello no sea estable.

Para poder desarrollar nuestro trabajo,  introduciremos en la ecuación los valores de x sumados con valor de Epsilon.

[pic 5]

Caso de [pic 6]

Sustituimos el valor de x por [pic 7] para poder estudiar lo anteriormente descrito. El resultado quedaría de la siguiente forma.

[pic 8]

Desarrollando obtendríamos:

[pic 9]

Ahora como queremos la parte lineal de la ecuación, quitamos el [pic 10] cuadrado y resultaría de esta forma:

[pic 11]

Resolución ecuación diferencial

Queremos conseguir una solución analítica de [pic 12] por lo que resolvemos mediante ecuaciones diferenciales. Transformando la ecuación anterior obtendríamos otra similar que sería:

[pic 13]

Como podemos observar tenemos una ecuación diferencial de primer orden. Para resolverla procedemos a calcular el polinomio característico de la misma.

[pic 14]

Resolviendo esto tendríamos dos soluciones para y.

[pic 15]

Observando dichos resultados podemos observar que los valores de y1 y y2 dependerán del valor de que demos a[pic 16]. Dependiendo el valor que tome [pic 17] la raíz cuadrada tendrá un valor positivo, negativo o cero que implicará directamente a que la ecuación diferencial tenga una solución real, imaginaria o doble.

Por este motivo se procederá a dividir en tres casos el desarrollo de la solución analitica, según sea la raíz positiva, negativa  o cero. Para ello operamos el radicando igualándolo a cero.

[pic 18]

De aquí obtenemos  que habrá tres intervalos según el valor de [pic 19], que [pic 20]>2, que [pic 21]<2 o que [pic 22]=2.

[pic 23]

A partir de ahora, resolveremos la ecuación diferencial dependiendo del valor de [pic 24] según los intervalos que hemos hallado anteriormente.

• β>2

Para [pic 25]>2 sabemos que tendremos dos soluciones reales ya que la raíz es positivo.

[pic 26]

Resolviendo esto tendríamos dos soluciones para y.

[pic 27]

Por lo tanto y quedaría.

[pic 28]

• Β=2

Para [pic 29]=2 sabemos que tenemos una única solución de la ecuación diferencial al anularse la raíz.

[pic 30]

Por lo tanto resulta

[pic 31]

• Β<2

Para [pic 32]<2 sabemos que tendremos dos soluciones las cuales serán imaginarias porque el resultado de la raíz es negativa.

[pic 33]

Resolviendo esto tendríamos dos soluciones para y.

[pic 34]

Vemos el resultado de la ecuación diferencial.

[pic 35]

Para reducir la complejidad en cuanto a notación de las ecuaciones, acordamos dicha igualdad en la definición:

...