Assume the current opposition is aphelic (Sun, Earth and Mars line up
perfectly, in that order, and Mars is at aphelion). Calculate the distance
between Mars and Earth in AU and km, as well as in light-minutes.
Separately calculate the perihelic opposition distance when Mars is at
perihelion.
That is, give another formula, not just a number! If I want you to give a
numerical answer, I usually ask to "evaluate" something, while "find" or
"derive" means analytical answer, and "calculate" means both analytical formula
and numerical value. I am forgetful and may depart from this convension, but in
general I want you to show the derivations with some comments on what you are
doing, not just pure symbolic math, and definitely not just numbers.
Also calculate the maximum distance in 1/2 of synodic period (January next year) when Mars will be in conjunction with the sun.
You are given the following data: $a_M = 1.524$ AU, 1 AU = 150 mln km, $c = 0.3$ mln km/s (speed of light), $e_M = 0.0934$, $a_E = 1$ AU (by definition), $e_E = 0.0167 \simeq 0$ (assume it's zero for this problem, since Mars's eccentricity, $e_M\gg e_E$, is much larger and more important).
SOLUTION
We did a sketch of the situation with red chalk marking position of Mars, talked about why Mars is a red planet, and then derived distance at opposition, $d$, in two cases (a,p = aphelion, perihelion). $$d_{a,p} = a_M(1\pm e_M) - a_E$$ Units (AU or other) check ok. Now we can substitute values: $d_{a,p} = \{0.666, 0.382\}\; \mbox{AU} = \{100, 57.2\}\; \mbox{mln km} = \{333, 191\}\; \mbox{light-seconds}.$
As you see, we kept 3 accurate digits in the answer only, since this many we also had in the input data. Never quote a long string of numbers from a calculator if that precision is not warranted by the precision of all inputs.
The last units (light-sec) are useful if you send a robotic rover to Mars and are driving it. You'll need to wait at least that may second for a reaction, and twice that to notice it after the signals are sent back to Earth.
Compute the ratio of brightness (i.e. flux of radiation received) and express as magnitude difference. As in the first problem, Earth's orbit can be assumed circular.
SOLUTION
As with all the problems involving illumination, brightness and (later in the
course) temperature calculations, we first mentally draw a sphere surrounding
the source of radiation (sun) wit radius equal to the distance to illuminated
body (Mars). The bolometric (all-wavelengths) flux of energy is solar
luminosity $L_\odot$ spread evenly on that imaginary sphere, so flux of solar
radiation at Mars, $F$, diring its aphelion/perihelion, is
$$F_{a,p} = \frac{L_\odot}{4\pi a_M^2 (1\pm e_M)^2}. $$
That needs to be multiplied by the cross section of Mars and by a constant
describig how well Mars scatters radiation (to be later received on Earth)
-- we call it albedo. But these things are all constant, the same at (a,p)
points, so they'll not be important for the ratio of fluxes below.
Next, light scattered off Mars has to get to Earth, so we mentally
distribute the scattered power on a big sphere with radius $d_{a,p}$
and area $4\pi d^2_{a,p}$, where expressions for $d^2_{a,p}$ have been
derived in the previous problem.
Finally, we divide the two fluxes received on Earth to get their ratio (many
constants cancel and will not be important for the result):
$$ \frac{I_p}{I_a} = \left[\frac{ (1+e_M)(a_M(1+e_M) - a_E) }
{ (1-e_M)(a_M(1- e_M)-a_E) }\right]^2 \,.$$
Units check out ok, and the ratio evaluates to 4.43, i.e. at the next
perihelic opposition Mars should be as much as 4.43x brighter than at the
present maximum (about -1.65 magnitudo in astronomical, logarithmic scale
$\Delta m = -2.5\log I_p/I_a$). This difference in brightness is totally due
to Mar's orbiotal eccentricity.
Time span from aphelion
around now to perihelion is half its orbital period or 687/2 days, but then
Mars won't be in opposition! Oppositions happen every 780 days, so you can do
the math of when we will see the next perihelic rather than aphelic
opposition.. or, if you want, you can ask google, exceptionally!
It's in Sept 2035. Go out on a clear night, find Mars in Gemini constellation
to the east of Orion, it's a visibly reddish 'star' as bright as the
nearby Sirius in Canes Venatici constellation (Hunting Dogs of hunter
Orion). Remember how it looks, to tell your friends in 2035 that you are
following Mars since its oppositions were bright but not SO bright.
We also started discussing the comet Dibiasky from the movie, but didn't get far. See below.
Above is the proof from the lecture that the polar-coordinate ($r,\theta$)
version of the equation of ellipse: $r(\theta) = a(1-e^2) /(1+e \cos \theta)$,
is identical with the more familiar Cartesian version of ellipse as a circle of
radius $a$ (equation $x^2+y^2=a^2$), contracted by a factor $b/a$ along
y-direction: $$\frac{x^2}{a^2} + \frac{y^2}{b^2} =1. $$
The proof furnishes the aspect ratio as a function of eccentricity $e$ in the
form $b = a \sqrt{1-e^2}$, which is equivalent to fact that the sum of distances
from two focal points of ellipse to every point on the curve is constant and
equal $2a$, as we have figured out in T1.
Assume inclination of the orbit $i=90^\circ$, aphelion at $r_a = 10^4$ AU, and perihelion at $r_p = 1 $ AU. We actually could assume a smaller perihelion distance, and the results would be practically identical, but it is easiest to take $r_p=1$ AU.
Use the laws of elliptic motion (see the introduction to assignment set A1). Derive $a$ and $e$ of the orbit, and perihelion velocity $v_p$, neglecting the presence of Earth. Then find the way to take into account an increase of the velocity by the time of impact, due to gravity of Earth. In the end, find the impact speed by considering the perpendicular velocities of the Earth $v_E = 30$ km/s, and $v_p$. Express the results in km/s. Use L-conservation or E-conservation law to find the answers.
SOLUTION
This blackboard
explains how to find $a$ and $e$ from $r_p$ and $e_a$. And these
explain how to find $v_p$. This can be done in a two-step way described in the text of the problem above. Formulae we find depend on $e$, which is $e=1-2\cdot 10^{-4}$, very close to $e=1$, and can be replaced with 1. (This is equivalent to saying that the comet came not from Oort cloud but from infinity, with initial zero speed. It's a very accurate approximation to truth. It also means that the extremely elongated elliptic trajectory inside the solar system is almost indistinguishable from parabola.)
The comet gains kinetic energy because it sinks gradually deeper into 2 potential wells, the main one of the sun's gravity and a smaller one of Earth's gravity. On contact with Earth, the depths (per unit mass of comet as test particle) are $-GM/a_E$ and $-GM_E/R_E$, correspondingly. You should check the units, they're ok. Here, $a_E = 1$ AU, and $R_E$ is the radius of the Earth. So, in inertial frame, and multiplying both sides of equation by 2, two times the comet's (specific) kinetic energy reaches value $$ v^2 = \frac{2G M_\odot}{a_E}+\frac{2GM_E}{R_E}.$$ We can phrase this as follows: in inertial frame, comet's speed at impact is the quadratic sum of the escape speeds from the sun and from the surface of the planet. Numerically, we found that the escape speed from the sun at 1 AU distance is 42.4 km/s, and I claim that the escape speed from Earth is 11.1 km/s (which can be verified either by plugging in the constants and using calculator, or by multiplying near-Earth circular speed 7.9 km/s by $\sqrt{2}$. Thus $$v = \sqrt{42.4^2+11.2^2}\;\mbox{km/s} = 43.8 \;\mbox{km/s}.$$ The speeds are in the inertial frame of the sun. Now we need to shift into the moving Earth's frame of reference, to obtain the impact speed. Subtracting two perpendicular vectors, one of the Earth's motion in inertial frame (at 30 km/s or more precisely 29.78 km/s), and the other of the vertical motion of the comet having speed $v$, will accomplish that. The result is this impact speed $$v_{imp} = \sqrt{42.4^2+11.1^2+29.78^2}\;\mbox{km/s} = 53.0\;\mbox{km/s}.$$
We will deal with the consequences of such a fast impact next time.
SOLUTION
We know the impact speed, it's 53 km/s. The mass can be estimated as volume ($4\pi /3) R^3 \sim 4R^3$, where $R=4.5e3$ m, so about 3.6e11 m$^3$) times the density 1.1e3 kg/m$^3$, comet's mass is thus about 4e14 kg. $$ E_K \approx 4e14 \; (53e3)^2/2 \sim 0.5e24\; {\mbox{J}}$$
This is $0.5e24 / 2e17 ~ 2.5 $ million Tsar Bomb energies. That's the number of thermeonuclear bombs each of which can evaporate or otherwise destroy a major city, and I don't know if there are even 1000 of them on our planet (there could be 2.5 million small towns and villages, though).
A quarter of a much smaller
Meteor Crater
near Winslow, Arizona, is seen here from a low-flying airplane
(click to blow up, pun intended)
The crater is circular and has ~1 km diameter. A stony-iron meteorite,
60 m across, fell 50k yr ago with an estimated impact speed 20 km/s.
In 2024, With one u/g student from DPES we hiked uphill to the rim
of the crater from the south. We did not descend into the crater. It was an
unofficial exploration, ignoring 'No tresspassing. Terrain monitored" signs.
(The crater is privately owned.) Estimate the kinetic energy in J, show that
it was equivalent to 100 Mt TNT, 2 or 3 times higher than the Tsar Bomba.
Additional question, appropriate for introductory Physics course: If an airplane flies just inside the crater's rim at constant altitude & constant speed $v=$ 60 m/s = 216 km/h = 117 kt, at what bank angle does it fly and how many g's does the passenger feel?
The equation of motion in rotating frame contains two gravitational attractions from massive bodies 1 and 2, centrifugal force and the Coriolis force. But Coriolis force in this problem is zero, since the test particle is at rest. We have shown that for any mass parameter $\mu = m_2/(m_1+m_2)$ the equilateral triangular position has sum of 3 accelerations equal to zero.
Look at the top of this page for the schematic drawing of the R3B situation
and the proof
which was also done by you on the blackboard
Next, the orbital period of ISS is roughly 90 minutes, so the visibility should last 2*0.0544*90 min = 10 minutes, if indeed we could see ISS from rise to setting time, which is hardly possible due to clouds, haze, trajectory not passing through zenith, and the unavoidable hills and obstacles. Notice that you gave me time which is 1/2 of the correct time. In practice, you might be lucky to see ISS for 5 minutes. There are sites that will tell you when and where to look, from a given city.
SOLUTION
Ans: 6.35 Earth radii from the center, 5.35 from surface (about 36000 km).
At the speed of light, you can send+receive the signal to+from a geostationary
communications satellite in about 0.2s, if it answers instantaneously.
So you can talk with between different, nearby, continents with a small delay.
Consider black-body particles (no scattering, just absorption) of physical density $2\times $ water (2000 kg/m$^3$), having the shape of a sphere of radius $s$. Derive the radiation pressure coefficient $\beta = |F_{rad}/F_{grav}|$.
SOLUTION
On the blackboard, I lost $c$, the speed of light, after the first few equalities..
it should be in every denominator!
$$ \beta = \frac{s_0}{s}$$
where
$$ s_0 = \frac{3}{16\pi G c\rho} \left(\frac{L_*}{M_*}\right).$$
Check the units. They are OK. Substituting the values of constants, we get
$s_0 = 0.28\, \mu$m (it's a fine dust size). Also, its less than a peak
wavelength of solar light $\lambda \approx 0.4\mu$m. So the geometrical optics
cannot really tell us what happens when $s=s_0$. Which doesn't mean that there
is something wrong with the formula. We simply should not use it for
$s < 0.5\,\mu$m. Within the domain of applicability of geom. optics,
radiation pressure rises to just above 0.5 of gravity.
Below the size $s\approx \lambda$ (wavelength of light), the wave
optics applies, and results in a certain decrease of radiation presssure
coefficient $\beta$, like in the log-log plot provided. But that's outside
of our scope here. Let's just say that we're not interested in
$s \lt \lambda \approx 0.5\mu$m.
Formula $\dot{E} = f \pi s^2$ is the usual
calculation of energy deposition rate, and $\dot{E} = \dot{p}\, c$ is the
relationship of photon's energy and momentum ($\dot{p} = m\,\dot{v}$).
And $\dot{p}$ is the same as force in Newtonian dynamics.
Don't yet pay attention to trajectories on the left, they're the solution to
the next problem.
Question: What is the eccentricity $e$ of the orbit that the debris particle
enters? (In other words, derive $e(\beta)$ relationship).
If you predict that large values of eccentricity are possible, it may be
relevant to the existence of quite extended debris disks around stars:
asteroids may collide within, say, 50 AU distance, but the dust disk can
extend to 1000 AU, if $e \sim 0.9$ or more.
SOLUTION
In the tutorial, we used energy integral to derive the result.
The proof goes like this:
The particle moves in two superposed force fields: one of gravity, the other
of radiation pressure. Both have the same radial dependence $r^{-2}$
on distance from the star, so there is a constant ratio $\beta$ between the
two. Potential energy being the sum of the two potentials, reads
$$-\frac{GM(1-\beta)}{r}.$$
Kinetic energy is initially $v_K^2/2 = (1/2)\,GM/r$, following from the
equality of centripetal and gravitational accelerations (radial force
balance). The total energy $E$ is constant along trajectory, since there is
no friction, and can be written at the initial point as
$$ E = -\frac{GM(1-\beta)}{r} +\frac{GM}{2r} = const.$$
For $\beta >0$, the inital point will be the closest point to the star,
a pericenter, where $r = a(1-e)$.
One the other hand, 2-body dynamics says that $E = -GM_*/(2a)$. In our case, to take into account an effective reduction of gravity by the opposing force of radiation pressure, we write $M_*$ (effective mass) as $M_*= (1-\beta) M $. So the second equation for $E$ is $E =-GM(1-\beta)/(2a)$. Comparing with the first expression for $E$, after some algebra, we obtain $$ e(\beta) = \frac{\beta}{1-\beta}.$$
Another way to derive the same result, even simpler, is to only consider
the constancy of specific angular momentum $L$. The addition of radiation
pressure force to the debris particle does not change its initial angular
momentum, because no non-radial forces are involved. Therefore, the initial
$L = \sqrt{GMr}$ that characterized both the parent and the debris bodies
before the onset of radiation effects, equals the angular momentum constant
at any point of the new trajectory:
$$L = \sqrt{GM_* a (1-e^2)}.$$
Notice that we did take into account radiation pressure by replacing
$M$ with $M_*= M(1-\beta)$.
Using $r = a(1-e)$ and cancelling common factor $GMa$, we get
$$ (1-\beta) (1-e^2) = 1-e,$$
or
$$ e(\beta) = \frac{\beta}{1-\beta}.$$
We notice that the debris will be an a bound, elliptic, orbit if $\beta < 1/2$. In particular, it will be on a circular orbit ($e=0$) if $\beta =0$. When $\beta=1/2$, debris particle enters a parabolic escape orbit. For larger $\beta$, as $\beta \to 1$, cf. the picture above (left side), the hyperbola gradually becomes a straight line. This is so, because $\beta =1$ is the case of zero net force. In that case, Galileo's law of inertia states that the path is a line, and the velocity is constant.
For $\beta > 1$, the orbit is hyperbolic but bet away from star.
As $\beta \to \infty$, the asymptotic trajectory simply shoots out
of the point of debris creation radially, away from the star.
SOLUTION
Summarizing our efforts, effective target area is $\pi b^2$, and parameter
of encounter $b$ is in a ratio with $R$ following from energy and angular
momentum conservation:
$$ \frac{b^2}{R^2} = 1 + \frac{2GM_\odot}{R v_\infty^2}. $$
We replaced $GM_\odot/R $ with $(R_{AU}/R)(GM_\odot/R_{AU}) = (1/8)(30$ km/s)$^2$
because at $R_{AU} = 1$ AU, Earth goes at the speed of whose square is
(30 km/s)$^2 = GM_\odot/R_{AU}$. That gives target area
$ b^2/R^2 = 1 + (2/8)30^2/28^2 \approx 1.28$, which isn't much; we have
honestly proven that the
gravitational enhancement of target area with
respect to the basic $\pi R^2$ is only by a factor 1.28.
The number of interstellar objects (let's call them comets) passing through
a sphere of radius $b$ per unit time is
$$ \dot{N} = n_c v_\infty (\pi b^2) = 1.28 \pi n_c v_\infty R^2$$
It does not matter from which directions those objects are coming, the
only important thing is their spatial density (number per unit volume) $n_c$
and their speed at infinity $v_\infty$. Even though every object comes
from a different direction, the cross section of an R-radius sphere that it
is aiming at is always perpendicular to the trajectory at infinity. We
turn all the incoming directions to be along x-axis in the figure on
the blackboard below (to which you astutely objected, thank you).
The number of events is just the product of $v_\infty n_c$ and the
target area, as if object were all coming from one direction.
The above formula allows us to estimate $n_c$. And the number $n$ of
comets ejected by one average star (i.e. each planetary system) is the
ratio of the number density of ejected comets to the number density of
stars, counted in a volume of the Galaxy big enough for good statistics:
$$ n = \frac{n_c}{n_*} = \frac{ \dot{N} }{ 1.28\pi {v_\infty}n_* R^2}$$
Units check out (no units on either side). We plug values w/o convertig them
to S.I. No calculator will be needed.
$$ n = \frac{ \mbox{pc}^3}{ 1.28 \cdot 5 \mbox{yr}\, v_\infty \;\pi\, 0.15\cdot
8^2 \mbox{AU}^2}$$
But $v_\infty$ is 28/30 of Earth's orbital speed, so (1yr)$ v_\infty =
(28/30) 2\pi$ AU = $\pi$ (28/15) AU, using the $2\pi$ AU circumference of
Earth's orbit. We will also need to convert pc to AU: pc/AU = 206265 ~ 2e5.
$$ n = \frac{ 8 \cdot 10^{15} \mbox{AU}^3}{ 1.28 \cdot 5 (28/15) \;\pi^2 0.15
\cdot 8^2 \mbox{AU}^3} \approx \frac{10^{15}}{ 5 (28/15) 1.5 (1.28\cdot 8)}
\approx \frac{10^{15}}{5 \cdot 28} \sim 10^{13}.$$
Notice how we did not do order-of-magnitude rounding initially, and waited
with it till the last step. Namely, it happens that you round quantities
like, say, 2.9 and 2.5 both to 1 in an intermediate result, losing one order
of magnitude in their product: 2.5*2.9 $ \sim 10$.
Here on the right is our brave attempt at quick arithmetic during the
tutorial (inexplicably, I somehow got the final order of magnitude right,
despite dropping one instance of $\pi$ in the upper-right corner, and the
final embarassing step 4/1200 = 1/4, which would actually be
$\sim 1/10$!):
The huge number we obtained is in a good agreement with the scenario of
planetary system formation presented in lectures, where a protoplanetary
disk is able to form up to $10^{14}$ planetesimals, but the giant planets
(using up most of them for construction of the cores) can easily eject a big
part the remaining planetesimals from their vicinity to the Galaxy.
If all our numerous assumptions are true (possible, if unlikely :-),
then about 1 out of 10 small bodies that form in a protoplanetary disk is
ejected to the Galaxy, based on the theoretical number of planetesimals,
and this exercise based on empirical knowledge. I haven't checked, but
think that this ratio calculated in dynamical simulations of planetary
formation is roughly the same.
Commentary:
A good intuition (that you can support by doing actual numerical, dynamical
calculations, if you spend a month learning how to do those) says that in
addition to those that ecaped, about ten percent of $10^{13}$ actually did
not manage to escape to, but stayed in the Oort cloud, barely bound
gravitationally to the sun. That would make Oort cloud initially consist of
$10^{12}$ or a trillion, i.e. thousand billion comets. You can just call
it an astronomical number. Or economical (there is one dollar of Canada's net
national debt per one comet, securely stored in the Oort cloud).
Not all survived in Oort cloud to this day: some hit the sun (or, less likely,
a planet) and some did escape after a stellar flyby in the solar neighborhood.
Yet most of them probably remained there, because astronomers counting the
arrival of new comets (from Oort cloud) think that, indeed, Oort cloud has
about a trillion of them! (Hence: 'securely stored').
So, the seemingly completely independent numbers of arrivals from the Oort
cloud and from the Galaxy, are not unrelated after all.
Everything seems
consistent and explainable by our theories, even though we can't confirm
every number observationally (e.g. we can't see Oort cloud comets in situ,
or any of the $10^{13}n_* \sim 10^{12}$ interstellar comets in every cubic
parsec, in situ). Moreover, the numbers may be of the same order of magnitude
in all planetary systems throughout the whole universe.
Solve the problem by computing the collisional time of the particles in the densest part of the disk, and comparing it with the age of the star (26 Myr).
Perpendicular optical thickness $\tau_\perp(r)$ can be approximated by the function $$ \tau_\perp = \tau_m \frac{2 x^2}{1+ x^4} $$ where $x = r/r_m$, $r_m =$ 100 AU, and $\tau_m$ is a nondimensional constant (to be computed). To understand that function better, find where and what maximum value does it attain.
SOLUTION
This picture contains a hint on what integral to write, to express
the total IR power of thermal radiation by dust, $I_{IR}$, though it's
missing the efficiency factor $1-A$ (fraction of intercepted light
that is absorbed), where $A=0.4$ denotes albedo (fraction scattered).
It's a good start, but we'll return to the full answers - prepare and
volunteer next time!
Your task was to show that this is a wrong explanation, by first (2a) deriving $v_\infty$ after ejection of a solar sail deliberately constructed to have radiation pressure coefficient $\beta = 1.5$ ... 3, around a star of sun's properties and unfurling the sail at the small distance $r_0=0.04$ AU, and then (2b) calculating how long it takes to travel at that speed the distance of $R_\odot = 8 $ kpc (sun-galactic center distance).
You managed to write (together) a proper equation for energy conservation, and we solved it and evaluated the answer: $$ v_\infty = \sqrt{2\beta -1}\,(r_0)=(\sqrt{2} ... \sqrt{5})\,30\, \mbox{km/s}\sqrt{\frac{1 \mbox{AU}}{r_0}} $$ Evaluation gives $v_\infty = 212$ ... 335 km/s. Being similar to or larger than the typical linear speed of rotation of disk galaxies (Milky Way rotates in the solar neighborhood at about 220 km/s), by appropriately directing the flight of the solar sail against the galactic rotation, you can cancel most of it, and let the apparatus plunge almost straight toward the center of the galaxy. Or direct it some other way, in any case the object will travel of for about $T_x =R_\odot/v_\infty$ until ir covers about one distance between the sun and the galactic center.
This time evaluates to about 37 Myr (millions of years).
So in a life of the Galaxy and the universe, there is plenty of time to
crisscross the Galaxy using a simple device of unfurled solar sail.
Our civilization is close to achieving the needed technology, I believe.
This outlines a problem that we discussed conceptually. Janus and Epimethus are
on orbits spaced by 50 km only, while their distance from Saturn's center is
about 151000 km. The satellites are actually larger than 50 km in size and would
collide, if they kept moving on circular orbits. But they don't and are saved
by a special kind of coorbital dance called a horseshoe orbit.
During a close (and slow) approach, they attract each other, which we showed i removes the angular momentum from the outer orbit and adds it to whoever (Janus or Epimethus) occupies the inner orbit at a time. The two orbits drift in opposite radial directions so fast that they switch outer to inner orbit or vice versa, and never come closer than distance $s$. Knowing the value of $s$ and the mass of Saturn, compute the synodic period of that horseshoe relative orbit! In other words, how often are the satellited reaching a minimum distance from each other?
SOLUTION hint
Use your knowledge of the synodic period, calculate it first without the physical interaction of the satellites, based in the given size of both initial orbits. Then shorten the time interval in proportion to how big a part of the full circle the horseshoe shown in the figure occupies.