Cálculo de un Eclipse Solar y Lunar. Ocultación y Tránsito/Tránsito/Teoría de un Tránsito Planetario por el Sol y Cálculo de los Elementos Besselianos

Capítulo Anterior

Contenidos

Capítulo Siguiente

Teoría de un Tránsito Planetario por el Sol y Cálculo de los Elementos Besselianos

Plano Fundamental o Principal de Referencia según Friedrich Wilhelm Bessel

Plano Fundamental o de Referencia según Friedrich Wilhelm Bessel

Ocurrencia de un Tránsito por el Sol

Habiendo hallado, por interpolación, el instante de la mínima distancia angular del centro del Planeta (Mercurio o Venus) con el centro del Sol en la Conjunción Inferior , se procede al cálculo de la distancia siendo perpendicular al Eje del "Cono de la Sombra" de Mercurio hacia el centro de la Tierra que corre por el Plano Fundamental o de Referencia, esto es, el valor de γ (gamma) en [Radios Terrestres], siendo positivo (+) hacia el Norte del centro de la Tierra y negativo (-) hacia el Sur:

γ = (((ΔTP * Seno(ΔPS)) / Seno((90 - (ΔPS))) / 6.378,14     (1)

Donde ΔPS es la mínima distancia angular en [°] entre el Planeta (Mercurio o Venus) y el Sol y ΔTP es la distancia en [kms] entre la Tierra y el Planeta (Mercurio o Venus) en ese mismo instante.

Luego se determina si ocurre o no un Tránsito por el Sol, según la siguiente condición:

Si -0,9972 < γ < 0,9972 entonces hay un Tránsito por el Sol. Para los tránsitos tangenciales al limbo solar se debe considerar la suma del semidiámetro del Planeta (Mercurio o Venus) al valor de γ.

Contactos Exteriores e Interiores del Planeta (Mercurio o Venus) con el Sol

Contactos Exteriores e Interiores del Planeta con el Sol

Entonces, en un Tránsito por el Sol:
Contactos Interiores del Planeta con el Sol si el observador está en F (DB) o en H (CA).
Contactos Exteriores del Planeta con el Sol si el observador está en E (CB) o en I (DA).
El Tránsito Medio ocurrirá si el observador está en G.
Vp es el vértice del "Cono Penumbral" del Planeta.
Vu es el vértice del "Cono de la Sombra (Umbra)" del Planeta.

Ecuaciones Fundamentales y Cálculos de los Elementos Besselianos. Ejemplos prácticos según el Tránsito de Mercurio por el Sol del 09.05.2016

Plano Fundamental o Principal de Referencia según Bessel

En la primer figura se describe el Plano Fundamental o Principal de Referencia según Bessel que pasa por el centro de la Tierra y siendo siempre perpendicular al Eje del "Cono de Sombra" de Mercurio. La línea OZ es paralela a tal Eje y tiene su origen también en el centro de la Tierra, y apunta hacia la esfera celeste, a una Coordenada Ecuatorial Geocéntrica: Ascensión Recta y Declinación del punto Z.

Con el fin de determinar las Coordenadas Rectangulares de Mercurio y del Sol, el eje de las x va desde el centro de la Tierra hacia el Punto Equinoccial Vernal (Marzo), el eje de las y a 90° del eje de las x y es positivo (+) hacia el punto del Ecuador cuya Ascensión Recta es igual a 90°. El eje de las z es positivo (+) si apunta hacia el hemisferio Norte Celeste.

Efemérides de Mercurio y del Sol

Sabiendo que la Conjunción Inferior Sol-Mercurio, en Ascensión Recta, ocurre a las 15:00:00 hs. (GMT = Greenwich Meridian Time) tomamos 9 horas para los cálculos respectivos. T₀ = 15 hs. es la hora central y anterior más cercana a tal conjunción, luego se realizan los cálculos para ±4 hs. a partir de esa T₀, es decir para las 11 hs., 12 hs., 13 hs., 14 hs., 15 hs., 16 hs., 17 hs., 18 hs. y 19 hs. (GMT).

Ambos astros están en Coordenadas Ecuatoriales Geocéntricas tomadas de las Efemérides diarias (00:00 GMT), y publicadas, por ejemplo, por The Astonomical Almanac donde las Ascensiones Rectas (α) de Mercurio y las del Sol están en el formato Hora, Minutos y Segundos, y las Declinaciones (δ) de Mercurio y las del Sol en el formato °, ' y ". También desde estas efemérides se toman la Paralaje Ecuatorial Horizontal de Mercurio π en [°] y la distancia r' Tierra-Sol en [U.A.].

Coordenadas Ecuatoriales Geocéntricas del Eje del "Cono de la Sombra" de Mercurio o del punto Z

Según la siguiente tabla se tienen las Constantes:

Calcular

b = Seno(π₀) / (r' * Seno(π))     (2)

luego g y G siendo esta última la distancia Mercurio-Sol (Centro a Centro) en [U.A.]

g = 1 - b     (3)

G = r' * g     (4)

después calcular la Ascensión Recta del punto Z o del Eje del "Cono de la Sombra" de Mercurio

a = α' - (b / ((1 - b) * Seno(δ'))) * Coseno(δ) * (α - α')     (5)

y la Declinación del punto Z o del Eje del "Cono de la Sombra" de Mercurio

d = δ' - (b / (1 - b)) * (δ - δ')     (6)

Coordenadas Rectangulares de Mercurio x,y,z

Calcular primero r la distancia Tierra-Mercurio (Centro a Centro) en [Radios Terrestres]

r = 1 / Seno(π)     (7)

luego las Coordenadas Rectangulares Geocéntricas de Mercurio x,y,z

x = r * Coseno(δ) * Seno(α - a)     (8)

y = r * Seno(δ - d) * Coseno((α - a) / 2)^2 + r * Seno(δ + d) * Seno((α - a) / 2)^2     (9)

z = r * Coseno(δ - d) * Coseno((α - a) / 2)^2 - r * Coseno(δ + d) * Seno((α - a) / 2)^2     (10)

Ángulo y Radio del "Cono de la Penumbra" de Mercurio (Contactos Exteriores)

Calcular el ángulo f₁ con vértice en Vp

f₁ = Aseno((Seno(H) + k * Seno(π₀)) / (r' * g))     (11)

luego c₁ u OVp que es la distancia desde el vértice del "Cono de la Penumbra" de Mercurio Vp hasta el Plano Fundamental o Principal de Referencia (O), por lo tanto siendo siempre positiva (+).

c₁ = z + k / Seno(f₁)     (12)

después el valor de i₁

i₁ = Tan(f₁)     (13)

y por último el Radio del "Cono de la Penumbra" de Mercurio l₁ (JO) sobre el Plano Fundamental o Principal de Referencia

l₁ = i₁ * c₁     (14)

Hora Sidérea Aparente y Ángulo Horario de Z en Greenwich (ω = 0)

Primero tomar de las Efemérides, por ejemplo del The Astonomical Almanac, la Hora Sidérea Aparente μ (00:00 GMT) y en [hms], luego calcularlas en [°] multiplicándolas por 15 (μ⁰). Después hallar el Ángulo Horario μ₁ siendo el ángulo horario del punto Z o Eje del "Cono de la Sombra" de Mercurio también en Greenwich (ω = 0)

μ₁ = μ⁰ - a     (15)

Coordenadas Rectangulares del Observador ξ,η,ζ

Calcular primero φ' en [°], ρ en [Radios Terrestres], θ en [°], B en [°] y A

φ' es la Latitud Geocéntrica según la siguiente fórmula

φ' = Atan((Seno(φ) * (1 - e^2)) / Coseno(φ))     (16)

donde φ es la Latitud Geográfica y e es la Excentricidad Terrestre según la tabla de las Constantes mencionada más arriba, o bien

e = (1 - b^2 / a^2)^0,5     (17)

siendo a el Radio Ecuatorial Terrestre y b el Radio Polar Terrestre, ambos ejes de la Elipsoide Terrestre también dados en la tabla de las Constantes mencionada más arriba.

ρ es la distancia en [Radios Terrestres] desde el Centro de la Tierra hasta la Latitud del Observador, siendo prácticamente constante en todo el tránsito

ρ = 0,99833132881 + 0,0007271 * Coseno(2 * φ) - 0,0000018 * Coseno(4 * φ)     (18)

θ = μ₁ - ω     (19)

donde ω es la Longitud del Observador hacia el Oeste (W) desde Greenwich (0° - 360°) y θ el Ángulo Horario del punto Z o del Eje del "Cono de la Sombra" de Mercurio.

B = Atan(Seno(φ') / (Coseno(φ') * Coseno(θ)))     (20)

A = ρ * Seno(φ') / Seno(B)     (21)

por último, calcular las Coordenadas Rectangulares del Observador ξ,η,ζ

ξ = ρ * Coseno(φ') * Seno(θ)     (22)

η = A * Seno(B - d)     (23)

ζ = A * Coseno(B - d)     (24)

donde ξ,η corresponden al Plano Fundamental o Principal de Referencia (equivalentes a las x,y de Mercurio).

Los valores de los cambios horarios para ξ,η,ζ están dados por

ξ' = μ' * Seno(1) * ρ * Coseno(φ') * Coseno(θ)     (25)

η' = μ' * Seno(1) * ξ * Seno(d) - d' * Seno(1) * ζ     (26)

ζ' = -μ' * Seno(1) * ξ * Coseno(d) + d' * Seno(1) * η     (27)

Los valores μ', d y d' los podemos encontrar en las tablas correspondientes (más arriba)

Ecuación Fundamental de los Tránsitos Planetarios por el Sol

Calcular la distancia Δ en [Radios Terrestres] medida desde el Observador hasta el Eje del "Cono de la Sombra" de Mercurio sobre el Plano Fundamental o Principal de Referencia, siendo la Ecuación Fundamental de los Tránsitos Planetarios por el Sol

Δ = ((x - ξ)^2 + (y - η)^2)^ 0,5    (28)

Otros Cálculos de Elementos Besselianos

Cálculo de: a₁', b' y c₁'

a₁' = -l₁' - μ' * Seno(1) * i₁ * x * Coseno(d)     (29)

b' = -y' + μ' * Seno(1) * x * Seno(d)     (30)

c₁' = x' + μ' * Seno(1) * (y * Seno(d) + i₁ * l₁ * Coseno(d))     (31)

Los valores μ', d, l₁, i₁, x, y, x' e y' los podemos encontrar en las tablas correspondientes (más arriba)

Cálculo de: E, e, F y f

E = Atan(b' / c₁')     (32)

e = b' / Seno(E)     (33)

F = Atan(d' / (μ' * Coseno(d))     (34)

f = (d' * Seno(1)) / Seno(F)     (35)

Los valores μ', d y d' los podemos encontrar en las tablas correspondientes (más arriba)

Capítulo Anterior

Contenidos

Capítulo Siguiente

Cálculo de un Tránsito

Capítulos

01	02	03
04	05	06

Notas de referencia