jueves, 8 de noviembre de 2012

Algebra vectorial con Python

Para acordarme de mis 49 años en actividades docentes, decidí tomar un curso de programación por Internet. 
En marzo de 1964, me sobraba un poco de tiempo en la UCR (había perdido una asignatura que me dejó renco en el plan de estudios), visité entonces a Manuel Enrique Castellón, por ese entonces Asesor de Matemáticas en el MEP y con el beneplácito de la directora del Liceo Rodrigo Facio (Lía Gómez de Brenes) di mis primeras lecciones de Matemática y Física en esa institución (no graduado; 8 lecciones x ¢ 22,50= ¢1800= $30, menos rebajos).  

En Física comenzábamos enseñando vectores, quizás solo un tratamiento de flechitas de diferente tamaño, para un lado y otro, trazadas con regla y transportador, apoyadas en el teorema de Pitágoras, de cosenos, de senos y la famosa ayudante de cálculo de esos tiempos, las Tablas Usuales. Los estudiantes le temían a este tema y recuerdo a un colega del Liceo de Costa Rica que le apodaron “vector” y cuando había un alumno nuevo, algo despistado, lo engatuzaban para que fuera a preguntarle algo, llamándole "don vector". En el LRF yo pasé a ser “coquito”, porque a la primera generación de ese liceo también los atormentaba la geometría y la física, y porque ya se veía que me iba a quedar sin cabello. 

Bueno el álgebra vectorial, independientemente de lo árida o agradable que se enseñe es útil para el manejo de cantidades que tienen magnitud y dirección, por ejemplo velocidades, aceleraciones y fuerzas de toda índole (de contacto, gravitacionales, electromagnéticas). Se requiere a veces encontrar la magnitud y dirección (ángulos) de vectores, sumarlos y restarlos. Algunos conceptos como el de trabajo, se definen en términos de un producto no direccionado de la fuerza y el desplazamiento (producto escalar) y otros como el torque, la cantidad de movimiento angular y la fuerza que un campo magnético ejerce sobre una carga en movimiento, se definen mediante un producto direccionado (producto vectorial).


Creo que esas operaciones y muchas otras debemos ser capaces de hacerlas con lápiz, papel y calculadora en cualquier lugar que estemos, son parte vital del entrenamiento de físicos, matemáticos e ingenieros. Esto no solo nos aportan la herramienta, sino además enriquece una línea de razonamiento aplicable en otras disciplinas y en la vida diaria.

El programa algebra_vectorial.py que le ofrezco aquí, se lo recomiendo para –verificar-  el resultado de sus operaciones con vectores, o cuando tenga alguna duda, luego de haber trabajado un buen rato en un problema de naturaleza vectorial. 
Quizás le ayude en décimo y undécimo año de Enseñanza Media y en la Física General universitaria. 
Si yo fuera usted, siempre lo tomaría como un apoyo, o último recurso, pues me sentiría orgulloso y satisfecho, si puedo maneja el álgebra vectorial con mis propios recursos mentales.
Con un poquito de conocimientos básicos sobre programación con Python, o cualquier otro lenguaje que prefiera, usted mismo lo puede codificar y mejorar. A lo mejor el mío tiene algunas “pulgas informáticas”, si las descubre por favor envíeme un comentario.
Recuerde además que con las operaciones básicas que tiene el programa, puedes hacer otras cosas usando la herramienta (o extendiendo usted mismo el programa), como:
  • Ángulo entre dos vectores, usando producto escalar (o vectorial).
  • Probar si dos vectores son paralelos o perpendiculares (con los productos).
  • Encontrar un vector unitario perpendicular a dos vectores dados (producto vectorial/magnitud).
Copie el código siguiente en el IDLE de Python, sálvelo, corra el módulo y úselo. Si quiere la lista de códigos (--.py) indíqueme una dirección de correo, por medio de un comentario y con gusto se lo envío.

El curso Learn to Program: The Fundamentals, de siete semanas, lo ofrece coursera. Lo dirigen de manera excelente, los profesores de University of Toronto, Jennifer Campbell y Paul Gries, se lo recomiendo.
Llegué casi hasta el final, (video clases, quices, examenes, proyectos), pero no lo terminé. Cometí la novatada de anotarme en tres más, incluyendo uno de Astronomía (Introduction to Astronomy), que espero iniciar el 27 de este mes.
Pero sí me enseñaron y aprendí las bases para continuar el estudio por mi propia cuenta, quizás con la colaboración de algún compañero de mi barrio.

print("ALGEBRA VECTORIAL (jav): Magnitud, ángulos, suma, resta y productos.")  
import math
print("")
print("Llame la función 'magnitud'.\nEscriba las coordenadas cartesianas [float]: positivos, negativos o cero;\ndos decimales.")     
print("Llame la función 'angulos'; 'A' es el ángulo con el eje-x y 'B' con el eje-z.")
print("Llame la función 'mas_menos'. Introdusca coordenadas.")
print("Llame la función 'productos'. Introduzca coordenadas.")

def magnitud(x,y,z):
    mg = math.sqrt(x**2 + y**2 + z**2)
    mag= round(mg,2)
    print("magnitud del vector", x,y,z)
    return mag

def angulos(x,y,z):
    Ar= math.atan2(y,x)
    A= math.degrees(Ar)
    Br= math.atan2(math.sqrt(x**2 + y**2), z)
    B= math.degrees(Br)
    print ("los ángulos A y B  en grados del vector", x,y,z, "son")
    return round(A,2), round (B,2)
          
def mas_menos(x1,y1,z1,x2,y2,z2):
    suma = x1 + x2, y1 + y2, z1 + z2
    resta = x1 - x2, y1 - y2, z1 - z2
    print ("la suma y la resta de los vectores", x1,y1,z1, "y", x2,y2,z2, "es:")
    return suma, resta

def producto(x1,y1,z1,x2,y2,z2):
    punto = x1*x2+y1*y2+z1*z2
    crus = round(y1*z2-y2*z1,3),round(x2*z1-x1*z2,3),round(x1*y2-x2*y1,3)
    mag1= math.sqrt(x1**2+y1**2+z1**2)
    mag2= math.sqrt(x2**2+y2**2+z2**2)
    magcrus= math.sqrt((y1*z2-y2*z1)**2+(x2*z1-x1*z2)**2+(x1*y2-x2*y1)**2)
    ang12r= math.acos(punto/(mag1*mag2))
    angulo12= math.degrees(ang12r)
    angle21r= math.asin(magcrus/(mag1*mag2))
    angle21= math.degrees(angle21r)
    print ("producto escalar", round(punto,2), "; producto vectorial", crus)
    print("ángulo V1->V2:",round(angulo12,2), round(angle21,2))
    print ("Si los valores son diferentes, analice un diagrama vectorial y decida.")

jueves, 1 de noviembre de 2012

Cometa y eclipse

El jueves 24 de octubre, durante la actividad para presentar el
Lunario de Cientec y conmemorar los 20 años de publicación ininterrumpida del mismo, mi buena amiga y colega de la Escuela de Física, en la Universidad de Costa Rica, Luz María Moya R., me hizo la siguiente pregunta:
Si yo sabía de un cometa que por los años 1947 a 1948, cuando ambos éramos chiquillos de escuela, que ella observó, desde--.


Bueno quiero agradecerle a “Luzma”, el compartir ese recuerdo conmigo y sin querer dejarme una agradable y entretenida tarea, con la que sigo aprendiendo más sobre fenómenos astronómicos.

Resulta que ahora con la tecnología informática y el acceso a bibliotecas virtuales, se puede hacer un poco de investigación histórico-científica muy interesante.
Primero le voy a echar por delante el corolario de esta entrada, que supongo lo sabemos, pero no lo tenemos tan a flor de piel, por si puede aplicarlo en Australia el 13 de noviembre, o luego:
“Durante un eclipse total de Sol, es posible ver o descubrir un cometa.”

Cometa eclipse
Ahora si desde el principio.
En 1948 no  fue la primera vez, que un cometa coincid con un eclipse, que, desde luego, es un evento raro.
Hay un reporte de observación desde Egipto, el
17 de mayo 1882, desde Chile el 16 de abril de 1893, desde Brasil, el 20 de mayo de 1947 y otro desde un sitio no especificado el 20 de julio de 1963.
Quizás recuerde la coincidencia algo reciente del
Cometa Hale-Bopp y el Eclipse Total de Sol del 9 de marzo de 1997.

El día primero de noviembre de 1948, hace hoy 64 años, ocurrió un Eclipse Total de Sol (número 19 de los 72 eclipses en  la serie Saros 142). Inició en Kenia, África y corrió por el Océano Índico hasta el sur de Australia, pero no fue de mucha duración, solo 115 segundos.
Desde luego, el eclipse y el cometa (en el momento de su descubrimiento) no podían verse desde Costa Rica, pues estaban situados para observadores en el hemisferio sur.
Observadores del eclipse en Nairobi, Kenia, descubren un cometa, a solo dos grados del Sol. Bastante brillante (m= -2) pues cuando se calculó su órbita se encontró que el perihelio había ocurrido el 27 de octubre.
Se le llamó Cometa Eclipse; O.S. 1948 XI, o C/1948 V1, lo que indica que es no periódico, o que la incertidumbre en su posible periodo es más de 200 años.

C/1948 V1
Su visibilidad durante crepúsculo matutino, mejoró a principios de diciembre, alcanzando magnitud cero y una cola de unos 30°. Ya para finales de diciembre había dejado de ser visible a simple vista. Dejó de verse hasta el 13 de abril de 1949. Fue reportado en la página 59-60 de la revista Sky & Telescope de enero de 1949.

Si introduce los elementos orbitales del cometa en un programa que genere efemérides, como Starry Night, encontrará que para el 16 de diciembre de 1948  a las 5:30 de la mañana sus coordenadas ecuatoriales (ascensión recta y declinación), coinciden con una región en la constelación Puppis, señalada en la figura 1.
P
erfectamente visible al suroeste por esa fecha en nuestro país. Es posible que Luz María haya visto al Cometa Eclipse de 1948 en una posición bastante cercana.



Referencias adicionales:

viernes, 26 de octubre de 2012

Telescopio observa un agujero negro gigante en el centro de la galaxia:

¿Cómo entender esto, si se dice que "de un agujero negro ni la luz puede salir"?
Bueno primero comencemos por caracterizar de una manera más precisa lo que es, o más bien lo que puede hacer un agujero negro (desde luego, hasta donde mi conocimiento básico me da), analizando algunos conceptos físicos simples que son aplicables a los agujeros negros:
  1. Velocidad de escape.
    Es la velocidad (ve)  que debe tener una partícula de masa m, para escapar de la atracción gravitacional de un cuerpo de masa M, si está a una distancia r del centro de éste.
    En Textos de física encuentra la derivación estricta de la fórmula de cálculo. De una manera simple digamos que la energía cinética de la partícula, requerida para escapar es:


  2. La energía potencial gravitatoria, no con referencia a la superficie del objeto de masa M, como se hace en el colegio, sino con referencia a una posición infinitamente alejada, donde la influencia de M es cero (Gravitación Universal) es:
     
  3. La energía mecánica total es la suma de las anteriores:
     
  4. Cuando la partícula realmente haya escapado, a una distancia infinita, ya no necesita moverse, su velocidad es cero, lo mismo que su energía potencial y mecánica. Entonces, por conservación de la energía:





    de donde se despeja la velocidad de escape:




  5. Observe que solo depende de la masa  M y de la distancia r (¡solo depende del agujero negro, para el caso de interés!)

    Usted puede entretenerse un rato calculado la velocidad de escape de una partícula que está en la atmósfera en un planeta, o en el Sol, ya sea desde la superficie, o desde una distancia r del centro. Solo necesita algunos datos de masa y radio de objetos del Sistema Solar.
    Para escapar de la superficie de la Tierra se necesita al menos
    11,2 km/s.
  6. La velocidad de escape de un fotón sería obviamente la velocidad de la luz c = 300 000 km/s.
    ¿Desde qué distancia (mínima) r del centro de un objeto de masa M podría escapar un fotón? Sería:

  7. Por ejemplo, para un agujero negro de 100 veces la masa del Sol (Ms= 1,99 x1030 kg), en unidades del Sistema Internacional ese valor es:


    ! unos 29 kilómetros!


  8. Esto sería el tamaño (radio máximo =  r) de ese agujero negro, o mejor digamos, el radio de una esfera (imaginaria) que define el llamado horizonte de eventos, afuera del cual, un fotón sí podría escapar (si se mueve en la dirección correcta, pero si está dentro de esa esfera, no podrá hacerlo.
    A eso se refiere la frase: “ni la luz puede escapar de un agujero negro.”
  9. Entonces, los agujeros negros no son puntos matemáticos, existe una distancia a la cual su influencia es fatal (el radio del horizonte de eventos), pero lo que esté fuera de ese horizonte, dependiendo lo que sea, más o menos está a salvo “y podría observarse”.
  10. El campo gravitatorio de un cuerpo de masa M  a una distancia r de su centro está dado por la ley de Newton de Gravitación Universal:

    Puede entretenerse un rato calculándolo para la superficie de cada planeta. Para la Tierra ese valor es el conocido 9,8 m/s2encontrará además que a una altitud de 2 radios terrestres encima de la Tierra, más o menos al principio de la estratosfera, es la cuarta parte!
  11. ¿Cuánto sería g para el citado agujero negro?

    mil quinientos cincuenta y tres millones de veces lo que usted está sintiendo ahora mismo, debido al campo gravitatorio de la Tierra.
    http://apod.nasa.gov/apod/ap060701.html

     
  12. Aún a una distancia de 10 veces el horizonte de eventos, ese campo solo ha decaído a una centésima parte y es capaz de producir aceleraciones increíblemente grandes a las partículas que caen atrapadas por el agujero negro, hasta puede despojar de gas y polvo la atmósfera de una estrella vecina.

    Cuando estas partículas aceleradas tienen carga eléctrica, emiten radiación electromagnética, (rayos x) que puede observarse mientras no estén dentro del horizonte de eventos.
    Esta es una de las maneras indirectas de observar un agujero negro, solo se requiere un telescopio de rayos x, que pueda enfocarlas y entregar una buena imagen. Parece que el telescopio NuSTAR puede hacerlo.
Referencias adicionales:

lunes, 15 de octubre de 2012

¡No es una caída libre, no es un salto al vacío!

La extraordinaria hazaña de Félix Baumgartner al romper la velocidad del sonido en el aire (343,2 m/s) es ciertamente una proeza, pues para alcanzar unos 1342 km/hora (372,8 m/s), se requiere caer en la atmósfera terrestre a lo largo de una gran distancia, tener un excelente sistema de frenado (el paracaídas principalmente) y de soporte vital para un ser humano, además de el diseño y construcción del globo aerostático, que con base en el Principio de Arquímedes, lo subió a 39045 km de altitud. Esto es solo la parte física o de ingeniería del asunto, queda todo el coraje, habilidad, valor y decisión para realizar el salto, condiciones que casi nadie tiene.

Si ustedes recuerdan su física elemental del colegio o primeros años de universidad, estarán de acuerdo conmigo  en que, ese salto fuera de lo común, no es una caída libre, un caso hipotético de caída, en la cual –no se toma en cuenta el  efecto de fuerza de rozamiento del aire
La caída del señor Baumgartner  es una caída real, afectada por todo lo que puede hacer la gravedad y la atmósfera.

Si fuera una caída libre, en los primeros 32,8 segundos se habría alcanzado la velocidad record reportada (1342 km/h) y bastarían solo 7,12 km para alcanzarla. Veamos la física del asunto:
  1. La fuerza que acelera a Félix hacia abajo es la fuerza de gravedad (su peso), pero como sabemos desde Galileo y Newton, la aceleración tiene un valor constante, dado por el campo gravitatorio de la Tierra. Solo que al ser la altura bastante considerable, vamos a tomar ese valor no como el conocido 9,8 m/s2, sino el promedio entre la gravedad a la altura máxima (unos 9,72 m/s2) y aquel, lo que da un valor de trabajo, para este análisis, de 9,76 m/s2.
  2. Entonces, si aplica dos de las relaciones más simples con que trabajó en décimo año


    Y t
    omando en cuenta que inicia el descenso desde el reposo (velocidad inicial cero), luego de que la cápsula de ascensión se ha detenido, usted puede probar fácilmente (con un poquito de álgebra) que bastaría un tiempo de caída libre de 38.2 s, durante los cuales desciende una distancia de 7,12 km.
     
  3. Ahora bien como el tiempo y la altitud reportadas son muy  diferentes; 4 minutos 20 segundos (=260 s) y 39,045 km, concluimos que en esta caída nada libre, la influencia de la resistencia del aire fue nada despreciable.
     
  4. Recuerde que la hipotética caída libre de los libros de física, no tiene velocidad límite. Si Félix hubiese caído libremente a lo largo de 39045 m habría alcanzado la sorprendente rapidez de 873 m/s. Pero esto no se puede en la atmósfera terrestre, debido al rozamaiento del aire.
    Además habría tardado (hasta el suelo), solo 89,4 s.
    Para esto último requiere aplicar la relación:
Como ve, esta gran hazaña, no es una caída libre, aunque así lo digan los encargados de la proeza, o lo repita algún periodista.
Me decía una amiga que observó conmigo la transmisión, “que no importa, eso se sobreentiende”, bueno si fuera así ¿para qué nos esforzamos en enseñar y en aprender? ?Será que una vez que ganamos un curso, le ordenamos al cerebro que haga "delete tema*. *".

Sitio de la imagen
Tampoco aunque lo repitan en la televisión y los periódios, -no es un salto al vacío-, a pesar de que puede sonar más heroico o quizás poético. 
La capa inferior de la atmósfera, la troposfera, tiene unos 18 km de altura. En ella es  donde está la mayor parte de los gases, es la más densa y donde ocurren todos los fenómenos meteorológicos comunes, evidentemente no es una región de vacío. Luego viene la estratosfera hasta unos 50 km de altitud, capa en la que inició el salto.
 

Sabemos que a nivel del mar, la presión atmosférica tiene su valor mayor (unos 101 kilopascales, equivalente al conocido 760 mm de mercurio). 
A 40 km de altura la presión ha disminuido a 18,7 kPa (=140 mm de Hg). La densidad del aire, responsable de la fuerza de rozamiento que ayuda a detener la caída de Félix, se comporta de manera similar.

Sitio de la imagen
Entonces esta caída en la atmosfera, está afectada por dos fuerzas; la simple fuerza de gravedad hacia abajo que contribuye al aumento de velocidad y la compleja fuerza de rozamiento con el aire dirigida hacia arriba. Compleja porque depende de varios factores, por ejemplo la forma, diseño y tamaño del traje y del paracaídas. 

Se puede suponer que la fuerza de rozamiento del aire depende del valor de la velocidad elevado al cuadrado, así que si ésta aumenta, también lo hace la fuerza de fricción, como bien lo sabe un motociclista.
Cuando se cae en un fluido como la atmósfera se alcanza una velocidad terminal, que ocurre cuando la fuerza de rozamiento del aire se iguala al peso, en esto se basa el diseño de los paracaídas.
La velocidad record del salto de Félix (372,8 m/s) debió ser medida (o calculada) posiblemente durante los primeros kilómetros de su caída (sin paracaídas), cuando el peso le ganaba al rozamiento. Esa pudo haber sido la velocidad terminal en esta  etapa, aún muy alta para tocar el suelo. 
Entonces se abre el paracaídas y con las nuevas condiciones se vuelve a alcanzar una segunda velocidad terminal que es la que permite el aterrizaje sin mucho riesgo.